ArcGIS + Python 批量裁剪、添加X/Y坐标脚本，arcgispython,前言前一段时间，同事

ArcGIS + Python 批量裁剪、添加X/Y坐标脚本，arcgispython,前言前一段时间，同事

前言

前一段时间，同事拿来的数据范围太大，用不了那么多（只需要一个乡镇的，结果拿来区县的），太多了加载也是问题。所以就让我给处理下。

由于文件较多，手动裁剪的话，我一个一个用ArcGIS工具箱中的工具进行裁剪，总体下来用了3个小时左右。后面由于需要再次缩小范围，这次再一个一个手动的话那就太坑了。我想到可以用Python写脚本来批量处理。

正文

在安装ArcGIS软件的时候是默认安装了Python的环境的，并且安装了arcpy库。但是，这个Python是2.7的，我以前学习的是用Python3.X的，我本来想着直接用3.X，但是在“import arcpy”的地方一直有问题，找了一些方法还是不行，最后只能用2.7的了。

下面是批量裁剪的 Python 脚本：

import arcpyimport globimport os arcpy.CheckOutExtension(‘Spatial‘) # 指定裁剪前的目录inws = r"D:\PythonClip\oldData" # 指定裁剪后的存放目录outws = r"D:\PythonClip\Newdata" #指定shp范围边界文件，即目标区域的边界mask = r"D:\PythonClip\panhuo.shp"cluster_tolerance="0.0000001 DecimalDegrees" #利用glob包，将inws下的所有shp文件读存放到shps 中shps = glob.glob(os.path.join(inws, "*.shp")) #循环shps 中的所有影像，进行按掩模提取操作for shp in shps:    outname = os.path.join(outws, os.path.basename(shp))  #指定输出文件的命名方式    print outname    try:        arcpy.Clip_analysis(shp,mask,outname)    except Exception as e:        print e.message

由于文件较多，所以是循环指定目录下的shp文件，一个一个裁剪并输出保存。

在裁剪后，又需要把文件的X/Y坐标添加到属性列表里。我想这个也用脚本来写吧，好没问题下面也顺利完成：

import arcpydef CalcXY():    try:        print "Set Env"        arcpy.env.workspace = "D:\\PythonClip\\Newdata"        print "Start AddField And Calculate centroid"        shps = arcpy.ListFiles("*.shp")        for shp in shps:            fieldList = arcpy.ListFields(shp,"","Double")            fieldNames = []            for field in fieldList:                fieldNames.append(field.name)            if "X" not in fieldNames:                arcpy.AddField_management(shp,"X","DOUBLE",20,4)                arcpy.AddField_management(shp,"Y","DOUBLE",20,4)                arcpy.CalculateField_management(shp, "X", "!SHAPE.CENTROID.X!","PYTHON_9.3")                arcpy.CalculateField_management(shp, "Y", "!SHAPE.CENTROID.Y!","PYTHON_9.3")                print shp        print "AddField And Calculate centroid End"    except Exception as e:        print e.message        CalcXY()CalcXY()

以上是添加X、Y的代码。

这些简单的脚本，在ArcGIS的帮助文档中都可以找到，而且每一个都有完整的例子。有一些Python基础的都可以实现。

ArcGIS + Python 批量裁剪、添加X/Y坐标脚本