0

这是我在这里的第一个问题，我会尽量说清楚。

我有两个 tiff 文件，其中一个是 20 x 20 km 区域和 2.93m 像素大小的 4 波段 PlanetScope 卫星图像。第二个 tiff 文件是一个单波段 Top10road 图像，其中包含像素所属道路类型的 0 到 9 标签，其中 0 是非道路。

这些文件可以通过以下链接下载： https ://www49.zippyshare.com/v/5pR3GPv3/file.html

RGBN 图像似乎有 5 条形状（6496、6723）和 uint16 类型。第五个波段只包含一个值，我不确定它包含什么。

我的目标是运行 CNN 和 U-net，通过语义分割从卫星图像中提取道路并比较结果。（这只是一个样本集，最终集将是 140 x 140 公里，但是我正在这个较小的子集上进行一些测试）

我无法查看 RGBN 图像。我已经尝试过 tifffile、PIL 并最终将每个波段导入为一个 numpy 数组，然后用 pyplot 将它们可视化。后者为每个波段单独生成了一个图像，但是我无法堆叠所有波段并查看 RGBN 图像。

以下代码是我尝试过的简化版本：

from matplotlib import pyplot as plt
import numpy
import rasterio
with rasterio.open('ede_subset_293_20180502_planetscope.tif') as src:
    r, g, b, n, x = src.read()

#works
plt.imshow(r, interpolation='nearest')
plt.show()

#does not work :(
rgbn = np.dstack((r,g,b,n))
plt.imshow(rgbn, interpolation='nearest')

rgb = np.dstack((r,g,b,n))
plt.imshow(rgb, interpolation='nearest')

预先感谢您的任何建议。

4

0 回答 0