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


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

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

             MapDocument类是地图文档(.mxd文件)对应的类。初始化参数是一个字符串,一般是.mxd文件的路径:

             mxd=arcpy.mapping.MapDocument(r"F:GeoDataChinaAreaChinaVector.mxd")

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

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

            

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

             lyr1=arcpy.mapping.Layer(r" F:GeoDataChinaAreaProvince.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})函数获取。

             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:GeoDataChinaAreaChinaVector.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:GeoDataChinaAreaProvince\"+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:GeoDataChinaAreaProvince")
    
    mxd=arcpy.mapping.MapDocument(r"F:GeoDataChinaAreaChinaVector.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:GeoDataChinaAreaProvince\"+paths[1]+".png")
    
    print 'ok'

    更多功能见ArcMap帮助文档Geoprocessing->ArcPy->Mapping Module。

  • 相关阅读:
    PowerTalk第一个版本儿控件
    PowerTalk第二个版本,支持(Msn回复信息)
    自然语言处理著作或期刊名称
    自然语言处理(NLP)网上资源整理 (转)
    TFIDF
    计算机科学及工程
    自然语言处理与计算语言学书籍汇总之一:国外书籍
    UVa 10696 f91
    缓存
    操作必须使用一个可更新的查询(转)
  • 原文地址:https://www.cnblogs.com/sunnyeveryday/p/6506347.html
Copyright © 2020-2023  润新知