Python+ArcPy 对栅格数据进行批量重投影

原创已于 2022-08-23 10:21:45 修改 · 3.7k 阅读

36 ·

CC 4.0 BY-SA版权

文章标签：

#python #arcgis #project management

于 2021-10-26 18:43:55 首次发布

python-GIS杂谈专栏收录该内容

10 篇文章

订阅专栏

本文介绍如何使用ArcPy模块在Python环境下批量将大量栅格数据投影转换到与指定Shapefile相同的投影。通过示例代码，展示了如何利用ArcGIS的ProjectRaster_management工具，实现对E盘data目录下所有tif文件的快速投影转换，并将结果保存到out目录。此方法适用于处理大量数据，提高工作效率。

部署运行你感兴趣的模型镜像

在对遥感图像进行投影时，常常遇到投影信息不一致的问题，因此需要对数据进行投影转换。

采用ArcGIS中Data Management --> Projection and Transformation --> Raster --> Project Raster可以对栅格文件进行逐一投影，但是当文件数量较大时，该操作并不是很方便。因此这里我们采用ArcPy将栅格数据批量转换到与shp文件相同的投影。

import os
import arcpy


arcpy.env.workspace = 'E:/data/'


def GetRaster(file_path, out_path):
    n=0
    files = os.listdir(file_path)
    for file in files:
        type = file.split(".")
        if type[-1] == "tif":
            print file
            raster = os.path.join(file_path,file)
            out_raster = os.path.join(out_path,file)
            projecttype = "Albers_Province.prj"      ## 目标投影信息的文件路径,可以是shapefile对应的prj文件
            arcpy.ProjectRaster_management(raster,out_raster, \
                                           projecttype, "NEAREST","#" , "", "","#")
            n=n+1
    print str(n)+" rasters are processed !!!"



if __name__ ==  "__main__":
    data_path = "./data/"    ## 栅格数据路径
    out_path = "./out/"      ## 输出文件路径
    GetRaster(data_path, out_path)

您可能感兴趣的与本文相关的镜像