使用GDAL在Python和Matplotlib中读取TIFF

Question

我正在尝试使用Python和MatPlotLib显示灰度TIFF文件，

到目前为止，我已经阅读了这个文件：

import scipy as N
import gdal
import sys
import matplotlib.pyplot as pyplot

try:
    tif = gdal.Open('filename.tif')
    tifArray = tif.ReadAsArray()
except:
    print 'The file does not exist.'
    sys.exit(0)

band1 = tif.GetRasterBand(1)
band2 = tif.GetRasterBand(2)
band3 = tif.GetRasterBand(3)

band1Array = band1.ReadAsArray()
band2Array = band2.ReadAsArray()
band3Array = band3.ReadAsArray()

但后来我不知道我还能做什么......我真是无能为力。 有人会帮我这个吗？

Answer 1

将文件处理成2个数组后，可以在matplotlib中使用任何绘制2D数组的函数，例如： cmap，imshow等。

以下是marbles example

的输出

img=mpimg.imread('MARBLES.TIF ')
imgplot = plt.imshow(img)

如果您只查看图像的band3，则会获得以下内容：

imgplot2 = plt.imshow(band3Array)
plt.show()

进一步了解MPL和2D阵列功能中的图像查看...