【无标题】

最新推荐文章于 2024-07-20 17:12:48 发布

巴比巴波

最新推荐文章于 2024-07-20 17:12:48 发布

阅读量92

点赞数 1

文章标签： python

本文链接：https://blog.csdn.net/qq_44920401/article/details/139098865

版权

gdal投影栅格

from osgeo import gdal
# 输入和输出路径
input_raster = r"E:\\5.tif"
output_raster = r"E:\\5_proTest.tif"


# 打开输入栅格数据集
input_ds = gdal.Open(input_raster)
input_proj = input_ds.GetProjection()
dataset = r'E:\\_pro.tif'
refer_ds = gdal.Open(dataset)
if input_ds:
    # 获取输入数据集的投影信息
    refer_proj = refer_ds.GetProjection()
    # 创建输出数据集
    driver = gdal.GetDriverByName("GTiff")
    output_ds = driver.Create(output_raster, input_ds.RasterXSize, input_ds.RasterYSize, 1, gdal.GDT_Int16)

    if output_ds:
        # 设置输出数据集的投影信息
        output_ds.SetProjection(refer_proj)
        output_ds.SetGeoTransform(refer_ds.GetGeoTransform())

        # 进行投影转换
        gdal.ReprojectImage(input_ds, output_ds, input_proj, refer_proj, gdal.GRA_NearestNeighbour)

        # 关闭数据集
        output_ds = None
        input_ds = None
        print("投影转换完成！")
    else:
        print("无法创建输出数据集。")
else:
    print("无法打开输入数据集。")