使用ArcPy制图模块实现网络电子地图打印.docx

电子地图的打印输出就是将电子地图转化为纸质地图。与纸质地图相比，电子地图最显著的特征是数据的存贮与显示相分离。由此产生电子地图的一系列新特点，如动态性、交互性、多媒体集成性、多尺度性等。但电子地图具有很强的设备依赖性，且不如纸质地图直观1。所以有些情况下将电子地图打印输出为纸质地图，既利用了电子地图制作简单、省时省力的优势，又具有了纸张使用方便的特点，逐渐成为纸质地图生产的一种新方法2 打印和导出电子地图是我们使用ArcMap常做的事，但在网络电子地图上要实现相同功能还是挺坎坷的，大家做法不一：（1）浏览器插件，常见的就是IE的ActiveX，优点是在浏览器端可以进行简单的制图调整，但缺点也要命，只支持IE、IE各版本难兼容、安全性限制所以此方法基本过时；（2）Web Service，发布Web服务调用ArcEngine组件，优点是采用服务的方式，通用性好，缺点就是需要部署ArcEngine，且专门开发。还有没有更好的方法呢？基于第2种思路，我们可以采用Server本身具备的GeoProcessing功能服务来实现，其中关键技术点就是使用ArcPy 的制图模块。系统客户端界面 ArcPy？没听过ArcPy说明你有点小out了作为ArcGIS 10新的Python库，ArcPy除了继承原有功能之外还增加了制图模块（ArcPy.Mapping）。ArcPy制图模块提供了对AO的高级封装，通过调用一系列方法实现修改地图范围、比例尺、符号化、图幅整饰等要素，实现自动化的地图打印或导出。接下来我们一步步了解下实现的具体过程：1、 在Python脚本中添加ArcPy引用：import arcpy2、 获取GP服务的输入参数（地图标题、地图范围、比例尺、模板样式、坐标参考等）：# 地图范围参数xMin = arcpy.GetParameterAsText(0)yMin = arcpy.GetParameterAsText(1)xMax = arcpy.GetParameterAsText(2)yMax = arcpy.GetParameterAsText(3)#坐标参考srValIn = arcpy.GetParameterAsText(4)srIn = arcpy.CreateObject(SpatialReference)srIn.loadFromString(srValIn)#地图比例尺mapScale = arcpy.GetParameterAsText(5)#地图模板mxd文档Layout = arcpy.GetParameterAsText(6)# 是否输出图层属性信息attributesFlag = arcpy.GetParameterAsText(7)if attributesFlag.lower() = true: layers_string = LAYERS_AND_ATTRIBUTESelse: layers_string = LAYERS_ONLY# 地图标题title = arcpy.GetParameterAsText(8)3、 打开选择的地图样式（Mxd模板），根据输入的参数修改其中的地图范围等信息：# 打开输入的地图模板mxd文档mxdName = os.path.join(LayoutRootPath, Layout) + .mxdmapDoc = arcpy.mapping.MapDocument(mxdName)# 指定地图范围和坐标参考dataFrames = arcpy.mapping.ListDataFrames(mapDoc)dataFrame. = dataFrames0dataFrameSR = dataFrame.spatialReference# 确保坐标参考一致pntInLL = arcpy.Point()pntInLL.X = xMinpntInLL.Y = yMinpntInUR = arcpy.Point()pntInUR.X = xMaxpntInUR.Y = yMaxpntGeomInLL = arcpy.PointGeometry(pntInLL,srIn)pntGeomInUR = arcpy.PointGeometry(pntInUR,srIn)pntGeomList = pntGeomList.append(pntGeomInLL)pntGeomList.append(pntGeomInUR)arcpy.CopyFeatures_management(pntGeomList, in_memory/extentPnts)rows = arcpy.SearchCursor(in_memory/extentPnts, , dataFrameSR)rows.reset()row = rows.next()newPntLL = row.shape.getPart(0)row = rows.next()newPntUR = row.shape.getPart(0)myExtent = dataFrame.extentarcpy.AddMessage(Old extent: %s %s %s %s % (myExtent.XMin, myExtent.YMin, myExtent.XMax, myExtent.YMax)myExtent.XMin = newPntLL.XmyExtent.YMin = newPntLL.YmyExtent.XMax = newPntUR.XmyExtent.YMax = newPntUR.YdataFrame.extent = myExtent# 打印制图内容nExt = dataFrame.extentnXMin = nExt.XMinnYMin = nExt.YMinnXMax = nExt.XMaxnYMax = nExt.YMaxnStr = %s %s %s %s % (nXMin, nYMin, nXMax, nYMax)myExtent = dataFrame.extentarcpy.AddMessage(New extent: %s %s %s %s % (myExtent.XMin, myExtent.YMin, myExtent.XMax, myExtent.YMax)#设置比例尺if mapScale != 0: dataFrame.scale = mapScale# 设置地图标题if title != : for elm in arcpy.mapping.ListLayoutElements(mapDoc, TEXT_ELEMENT): if elm.text = Map Title: elm.text = title4、 用导出方法，输出PDF或MXD文档到临时目录，代码如下：arcpy.AddMessage( Saving as PDF)# 输出到临时目录scratch = arcpy.env.scratchWorkspace#arcpy.AddMessage(scratch)#if not scratch: scratch = C:/TempOutputFile = os.path.join(scratch,Map.pdf)arcpy.mapping.ExportToPDF(mapDoc,OutputFile, layers_attributes=layers_string)5、 把输出的PDF或MXD文档的URL作为GP服务的输出参数，传给客户端，进行下载，代码如下：# 得到执