Python 中的 GDAL 库可以用于地理数据处理,其中 warp 功能提供了重投影功能,可以将数据从一个投影系统转换到另一个投影系统。
使用 gdal.Warp() 函数可以实现对数据的重投影,该函数的语法如下:
gdal.Warp(destNameOrDestDS, srcDSOrSrcDSTab, options = options)
例如,以下代码演示了如何将一个 GeoTIFF 文件从 WGS 84 投影系统转换到 UTM 坐标系统:
from osgeo import gdal
读取 GeoTIFF 文件
src_ds = gdal.Open('input.tif')
设置目标文件名称
dst_file = 'output.tif'
设置重投影选项
options = gdal.WarpOptions(dstSRS='EPSG:32632', format='GTiff')
进行重投影
gdal.Warp(dst_file, src_ds, options=options)