使用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等.

这是大理石示例的输出

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

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

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

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