1、准备 安装python包
apt-get update
apt-get install -y python3-gdal
pip install gdal
pip install pyproj
2、下载高程tif数据
全国各省12.5m高程数据-免费下载-资源下载-数字地球开放平台
注册下,下载他的网盘下载工具,然后选取要使用的区域下载。我下载的宁夏的。
3、python代码
以下的代码是专门针对以上的高程数据写的代码。不同的高程数据可能需要不同的投影计算方法。
使用Pyproj库进行投影转换:
通过pyproj.Proj方法,确认源投影(EPSG:4326)和目标投影(EPSG:32648)。
将输入的经纬度数据从WGS 84转换为UTM 48N坐标系统。
import gdal
import osr
from pyproj import Proj, transform
def print_geo_transform_info(transform):
print(f"Geotransform:")
print(f" Origin: ({transform[0]}, {transform[3]})")
print(f