利用Arcpy发布地图服务，制作切片

以你之姓@ 2022-07-24 09:10 474阅读 0赞

因为工作需求，我需要把管理员上传的图片进行切片，供用户下载切片离线浏览。后台切片功能主要分为以下几个步骤：

1. 制作地图文档（*.mxd）;
    2. 发布地图文档;
    3. 制作服务器缓存;
    4. 生成切片；
    5. 打包成zip

ArcGIS版本是10.2，下面这些代码应该在10.0上无法运行。生成mxd文档时不支持栅格数据，我的图片大都是jpeg格式的，所以最终我更换了ArcGIS版本，妥协了…

## Arcpy制作地图文档 ##

我没有找打直接生成地图文档的示例代码，因此我使用了SaveAs的方法，就是利用已有的空白文档（我称之为“模板”），加载上jpg图片后再“另存为”。这样的一个好处就是我可以直接在模板中定义地理坐标系，后面再发布为地图服务时会用到的，否则报错。

# 创建map document
        def CreateMxd(imagepath,mxdpath):
            dirname=os.path.dirname(imagepath)
            imagename=os.path.basename(imagepath)
            dotindex=imagename.index('.')
            name=imagename[0:dotindex]
            new_mxd=os.path.abspath(dirname+"/"+name+".mxd")
            rasterLayer="raster"
    
            temp_mxd = arcpy.mapping.MapDocument(mxdpath)
            df=arcpy.mapping.ListDataFrames(temp_mxd,"Layers")[0]
            arcpy.MakeRasterLayer_management(imagepath,rasterLayer,"","","")
            addLayer=arcpy.mapping.Layer(rasterLayer)
            arcpy.mapping.AddLayer(df,addLayer,"TOP")
            temp_mxd.saveACopy(new_mxd)
            del temp_mxd
            return new_mxd

这个函数的参数是“图片的位置”和“模板的位置”，然后将图片的名称提取出来作为新mxd文件的名称。

需要注意的是，直接添加jpg格式的图层在ArcGIS10.2中是允许的，但生成的mxd文档里面的图层样式是采用“拉伸（Stretched）”方式展示的，图片的样式直接成为灰色调了，而我需要的是“RGB合成（RGB Composite）”。因此我使用了这种“先创建临时栅格图层，再添加到mxd”的方式。  
　　![使用RGB合成][RGB]

## 发布地图文档 ##

我们都知道，在ArcGIS10.1之后都是先将mxd文档转为sd文档，再进行发布的。这样做也是有好处的（听说是），不再赘述。官网文档给出两种方式发布地图，一是需要ArcGIS Server服务器名、用户名和密码方式验证登录，二是采用arcgis server连接文件ags的方式直接验证。我采用了后者，比如我的连接文件位置在“C:\\Users\\Administrator\\AppData\\Roaming\\ESRI\\Desktop10.2\\ArcCatalog\\arcgis on WIN-20150327EEH\_6080 (admin).ags”。直接引用这个地址作为参数传入下面这个函数就行。

# 发布服务
        def PublishService(mxdpath,agspath):
            new_mxd = arcpy.mapping.MapDocument(mxdpath)
            dirname=os.path.dirname(mxdpath)
            mxdname=os.path.basename(mxdpath)
            dotindex=mxdname.index('.')
            servicename=mxdname[0:dotindex]
    
            sddraft = os.path.abspath(servicename + '.sddraft')
            sd = os.path.abspath(servicename + '.sd')
            if os.path.exists(sd):
            os.remove(sd)
            #创建服务定义草稿draft
            arcpy.CreateImageSDDraft(new_mxd, sddraft, servicename, 'ARCGIS_SERVER', agspath,False,None, "Ortho Images", "ortho images,image service")
    
            #分析草稿draft
            analysis = arcpy.mapping.AnalyzeForSD(sddraft)
            #打印分析结果
            print "The following information was returned during analysis of the MXD:"
            for key in ('messages', 'warnings', 'errors'):
            print '----' + key.upper() + '---'
            vars = analysis[key]
            for ((message, code), layerlist) in vars.iteritems():
            print '', message, ' (CODE %i)' % code
            print ' applies to:',
            for layer in layerlist:
            print layer.name,
            print
    
            # 发送草稿至服务器
            if analysis['errors'] == {}:
            # Execute StageService. This creates the service definition.
            arcpy.StageService_server(sddraft, sd)
            # Execute UploadServiceDefinition. This uploads the service definition and publishes the service.
            arcpy.UploadServiceDefinition_server(sd, agspath)
            print "Service successfully published"
            else: 
            print "Service could not be published because errors were found during analysis."
    
            print arcpy.GetMessages()
            return agspath.replace(".ags","/"+servicename+".MapServer")

同样两个参数，这个函数需要mxd文档路径和ags文件路径作为参数。中间需要先生成服务定义文件（*.sd）和服务定义草稿文件（*.sddraft）,然后对sddraft文件进行分析，若中间没有出现错误，现将其发送至服务器。然后再将sd文件，通过ags文件发送至服务器，更新服务。

这段代码最后我返回了刚刚发布的服务地址，后面将会用到。

## 制作服务器缓存 ##

注意哟，这里是定义服务器缓存，不是生成切片地方。这里是定义，下一步才是生成。就好比“打开存有《名侦探柯南》的云盘文件夹”，下一步才是“点击播放按钮”。

#制作地图服务器缓存
        def CreateCache(inputService,cachepath):
            # List of input variables for map service properties
            tilingSchemeType = "NEW"
            scalesType = "CUSTOM"
            numOfScales = "5"#此参数无效
            scales = [500000000,250000000,125000000,64000000,32000000]
            dotsPerInch = "96"
            tileOrigin = "0 0"
            tileSize = "256 x 256"
            cacheTileFormat = "PNG8"
            tileCompressionQuality = ""
            storageFormat = "COMPACT"
            predefinedTilingScheme = ""
    
            try:
                starttime = time.clock()
                result = arcpy.CreateMapServerCache_server(inputService,cachepath,tilingSchemeType, scalesType, numOfScales, dotsPerInch,tileSize, predefinedTilingScheme,tileOrigin, scales,cacheTileFormat,tileCompressionQuality,storageFormat)
    
                # print messages to a file
                while result.status < 4:
                    time.sleep(0.2)
                resultValue = result.getMessages()
                print "completed " + str(resultValue)
    
            except Exception, e:
                # If an error occurred, print line number and error message
                tb = sys.exc_info()[2]
                print "Failed at step 1 \n" "Line %i" % tb.tb_lineno
                print e.message
    
            print "Executed creation of Map server Cache schema "

这里主要用了`CreateMapServerCache_server`这个函数，具体函数的意义这里就不再多说了，完全可以自己去查看ArcGIS文档。我是使用了自定义的方式切片，可以采用STANDARD方式切片，然后设置numOfScales的值就可以了。

## 生成切片 ##

这里才是生成切片的地方，用到的是Arcpy的`ManageMapServerCacheTiles_server`函数，同样具体的参数还是需要自己去查官方文档。帮助文档上面说参数`waitForJobCompletion`的值为“WAIT”时，可以从其他地方查到具体的切片进度，这里我还没去做，看以后的需求吧。

#生成瓦片
        def CreateTiles(inputService):
            scales = ""
            numOfCachingServiceInstances = 2
            updateMode = "RECREATE_ALL_TILES"
            areaOfInterest = ""
            waitForJobCompletion = "WAIT"
            updateExtents = ""
    
            try:
                result = arcpy.ManageMapServerCacheTiles_server(inputService, scales,updateMode,numOfCachingServiceInstances,areaOfInterest, updateExtents,waitForJobCompletion)
                #print messages to a file
                while result.status < 4:
                    time.sleep(0.2)
                resultValue = result.getMessages()
                print "completed " + str(resultValue)
                print "Created cache tiles for given schema successfully"
            except Exception, e:
                # If an error occurred, print line number and error message
                tb = sys.exc_info()[2]
                print "Failed at step 1 \n" "Line %i" % tb.tb_lineno
                print e.message
    
            print "Created Map server Cache Tiles "

这个函数需要上一步返回的参数，就是那个以“*\**.MapServer”结尾的全路径。

## 压缩切片文件 ##

将切片打包主要是为了便于用户下载，在线预览的话就不需要这一步了。下面的代码是参考了[linda1000][]的博文，有兴趣的可以自己前去看看。

这里的两个参数分别是要压缩的切片路径，然后是压缩后的输出路径（注意这里要以\*.zip结尾，包括压缩文件名）.

# MakeZipFile
        def MakeZipFile(filepath,zippath):
            filelist = []
            #Check input ...
            fulldirname = os.path.abspath(filepath)
            fullzipfilename = os.path.abspath(zippath)
            print "Start to zip %s to %s ..." % (fulldirname, fullzipfilename)
            if not os.path.exists(fulldirname):
                print "Dir/File %s is not exist" % fulldirname
                return
            if os.path.isdir(fullzipfilename):
                tmpbasename = os.path.basename(filepath)
                fullzipfilename = os.path.normpath(os.path.join(fullzipfilename, tmpbasename))
    
            #Get file(s) to zip ...
            if os.path.isfile(filepath):
                filelist.append(filepath)
                filepath = os.path.dirname(filepath)
            else:
                #get all file in directory
                for root, dirlist, files in os.walk(filepath):
                    for filename in files:
                        filelist.append(os.path.join(root,filename))
    
            #Start to zip file ...
            destZip = zipfile.ZipFile(fullzipfilename, "w")
            for eachfile in filelist:
                destfile = eachfile[len(filepath):]
                print "Zip file %s..." % destfile
                destZip.write(eachfile, destfile)
            destZip.close()
            print "Zip folder succeed!"

## 调用 ##

这段代码我是使用Java runtime直接调用的，可以使用这种方式一套流程走下来。但也有不大合适的地方，比如瓦片被删掉了，我只想生成一下切片等等，这样就不合适了。

#********************************************
        # 主 函 数 *
        #********************************************
        # 程序运行所需参数：
        # 1.图片所在位置（filepath+filename）
        # 2.mxd模板位置（filepath+filename）
        # 3.arcgis server连接文件（ags）位置
        # 4.生成切片缓存位置（out_path）
        # 5.生成压缩文件位置（out_path+zipname）
        if __name__ == '__main__':
            #验证路径是否正确
            imagepath=sys.argv[1]
            mxdpath=sys.argv[2]
            agspath=sys.argv[3]
            cachepath=sys.argv[4]
            zippath=sys.argv[5]
            if not os.path.isfile(imagepath):
                print "图片文件不存在"
                sys.exit()
            if not os.path.isfile(mxdpath):
                print "mxd模板不存在"
                sys.exit()
            if not os.path.isfile(agspath):
                print "arcgis server连接文件不存在"
                sys.exit()
            if not os.path.isdir(cachepath):
                print "缓存目录不存在"
                sys.exit()
            # 创建mxd
            new_mxd=CreateMxd(imagepath,mxdpath)
            print "Finished CreateMxd"
            # 发布服务
            inputService=PublishService(new_mxd,agspath)
            print "Finished PublishService"
            # 制作服务器缓存
            CreateCache(inputService,cachepath)
            print "Finished CreateCache"
            # 生成瓦片
            CreateTiles(inputService)
            print "Finished CreateTiles"
            # 压缩文件
            tcachepath=os.path.abspath(cachepath+"/"+new_mxd[new_mxd.rindex("\\")+1:new_mxd.rindex(".")])
            print tcachepath
            MakeZipFile(tcachepath,zippath)
            print "Finished MakeZipFile"

## 更新 ##

使用过程中，发现无法使得切出更小比例尺的切片，最小的切片也是１:500000000。但当用户上传的图片过大时，切出的图片无法在视图内完整的看到全貌。这显然不合理，最后终于找到了解决方法。草稿文件上传至服务器前对其进行修改，sddraft文件是指上是一个符号xml标准的文件。使用python直接可以进行操作，将里面的minScale改为更大

#修改草稿draft
        # read sddraft xml
        doc = DOM.parse(sddraft)
    
        # turn on caching in the configuration properties
        configProps = doc.getElementsByTagName('ConfigurationProperties')[0]
        propArray = configProps.firstChild
        propSets = propArray.childNodes
        for propSet in propSets:
            keyValues = propSet.childNodes
            for keyValue in keyValues:
                if keyValue.tagName == 'Key':
                    if keyValue.firstChild.data == "minScale":
                        # turn on caching
                        keyValue.nextSibling.firstChild.data = "32000000000"
    
        # output to a new sddraft
        if os.path.exists(sddraft): os.remove(sddraft)
        f = open(sddraft, 'w')
        doc.writexml( f )
        f.close()

[下载源码][Link 1]

[RGB]: /images/20220722/70a6a6103dab4041b75e8f1d3740cd89.png
[linda1000]: http://blog.csdn.net/linda1000/article/details/10432133
[Link 1]: http://download.csdn.net/detail/dahongdahong/9763241