[ARCGIS]使用ARCPY+PYTHON+FME批量出图

应用场景：

最近在做某区退耕还林验收项目，要求按每个村出比例尺1：10000的JPG格式的图片，由于村比较多（接近400个，很多村在1：10000比例尺下无法显示完整，需要出多张图片），所以考虑批量完成该工作。

思路：

1.设定好MXD模板布局视图，调整好所有图例，将标题设置为变量，可以根据某个村的名字自动变化；

2.一个村无法只导出一张图片的情况，需要生成辅助图层分幅，分幅使用FME，生成图形的外接矩形，计算矩形长宽以及长宽/出图的幅宽，向上取整为倍数，使用切片工具按照两个倍数值进行切片，然后使用裁剪工具对原图形进行裁剪，只取出inside部分；

分幅之后大约有1000个图幅，后续再用造林小班和分幅图叠加，筛选出有叠加部分并排序，生成唯一字段。

3.在进行2步骤的时候，因为矩形之间有交叉部分必须要单独对每个图形进行处理，因此可以使用FME WORKSPACERUNER转换器循环执行，写入GDB文件中；

4.将GDB所有属性值导出为EXCEL文件待用；

5.根据4生成的EXCEL文件使用PYTHON生成若干个（数量和GDB文件中的矢量图形数量相等，名称和GDB文件中唯一值对应）MXD文档，待用；

6.使用ARCPY选中要素，缩放至选中要素，清除所选要素（注意这里有坑），更改标题变量值；

7.使用ARCPY导出JPEG图片。

具体步骤：

1 MXD模板制作

MXD模板制作有两个要点。

第一个是要对需要发生变化的文本、图例赋上唯一名称。右键选中文本或图例，点击“属性”→“大小和位置”→“元素名称”，自己按照需求添加。

第二个是要固定比例尺，因为不固定比例尺，在后续使用ARCPY缩放至选中要素功能时会自动缩放到不同的比例尺（每个村大小不同导致，也是因为如此需要分幅）。右键点击“图层”→“数据库”，设定“范围”为“固定比例”，比例调整为所需大小，此处我设置为1：10000.

2 行政村分幅

下图为FME模板

3 使用PYTHON生成MXD文档

考虑到需要对过程文件的记录，选择先生成MXD文档的方式（也可以在调整图形时直接生成），因为习惯用PYTHON3，所以这一步骤还是用的PYTHON3写的，比较简单，以下是代码。

# -*- coding: utf-8 -*-

import shutil,sys
from openpyxl import load_workbook
old_file_path=r"F:\arcgis\2018.mxd"

if __name__=='__main__':
    workbook=load_workbook(r'F:\arcgis\TEST.xlsx')
    sheet=workbook.get_sheet_by_name("CJ")
    rows=sheet.rows
    cols=sheet.columns
    max_rows=sheet.max_row
    max_columns=sheet.max_column
    Lfile=[]
    for i in sheet["H"]:
        if i.value!='IID':
            Lfile.append(i.value)
    print(Lfile)
    for i in Lfile:
        new_file_path=r"F:\arcgis\{0}.mxd".format(i)
        shutil.copyfile(old_file_path,new_file_path)

4 使用ARCPY生成重新生成MXD文档

必须使用arcgis自带的python，我使用的是arcgis10.3，自带python2.7,以下是代码。

# -*- coding: utf-8 -*-
import arcpy,os,sys

reload(sys)
sys.setdefaultencoding( "utf-8" )

def resetmxd(iid):
    filename="{0}.mxd".format(iid)
    mxd = arcpy.mapping.MapDocument(filename)

    lyr = arcpy.mapping.ListLayers(mxd)[0]  # 需要改动的图层，这里是选择要素的图层
    lyr1 = arcpy.mapping.ListLayers(mxd)[4]  # 需要缩放高亮显示的图层，可选

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

    rows = arcpy.SearchCursor(lyr)   #行数
    for row in rows:
        if row.getValue("IID")==iid:
            title=u'2018年度{0}退耕还林'.format(row.getValue("title"))

    query = "IID='{0}'".format(iid)      #检索条件
    print query
    arcpy.SelectLayerByAttribute_management(lyr, "NEW_SELECTION", query)  # 选中一个新的要素

    df.zoomToSelectedFeatures()  # 缩放至选中要素
    el = arcpy.mapping.ListLayoutElements(mxd, "TEXT_ELEMENT")[-1]
    el.text = title  # 更新的是文本，不是教程上的元素名称
    arcpy.SelectLayerByAttribute_management(lyr, "CLEAR_SELECTION")  # 清除所选要素
    #
    # arcpy.SelectLayerByAttribute_management(lyr1, "NEW_SELECTION", "FID=8")
    # # 选择是否高亮显示

    arcpy.mapping.MapDocument.save(mxd)  # 重新保存mxd

    arcpy.RefreshActiveView()     #刷新视图
    arcpy.RefreshTOC()
    del mxd

if __name__=='__main__':
    L=['5002301002011', '5002301002001', '5002301002041', '5002301002042', '5002301002031', '5002301002032', '5002301002051', '5002301002052', '5002301002053']
    for i in L:
        resetmxd(i)

因为还在测试阶段，所以上述都是测试数据，比较乱，到时候拿到正式的小班数据后会对文档进行整理。

这里有个两个选择，第一高亮显示行政区（不能高亮显示分幅的行政区，看起来会很奇怪），第二不高亮显示行政区，这个按照自己的需求决定，所以这里有个坑：清除所选要素（不高亮显示）的方法还是SelectLayerByAttribute_management（）不过第二个参数值需要改为“CLEAR_SELECTION”。

5 使用ARCPY出图

# -*- coding: utf-8 -*-
import arcpy,os,sys

reload(sys)
sys.setdefaultencoding( "utf-8" )

if __name__=='__main__':
    L=['5002301002011', '5002301002001', '5002301002041', '5002301002042', '5002301002031', '5002301002032', '5002301002051', '5002301002052', '5002301002053']
    for iid in L:
        filename = "{0}.mxd".format(iid)
        mxd = arcpy.mapping.MapDocument(filename)
        new_path='{0}.jpg'.format(iid)
        arcpy.mapping.ExportToJPEG(mxd,new_path,resolution=300)

 

ExportToJPEG方法参数参考文档：http://desktop.arcgis.com/zh-cn/arcmap/10.3/analyze/arcpy-mapping/exporttojpeg.htm

先记录大概的思路和过程，等项目完成再完善此文档，欢迎友好讨论！