坐标系转换 tif_WGS84 to UTM

Python实现tif_WGS84到UTM坐标系转换
import rasterio
from rasterio.warp import calculate_default_transform, reproject, Resampling
from pyproj import CRS

# 输入输出路径
src_path = r"F:\EGM08_gravity_data\EGM08_gravity_data.tif"
dst_path = r"F:\EGM08_gravity_data\EGM08_gravity_data_utm.tif"

# 定义自定义横向墨卡托投影(中央经线-45°,标准纬线70°N)
dst_crs = CRS.from_proj4("""
    +proj=tmerc +ellps=WGS84 +lat_0=90 +lon_0=-45 +k=0.994 
    +x_0=2000000 +y_0=2000000 +units=m +no_defs
""")

with rasterio.open(src_path) as src:
    # 计算转换参数(优化极区算法)
    transform, width, height = calculate_default_transform(
        src.crs, dst_crs,
        src.width, src.height,
        *src.bounds,
        densify_pts=100  # 加密计算点避免北极形变[5](@ref)
    )

    # 配置输出元数据
    kwargs = src.meta.copy()
    kwargs.update({
        'crs': dst_crs,
        'transform': transform,
        'width': width,
        'height': height,
        'compress': 'lzw'  # 启用LZW压缩减少文件体积
    })

    # 执行重投影
    with rasterio.open(dst_path, 'w', ** kwargs) as dst:
        for i in range(1, src.count + 1):
            reproject(
                source=rasterio.band(src, i),
                destination=rasterio.band(dst, i),
                src_transform=src.transform,
                src_crs=src.crs,
                dst_transform=transform,
                dst_crs=dst_crs,
                resampling=Resampling.cubic_spline  # 三次样条插值保持重力场连续性[9](@ref)
            )
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值