python使用arcpy.mapping模块批量出图_Python

python使用arcpy.mapping模块批量出图

2020-09-22 09:49乐水悠悠 Python

出图是项目里常见的任务，有的项目甚至会要上百张图片，所以批量出土工具很有必要，这篇文章主要介绍了python使用arcpy.mapping模块批量出图，有兴趣的可以了解一下。

出图是项目里常见的任务，有的项目甚至会要上百张图片，所以批量出土工具很有必要。arcpy.mapping就是ArcGIS里的出图模块，能快速完成一个出图工具。

arcpy.mapping模块里常用的类有MapDocument、DataFrame、Layer、DataDrivenPages和TextElement。

MapDocument类是地图文档（.mxd文件）对应的类。初始化参数是一个字符串，一般是.mxd文件的路径：

1	`mxd=arcpy.mapping.MapDocument(r"F:\GeoData\ChinaArea\ChinaVector.mxd")`

DataFrame类用于操作地图内的Data Frame（即下图的Layers），能够控制地图的范围、比例尺等。用arcpy.mapping.ListDataFrames(map_document, {wildcard})函数获取。

1	`df=` `arcpy.mapping.ListDataFrames(mxd)[0]`

Layer类用于操作具体的图层。能够控制图斑的样式、可见性等。可以用.lyr文件的路径初始化，也可以通过arcpy.mapping.ListLayers(map_document_or_layer, {wildcard}, {data_frame})函数获取。

									lyr1=arcpy.mapping.Layer(r" F:\GeoData\ChinaArea\Province.lyr")

									df.addLayer(lyr1)

									lyr2=arcpy.mapping.ListLayer(mxd,"",df)[0]

DataDrivenPages类需要配合ArcMap中的Data Driven Pages工具使用。用于一个矢量文件内的全部或部分图斑每个出一张图的情况。

TextElement类用于操作地图上的文字，比如图名、页数。通过arcpy.mapping.ListLayoutElements (map_document, {element_type}, {wildcard})函数获取。

1	`txtElm=arcpy.mapping.ListLayoutElements(mxd,"TEXT_ELEMENT")[0]`

常见的出图模式有两种：一个矢量文件里每个图斑出一张图，一个文件夹下每个矢量文件出一张图。

每个图斑出一张图：

这种情况有Data Driven Pages工具配合最好。打开ArcMap的Customize->Toolbars->Data Driven Pages，设置好图层、名称字段、排序字段、显示范围和比例尺，保存地图。

									# coding:utf-8

									import arcpy

									mxd=arcpy.mapping.MapDocument(r"F:\GeoData\ChinaArea\ChinaVector.mxd")

									for pageNum in range(1,mxd.dataDrivenPages.pageCount):

									     mxd.dataDrivenPages.currentPageID=pageNum

									     mapName=mxd.dataDrivenPages.pageRow.getValue(mxd.dataDrivenPages.pageNameField.name)

									     print mapName

									     arcpy.mapping.ExportToPNG(mxd,r"F:\GeoData\ChinaArea\Province\\"+mapName+".png")

									print 'ok'

一个文件夹下的每个矢量文件出一张图：

									# coding:utf-8

									import arcpy

									import os

									def GetShpfiles(shpdir):

									     shpfiles=[]

									     allfiles=os.listdir(shpdir)

									     for file in allfiles:

									          if os.path.isfile(file):

									              if file.endswith('.shp'):

									                   shpfiles.append(file)

									          else:

									              shpfiles.extend(GetShpfiles(file))

									     return shpfiles

									allshps=GetShpfiles(r"F:\GeoData\ChinaArea\Province")

									mxd=arcpy.mapping.MapDocument(r"F:\GeoData\ChinaArea\ChinaVector.mxd")

									lyr=arcpy.mapping.ListLayer(mxd)[0]

									for shp in allshps:

									     paths=os.path.split(shp)

									     print paths[1]

									     lyr.replaceDataSource(paths[0],"SHAPEFILE_WORKSPACE",paths[1])

									     arcpy.mapping.ExportToPNG(mxd,r"F:\GeoData\ChinaArea\Province\\"+paths[1]+".png")

									print 'ok'