ArcGIS Python批量裁剪、批量合并，MXD转图片PDF等源码|代码|sys|rows|视频文件|python|arcgis|import

分享至

10.4.3ArcGIS Python其他例子 10.4.3.1 计算面层的面积首先写好脚本如下：

import arcpyfrom arcpy import envimport os

import sys

fc= arcpy.GetParameterAsText(0)fieldname= arcpy.GetParameterAsText(1)shapeName = arcpy.Describe(fc).shapeFieldName

rows = arcpy.UpdateCursor(fc)

i=1;for row in rows:#记录循环feat = row.getValue(shapeName)row.setValue(fieldname,feat.area)arcpy.AddMessage("No:"+str(i)+":"+str(feat.area))rows.updateRow(row)

i=i+1;

保存成py文件，新建一个工具箱，工具箱右键添加->脚本，找到对应py，设置参数，如下：

面积字段获取自面图层，面积字段是双精度，过滤器是

如果不导入脚本，必须将对应py文件到拷贝到对应位置。否则不能使用。工具在：chp10\python\其他.tbx\UpdateCursor(计算面积)，可以直接运行，右键编辑查看代码。
10.4.3.2 更新字段值为顺序号

在：chp10\python\其他.tbx\更新字段为值顺序号，可以直接运行，右键编辑查看代码

代码如下：

import arcpyfrom arcpy import envimport os

import sys

fc= arcpy.GetParameterAsText(0)fieldname= arcpy.GetParameterAsText(1)k= arcpy.GetParameter(2) ##开始值

rows = arcpy.UpdateCursor(fc)

i=kfor row in rows:row.setValue(fieldname,i)arcpy.AddMessage("No:"+str(i))rows.updateRow(row)i=i+1;如果需要进度条代码如下：

import arcpyfrom arcpy import envimport os

import sysdefaultencoding = 'utf-8'if sys.getdefaultencoding() != defaultencoding:reload(sys)sys.setdefaultencoding(defaultencoding)

fc= arcpy.GetParameterAsText(0)fieldname= arcpy.GetParameterAsText(1)k= arcpy.GetParameter(2)num=arcpy.GetCount_management(fc)n=int(str(num))arcpy.SetProgressor("step", "更新:"+fc+",Field="+fieldname,0,n,1)rows = arcpy.UpdateCursor(fc)

i=k

for row in rows:arcpy.SetProgressorLabel("正在等待....")row.setValue(fieldname,i)arcpy.AddMessage("当前:"+str(i))rows.updateRow(row)arcpy.SetProgressorPosition()i=i+1arcpy.ResetProgressor()del row

作用：更新一个表（含要素）一个字段值为顺序号

参数如下，字段参数获取自图层

运行界面

10.4.3.3 矢量批量裁剪的代码

在：chp10\python\矢量批量裁剪.tbx\矢量批量裁剪，可以直接运行，右键编辑查看代码

代码如下：

# 1.py

# Created on: 星期日一月 10 2018 11:02:13 上午

# (generated by ArcGIS/ModelBuilder)

# Usage: 矢量图批量切割， by 闫磊 4个参数

# 原始数据是图层，可以多选

# 切割工具是是接幅表或者行政

# 字段是输出mdb名称

# 输出路径

# Create Geoprocessing Object

import sys, os, string

import arcpy

from arcpy import env

defaultencoding = 'utf-8'

if sys.getdefaultencoding() != defaultencoding:

reload(sys)

sys.setdefaultencoding(defaultencoding)

arcpy.env.overwriteOutput = True

inworkspace = arcpy.GetParameterAsText(0)

arcpy.AddMessage("输入数据="+inworkspace)

clipshp = arcpy.GetParameterAsText(1)

arcpy.AddMessage("裁剪=clipshp"+clipshp)

fieldname= arcpy.GetParameterAsText(2)

arcpy.AddMessage("字段=fieldname"+fieldname)

outworkspace = arcpy.GetParameterAsText(3)

arcpy.AddMessage("输出="+outworkspace)

mdbbool = arcpy.GetParameterAsText(4)

arcpy.AddMessage("是否mdb="+mdbbool)

desc = arcpy.Describe(clipshp)

filepath=desc.CatalogPath

p=filepath.find(".mdb")

ftype="String"

for field in desc.fields:

if field.Name ==fieldname:

ftype=field.Type

break

arcpy.AddMessage(u"默认地理数据库:"+arcpy.env.scratchWorkspace)

jfb_Select=arcpy.env.scratchWorkspace+"\yl999"#不能c:\要c:\\或者 c:/

rows = arcpy.SearchCursor(clipshp)

#arcpy.AddMessage(u"5=执行到这里")

row = rows.next()

#arcpy.AddMessage(u"6=执行到这里")

while row:

#arcpy.AddMessage(u"7=执行到这里")

fieldvalue =""+ str(row.getValue(fieldname))

#arcpy.AddMessage(u"值fieldvalue="+fieldvalue)

if p>0: #mdb

Expression="["+fieldname +"]="

else:

Expression="\""+fieldname +"\"="

#arcpy.AddMessage(u"表达式Expression1="+Expression)

if ftype=="String":

Expression=Expression+"'"+fieldvalue+"'"

else:

Expression=Expression+fieldvalue

#arcpy.AddMessage(u"Expression2="+Expression)

arcpy.Select_analysis(clipshp, jfb_Select,Expression)

#arcpy.AddMessage(u"6=clipshp"+clipshp)

out_mdb=""

#arcpy.AddMessage("======================================================out_mdb"+out_mdb)

if mdbbool=="true":

out_mdb=outworkspace + "\\"+fieldvalue+".mdb" #os.path.basename(dataset)

else:

out_mdb=outworkspace + "\\"+fieldvalue+".gdb"

arcpy.AddMessage(u"out_mdb"+out_mdb)

if not arcpy.Exists(out_mdb):

if mdbbool=="true":

arcpy.CreatePersonalGDB_management(os.path.dirname(out_mdb),os.path.basename(out_mdb))

else:

arcpy.CreateFileGDB_management(os.path.dirname(out_mdb),os.path.basename(out_mdb))

mydatasets= string.split(inworkspace,";")

for dataset in mydatasets:

try:

mylayer=os.path.basename(dataset)

arcpy.AddMessage(u"clip:"+dataset+" to "+out_mdb+"\\"+ mylayer)

mylayer=mylayer.replace("(","")

mylayer=mylayer.replace(")","")

arcpy.Clip_analysis(dataset, jfb_Select,out_mdb+"\\"+ mylayer, "")

except Exception, ErrorDesc:

#If an error set output boolean parameter "Error" to True.

arcpy.AddError(str(ErrorDesc))

row = rows.next()

if arcpy.Exists(jfb_Select):

arcpy.Delete_management(jfb_Select)

作用：使用一个矢量图层批量裁剪多个矢量数据，字段值是裁剪后数据库名称

图10-53 批量裁剪参数设置界面

界面如下：

图10-54 批量裁剪运行界面

是否mdb，勾上是mdb,不勾是gdb

10.4.3.4 矢量批量合并

在：chp10\python\批量合并.tbx\数据批量合并下，可以直接运行，右键编辑可以查看代码

代码如下：

import sys

import arcpy

import string

try:

workspace =arcpy.GetParameterAsText(0) #'C:\Users\Administrator\Desktop\\cc'

outdb =arcpy.GetParameterAsText(1) #'C:\Users\Administrator\Desktop\\lutian.mdb'

arcpy.env.workspace = workspace

arcpy.AddMessage("outdb:"+outdb)

files = arcpy.ListWorkspaces("","")

for File in files:

arcpy.AddMessage("File:"+File)

arcpy.env.workspace = outdb

fcs = arcpy.ListFeatureClasses()

for fc in fcs:

arcpy.AddMessage("fc:"+fc)

if arcpy.Exists(File + "\\" + fc):

arcpy.Append_management([ File + "\\" + fc], outdb + "\\" + fc,"NO_TEST","","")

else:

arcpy.AddMessage("not exists:"+File + "\\" + fc)

fcs = arcpy.ListTables()

for fc in fcs:

arcpy.AddMessage("fc:"+fc)

if arcpy.Exists(File + "\\" + fc):

arcpy.Append_management([File + "\\" + fc], outdb + "\\" + fc,"NO_TEST","","")

else:

arcpy.AddMessage("not exists:"+File + "\\" + fc)

dss = arcpy.ListDatasets()

for ds in dss:

arcpy.AddMessage("ds:"+ds)

arcpy.env.workspace = outdb+"\\"+ds

fcs1 = arcpy.ListFeatureClasses()

for fc1 in fcs1:

arcpy.AddMessage("fc1:"+fc1)

if arcpy.Exists(File + "\\" + ds + "\\" + fc1):

arcpy.Append_management([File + "\\" + ds + "\\" + fc1], outdb + "\\" + ds + "\\" + fc1,"NO_TEST","","")

else:

arcpy.AddMessage("not exists:"+File + "\\" + ds + "\\" + fc1)

except arcpy.ExecuteError:

arcpy.AddWarning(arcpy.GetMessages())

参数如下，工作空间可以是地理数据库（文件地理数据或者个人地理数据），也是可以一个文件夹，这里建议是数据库

图10-55 批量合并参数设置界面

运行界面如下：

图10-56 批量合并运行界面

10.4.3.5 影像批量裁剪

在：chp10\python\影像切割.tbx\影像切割，可以直接运行，右键编辑可以查看程序源代码

源代码如下：

import sys, os, string,types

import arcpy

from arcpy import env

arcpy.env.overwriteOutput = True

oldraster = arcpy.GetParameterAsText(0)

arcpy.AddMessage("1oldraster="+oldraster)

clipshp = arcpy.GetParameterAsText(1)

arcpy.AddMessage("2clipshp="+clipshp)

fieldname= arcpy.GetParameterAsText(2)

arcpy.AddMessage("3fieldname="+fieldname)

outworkspace= arcpy.GetParameterAsText(3)

arcpy.AddMessage("4="+outworkspace)

arcpy.CheckOutExtension("spatial")

rows = arcpy.SearchCursor(clipshp)

jfb_Select=outworkspace+"/temp.shp" #不能c:\要c:\\或者c:/

for row in rows:

try:

b=1

value=row.getValue(fieldname)

#gp.AddMessage("value="+value)

if (type(value) is types.IntType):

fieldvalue = str(value)

b=2

elif (type(value) is types.StringType): #是否string类型

fieldvalue = value

else:

fieldvalue = str(value)

arcpy.AddMessage("fieldvalue="+fieldvalue)

if b==2:

Expression="\""+fieldname +"\" ="+fieldvalue+""

else:

Expression="\""+fieldname +"\" ='"+fieldvalue+"'"

arcpy.AddMessage("Expression="+Expression+",jfb_Select="+jfb_Select+",clipshp="+clipshp)

arcpy.Select_analysis(clipshp, jfb_Select, Expression)

out_raster =outworkspace+"/"+fieldvalue+".tif"

arcpy.gp.ExtractByMask_sa(oldraster, jfb_Select, out_raster)

except Exception, ErrorDesc:

#If an error set output boolean parameter "Error" to True.

arcpy.AddError(str(ErrorDesc))

if arcpy.Exists(jfb_Select):

arcpy.Delete_management(jfb_Select)

作用：使用一个矢量数据，批量裁剪一个影像，矢量字段值是裁剪后影像的数据名，格式为tif

参数：

图10-57 影像批量裁剪参数设置界面

运行界面如下：

图10-58 影像批量裁剪运行界面

10.4.3.6 批量压缩数据库

在chp10\python\压缩数据.tbx\批量压缩数据库，可以直接运行，右键编辑可以查看程序源代码

源代码如下：

import arcpy

import sys

import os

from os.path import join, getsize

reload(sys)

sys.setdefaultencoding('utf8')

inpath = arcpy.GetParameterAsText(0)

for root, dirs, files in os.walk(inpath):

arcpy.env.workspace = root

workspaces = arcpy.ListWorkspaces("*", "Access")

for workspace in workspaces:

try:

arcpy.AddMessage("正在压缩mdb: " + str(workspace))

arcpy.Compact_management(workspace)

except Exception, ErrorDesc:

arcpy.AddError(str(ErrorDesc))

workspaces = arcpy.ListWorkspaces("*", "FileGDB")

for workspace in workspaces:

try:

arcpy.AddMessage("正在压缩gdb: " + str(workspace))

arcpy.Compact_management(workspace)

except Exception, ErrorDesc:

arcpy.AddError(str(ErrorDesc))

作用；批量压缩一个文件夹（含子文件夹）下所有的mdb和gdb

参数：

图10-59 批量压缩数据库参数设置界面

运行界面：

图10-60 批量压缩数据库运行界面

10.4.3.7 mxd文档批量导出图片

在：chp10\python\批量mxd转图片.tbx\mxd转图片，可直接运行，右键编辑可以查看程序源代码

源代码如下：

import arcpy

import os

import sys

indata = arcpy.GetParameterAsText(0)

outdata = arcpy.GetParameterAsText(1)

EXT = arcpy.GetParameterAsText(2)

dpi = arcpy.GetParameter(3)

EXT=EXT.upper()

for root, dirs, files in os.walk(indata): #会遍历所有文件包括子文件里的数据

arcpy.AddMessage(u"目录:" + root)

rootList=root.split('.')

n=len(rootList)

if (n<2): #不考虑进入 gdb 目录下循环

for name in files:

#arcpy.AddMessage(u"name:" + name)

mystr=name.upper()

if (mystr.endswith('.MXD')): #找到扩展名.mxd 的数据

try:

fs = os.path.join(root,name)

mxd = arcpy.mapping.MapDocument(fs)

arcpy.AddMessage(u" 保存了 :" + outdata+"/" +

name.split('.')[0]+"."+EXT)

if EXT=="JPG":

arcpy.mapping.ExportToJPEG(mxd,outdata +"/"+name.split('.')[0],resolution=dpi)

elif EXT=="TIF":

arcpy.mapping.ExportToTIFF(mxd,outdata +"/"+name.split('.')[0],resolution=dpi)

elif EXT=="PDF":

arcpy.mapping.ExportToPDF(mxd,outdata +"/"+name.split('.')[0],resolution=dpi)

elif EXT=="PNG":

arcpy.mapping.ExportToPNG(mxd,outdata +"/"+name.split('.')[0],resolution=dpi)

elif EXT=="EPS":

arcpy.mapping.ExportToEPS(mxd,outdata +"/"+name.split('.')[0],resolution=dpi)

elif EXT=="EMF":

arcpy.mapping.ExportToEMF(mxd,outdata +"/"+name.split('.')[0],resolution=dpi)

del mxd

except Exception, ErrorDesc:

arcpy.AddError(str(ErrorDesc))

作用；批量一个文件夹下所有的mxd导出图片，可以是PDF，JPG，PNG，TIF，EPS，EMF

参数：

图10-61 mxd转图片参数设置

扩展名是值列表

图10-62 值列表输入界面

运行界面如下：

图10-63 mxd转图片运行界面

文章来源：GISAI

如有侵权请联系删除 cehui8@qq.com

行业精选

网课测绘师飞手证 cors账号

网课：

测绘师：

飞手证：

航测：

硬件：

技师：

实战班：

其他：

自助查询

页面底部→关注中测网

对话框发送【关键词】获取精选内容

【矢量地图】各省 /市 /县 shp矢量数据

【ZBBZ】 CAD /CASS坐标、标注插件

【规范】测量规范、常用Excel表格

【CASS】快捷键、高手技巧、操作视频

【RTK】南方/华测各类RTK操作视频

【GIS】 ArcGIS /MapGIS技巧、快捷键

【航测】无人机作业流程、规范、案例

【资质】测绘资质查询 /申办、甲级名单

【定额】测绘生产成本定额、勘察收费

【招标】测绘地信招标信息 /公告

【地图】奥维地图、谷歌地图下载/安装

【图片】逆天测绘图片、测绘日常实拍图

特别声明：以上内容(如有图片或视频亦包括在内)为自媒体平台“网易号”用户上传并发布，本平台仅提供信息存储服务。

Notice: The content above (including the pictures and videos if any) is uploaded and posted by a user of NetEase Hao, which is a social media platform and only provides information storage services.