teratail header banner
teratail header banner
質問するログイン新規登録

回答編集履歴

1

追記

2018/11/13 10:43

投稿

fana
fana

スコア12229

answer CHANGED
@@ -14,4 +14,81 @@
14
14
  > ヒストグラムの作成方法がわかりません.
15
15
 
16
16
  度数分布を求めるだけですから,難しい話ではありません。
17
- 画像のy方向サイズ分の配列でも用意して(要素は全0に初期化しておき),画像の全画素を走査して,エッジ点のy座標に対応する配列要素をインクリメントすることを繰り返せばよいでしょう.
17
+ 画像のy方向サイズ分の配列でも用意して(要素は全0に初期化しておき),画像の全画素を走査して,エッジ点のy座標に対応する配列要素をインクリメントすることを繰り返せばよいでしょう.
18
+
19
+ ---
20
+
21
+ とりあえず,ものすごーく雑に実装してみました.
22
+ (当方,pythonがダメなので,C++ですが)
23
+
24
+ ```ここに言語を入力
25
+ int main(void)
26
+ {
27
+ //画像読込
28
+ cv::Mat Src = cv::imread( "SrcImg.png", CV_LOAD_IMAGE_GRAYSCALE );
29
+ if( Src.empty() )return 0;
30
+ cv::imshow( "Src", Src );
31
+
32
+ //ヒストグラム作成
33
+ std::vector< unsigned int > Hist( Src.rows, 0 ); //y方向サイズ分の配列.全要素0で初期化.
34
+ for( int y=1; y+1<Src.rows; ++y )
35
+ {
36
+ const unsigned char *pU = Src.ptr<unsigned char>( y-1 );
37
+ const unsigned char *pC = Src.ptr<unsigned char>( y );
38
+ const unsigned char *pD = Src.ptr<unsigned char>( y+1 );
39
+ for( int x=1; x+1<Src.cols; ++x )
40
+ {
41
+ int Dx = std::max( abs( pC[x]-pC[x-1] ), abs( pC[x]-pC[x+1] ) );
42
+ int Dy = std::max( abs( pC[x]-pU[x] ), abs( pC[x]-pD[x] ) );
43
+ if( Dy>32 && Dy > Dx*4 ) //※エッジ点判定.雑すぎるが.
44
+ { ++Hist[y]; }
45
+ }
46
+ }
47
+ unsigned int MaxVal = *std::max_element( Hist.begin(), Hist.end() ); //ヒストグラムの最大値
48
+
49
+ {//表示:ヒストグラム画像
50
+ cv::Mat HistImg = cv::Mat::zeros( Src.rows, 100, CV_8U );
51
+ for( int y=0; y<Src.rows; ++y ){ cv::line( HistImg, cv::Point(0,y), cv::Point(100,y), cv::Scalar( 255.0 * Hist[y] / MaxVal ) ); }
52
+ cv::imshow( "Hist", HistImg );
53
+ }
54
+
55
+ //てきとー処理でヒストグラムの山(=線のy座標)を見つける
56
+ std::vector< int > LineY; //ここに見つけた山のy座標群を格納
57
+ {
58
+ unsigned int Thresh = MaxVal * 0.75; //閾値
59
+ int RiseY = -1; //閾値を越えた箇所記憶用.越えてないときは-1
60
+ for( int y=0; y<Hist.size(); ++y )
61
+ {
62
+ if( RiseY < 0 )
63
+ {
64
+ if( Hist[y]>=Thresh )
65
+ { RiseY = y; }
66
+ }
67
+ else
68
+ {
69
+ if( Hist[y]<Thresh )
70
+ {
71
+ LineY.push_back( (RiseY + y)/2 );
72
+ RiseY = -1;
73
+ }
74
+ }
75
+ }
76
+ }
77
+ {//表示:線検出結果
78
+ cv::Mat Result;
79
+ cv::cvtColor( Src, Result, CV_GRAY2BGR );
80
+ for( int y : LineY ){ cv::line( Result, cv::Point(0,y), cv::Point(Result.cols,y), cv::Scalar( 0,0,255 ) ); }
81
+ cv::imshow( "Lines", Result );
82
+ }
83
+
84
+ //
85
+ cv::waitKey();
86
+ return 0;
87
+ }
88
+ ```
89
+
90
+ 結果の絵.
91
+ 左から,原画,ヒストグラムを可視化したもの,線の位置を赤色で書いたもの.
92
+ (原画はてきとーにペイントで手書きしました.
93
+  線だけだと微秒なので,雑音として,小学校のときにこんな記号習ったような…?というのを何となく散りばめました.)
94
+ ![イメージ説明](7dd2d2bd34c6ccd9fd52d94dfba55252.png)