下面是一个示例代码,实现从 GeoTIFF 影像中读取数据并转换为 UTM 坐标系:
import gdalfrom pyproj import Proj, transform
# 读取 GeoTIFF 影像
dataset = gdal.Open("image.tif")
# 获取影像的投影坐标系
wgs84 = Proj(proj='latlong', datum='WGS84')
utm_zone = dataset.GetUTMZone()
utm = Proj(proj='utm', zone=utm_zone, datum='WGS84')
# 读取影像的像素坐标