使用ArcPy对栅格数据进行批量投影坐标转换

最新推荐文章于 2024-03-26 21:27:44 发布

weixin_33805992

最新推荐文章于 2024-03-26 21:27:44 发布

阅读量2.8k

点赞数 3

文章标签： python

原文链接：https://my.oschina.net/JasonZhang/blog/667855

版权

2019独角兽企业重金招聘Python工程师标准>>>

# -*- coding:utf-8 -*- 
##========================
##批量栅格数据投影转换


import arcpy,os,os.path,init

def projectRaster(rootPath):
    try:
        
        ##arcpy工作目录
        root_path = rootPath
        arcpy.env.workspace = root_path

        ##待处理文件所在目录(相对于根目录)
        input_path = "tif-source"
        output_path = "tif-target"

        ##源坐标系 "CGCS2000_3_Degree_GK_CM_123E" 
        sourceSR = arcpy.SpatialReference("CGCS2000 3 Degree GK CM 123E")
        ##目标坐标系(WGS 1984 Web Mercator Auxiliary Sphere)
        targetSR = arcpy.SpatialReference("WGS 1984 Web Mercator (auxiliary sphere)")

        ##遍历目录，查找栅格数据
        files = os.listdir(root_path+os.sep+input_path)
        for f in files:
            if os.path.splitext(f)[1].upper() == ".TIF":
                fileName = os.path.splitext(f)[0] + ".tif"
                in_dataset = input_path + os.sep + fileName
                out_dataset = output_path + os.sep + fileName

                print "begin project "+in_dataset+" from: " +sourceSR.name+" to: "+targetSR.name
                
                arcpy.ProjectRaster_management(in_dataset, out_dataset, targetSR, "NEAREST",\
                                       "#", "#", "#",sourceSR)

        print "project success!"
        
    except arcpy.ExecuteError:
        print "Project Raster example failed."
        print arcpy.GetMessages()

################################################
if __name__ == '__main__':

    #指定处理文件根目录
    root_path = r"E:\gis"
    projectRaster(root_path)

转载于:https://my.oschina.net/JasonZhang/blog/667855