在 PIL 中打开 GeoTIFF 问题答案

【问题标题】：GeoTIFF issue with opening in PIL在 PIL 中打开 GeoTIFF 问题
【发布时间】：2016-10-18 14:49:31
【问题描述】：

每次我在 python（尝试 PIL、matplotlib、scipy、openCV）中打开正射影像的 GeoTIFF 图像时，图像都会出错。一些角落正在被裁剪，但图像保持其原始形状。如果我在 Photoshop 中手动将 tif 转换为 png 并加载它，它确实可以正常工作。所以看起来 PIL 在处理 tif 文件时遇到了一些问题，其中的对象没有填满整个画布。有人有解决这个问题的办法吗？

部分原图：

打开后：

【问题讨论】：

标签： python image scipy tiff geotiff

【解决方案1】：

如果你把你正在使用的图的链接放上去（如果它是免费的），那就太好了。我从here 下载了一个示例GeoTIFF 图像，并使用gdal 打开它。

geotiff.ReadAsArray() 的形状是(3, 1024, 2048) 所以我把它转换成(1024, 2048, 3) (RGB) 并用imshow 打开它：

import gdal
gdal.UseExceptions()
import matplotlib.pyplot as plt
import numpy as np

geotiff = gdal.Open('/home/vafanda/Downloads/test.tif')
geotiff_arr= geotiff.ReadAsArray()
print np.shape(geotiff_arr)
geotiff_shifted = np.rollaxis(geotiff_arr,0,3)
print "Dimension converted to: "
print np.shape(geotiff_shifted)
plt.imshow(geotiff_shifted)
plt.show()

结果：

【讨论】：

@seanzand 没问题。仅供参考，如果您对答案感到满意，您可以接受:)
@seanzand 不，我的意思是您可以接受答案中的解决方案对您有用。以this为例，如果你点击一个答案的复选标记，它将被“接受”