python 读取geotiff_使用geotiff中的gdal python读取高程

如果您希望将所有高程值读取为numpy数组，则通常会执行以下操作：

from osgeo import gdal

gdal.UseExceptions()

ds = gdal.Open('test_data.tif')

band = ds.GetRasterBand(1)

elevation = band.ReadAsArray()

print elevation.shape

print elevation

elevation 将是一个2D numpy数组 . 如果您想快速绘制值，可以使用 matplotlib ：

import matplotlib.pyplot as plt

plt.imshow(elevation, cmap='gist_earth')

plt.show()

如果你想看到一个带有正确的* x，y坐标的图，你会做类似的事情：

nrows, ncols = elevation.shape

# I'm making the assumption that the image isn't rotated/skewed/etc.

# This is not the correct method in general, but let's ignore that for now

# If dxdy or dydx aren't 0, then this will be incorrect

x0, dx, dxdy, y0, dydx, dy = ds.GetGeoTransform()

x1 = x0 + dx * ncols

y1 = y0 + dy * nrows

plt.imshow(elevation, cmap='gist_earth', extent=[x0, x1, y1, y0])

plt.show()