python - 从Geotiff使用Gdal python读取高程

我正在使用GDAL加载geotiff文件。我设法读取坐标X，Y，但未读取海拔。

以前有人处理过类似案件吗？

问候，

最佳答案

如果您希望将所有高程值读取到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()