Arcpy批量投影shp

最新推荐文章于 2024-06-22 15:00:06 发布

蒋户川柯基

最新推荐文章于 2024-06-22 15:00:06 发布

阅读量331

点赞数

文章标签： python

本文链接：https://blog.youkuaiyun.com/m0_46663205/article/details/129747656

版权

该脚本在D:/Desktop/diantongjishp目录下工作，通过arcpy模块遍历所有.shp文件，将它们投影到WGS1984UTMZone47N（32647）坐标系，并保存到diantongjishppro目录下，每个文件名后添加_pro作为标识。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

# -*- coding: utf-8 -*-
#批量投影

import arcpy
import os

# Set environment settings
arcpy.env.workspace = "D:/Desktop/diantongjishp"
inspace = "D:\\Desktop\\diantongjishp\\"

list_shapefiles = arcpy.ListFiles("*shp")
#print list_shapefiles

outspace = "D:\\Desktop\\diantongjishppro\\"

outCS = arcpy.SpatialReference(32647)
            
for shp_name in list_shapefiles:

    print shp_name

    proshp = inspace + shp_name
    
    #投影后输出文件
    outname = outspace + shp_name.split('.')[0] + '_pro'+'.shp'
    
    #批量投影
    arcpy.Project_management(proshp, outname, outCS)
    print ("         ---- OK！ ----         ")
print("   ---   All project is OK!   ---   ")