在python中绘制多光谱卫星tiff图像

时间:2017-03-06 00:36:58

标签: python plot gis tiff

我试图使用tifffile lib在python中验证/绘制卫星图像。

对于有三个频道的RGB照片,我没有遇到问题,可以按预期查看它们但是当我尝试8频道SWIR图像时,事情似乎有所破坏:

img_3 = np.transpose(tiff.imread("imageRGB.tif"),(1,2,0))
img_a = np.transpose(tiff.imread("imageA.tif"),(1,2,0))

img_3.shape: (3349L, 3396L, 3L)
img_a.shape: (3349L, 3396L, 8L)


tifffile.imshow(img_a) 

enter image description here

虽然tifffile.imshow(img_3)在一个漂亮的矩形中显示了预期的图像。 我可以绘制(部分)img_a的唯一方法是,如果我选择一个特定的频道,如:

tifffile.imshow(img_a[:,:,1]) 

我缺少什么?

1 个答案:

答案 0 :(得分:0)

specteral维度应位于第一个轴上,tiff.imshow()才能正常工作。

tiff.imshow(np.rollaxis(img_a,2,0))可以解决问题。