質問編集履歴
1
全コードを追加
title
CHANGED
|
File without changes
|
body
CHANGED
|
@@ -12,7 +12,7 @@
|
|
|
12
12
|
> error: OpenCV(3.4.2) C:(...): error: (-215:Assertion failed) mv[i].size == mv[0].size && mv[i].depth() == depth in function 'cv::merge'
|
|
13
13
|
|
|
14
14
|
rgbをmergeする時は問題ないのですが、ycrcbをmergeする時にエラーが出ます。
|
|
15
|
-
また、red, green, blue, y, cr,cbのサイズ、クラス
|
|
15
|
+
また、以下はred, green, blue, y, cr,cbのサイズ、クラスです。
|
|
16
16
|
> red: (340, 510) <class 'numpy.ndarray'>
|
|
17
17
|
green: (340, 510) <class 'numpy.ndarray'>
|
|
18
18
|
blue: (340, 510) <class 'numpy.ndarray'>
|
|
@@ -20,4 +20,101 @@
|
|
|
20
20
|
cr: (340, 510) <class 'numpy.ndarray'>
|
|
21
21
|
cb: (340, 510) <class 'numpy.ndarray'>
|
|
22
22
|
|
|
23
|
-
配列のサイズもクラスあっていて3チャンネルでmergeしてるので、問題ないと思うのですが、なぜエラーが出るのでしょうか？教えていただきたいです。
|
|
23
|
+
配列のサイズもクラスあっていて3チャンネルでmergeしてるので、問題ないと思うのですが、なぜエラーが出るのでしょうか？教えていただきたいです。
|
|
24
|
+
|
|
25
|
+
[追記] 以下は全コードです。
|
|
26
|
+
```python
|
|
27
|
+
import cv2
|
|
28
|
+
import numpy as np
|
|
29
|
+
from matplotlib import pyplot as plt
|
|
30
|
+
|
|
31
|
+
# カラー画像の対応 (輝度だけヒストグラム平坦化)
|
|
32
|
+
# 輝度と色味を分ける色空間 YCbCr (Y:輝度、Cb:青の色差、Cr:赤の色差) を利用するとよいです。
|
|
33
|
+
# 色味の問題は、RGB を一旦この YCbCr に変換して Y だけヒストグラム平坦化すれば解決します。
|
|
34
|
+
|
|
35
|
+
# Image画像をRGBL値に分解するメソッド
|
|
36
|
+
def getRGBLFromImage(img):
|
|
37
|
+
rgb = cv2.cvtColor(img, cv2.COLOR_BGR2RGB)
|
|
38
|
+
red, green, blue = rgb[:,:,0], rgb[:,:,1], rgb[:,:,2]
|
|
39
|
+
ycrcb = cv2.cvtColor(img, cv2.COLOR_BGR2YCrCb)
|
|
40
|
+
y,cr,cb = cv2.split(ycrcb)
|
|
41
|
+
return red, green, blue, y, cb, cr
|
|
42
|
+
|
|
43
|
+
# RGBL値のヒストグラムを個別にとるメソッド
|
|
44
|
+
def createHistogram(src,ax,col):
|
|
45
|
+
X = np.arange(256)
|
|
46
|
+
hist, bins = np.histogram(src.ravel(),bins=256,range=(0,255))
|
|
47
|
+
ax.bar(X,hist,color=col)
|
|
48
|
+
ax.set_xlim(0,255)
|
|
49
|
+
ax.set_xlabel("Pixel value", fontsize=10)
|
|
50
|
+
ax.set_ylabel("Number of pixcels", fontsize=10)
|
|
51
|
+
ax.set_title("Image(" + col + ") Histogram")
|
|
52
|
+
ax.grid()
|
|
53
|
+
return hist
|
|
54
|
+
|
|
55
|
+
#Y、Cr、Cb成分からImage画像を生成するメソッド
|
|
56
|
+
def createImageFromYCrCb(new_y,cr,cb):
|
|
57
|
+
print(new_y.shape, cr.shape, cb.shape)
|
|
58
|
+
a = cv2.merge((red, green, blue))
|
|
59
|
+
ycrcb = cv2.merge((new_y,cr,cb))
|
|
60
|
+
rgb = cv2.cvtColor(ycrcb, cv2.COLOR_YCrCb2BGR)
|
|
61
|
+
cv2.imwrite("result.jpg", rgb)
|
|
62
|
+
|
|
63
|
+
# ヒストグラムを平坦化するメソッド
|
|
64
|
+
def equalizeHistogram(src_hist, src):
|
|
65
|
+
cdf = src_hist.cumsum() # np.cumsum():要素を足し合わせて配列として出力。
|
|
66
|
+
cdf = np.array(cdf*255/(w*h))
|
|
67
|
+
print('cdf.shape',cdf.shape)
|
|
68
|
+
new_y = np.zeros((h,w))
|
|
69
|
+
for i in range(h):
|
|
70
|
+
for j in range(w):
|
|
71
|
+
new_y[i][j] = int(cdf[y[i][j]])
|
|
72
|
+
return new_y
|
|
73
|
+
|
|
74
|
+
# 画像ファイルを読み込む
|
|
75
|
+
#img = cv2.imread("image/apples.jpg")
|
|
76
|
+
#img = cv2.imread("image/castle.jpg")
|
|
77
|
+
#img = cv2.imread("image/forest.jpg")
|
|
78
|
+
|
|
79
|
+
h,w,_ = img.shape
|
|
80
|
+
|
|
81
|
+
# Image画像をR,G,B,Lに分解する
|
|
82
|
+
red, green, blue, y, cr, cb = getRGBLFromImage(img);
|
|
83
|
+
print('red.shape:',red.shape)
|
|
84
|
+
print('green.shape:',green.shape)
|
|
85
|
+
print('blue.shape:',blue.shape)
|
|
86
|
+
print('y.shape:',y.shape)
|
|
87
|
+
print('cr.shape:',cr.shape)
|
|
88
|
+
print('cb.shape:',cb.shape)
|
|
89
|
+
|
|
90
|
+
# R,G,B,Lそれぞれのヒストグラムを求め表示する
|
|
91
|
+
fig = plt.figure(figsize=(10,15),dpi=100)
|
|
92
|
+
ax1 = fig.add_subplot(2,2,1)
|
|
93
|
+
ax2 = fig.add_subplot(2,2,2)
|
|
94
|
+
ax3 = fig.add_subplot(2,2,3)
|
|
95
|
+
ax4 = fig.add_subplot(2,2,4)
|
|
96
|
+
createHistogram(red, ax1, 'red')
|
|
97
|
+
createHistogram(green, ax2, 'green')
|
|
98
|
+
createHistogram(blue, ax3, 'blue')
|
|
99
|
+
y_hist = createHistogram(y, ax4, 'y')
|
|
100
|
+
|
|
101
|
+
#ヒストグラム平坦化を行う
|
|
102
|
+
print('red.type:',type(red))
|
|
103
|
+
print('green.type:',type(green))
|
|
104
|
+
print('blue.type:',type(blue))
|
|
105
|
+
print('y.type:',type(y))
|
|
106
|
+
print('cr.type:',type(cr))
|
|
107
|
+
print('cb.type:',type(cb))
|
|
108
|
+
new_y = equalizeHistogram(y_hist, y)
|
|
109
|
+
|
|
110
|
+
# Y,Cr,Cb成分からImage画像を生成する
|
|
111
|
+
dst=createImageFromYCrCb(new_y,cr,cb);
|
|
112
|
+
cv2.imshow('result',dst)
|
|
113
|
+
cv2.waitKey(0)
|
|
114
|
+
|
|
115
|
+
# R,G,B,Lそれぞれのヒストグラムを求め表示する
|
|
116
|
+
createHistogram(red, ax1, 'red')
|
|
117
|
+
createHistogram(green, ax2, 'green')
|
|
118
|
+
createHistogram(blue, ax3, 'blue')
|
|
119
|
+
createHistogram(y, ax4, 'yellow')
|
|
120
|
+
```
|