使用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中绘制二维数组的任何函数，例如cmap、imshow等

这是带有marbles example的输出

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

如果只查看图像的第3条，则会得到以下结果：

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

在MPL和2D阵列函数中进一步研究图像查看。。。