用ArcGIS的Python腳本計算采煤沉陷區(qū)土地損毀程度

朱寶才，王文龍，梁文俊

(1.山西農(nóng)業(yè)大學(xué) 林學(xué)院，山西 太谷 030801； 2.西北農(nóng)林科技大學(xué) 水土保持研究所黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室，陜西 楊凌 712100； 3.中國科學(xué)院 水利部水土保持研究所，陜西 楊凌 712100)

?

用ArcGIS的Python腳本計算采煤沉陷區(qū)土地損毀程度

朱寶才1,2，王文龍2,3*，梁文俊1

(1.山西農(nóng)業(yè)大學(xué) 林學(xué)院，山西 太谷 030801； 2.西北農(nóng)林科技大學(xué) 水土保持研究所黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室，陜西 楊凌 712100； 3.中國科學(xué)院 水利部水土保持研究所，陜西 楊凌 712100)

摘要：為計算井工煤礦采煤沉陷區(qū)土地損毀程度，利用ArcGIS的Python腳本語言編制軟件，實現(xiàn)通過輸入相同坐標系的下沉等值線、附加傾斜等值線及水平變形等值線3個DWG格式文件和土地利用現(xiàn)狀SHP格式文件，直接輸出土地損毀程度的SHP格式文件，從而獲得沉陷區(qū)土地損毀程度預(yù)測圖。通過對三元吉祥煤礦的實例分析，該軟件亦可直接用于其他井工煤礦開采沉陷區(qū)土地損毀程度計算，且方便基層技術(shù)人員使用和推廣。

關(guān)鍵詞：采煤沉陷區(qū)；土地損毀；預(yù)測；GIS；Python

2011年5月，國土資源部制定《土地復(fù)墾方案編制規(guī)程》，共分7個部分，其中第3部分：井工煤礦(TD/T 1031.3 —2011)附錄A 給出井工煤礦開采土地損毀預(yù)測方法，已有成熟軟件[1～8]計算下沉、傾斜變形和水平變形；附錄 B 給出采煤沉陷土地損毀程度分級參考標準，尚無成熟軟件進行土地損毀程度分級計算，要在GIS軟件中結(jié)合土地利用現(xiàn)狀圖手動分析，效率較低，且容易出錯。為快速準確地計算采煤沉陷土地損毀程度，本研究以三元吉祥煤礦為例根據(jù)土地復(fù)墾方案編制規(guī)程的采煤沉陷土地損毀程度分級參考標準，基于井工煤礦開采土地損毀預(yù)測和土地利用現(xiàn)狀數(shù)據(jù)，用ArcGIS 腳本語言Python編制“采煤沉陷土地損毀程度分級計算軟件V1.0”(已登記計算機軟件著作權(quán)，證書號：軟著登字第1086259號)實現(xiàn)采煤沉陷區(qū)土地損毀程度計算。

1研究區(qū)概況與研究方法

1.1研究區(qū)概況

三元吉祥煤礦地理坐標為東經(jīng) 113°03′09″～113°05′11″，北緯 36°22′07″～36°22′56″，海拔最高1 029.92 m，最低899.40 m，位于山西省長治市郊區(qū)黃碾鎮(zhèn)，太行山中段西側(cè)，長治盆地東部，總體北高南低，廣為第四系黃土覆蓋，北部有少許基巖出露，溝谷較發(fā)育，無常年性河流，雨季洪水匯入濁漳河。屬于暖溫帶半濕潤大陸性季風(fēng)氣候，四季分明，雨熱同季，溫和適中。據(jù)長治市氣象局1961—2010年資料統(tǒng)計，年均氣溫 9.1 ℃，最高氣溫37.6 ℃，最低氣溫為-23.0 ℃。年均降雨量595 mm，最大降雨量832.9 mm(1971年)，最小降雨量340.9 mm(1965年)，雨季多集中在七、八2月，年均蒸發(fā)量1 558 mm，年平均無霜期170 d，冰凍期為11月至次年3月，最大凍土深度75 cm。一般冬季多西北風(fēng)，夏季多東南風(fēng)，最大風(fēng)力10級。土壤類型為褐土性壤土，農(nóng)作物主要是玉米和小麥，喬木主要是油松和刺槐，灌木以荊條、黃刺玫和胡枝子等為主，草本以白羊草草叢和蒿類草叢為主。

1.2數(shù)據(jù)來源

土地利用現(xiàn)狀圖數(shù)據(jù)源自第二次全國土地調(diào)查數(shù)據(jù)和2011年土地利用變更數(shù)據(jù)，為3度帶高斯投影超圖格式文件，在超圖中轉(zhuǎn)換文件格式為SHP格式，并重命名為LandUse_3.shp。

井工煤礦開采土地損毀預(yù)測數(shù)據(jù)采用概率積分法，通過輸入與土地利用現(xiàn)狀數(shù)據(jù)坐標系統(tǒng)一致的地形圖、煤層底板等高線和工作面平面位置等數(shù)據(jù)，應(yīng)用中國礦業(yè)大學(xué)“任意形多工作面多線段開采沉陷預(yù)計系統(tǒng)(MSPS)[1]”獲得，包含10 mm下沉等值線和起點與間距均為0.5 m的下沉等值線文件命名為XiaChen.DWG。包含起點為0 mm·m-1，間距為2 mm·m-1的附加傾斜等值線文件命名為QingXie.DWG，包含起點為0 mm·m-1，間距為1 mm·m-1的水平變形等值線文件命名為ShuiPing.DWG。

以上數(shù)據(jù)皆存儲于c:data中。

1.3數(shù)據(jù)處理

依據(jù)《土地復(fù)墾方案編制規(guī)程—第3部分：井工煤礦》推薦的“采煤沉陷土地損毀程度分級參考標準”，應(yīng)用極限條件法確定各地類各等級土地損毀范圍，再按原土地利用現(xiàn)狀圖、輕度以上損毀范圍，中度以上損毀范圍和重度損毀范圍自下而上的順序疊加即得土地損毀程度圖，其中：

輕度以上損毀范圍：依據(jù)國家煤炭工業(yè)局2000年制定的《建筑物、水體、鐵路及主要井巷煤柱留設(shè)與壓煤開采規(guī)程》確定10 mm下沉等值線包含區(qū)域。

水田中度以上損毀范圍：{下沉>1.0 m區(qū)域}∪{水平變形>3.0 mm·m-1區(qū)域}∪{附加傾斜>4.0 mm·m-1區(qū)域}；水澆地中度以上損毀范圍：{下沉>1.5 m區(qū)域}∪{水平變形>4.0 mm·m-1區(qū)域}∪{附加傾斜>6.0mm·m-1區(qū)域}；旱地中度以上損毀范圍：{下沉>2.0 m區(qū)域}∪{水平變形>8.0 mm·m-1區(qū)域}∪{ 附加傾斜20.0 mm·m-1區(qū)域}；林地、草地中度以上損毀范圍：{下沉>2.0 m區(qū)域}∪{水平變形>8.0 mm·m-1區(qū)域}∪{ 附加傾斜>20.0 mm·m-1區(qū)域}。

水田重度損毀范圍：{下沉>2.0 m區(qū)域}∪{水平變形>6.0 mm·m-1區(qū)域}∪{附加傾斜>10.0 mm·m-1區(qū)域}；水澆地中度損毀范圍：{下沉>3.0 m區(qū)域}∪{水平變形>8.0 mm·m-1區(qū)域}∪{ 附加傾斜>12.0 mm·m-1區(qū)域}；旱地中度損毀范圍：{下沉>5.0 m區(qū)域}∪{水平變形>16.0 mm·m-1區(qū)域}∪{附加傾斜>40.0 mm·m-1區(qū)域}；林地、草地中度損毀范圍：{下沉>6.0 m區(qū)域}∪{水平變形>20.0 mm·m-1區(qū)域}∪{附加傾斜>50.0 mm·m-1區(qū)域}

1.4Python實現(xiàn)

Python是開源編程語言，ArcGIS9.0開始引入Python。封裝大部分ArcGIS 空間分析工具于Arcgisscripting 站點包中，從ArcGIS10.0 開始，站點包更名為Arcpy，Python腳本可獨立于ArcGIS運行[9～14]。

1.導(dǎo)入ArcPy。

ArcPy提供所有地理處理工具使用 Python 語言操作的接口，并為處理和查詢 GIS 數(shù)據(jù)提供函數(shù)和類，如env 類包含所有地理處理環(huán)境。通過導(dǎo)入env類，設(shè)置工作空間為c:/data文件夾，代碼為：

import arcpy

from arcpy import env

env.workspace = “C:/data”

2.創(chuàng)建個人地理數(shù)據(jù)庫存放臨時文件

在C:/data創(chuàng)建名為Temp.gdb的個人地理數(shù)據(jù)庫存放地理處理的臨時文件，方便用后統(tǒng)一刪除，代碼為：

arcpy.CreateFileGDB_management(env.workspace, “Temp.gdb”)

3.分別轉(zhuǎn)換CAD格式的下沉等值線文件XiaChen.DWG，水平變形等值線文件ShuiPing.DWG和附加傾斜等值線文件QingXie.DWG為地理數(shù)據(jù)庫Temp.gdb的要素類XiaChen，ShuiPing和QingXie。以轉(zhuǎn)換CAD格式下沉等值線為地理數(shù)據(jù)庫為例，代碼如下：

arcpy.CADToGeodatabase_conversion(“XiaChen.DWG”, “Temp.gdb”, “XiaChen”, “1000”)

拷貝土地利用現(xiàn)狀LandUse_3.shp為Temp.gdb的要素類LandUse_3。代碼為：

arcpy.CopyFeatures_management(“LandUse_3.shp”,“Temp.gdb/LandUse_3”)

4.用要素類XiaChen，ShuiPing和QingXie等值線創(chuàng)建臨時圖層，擬存放臨界下沉等值線(值為10, 1000, 1500, 2000, 3000, 5000和6000)、臨界水平變形等值線(值為3, 4, 6, 8, 16, 20)和臨界附加傾斜等值線(值為4,6,10,12,20,40,50)，圖層名命名采取“要素類名”+“_”(下劃線)+“臨界值”的方法，擬存放下沉值為10的等值線圖層命名為“XiaChen_10”；擬存放水平變形值為3的圖層命名為“ShuiPing_3”；擬存放附加傾斜值為4的圖層命名為“QingXie_4”。其代碼示例如下：

arcpy.MakeFeatureLayer_management(“Temp.gdb/XiaChen/Polyline”, “XiaChen_10”)

arcpy.MakeFeatureLayer_management(“Temp.gdb/ShuiPing/Polyline1”, “ShuiPing_3”)

arcpy.MakeFeatureLayer_management(“Temp.gdb/QingXie/Polyline2”, “QingXie_4”)

用要素類LandUse_3創(chuàng)建臨時圖層，擬存放水田、水澆地、旱地，林地和草地及全部地類的圖層分別命名為LandUse_011、LandUse_012、LandUse_013和LandUse_03_04和LandUse。

5.通過屬性“Elevation”選擇臨界下沉等值線、臨界水平變形等值線和臨界附加傾斜等值線，圖層名中包含“臨界值”即為臨界等值線值，XiaChen_10選擇的屬性即為“Elevation”=10。因水平變形和附加傾斜值用正負號表示變形方向，故對其等值線值取絕對值。其代碼示例如下：

arcpy.SelectLayerByAttribute_management(“XiaChen_10”, “NEW_SELECTION”, ‘“Elevation”= 10’)

arcpy.SelectLayerByAttribute_management(“ShuiPing_3”, “NEW_SELECTION”, ‘Abs(“Elevation”)=3’)

arcpy.SelectLayerByAttribute_management(“QingXie_4”, “NEW_SELECTION”, ‘Abs(“Elevation”) = 4’)

通過地類編碼屬性“DLBM”選擇地類，011為水田，012為水澆地，013為旱地，031為有林地，以上各類通過新選擇，再將032、033、041、042和043添加選擇到LandUse_03_04圖層中。其代碼示例如下：

arcpy.SelectLayerByAttribute_management(“LandUse_03_04”, “NEW_SELECTION”, ‘CAST(“DLBM” AS INTEGER) = 031’)

arcpy.SelectLayerByAttribute_management(“LandUse_03_04”, “ADD_TO_SELECTION”, ‘CAST(“DLBM” AS INTEGER) = 032’)

6.將臨界等值線轉(zhuǎn)換為臨界面，輸入文件為已通過對應(yīng)屬性選擇的臨界等值線圖層，輸出文件命名為“polygon_”+原臨界等值線文件名，如XiaChen_10等值線轉(zhuǎn)換為面后，命名為“polygon_XiaChen_10”，其代碼示例如下：

arcpy.FeatureToPolygon_management(“XiaChen_10”, “Temp.gdb/polygon_XiaChen_10”, “”, “NO_ATTRIBUTES”, “”)

7.確定輕度及輕度以上損毀范圍。在polygon_XiaChen_10屬性表中添加損毀等級字段“SHDJ”，字段類型為SHORT，在字段中可存儲9位。并將“SHDJ”字段統(tǒng)一賦值為“1”，即輕度損毀。代碼實現(xiàn)如下：

arcpy.AddField_management(“Temp.gdb/polygon_XiaChen_10”, “SHDJ”, “SHORT”, 9, “”, “”,“SHDJ”, “NULLABLE”)

arcpy.CalculateField_management(“Temp.gdb/polygon_XiaChen_10”, “SHDJ”, “1”, “”)

8.分地類確定中度及中度以上和重度可能損毀范圍。依據(jù)下沉、水平變形、附加傾斜臨界值聯(lián)合確定水田、水澆地、旱地、林草地中度及中度以上可能損毀范圍，要素類分別命名為Land011_Class_2，Land012_Class_2，Land013_Class_2和Land03_04_Class_2。水田、水澆地、旱地、林草地重度可能損毀范圍要素類分別命名為Land011_Class_3，Land012_Class_3， Land013_Class_3和Land03_04_Class_3。其實現(xiàn)代碼示例如下：

arcpy.Union_analysis ([“Temp.gdb/polygon_XiaChen_1000”,“Temp.gdb/polygon_ShuiPing_3”, “Temp.gdb/polygon_QingXie_4”], “Temp.gdb/Land011_Class_2”,“ONLY_FID”,“”)

在以上各要素屬性表中添加損毀等級字段“SHDJ”，設(shè)置同前，并將損毀等級為中度及中度以上要素(即名稱中有Class_2)的“SHDJ”字段統(tǒng)一賦值為“2”，損毀等級為重度要素(即名稱中有Class_3)的“SHDJ”字段統(tǒng)一賦值為“3”。

9.融合，為使制圖中同一地塊不被分割成具有相同屬性的多個地塊，將以上要素按“SHDJ”字段進行融合，融合后的要素名在原要素名上加“_dissolve”。其代碼示例如下：

arcpy.Dissolve_management(“Temp.gdb/Land011_Class_2”, “Temp.gdb/Land011_Class_2_dissolve”, “SHDJ”, “”, “”, “”)

10.分地類確定中度和重度可能損毀范圍。更新要素類屬性，用水田重度可能損毀范圍Land011_Class_3_dissolve更新中度及中度以上可能損毀范圍Land011_Class_2_dissolve，得到水田中度和重度可能損毀范圍Land011_Class_2_3，代碼實現(xiàn)如下：

arcpy.Update_analysis(“Temp.gdb/ Land011_Class_2_dissolve”, “Temp.gdb/Land011_Class_3_dissolve”,“Temp.gdb/Land011_Class_2_3”)

同樣方法，分別得到水澆地、旱地、林草地的中度和重度可能損毀范圍要素類Land012_Class_2_3，Land013_Class_2_3和Land03_04_Class_2_3。

11.標識(Identity)要素類，用輕度及輕度以上損毀范圍要素Temp.gdb/polygon_XiaChen_10標識土地利用Temp.gdb/LandUse，得到輕度及輕度以上損毀土地Temp.gdb/LandUse_Damage_1，實現(xiàn)代碼如下：

arcpy.Identity_analysis (“Temp.gdb/LandUse”, “Temp.gdb/polygon_XiaChen_10”, “Temp.gdb/LandUse_Damage_1”)

用水田中度和重度可能損毀范圍Land011_Class_2_3標識水田土地利用LandUse_011，得到水田中度和重度損毀范圍LandUse011_Damage_2_3。實現(xiàn)代碼如下：

arcpy.Identity_analysis (“LandUse_011”, “Temp.gdb/Land011_Class_2_3”, “Temp.gdb/LandUse011_Damage_2_3”)

同樣方法分別得到水澆地、旱地、林草地的中度和重度土地損毀要素類Land012_ Damage _2_3，Land013_ Damage _2_3和Land03_04_ Damage _2_3。

12.逐地類更新(Update)，用水田的中度和重度土地損毀要素類LandUse011_Damage_2_3更新輕度及輕度以上損毀土地LandUse_Damage_1，再用水澆地的中度和重度土地損毀要素類Land012_ Damage _2_3更新其結(jié)果，依次Land013_ Damage _2_3和Land03_04_ Damage _2_3更新前一更新結(jié)果，實現(xiàn)代碼示例如下：

arcpy.Update_analysis (“Temp.gdb/LandUse_Damage_1”, “Temp.gdb/LandUse011_Damage_2_3”, “Temp.gdb/LandUse_Damage_011OK”)

arcpy.Update_analysis (“Temp.gdb/LandUse_Damage_013OK ”, “Temp.gdb/LandUse03_04_Damage_2_3”, “LandUse_Damage”)

將上述代碼保存為腳本文件Damage.py。

1.5軟件運行

在ArcGIS10.2可運行的狀態(tài)下，用Python2.7 打開本腳本程序文件Damage.py，按F5鍵運行腳本。

為獲得數(shù)據(jù)處理所用時間，在開始運行軟件時獲取系統(tǒng)時間，單位是秒，數(shù)據(jù)處理結(jié)束后再次獲取系統(tǒng)時間，兩者之差即為數(shù)據(jù)處理總時間。本次數(shù)據(jù)處理在Intel(R) Core(TM) i5-4210Y CPU@1.50 GHz四核處理器，4.00 GB內(nèi)存，128 GB固態(tài)硬盤的條件下運行240 s。

計算結(jié)果土地損毀程度文件LandUse_Damage.shp存儲于C:data文件夾中，中間過程數(shù)據(jù)存儲于Temp.gdb地理數(shù)據(jù)庫中，LandUse_Damage.shp文件相對于LandUse_3.shp增加了三個字段，分別是SHDJ(損毀等級)、Shape_Leng(周長)和Shape_Area(面積)，SHDJ字段所列數(shù)字即為該小班采煤沉陷土地損毀程度，其中0代表無損毀；1代表輕度損毀；2代表中度損毀；3代表重度損毀。Shape_Area字段為該小班面積。

2結(jié)果與分析

基于“采煤沉陷土地損毀程度分級計算軟件V1.0”計算的三元吉祥煤礦土地損毀程度預(yù)測面積統(tǒng)計見表1，由于篇幅所限，未附土地損毀程度預(yù)測圖。各地類輕度損毀總面積為361.47 hm2，占損毀總面積的63.20%，重度損毀總面積123.42 hm2，占損毀總面積的21.58%，可見三元吉祥煤礦開采沉陷區(qū)土地損毀以輕度為主。損毀的耕地全部為旱地，面積為257.88 hm2，占總損毀土地面積的45.09%，其中輕度損毀占23.85%，中度損毀占7.80%，重度損毀占13.44%。損毀林地包括有林地和其他林地，面積為121.74 hm2，占總損毀土地面積的21.28%，其中輕度損毀占11.79，中度損毀占5.17%，重度損毀占4.33%。損毀草地全部為其他草地，面積為66.32 hm2，占總損毀土地面積的17.67%，其中輕度損毀占11.60%，中度損毀占2.26%，重度損毀占3.81%。耕地林地草地共損毀480.66 hm2，占損毀總面積的84.04%，主要原因是礦區(qū)內(nèi)以耕地林地草地為主。

表1三元吉祥煤礦采煤沉陷區(qū)土地損毀程度預(yù)測面積統(tǒng)計表/hm2

Table 1Sanyuan Jixiang coal mining subsidence land damage statistical prediction of surface area

地類名稱Landuse輕度損毀Lightdamage中度損毀Mediumdamage重度損毀Severedamage總計Total旱地136.4344.6076.85257.88有林地31.2117.4912.5161.21其他林地36.2012.0812.2560.53其他草地66.3212.9021.82101.04鐵路用地3.140.000.003.14農(nóng)村道路0.300.000.000.30設(shè)施農(nóng)用地1.260.000.001.26沼澤地0.710.000.000.71裸地6.050.000.006.05城市1.420.000.001.42村莊47.260.000.0047.26采礦用地30.830.000.0030.83風(fēng)景名勝及特殊用地0.340.000.000.34總計361.4787.07123.42571.96

基于“采煤沉陷土地損毀程度分級計算軟件V1.0”計算的三元吉祥煤礦土地損毀程度預(yù)測小班731個，基于手動ArcGIS預(yù)測方法，獲得土地損毀程度小班也是731個，將兩種方法獲得的土地損毀程度小班分別按面積從大到小排序，逐一對比，發(fā)現(xiàn)對應(yīng)地類完全相同，對應(yīng)面積差小于0.001 m2，相對誤差小于0.1%。

3結(jié)論和討論

基于“采煤沉陷土地損毀程度分級計算軟件V1.0”進行煤礦土地損毀程度預(yù)測可達到手動ArcGIS分析相同的精度。利用“采煤沉陷土地損毀程度分級計算軟件V1.0”能快速地計算采煤沉陷區(qū)土地損毀程度，可直接用于其他井工煤礦開采沉陷區(qū)土地損毀程度計算，為土地復(fù)墾方案，生態(tài)環(huán)境恢復(fù)治理方案及水土保持方案等提供準確的土地損毀程度數(shù)據(jù)。

“采煤沉陷土地損毀程度分級計算軟件V1.0” 處理時間大量縮短，技術(shù)要求顯著降低，可不打開ArcGIS，直接在Python2.7中運行，運行前只需把土地利用現(xiàn)狀、下沉等值線、附加傾斜等值線和水平變形等值線文件按指定格式和名稱存儲于指定位置即可，方便簡單，降低了采煤沉陷區(qū)土地損毀程度計算門檻，方便基層技術(shù)人員使用和推廣。

參考文獻

[1]蔡懷恩. 開采沉陷預(yù)計的方法及發(fā)展趨勢[J]. 露天采礦技術(shù), 2007, 22(4):43-44.

[2]康建榮, 王金莊, 溫澤民. 任意形多工作面多線段開采沉陷預(yù)計系統(tǒng)(MSPS)[J]. 礦山測量, 2000, 28(1):24-27.

[3]李明, 李琰慶, 王紅梅. 任意形狀工作面開采沉陷預(yù)計系統(tǒng)開發(fā)[J]. 礦業(yè)安全與環(huán)保, 2008, 35(5):18-21.

[4]李培現(xiàn), 譚志祥, 齊公玉, 等. 基于MATLAB的開采沉陷預(yù)計系統(tǒng)[J]. 中國礦業(yè), 2008,17(11):72-76.

[5]劉虎, 吳侃, 陳冉麗. 基于AutoCAD的開采沉陷預(yù)計分析軟件研究[J]. 礦業(yè)安全與環(huán)保, 2009, 36(2):22-23.

[6]謝和平, 周宏偉, 王金安, 等. FLAC在煤礦開采沉陷預(yù)測中的應(yīng)用及對比分析[J]. 巖石力學(xué)與工程學(xué)報, 1999, 18(4):29-33.

[7]張慶松, 高延法, 劉松玉, 等. 基于粗集與神經(jīng)網(wǎng)絡(luò)相結(jié)合的巖移影響因素分析與開采沉陷預(yù)計方法研究[J]. 煤炭學(xué)報, 2004, 29(1):22-25.

[8]柴華彬, 鄒友峰, 袁占良, 等. SuperMap系統(tǒng)在開采沉陷預(yù)計分析中的應(yīng)用[J]. 西安科技大學(xué)學(xué)報, 2005,25(1):52-56.

[9]方圣輝, 張玉賢, 佃袁勇, 等. 基于Python的ArcGIS地理數(shù)據(jù)批處理[J]. 測繪與空間地理信息, 2015, 38(1):1-2.

[10]焦洋, 鄧鑫, 李勝才. 基于Python的ArcGIS空間數(shù)據(jù)格式批處理轉(zhuǎn)換工具開發(fā)[J]. 現(xiàn)代測繪, 2012, 35(3):54-55.

[11]李強, 白建榮, 李振林, 等. 基于Python的數(shù)據(jù)批處理技術(shù)探討及實現(xiàn)[J]. 地理空間信息, 2015, 13(2):54-56.

[12]彭海波, 向洪普. 基于Python的空間數(shù)據(jù)批量處理方法[J]. 測繪與空間地理信息, 2011, 34(4):81-82.

[13]邵保華, 田學(xué)志. 談Python在Arcgis地理處理中的應(yīng)用[J]. 林業(yè)勘查設(shè)計, 2012, 41(2):99-100.

[14]朱海金, 翁偉琳. 基于python數(shù)字高程模型地形數(shù)據(jù)批量提取[J]. 地理空間信息, 2013, 11(5):136-138.

(編輯：馬榮博)

Calculate land damage degree in the coal mining subsidence area by Python, ArcGIS script language

Zhu Baocai1, 2, Wang Wenlong2, 3*， Liang Wenjun1

(1.CollegeForestry,ShanxiAgriculturalUniversity,Taigu030801,China; 2.StateKeyLaboratoryofSoilandErosionandDrylandFarmingontheLoessPlateau,InstituteofSoilandWaterConservation,NorthwestA&FUniversity,Yangling712100,China; 3.InstituteofSoilandWaterConservation,ChineseAcademyofSciencesandMinistryofWaterResources,Yangling712100,China)

Abstract:For calculating land damage degree in the coal mining subsidence area, the Python that attached to ArcGIS, a script language compilation software, were used to predict the land damage degree area of every land use type. Land damaged degree result characterized as the SHP file format can be outputted directly and then the predicted figures of land damaged can be obtained, through the way that the subsidence contour, additional tilting contour and horizontal deformation contour files with DWG format and the land use file with SHP format were inputted altogether in the Python software. The satisfactory result was obtained from SanYuanJiXiang Coal Mine. Therefore, the software can be applied directly to calculate land damaged degree in other coal mining subsidence area. It is promoted easily and convenient for the application of basic technician.

Key words:Coal mining subsidence area; Land damage; Prediction; GIS; Python

中圖分類號：F301.24

文獻標識碼：A

文章編號：1671-8151(2016)05-0371-06

基金項目：國家自然科學(xué)基金項目(31500523，41501201，41571275)，中國科學(xué)院西部行動計劃項目(KZCX2-XB3-13-4)

作者簡介：朱寶才(1978-)，男(漢)，內(nèi)蒙古巴林左旗人，講師，博士研究生，研究方向：土壤侵蝕與水土保持*通訊作者：王文龍，博士，研究員，博士生導(dǎo)師。Tel：13669222856；E-mail：wlwang@nwsuaf.edu.cn

收稿日期：2015-12-18 修回日期：2016-01-28