栅格数据批量投影

---

写在前面

对地理学及相关领域的同学来说，基本都是要涉及处理栅格数据的（这里以tif格式的数据为例）。有时候可能会需要对一些相同类型的栅格数据进行批量投影（数据本身是地理坐标系or其他类型的投影坐标系）。
这时候用代码来操作无疑是个较为快速的选择。

---

提示：以下是本篇文章正文内容，下面案例可供参考

一、关于EPSG代码

不同的EPSG代码可以对应不同的投影坐标系，具体介绍可以参考【EPSG码与常见的我国投影坐标系对应关系】这篇博文，可以在epsg.io网站上直接查询，或者根据上面博文中提供的方法直接计算，都很方便。

网站首页如下

二、案例代码

代码如下：

import os
from osgeo import gdal

# 指定待投影的文件夹路径
input_folder = r'D:\work'#示例路径（绝对路径）
# 指定投影后保存的文件夹路径
output_folder = r'E:\Tmp'
# 指定目标投影坐标系
target_crs = 'EPSG:4549'  # CGCS2000_3_Degree_GK_CM_90E的EPSG代码【这里根据自己的需求进行修改🙂】

# 获取文件夹下所有GeoTIFF文件
raster_files = [f for f in os.listdir(input_folder) if f.endswith('.tif')]

# 批量进行栅格投影
for file in raster_files:
    input_path = os.path.join(input_folder, file)
    output_path = os.path.join(output_folder, file)

    # 执行投影转换
    gdal.Warp(output_path, input_path, dstSRS=target_crs,xRes=30.066966, yRes=30.066966)#指定输出栅格数据的分辨率（即像元大小）
   # ⚠️这里记得设置一下，因为我发现如果不设置，后面投影完成之后分辨率会变，到时候再进行重采样比较麻烦🌚

print("栅格数据投影完成。")

---

提示：

可以在ArcGIS中查看栅格的像元分辨率

总结

个人感觉是一个简单但算是经常会用到的工具，希望能对大家有所帮助！