teratail
回答率
85.36%
teratail header banner
teratail header banner
質問するログイン新規登録
トップOpenCVに関する質問２つの値を元に画像上にプロットする　python opencv

編集履歴

回答編集履歴

2

たぶんOK

2020/01/05 09:37

投稿

退会済みユーザー
test CHANGED
@@ -1,4 +1,336 @@
1
- **2020-01-05追補**
1
+ **2020-01-05追補2**
2
+
3
+ ![上--(重心X,矩形Y)、下--根拠](dbedcf3c1baf96774d25eb95f4549952.png)
4
+
5
+ 上 (重心X,矩形Y)を青丸で表示
6
+
7
+ 下 根拠を図示
8
+
9
+
10
+
11
+ ```Python3
12
+
13
+ import cv2
14
+
15
+ import numpy as np
16
+
17
+ from matplotlib import pyplot as plt
18
+
19
+
20
+
21
+
22
+
23
+
24
+
25
+ src= cv2.imread("kuro1.png")
26
+
27
+
28
+
29
+
30
+
31
+
32
+
33
+ #ユニフォーム抽出
34
+
35
+ dst1 = cv2.inRange(src, (3,3,3), (10,8,10))
36
+
37
+ dst2 = cv2.blur(dst1, ksize=(5,5))
38
+
39
+ ret,thresh1 = cv2.threshold(dst2,127,255,cv2.THRESH_BINARY)
40
+
41
+
42
+
43
+ # ラベリング処理（x座標）
44
+
45
+ label = cv2.connectedComponentsWithStats(thresh1)
46
+
47
+
48
+
49
+ stats = label[2]
50
+
51
+ area = stats[:, cv2.CC_STAT_WIDTH] * stats[:, cv2.CC_STAT_HEIGHT]
52
+
53
+ top2_idx = area.argsort()[-3:-1] # 2番目、3番めに面積が大きいラベル (1番目は背景なので除く)
54
+
55
+
56
+
57
+ # オブジェクト情報を項目別に抽出
58
+
59
+ data = label[2]
60
+
61
+ center = label[3]
62
+
63
+
64
+
65
+ # ラベリング結果書き出し用に二値画像をカラー変換
66
+
67
+ color_src = cv2.cvtColor(thresh1, cv2.COLOR_GRAY2BGR)
68
+
69
+
70
+
71
+ # オブジェクト情報を利用してラベリング結果を表示
72
+
73
+ lst_pos = []
74
+
75
+ for i in top2_idx:
76
+
77
+ # 各オブジェクトの外接矩形を赤枠で表示
78
+
79
+ x0 = data[i][0]
80
+
81
+ y0 = data[i][1]
82
+
83
+ x1 = data[i][0] + data[i][2]
84
+
85
+ y1 = data[i][1] + data[i][3]
86
+
87
+ cv2.rectangle(color_src, (x0, y0), (x1, y1), (0, 0, 255))
88
+
89
+
90
+
91
+ # 各オブジェクトの重心座標をに黄文字で表示
92
+
93
+ cv2.putText(color_src, "X: " + str(int(center[i][0])), (x1 - 10, y1 + 15), cv2.FONT_HERSHEY_PLAIN, 1, (0, 255, 255))
94
+
95
+
96
+
97
+ # 追加コード: 重心座標を書き出す
98
+
99
+ # img = thresh1[y0:y1,x0:x1]
100
+
101
+ # lst_pos.append((y0,y1,x0,x1,int(center[i][0])))
102
+
103
+ lst_pos.append((int(center[i][0]),int(center[i][1])))
104
+
105
+
106
+
107
+
108
+
109
+ # x座標の位置
110
+
111
+ cv2.imshow('x',color_src)
112
+
113
+
114
+
115
+ #「背景」となる画像の取り込み（グレースケール）
116
+
117
+ img_src01 = cv2.imread("zero1.png", 0)
118
+
119
+
120
+
121
+ #「差分」をもった画像の取り込み（グレースケール）
122
+
123
+ img_src02 = cv2.imread("kuro1.png", 0)
124
+
125
+
126
+
127
+ #「背景差分」計算用オブジェクトの作成　
128
+
129
+ # bgObj = cv2.bgsegm.createBackgroundSubtractorMOG()
130
+
131
+ # 修正点: ただただ手元のOpenCVでエラーが出ないようにするための目的で変更。
132
+
133
+ bgObj = cv2.createBackgroundSubtractorMOG2()
134
+
135
+
136
+
137
+ #差分となっている「前景領域」に対してマスクをかける
138
+
139
+ fgmask = bgObj.apply(img_src01)
140
+
141
+ fgmask = bgObj.apply(img_src02)
142
+
143
+
144
+
145
+ # ラベリング処理(y座標)
146
+
147
+ label2 = cv2.connectedComponentsWithStats(fgmask)
148
+
149
+
150
+
151
+ ### 追加した部分1
152
+
153
+ stats2 = label2[2]
154
+
155
+ area2 = stats2[:, cv2.CC_STAT_WIDTH] * stats2[:, cv2.CC_STAT_HEIGHT]
156
+
157
+ top2_idx2 = area2.argsort()[-3:-1] # 2番目、3番めに面積が大きいラベル (1番目は背景なので除く)
158
+
159
+
160
+
161
+ # オブジェクト情報を項目別に抽出
162
+
163
+ data2 = label2[2]
164
+
165
+ center2 = label2[3]
166
+
167
+
168
+
169
+ # ラベリング結果書き出し用に二値画像をカラー変換
170
+
171
+ color_src2 = cv2.cvtColor(fgmask, cv2.COLOR_GRAY2BGR)
172
+
173
+
174
+
175
+ # 追加コード: 適当に追加
176
+
177
+ img_col = cv2.imread("zero1.png")
178
+
179
+
180
+
181
+ # オブジェクト情報を利用してラベリング結果を表示
182
+
183
+ for i in top2_idx2:
184
+
185
+ # 各オブジェクトの外接矩形を赤枠で表示
186
+
187
+ x0 = data2[i][0]
188
+
189
+ y0 = data2[i][1]
190
+
191
+ x1 = data2[i][0] + data2[i][2]
192
+
193
+ y1 = data2[i][1] + data2[i][3]
194
+
195
+ cv2.rectangle(color_src2, (x0, y0), (x1, y1), (0, 0, 255))
196
+
197
+
198
+
199
+ # 各オブジェクトの重心座標をに黄文字で表示
200
+
201
+ cv2.putText(color_src2, "Y: " + str(int(y1)), (x1 - 10, y1 + 30), cv2.FONT_HERSHEY_PLAIN, 1, (0, 255, 255))
202
+
203
+
204
+
205
+ ##############################################
206
+
207
+ # 内外判定
208
+
209
+ poly = np.array([[int(x0),int(y0)], [int(x0),int(y1)], [int(x1),int(y1)], [int(x1),int(y0)]])
210
+
211
+ for j in range(len(lst_pos)):
212
+
213
+
214
+
215
+ if cv2.pointPolygonTest( poly, lst_pos[j], False)== 1:
216
+
217
+ # print("in")
218
+
219
+ # 青で描画
220
+
221
+ img_col = cv2.circle(img_col,(lst_pos[j][0],y1), 10, (255,0,0), 2)
222
+
223
+ # else:
224
+
225
+ # print("out")
226
+
227
+ ##############################################
228
+
229
+
230
+
231
+ cv2.imshow('img_col',img_col)
232
+
233
+
234
+
235
+ cv2.imshow('y',color_src2)
236
+
237
+ ##
238
+
239
+ ##
240
+
241
+ ### この下からコートの画像に値を点で表したい！！！！！！！！！！！！！！！
242
+
243
+ ##
244
+
245
+ ##midori=cv2.imread("coat2.png")
246
+
247
+ ##
248
+
249
+ ### ！！！！！！！！！！！新たにここを追加！！！！！！！！！！
250
+
251
+ ### この座標部分をどう直したら良いかがわかりません！！！！！！！！！！！
252
+
253
+ ##
254
+
255
+ ##
256
+
257
+ ### 変更点: coatが未定義のため適当に変更
258
+
259
+ ### color_src3 = cv2.cvtColor(coat, cv2.COLOR_GRAY2BGR)
260
+
261
+ ### coat3=cv2.circle(color_src3,(int(center[i][0]),int(top2_idx2(y1))), 63, (0,0,255), -1)
262
+
263
+ ##
264
+
265
+ ### 変更点: 画像を適当に読み込み
266
+
267
+ ### img_col = cv2.imread("zero1.png")
268
+
269
+ ##img_col2 = img_col.copy()
270
+
271
+ ##
272
+
273
+ ### 変更点: Centerを赤で描画
274
+
275
+ ##for i in range(len(center)):
276
+
277
+ ## img_col2 = cv2.circle(img_col2,(int(center[i][0]),int(center[i][1])), 3, (0,0,255), -1)
278
+
279
+ ##
280
+
281
+ ##
282
+
283
+ ### 変更点: 矩形の頂点を緑で描画
284
+
285
+ ##for i in top2_idx2:
286
+
287
+ ## x0 = data2[i][0]
288
+
289
+ ## y0 = data2[i][1]
290
+
291
+ ## x1 = data2[i][0] + data2[i][2]
292
+
293
+ ## y1 = data2[i][1] + data2[i][3]
294
+
295
+ ## img_col2 = cv2.circle(img_col2,(int(x0),int(y0)), 3, (0,255,0), -1)
296
+
297
+ ## img_col2 = cv2.circle(img_col2,(int(x0),int(y1)), 3, (0,255,0), -1)
298
+
299
+ ## img_col2 = cv2.circle(img_col2,(int(x1),int(y0)), 3, (0,255,0), -1)
300
+
301
+ ## img_col2 = cv2.circle(img_col2,(int(x1),int(y1)), 3, (0,255,0), -1)
302
+
303
+ ##
304
+
305
+ ### cv2.imshow("img_col",img_col)
306
+
307
+ ##
308
+
309
+ ### 変更点: Centerのデバグ用
310
+
311
+ ##img_merged = cv2.addWeighted(img_col2, 0.8, color_src, 0.2, 0)
312
+
313
+ ##cv2.imshow("img_merged",img_merged)
314
+
315
+
316
+
317
+ cv2.imwrite("img_col2_merged.png",np.vstack((img_col,img_merged)))
318
+
319
+ cv2.waitKey(0)
320
+
321
+ cv2.destroyAllWindows()
322
+
323
+
324
+
325
+ ```
326
+
327
+
328
+
329
+ ---
330
+
331
+
332
+
333
+ **2020-01-05追補1**
2
334
 
3
335
  > この画像にプログラムに書かれているラベリング処理を２回した時に求めた
4
336
 

1

追補

2020/01/05 09:36

投稿

退会済みユーザー
test CHANGED
@@ -1,3 +1,305 @@
1
+ **2020-01-05追補**
2
+
3
+ > この画像にプログラムに書かれているラベリング処理を２回した時に求めた
4
+
5
+ > x座標とy座標の値を赤い点として表示させたいです
6
+
7
+ >
8
+
9
+ > 赤い点と書いてありますがわかりやすくばつ印などでも良いです
10
+
11
+
12
+
13
+ プログラム全体の意図が良く掴めていませんが、「座標をプロットする」という点において、希望箇所（！！！の辺り）に対応コードを追加したため、これで解決できると思います。
14
+
15
+
16
+
17
+ また、赤い点、とありますが、Centerと矩形と両方とも位置を示すラベルと思いますので、色を分けて表示させました。
18
+
19
+ ※OpenCVで描画する際にバツより丸の方が実装しやすいため丸で実装しています。
20
+
21
+
22
+
23
+ ![座標追加イメージ](2ef58c920562a140c496e70be5090202.png)
24
+
25
+ 赤: Center
26
+
27
+ 緑: 矩形の頂点
28
+
29
+
30
+
31
+ コメントの
32
+
33
+ > `for i in top2_idx2`の部分を
34
+
35
+ > `j`に変えてみることで解決したりするのでしょうか？
36
+
37
+ については、今回は二重ループの出番がないので使用していません。
38
+
39
+ 二重ループで`i`とは別の変数を使いたいときに`j`の出番があります。
40
+
41
+
42
+
43
+ 今回は！！！の部分にfor文を追加しましたが、ラベル作成時にこの処理もまとめてやってしまった方がfor文が少なくなりますので、より早くそして、バグの発生する余地の少ないコードになると思います。
44
+
45
+
46
+
47
+ ```Python3
48
+
49
+ import cv2
50
+
51
+ import numpy as np
52
+
53
+ from matplotlib import pyplot as plt
54
+
55
+
56
+
57
+ src= cv2.imread("kuro1.png")
58
+
59
+
60
+
61
+ #ユニフォーム抽出
62
+
63
+ dst1 = cv2.inRange(src, (3,3,3), (10,8,10))
64
+
65
+ dst2 = cv2.blur(dst1, ksize=(5,5))
66
+
67
+ ret,thresh1 = cv2.threshold(dst2,127,255,cv2.THRESH_BINARY)
68
+
69
+
70
+
71
+ # ラベリング処理（x座標）
72
+
73
+ label = cv2.connectedComponentsWithStats(thresh1)
74
+
75
+
76
+
77
+ stats = label[2]
78
+
79
+ area = stats[:, cv2.CC_STAT_WIDTH] * stats[:, cv2.CC_STAT_HEIGHT]
80
+
81
+ top2_idx = area.argsort()[-3:-1] # 2番目、3番めに面積が大きいラベル (1番目は背景なので除く)
82
+
83
+
84
+
85
+ # オブジェクト情報を項目別に抽出
86
+
87
+ data = label[2]
88
+
89
+ center = label[3]
90
+
91
+
92
+
93
+ # ラベリング結果書き出し用に二値画像をカラー変換
94
+
95
+ color_src = cv2.cvtColor(thresh1, cv2.COLOR_GRAY2BGR)
96
+
97
+
98
+
99
+ # オブジェクト情報を利用してラベリング結果を表示
100
+
101
+ for i in top2_idx:
102
+
103
+ # 各オブジェクトの外接矩形を赤枠で表示
104
+
105
+ x0 = data[i][0]
106
+
107
+ y0 = data[i][1]
108
+
109
+ x1 = data[i][0] + data[i][2]
110
+
111
+ y1 = data[i][1] + data[i][3]
112
+
113
+ cv2.rectangle(color_src, (x0, y0), (x1, y1), (0, 0, 255))
114
+
115
+
116
+
117
+ # 各オブジェクトの重心座標をに黄文字で表示
118
+
119
+ cv2.putText(color_src, "X: " + str(int(center[i][0])), (x1 - 10, y1 + 15), cv2.FONT_HERSHEY_PLAIN, 1, (0, 255, 255))
120
+
121
+
122
+
123
+ cv2.imshow('x',color_src)
124
+
125
+
126
+
127
+ #「背景」となる画像の取り込み（グレースケール）
128
+
129
+ img_src01 = cv2.imread("zero1.png", 0)
130
+
131
+
132
+
133
+ #「差分」をもった画像の取り込み（グレースケール）
134
+
135
+ img_src02 = cv2.imread("kuro1.png", 0)
136
+
137
+
138
+
139
+ #「背景差分」計算用オブジェクトの作成　
140
+
141
+ # bgObj = cv2.bgsegm.createBackgroundSubtractorMOG()
142
+
143
+ # 修正点: ただただ手元のOpenCVでエラーが出ないようにするための目的で変更。
144
+
145
+ bgObj = cv2.createBackgroundSubtractorMOG2()
146
+
147
+
148
+
149
+ #差分となっている「前景領域」に対してマスクをかける
150
+
151
+ fgmask = bgObj.apply(img_src01)
152
+
153
+ fgmask = bgObj.apply(img_src02)
154
+
155
+
156
+
157
+ # ラベリング処理(y座標)
158
+
159
+ label2 = cv2.connectedComponentsWithStats(fgmask)
160
+
161
+
162
+
163
+ ### 追加した部分1
164
+
165
+ stats2 = label2[2]
166
+
167
+ area2 = stats2[:, cv2.CC_STAT_WIDTH] * stats2[:, cv2.CC_STAT_HEIGHT]
168
+
169
+ top2_idx2 = area2.argsort()[-3:-1] # 2番目、3番めに面積が大きいラベル (1番目は背景なので除く)
170
+
171
+
172
+
173
+ # オブジェクト情報を項目別に抽出
174
+
175
+ data2 = label2[2]
176
+
177
+ center2 = label2[3]
178
+
179
+
180
+
181
+ # ラベリング結果書き出し用に二値画像をカラー変換
182
+
183
+ color_src2 = cv2.cvtColor(fgmask, cv2.COLOR_GRAY2BGR)
184
+
185
+
186
+
187
+ # オブジェクト情報を利用してラベリング結果を表示
188
+
189
+ for i in top2_idx2:
190
+
191
+ # 各オブジェクトの外接矩形を赤枠で表示
192
+
193
+ x0 = data2[i][0]
194
+
195
+ y0 = data2[i][1]
196
+
197
+ x1 = data2[i][0] + data2[i][2]
198
+
199
+ y1 = data2[i][1] + data2[i][3]
200
+
201
+ cv2.rectangle(color_src2, (x0, y0), (x1, y1), (0, 0, 255))
202
+
203
+
204
+
205
+ # 各オブジェクトの重心座標をに黄文字で表示
206
+
207
+ cv2.putText(color_src2, "Y: " + str(int(y1)), (x1 - 10, y1 + 30), cv2.FONT_HERSHEY_PLAIN, 1, (0, 255, 255))
208
+
209
+
210
+
211
+ cv2.imshow('y',color_src2)
212
+
213
+
214
+
215
+
216
+
217
+ # この下からコートの画像に値を点で表したい！！！！！！！！！！！！！！！
218
+
219
+
220
+
221
+ midori=cv2.imread("coat2.png")
222
+
223
+
224
+
225
+ # ！！！！！！！！！！！新たにここを追加！！！！！！！！！！
226
+
227
+ # この座標部分をどう直したら良いかがわかりません！！！！！！！！！！！
228
+
229
+
230
+
231
+
232
+
233
+ # 変更点: coatが未定義のため適当に変更
234
+
235
+ # color_src3 = cv2.cvtColor(coat, cv2.COLOR_GRAY2BGR)
236
+
237
+ # coat3=cv2.circle(color_src3,(int(center[i][0]),int(top2_idx2(y1))), 63, (0,0,255), -1)
238
+
239
+
240
+
241
+ # 変更点: 画像を適当に読み込み
242
+
243
+ img_col = cv2.imread("zero1.png")
244
+
245
+
246
+
247
+ # 変更点: Centerを赤で描画
248
+
249
+ for i in range(len(center)):
250
+
251
+ img_col = cv2.circle(img_col,(int(center[i][0]),int(center[i][1])), 3, (0,0,255), -1)
252
+
253
+
254
+
255
+
256
+
257
+ # 変更点: 矩形の頂点を緑で描画
258
+
259
+ for i in top2_idx2:
260
+
261
+ x0 = data2[i][0]
262
+
263
+ y0 = data2[i][1]
264
+
265
+ x1 = data2[i][0] + data2[i][2]
266
+
267
+ y1 = data2[i][1] + data2[i][3]
268
+
269
+ img_col = cv2.circle(img_col,(int(x0),int(y0)), 3, (0,255,0), -1)
270
+
271
+ img_col = cv2.circle(img_col,(int(x0),int(y1)), 3, (0,255,0), -1)
272
+
273
+ img_col = cv2.circle(img_col,(int(x1),int(y0)), 3, (0,255,0), -1)
274
+
275
+ img_col = cv2.circle(img_col,(int(x1),int(y1)), 3, (0,255,0), -1)
276
+
277
+
278
+
279
+ cv2.imshow("img_col",img_col)
280
+
281
+
282
+
283
+ # 変更点: Centerのデバグ用
284
+
285
+ # img_merged = cv2.addWeighted(img_col, 0.8, color_src, 0.2, 0)
286
+
287
+ # cv2.imshow("img_merged",img_merged)
288
+
289
+
290
+
291
+ cv2.waitKey(0)
292
+
293
+ cv2.destroyAllWindows()
294
+
295
+ ```
296
+
297
+
298
+
299
+ ---
300
+
301
+
302
+
1
303
  > 上記の２箇所のラベリング処理の部分のputtextでx座標とy座標を渡しているので
2
304
 
3
305
  > その値を新たな　zero1.pngの画像に点として表示したいです