ArcPython地理處理大殺器

ArcPython，地理處理大殺器Esri(中國)信息技術(shù)西安分公司嚴(yán)瑞、范秦光

大家好，我是Esri西安辦嚴(yán)瑞，很榮幸能夠在2016年第一堂公開課和大家見面，在這里給大家拜一個(gè)晚年，祝各位新的一年身體健康、萬事如意！需要跟大家說一聲抱歉，由于我嗓子發(fā)炎，聲音嘶啞，本堂課我將會(huì)用旁邊文字替代音頻的形式來完成，小伙伴們

另：本堂課是以實(shí)戰(zhàn)練習(xí)為主的課，所以更多的內(nèi)容是以實(shí)際操作為主，大家只要看明白我的練習(xí)操作就可以，課后會(huì)將所有源碼和數(shù)據(jù)提供給大家，是全部哦！希望各位小伙伴們能夠理解！順便可以體驗(yàn)一下Esri中國史上第一堂”無聲公開課“哦！1ArcPython在ArcGIS中的用法2與ArcGISDesktop實(shí)戰(zhàn)篇3與ArcGISServer實(shí)戰(zhàn)篇4總結(jié)注：本講座針對有一定python和arcpy根底的朋友~第三局部內(nèi)容我邀請了我的同事范秦光同學(xué)為大家講解，他在server與python的結(jié)合應(yīng)用中具有豐富的經(jīng)驗(yàn)，這一局部內(nèi)容大概有十五分鐘的時(shí)間，在這里感謝范范！從未了解過ArcGIS中Python的朋友請先看….優(yōu)酷搜索關(guān)鍵詞“ArcGISPython”……1.現(xiàn)有工具實(shí)現(xiàn)不了，可以用python2.流程化需要時(shí)，可以使用python3.沒有AE等二次開發(fā)環(huán)境4.其他特殊場景什么情況下使用ArcPython？優(yōu)點(diǎn)：便捷、輕量級ArcPython在ArcGIS中的用法ArcGIS中python的幾種用法1.字段計(jì)算器2.制作成gp工具3.PythonwindowArcGIS中python的幾種用法4.PythonAdd-in5.在ArcGISServer中使用python與Desktop實(shí)戰(zhàn)篇為大家重磅打造典型應(yīng)用案例22個(gè)1.自增計(jì)算(字段計(jì)算器)total=0defaccumulate(increment):globaltotaliftotal:total+=incrementelse:total=incrementreturntotal為ID字段值設(shè)置自增2.記錄地塊權(quán)屬人變更的要素(字段計(jì)算器)判斷土地權(quán)屬人是否變化，不變那么Mark值為1，否那么為02.記錄地塊權(quán)屬人變更的要素(字段計(jì)算器)權(quán)屬人沒有變更權(quán)屬人變更將存儲(chǔ)在字段值中類似379.826的數(shù)據(jù)處理成K379+8263.樁號數(shù)據(jù)處理(字段計(jì)算器)第一步：!樁號!.replace(".","+")第二步：“K“+!樁號!交通行業(yè)通常需要將道路的里程數(shù)據(jù)存儲(chǔ)成K389+34的形式。1.自增計(jì)算2.記錄地塊權(quán)屬人變更的要素3.樁號數(shù)據(jù)處理操作演示：importarcpyfrom

arcpyimport

envenv.workspace=r'C:\Users\yanrui\Documents\ArcGIS\Default.gdb'fcs=arcpy.ListFeatureClasses()for

fcin

fcs:

arcpy.Delete_management(fc)4.刪除Default.gdb中的所有要素類、表、柵格批量刪除要素類：rasters=arcpy.ListRasters()for

rasterin

rasters:

arcpy.Delete_management(raster)批量刪除柵格：批量刪除表：tables=arcpy.ListTables()for

tablein

tables:

arcpy.Delete_management(table)ArcGISDesktop安裝好后環(huán)境變量中默認(rèn)輸出路徑是在C盤…的Default.gdb中，我們通常需要清理這個(gè)文件夾。但Toolbox沒有批量清理，這里可以使用python5.空間隨機(jī)抽取假設(shè)干數(shù)這個(gè)案例是幫助一個(gè)客戶從一個(gè)要素類中隨機(jī)抽取假設(shè)干個(gè)數(shù)據(jù)。然后再輸出。6.地震目錄自動(dòng)空間化第一步：Eclipse下載PyDev插件第二步：下載cx_Oracle.pyd包第三步：代碼實(shí)現(xiàn)第四步：Windows設(shè)置方案任務(wù)本練習(xí)在eclipse環(huán)境中完成(Eclipse環(huán)境中)將存儲(chǔ)在oracle中的地震數(shù)據(jù)〔包含x,y，時(shí)間，地點(diǎn)〕動(dòng)態(tài)的插入到矢量數(shù)據(jù)中。這個(gè)腳本是在eclipse環(huán)境中完成，為小伙伴們提供了文檔操作說明文檔。!type!+'%d'%!ID!7.相同字段值排序?qū)ν活惙N植地?cái)?shù)據(jù)進(jìn)行排序，原始字段是type，更新后是code字段。8.插入點(diǎn)(制作成GP工具、發(fā)布效勞、js頁面調(diào)用)使用arcpy結(jié)合python，給一個(gè)點(diǎn)圖層插入數(shù)據(jù)，制作成gp工具。8.插入點(diǎn)(制作成GP工具、發(fā)布效勞、js頁面調(diào)用)詳細(xì)操作步驟請參見操作手冊：工具發(fā)布成效勞，以及js頁面調(diào)用過程都在操作文檔中，文檔會(huì)提供給小伙伴們。4.刪除Default.gdb中的所有要素類、表、柵格5.空間隨機(jī)抽取假設(shè)干數(shù)6.地震目錄自動(dòng)空間化7.相同字段值排序8.插入點(diǎn)操作演示：9.字段拆分9.字段拆分10.讀取Excell數(shù)據(jù)生成多邊形11.讀取txt生成多邊形12.圖層批量合并9.字段拆分10.讀取Excell數(shù)據(jù)生成多邊形11.讀取txt生成多邊形12.圖層批量合并操作演示：13.shp數(shù)據(jù)至mdb14.打印線段折點(diǎn)坐標(biāo)importarcpymxd=arcpy.mapping.MapDocument(r“C:\Project\Project.mxd")mxd.findAndReplaceWorkspacePaths(r"C:\Project\Data",r"C:\Project\Data2")mxd.saveACopy(r"C:\Project\Project2.mxd")delmxd15.修復(fù)數(shù)據(jù)鏈接數(shù)據(jù)源批量替換：打印數(shù)據(jù)連接損壞的圖層：import

arcpymapdoc=arcpy.mapping.MapDocument("CURRENT")brokenlist=arcpy.mapping.ListBrokenDataSources(mapdoc)for

lyrin

brokenlist: print

del

mapdoc16.修改柵格數(shù)據(jù)分辨率ApplyEnvironment函數(shù)：僅僅限于Cellsize,

CurrentWorkspace,

Extent,

Mask,OutputCoordinateSystem,ScratchWorkspace,和SnapRaster等操作17.批量出圖18.坡度分級提取并計(jì)算面積13.Shp數(shù)據(jù)至mdb14.打印折線節(jié)點(diǎn)15.略16.修改柵格數(shù)據(jù)分辨率17.批量出圖18.坡度分級提取計(jì)算面積操作演示：下載PythonAdd-InWizard:

19.縮放至選中的features(PythonAdd-In)詳細(xì)操作步驟請參見操作手冊：pythonadd-in插件支持的最小版本是10.1，默認(rèn)也是10.1，我機(jī)器上安裝的是ArcGISDesktop10.3，所以需要修改config.xml里面的版本號，否那么在安裝插件時(shí)會(huì)報(bào)“NoGUIcomponentsfoundinthisAdd-in.Add-Inversiondoesnotmatch”：注意：20.柵格數(shù)據(jù)批量拷貝至對應(yīng)文件夾21.獲取多邊形頂點(diǎn)/DownloadfreeGISgeoproccessingtools,learnaboutESRI'sArcGISandArcPy,getfieldcalculatortips,andmore……22.制圖腳本樣例19.縮放至選中的features20.柵格數(shù)據(jù)批量拷貝至對應(yīng)文件夾21.獲取多邊形頂點(diǎn)20.制圖腳本樣例操作演示：下面的內(nèi)容，有請范秦光范范同學(xué)為大家講解~~~與Server實(shí)戰(zhàn)篇ArcGISServer命令行實(shí)用程序使用ArcGISRESTAPI編寫Python腳本利用ArcPy編寫Python腳本實(shí)現(xiàn)效勞的發(fā)布位于ArcGISforServer安裝目錄tools/admin目錄下

提供一套管理ArcGISServer的腳本命令行工具

1、ArcGISServer命令行實(shí)用程序

1、ArcGISServer命令行實(shí)用程序可以進(jìn)行以下操作啟動(dòng)和停止效勞從集群添加和移除計(jì)算機(jī)從效勞定義(.sd)文件發(fā)布效勞備份或復(fù)原站點(diǎn)配置定義新地圖緩存構(gòu)建和刪除地圖緩存在松散和緊湊格式之間轉(zhuǎn)換地圖緩存實(shí)戰(zhàn)演練：通過Terminal終端調(diào)用使用

1、ArcGISServer命令行實(shí)用程序?qū)崙?zhàn)案例:使用Python調(diào)用ArcGISAdminAPI修改ArcGISServer參數(shù)客戶需求介紹：在客戶的業(yè)務(wù)中，對ArcGISforServer的安裝全部是通過腳本自動(dòng)一鍵執(zhí)行，因?yàn)榭蛻羰褂玫腉P效勞以及地圖效勞數(shù)量多高并發(fā)并且要求7X24小時(shí)運(yùn)行，因此通過Python調(diào)用ArcGISAdminAPI構(gòu)造可執(zhí)行腳本，ArcGISServer安裝完成后，調(diào)用該腳本，調(diào)整ArcGISforServer的默認(rèn)配置。2、使用ArcGISRESTAPI編寫Python腳本實(shí)戰(zhàn)案例:使用Python調(diào)用ArcGISAdminAPI修改ArcGISServer參數(shù)通過Python構(gòu)筑HTTPrequests，將Appservermaximumheapsize和SOCmaximumheapsize的值分別修改為512和128。2、使用ArcGISRESTAPI編寫Python腳本實(shí)戰(zhàn)案例:使用Python調(diào)用ArcGISAdminAPI修改ArcGISServer參數(shù)

步驟1運(yùn)行腳本時(shí)需要有管理員權(quán)限，必須提供一個(gè)具有管理員權(quán)限的用戶名和密碼：

步驟2在腳本中獲取和使用令牌步驟3配置參數(shù)，構(gòu)筑請求2、使用ArcGISRESTAPI編寫Python腳本案例1、使用Python發(fā)布構(gòu)筑GP效勞客戶需求介紹：大量的圖片〔.tif.jpg〕并且?guī)в杏涗涀鴺?biāo)信息的文件〔.tfw.jpgw〕，需要將這些文件導(dǎo)入到SDE數(shù)據(jù)庫中，作為背景底圖使用。

3、利用ArcPy編寫Python腳本實(shí)現(xiàn)效勞的發(fā)布Python、ArcPyAPIArcCatalog\ArcToolboxArcMap(ShareAsService)或ArcPy發(fā)布SD設(shè)計(jì)和制作腳本制作腳本工具調(diào)試運(yùn)行、發(fā)布效勞設(shè)計(jì)和制作腳本

步驟1定義參數(shù)

步驟2將圖片作為柵格數(shù)據(jù)集導(dǎo)入到臨時(shí)工作空間中、并定義空間參考步驟3配進(jìn)行投影變換，輸出路徑設(shè)置為SDE數(shù)據(jù)庫3、利用ArcPy編寫Python腳本實(shí)現(xiàn)效勞的發(fā)布制作腳本工具3、利用ArcPy編寫Python腳本實(shí)現(xiàn)效勞的發(fā)布調(diào)試運(yùn)行、發(fā)布效勞3、利用ArcPy編寫Python腳本實(shí)現(xiàn)效勞的發(fā)布案例2、使用Python發(fā)布地圖效勞客戶需求介紹：在客戶的業(yè)務(wù)中，預(yù)先定義好地圖文檔〔.mxd〕,根據(jù)業(yè)務(wù)數(shù)據(jù)變化動(dòng)態(tài)替換其中的圖層數(shù)據(jù)，然后將其發(fā)布成地圖效勞。準(zhǔn)備地圖文檔構(gòu)思設(shè)計(jì)制作發(fā)布地圖效勞的GP腳本定義參數(shù)添加?xùn)鸥駡D層生成草稿〔sddraft〕生成效勞定義〔sd〕發(fā)布