您好,登錄后才能下訂單哦!
這篇文章給大家分享的是有關(guān)python中arcpy的示例分析的內(nèi)容。小編覺(jué)得挺實(shí)用的,因此分享給大家做個(gè)參考,一起跟隨小編過(guò)來(lái)看看吧。
有多個(gè)文件地理數(shù)據(jù)庫(kù)(gdb),數(shù)據(jù)庫(kù)內(nèi)有多個(gè)面要素類圖層,每個(gè)圖層不能有自重疊,也不能和其他圖層重疊。所以,需要為每個(gè)文件地理數(shù)據(jù)庫(kù)(gdb)進(jìn)行拓?fù)錂z查。
1.在輸出文件夾下新建與輸入文件夾下同名的gdb。
2.在gdb內(nèi)新建要素類數(shù)據(jù)集。
3.在數(shù)據(jù)集內(nèi)創(chuàng)建拓?fù)洹?/p>
4.向拓?fù)渲刑砑右仡悺?/p>
5.添加拓?fù)湟?guī)則。
6.拓?fù)潋?yàn)證。
7.導(dǎo)出拓?fù)溴e(cuò)誤。
將需要拓?fù)涞膅db放入同一個(gè)文件夾下,用python代碼串連起來(lái)利用循環(huán)遍歷gdb便可以進(jìn)行批量拓?fù)淞恕?/p>
# -*- coding: gbk -*- import os import arcpy import sys reload(sys) sys.setdefaultencoding('utf-8') arcpy.AddMessage(sys.getdefaultencoding()) arcpy.env.XYResolution = "0.0001 Meters" #設(shè)置XY分辨率 arcpy.env.XYTolerance = "0.001 Meters" #設(shè)置XY容差 def checkTopology(in_path, out_path): cout = 0 fail = 0 warning = 0 fail_list = [] warning_list = [] arcpy.env.workspace = in_path workspaces = arcpy.ListWorkspaces("*", "ALL") for workspace in workspaces: arcpy.AddMessage("=" * 60) cout += 1 arcpy.AddMessage(" (" + str(cout) + ") " + os.path.basename(workspace)) arcpy.AddMessage(" ") try: # 將每個(gè)gdb設(shè)為工作區(qū) arcpy.env.workspace = workspace fc_list = arcpy.ListFeatureClasses() #列出gdb內(nèi)的要素類圖層 in_fc_path_list = [] for fc in fc_list: in_fc_path_list.append(os.path.join(workspace,fc)) # 在輸出路徑out_path下新建gdb gdb_name1 = os.path.basename(workspace) result_gdb = os.path.join(out_path, gdb_name1) if not arcpy.Exists(result_gdb): arcpy.AddMessage(" 新建 " + result_gdb + " ...") arcpy.CreateFileGDB_management(out_path, gdb_name1) else: arcpy.AddWarning(result_gdb + " 已存在!") warning_list.append(result_gdb) warning += 1 continue arcpy.AddMessage(" 在 " + result_gdb + " 內(nèi)新建要素?cái)?shù)據(jù)集...") dataset_name = "dataset" dataset_path = os.path.join(result_gdb,dataset_name) # 引用包含要應(yīng)用的空間參考的要素類或要素?cái)?shù)據(jù)集, # 將in_fc_path_list[0]的坐標(biāo)作為要素?cái)?shù)據(jù)集的坐標(biāo),所以gdb內(nèi)的所有要素類應(yīng)該是統(tǒng)一坐標(biāo) arcpy.CreateFeatureDataset_management(result_gdb, dataset_name, in_fc_path_list[0]) #導(dǎo)入要素類到數(shù)據(jù)集 arcpy.AddMessage(" 向 " + dataset_path + " 導(dǎo)入要素類圖層...") arcpy.FeatureClassToGeodatabase_conversion(in_fc_path_list,dataset_path) #在dataset數(shù)據(jù)內(nèi)創(chuàng)建拓?fù)? arcpy.AddMessage(" 創(chuàng)建拓?fù)?..") topology_name = "topology" topology_path = os.path.join(dataset_path, topology_name) arcpy.CreateTopology_management(dataset_path, topology_name) # 將導(dǎo)入dataset的要素類添加到拓?fù)渲? arcpy.AddMessage(" 向拓?fù)渲刑砑右仡?..") dataset_fc_path_lsit = [] #數(shù)據(jù)集中的要素類絕對(duì)路徑列表 for ifc_name in fc_list: # 拼接數(shù)據(jù)集中的要素類絕對(duì)路徑 dataset_fc_path = os.path.join(dataset_path,ifc_name) dataset_fc_path_lsit.append(dataset_fc_path) arcpy.AddFeatureClassToTopology_management(topology_path, dataset_fc_path, "1", "1") arcpy.AddMessage(" 添加拓?fù)湟?guī)則...") for i in range(len(dataset_fc_path_lsit)): fc_path2 = dataset_fc_path_lsit[i] #規(guī)則1:不能重疊 arcpy.AddRuleToTopology_management(topology_path, "Must Not Overlap (Area)", fc_path2, "", "", "") for j in range(i + 1, len(dataset_fc_path_lsit)): fc_path3 = dataset_fc_path_lsit[j] # 規(guī)則2:不能重疊與其他要素類重疊 arcpy.AddRuleToTopology_management(topology_path, "Must Not Overlap With (Area-Area)", fc_path2, "", fc_path3, "") arcpy.AddMessage(" 拓?fù)潋?yàn)證...") arcpy.ValidateTopology_management(topology_path, "Full_Extent") arcpy.AddMessage(" 導(dǎo)出拓?fù)溴e(cuò)誤...") arcpy.ExportTopologyErrors_management(topology_path, dataset_path, "topoError") except Exception as e: arcpy.AddError(e.message) arcpy.AddError(" " + os.path.basename(workspace) + " 失敗!") fail_list.append(os.path.basename(workspace)) fail += 1 arcpy.AddMessage('+' * 60) arcpy.AddMessage(u" 成功:" + str(cout - fail - warning) + u" 個(gè)!") if warning > 0: arcpy.AddWarning(u" 警告:" + str(warning) + u" 個(gè)! 如下:") arcpy.AddWarning("####" + '*' * 20) for fff in warning_list: arcpy.AddWarning(" " + fff) arcpy.AddWarning("####" + '*' * 20) if fail > 0: arcpy.AddError(u" 失?。?quot; + str(fail) + u" 個(gè)! 如下:") arcpy.AddError("####" + '*' * 20) for ff in fail_list: arcpy.AddError(" " + ff) arcpy.AddError("####" + '*' * 20) arcpy.AddMessage('+' * 60) if __name__ == "__main__": in_path_globle = arcpy.GetParameterAsText(0) out_path_globe = arcpy.GetParameterAsText(1) checkTopology(in_path_globle,out_path_globe)
感謝各位的閱讀!關(guān)于“python中arcpy的示例分析”這篇文章就分享到這里了,希望以上內(nèi)容可以對(duì)大家有一定的幫助,讓大家可以學(xué)到更多知識(shí),如果覺(jué)得文章不錯(cuò),可以把它分享出去讓更多的人看到吧!
免責(zé)聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng),如果涉及侵權(quán)請(qǐng)聯(lián)系站長(zhǎng)郵箱:is@yisu.com進(jìn)行舉報(bào),并提供相關(guān)證據(jù),一經(jīng)查實(shí),將立刻刪除涉嫌侵權(quán)內(nèi)容。