python - 创建 RGB 复合 SAR 图像

Question

我是 Python 编程的新手，我需要你的帮助。在发布之前，我总是先对我的问题进行研究。

我有 tiff 格式的 SAR 双偏振图像(2^16 灰度值)。在这个 tiff 图像中有两个波段。第一个波段 (HH_band) 是水平极化 channel ，第二个波段 (HV_band) 是垂直极化 channel 。我想创建一个 RGB 合成图像。为此，我需要按如下方式分层堆叠两个 channel :

1. 获取第一个波段(HH_band)
2. 获取第二个频段(HV_band)
3. 获取比率 (HH_band/HV_band)

我知道有很多人发布过类似的帖子(自然色的 RGB 合成图像)。我尝试使用 openCV 库中的 cv2.merge 或 cv2.split 但没有用。我认为用 Python 创建 SAR RGB 图像会相对容易(正如我看到的一些关于创建 LANDSAT RGB 图像的帖子)，但我陷入了困境。

如果有任何帮助，我将不胜感激。

Answer 1

这是一种以编程方式完成乐队组合的可能方法:

import numpy as np

tif = io.imread('dual_polarization_image.tif')

band = {'HH': 0, 'HV': 1}

r = tif[:, :, band['HH']]
g = tif[:, :, band['HV']]

hh = r.astype(np.float64)
hv = g.astype(np.float64)

b = np.divide(hh, hv, out=np.zeros_like(hh), where=hv!=0)

rgb = np.dstack((r, g, b.astype(np.uint16)))

备注:

• 可以通过简单地重新定义字典 ​​band 的值来处理 TIFF 图像中波段的不同排列。
• 在计算波段比之前，需要将数据转换为 np.float64。
• 我利用了 universal functions 的 where 选项以避免出现零除法警告。
• 为了使合成成为可能，波段比率(蓝色 channel )必须转换回与原始波段(红色和绿色 channel )相同的类型(即 np.uint16) ).