RGB画像の作成方法

RGB画像

RGB画像を作る時は三つのバンドを利用する。
それぞれのバンドを青、緑、赤に着色して合成カラー画像を作成する。
ここではバイナリデータ形式の大気上端反射率を使用している。
Sentinel2Aにおいて私が使ったバンドは下記の通りである。

VN02を青の波長帯
VN03を緑の波長帯
VN04を赤の波長帯
VN08を近赤外の波長帯
VN11,SW03なども利用した。

#ture color　本来の色が見える
rgb_3d = np.zeros((nxh,nyh,3), dtype=np.uint8)
rgb_3d[:, :, 0]=brefv04[:, :] #Red
rgb_3d[:, :, 1]=brefv03[:, :] #Green
rgb_3d[:, :, 2]=brefv02[:, :] #Blue
img = Image.fromarray(rgb_3d, 'RGB')
# natural color　植物の特徴的に示すことができる
rgb_3d = np.zeros((nxh,nyh,3), dtype=np.uint8)
rgb_3d[:, :, 0]=brefv04[:, :] #Red
rgb_3d[:, :, 1]=brefv08[:, :] #Green
rgb_3d[:, :, 2]=brefv03[:, :] #Blue
img = Image.fromarray(rgb_3d, 'RGB')
#false color　
rgb_3d = np.zeros((nxh,nyh,3), dtype=np.uint8)
rgb_3d[:, :, 0]=brefv08[:, :] #Red
rgb_3d[:, :, 1]=brefv04[:, :] #Green
rgb_3d[:, :, 2]=brefv03[:, :] #Blue
img = Image.fromarray(rgb_3d, 'RGB')
#よく利用したパターン（色が鮮やかに出る）
rgb_3d = np.zeros((nxh,nyh,3), dtype=np.uint8)
rgb_3d[:, :, 0]=brefsw03[:, :] #Red
rgb_3d[:, :, 1]=brefv11[:, :] #Green
rgb_3d[:, :, 2]=brefv08[:, :] #Blue
img = Image.fromarray(rgb_3d, 'RGB')

plt.imshow(img, origin='upper')
plt.show()

今日はここまで！