使用python语言,编写一段代码,通过gadl读取tif影像数据并将wgs84坐标进行UTM投影
时间: 2023-02-27 15:29:10 浏览: 190
下面是一个示例代码,实现从 GeoTIFF 影像中读取数据并转换为 UTM 坐标系:
```
import gdal
from 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')
# 读取影像的像素坐标值
cols = dataset.RasterXSize
rows = dataset.RasterYSize
geotransform = dataset.GetGeoTransform()
# 将像素坐标转换为经纬度坐标
for row in range(rows):
for col in range(cols):
x = geotransform[0] + col * geotransform[1] + row * geotransform[2]
y = geotransform[3] + col * geotransform[4] + row * geotransform[5]
lon, lat = transform(utm, wgs84, x, y)
# 在此处可以进行其他操作,例如存储转换后的坐标值
```
在代码中,我们使用了 GDAL 和 pyproj 库。GDAL 是一个开源的数字地理信息系统(GIS)库,可以读取和处理多种格式的地理数据;pyproj 库则提供了许多投影转换功能。
通过这段代码,您可以将 WGS84 坐标系的经纬度坐标转换为 UTM 坐标系。
阅读全文