Python 如何将对影像建立高程格网

Python 可以使用 GDAL 库来处理影像和高程数据。具体步骤如下：

1. 导入 GDAL 库

from osgeo import gdal

2. 打开影像文件

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

3. 获取影像的地理参考信息

geotransform = ds.GetGeoTransform()

4. 获取影像的投影信息

projection = ds.GetProjection()

5. 读取影像数据

band = ds.GetRasterBand(1)
data = band.ReadAsArray()

6. 建立高程格网

import numpy as np
from scipy.interpolate import griddata

# 获取影像的坐标信息
x_size = ds.RasterXSize
y_size = ds.RasterYSize
x = np.linspace(geotransform[0], geotransform[0] + geotransform[1] * x_size, x_size)
y = np.linspace(geotransform[3], geotransform[3] + geotransform[5] * y_size, y_size)
xx, yy = np.meshgrid(x, y)

# 获取影像的高程信息
z = data.flatten()

# 建立高程格网
xi = np.linspace(x.min(), x.max(), 1000)
yi = np.linspace(y.min(), y.max(), 1000)
xi, yi = np.meshgrid(xi, yi)
zi = griddata((x, y), z, (xi, yi), method='cubic')

7. 可以将高程格网保存为 GeoTIFF 文件

driver = gdal.GetDriverByName('GTiff')
out_ds = driver.Create('elevation.tif', 1000, 1000, 1, gdal.GDT_Float32)
out_ds.SetGeoTransform((xi.min(), (xi.max()-xi.min())/1000, 0, yi.min(), 0, (yi.max()-yi.min())/1000))
out_ds.SetProjection(projection)
out_ds.GetRasterBand(1).WriteArray(zi)
out_ds.FlushCache()