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

Posted 漠里

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了ArcGIS + Python 批量裁剪添加X/Y坐标脚本相关的知识,希望对你有一定的参考价值。

前言


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

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

正文

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

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

import arcpy
import glob
import os
 
arcpy.CheckOutExtension(Spatial)
 
# 指定裁剪前的目录
inws = r"D:PythonClipoldData"
 
# 指定裁剪后的存放目录
outws = r"D:PythonClipNewdata"
 
#指定shp范围边界文件,即目标区域的边界
mask = r"D:PythonClippanhuo.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 arcpy

def 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坐标脚本的主要内容,如果未能解决你的问题,请参考以下文章

分享一个ArcGIS批量裁剪批量投影批量提取值至点的工具

分享一个ArcGIS批量裁剪批量投影批量提取值至点的工具

arcgis模型构建器使用(二):按要素属性批量裁剪栅格数据

如何利用arcgis批量裁切矢量数据

ArcGIS风暴ArcGIS 10.6创建规则格网(渔网fishnet)图文经典详解

Arcgis之矢量数据裁剪