基于ArcPy的河南省小流域劃分輔助工具集設(shè)計(jì)與實(shí)現(xiàn)

王松倫 丁立 潘子軒 賈智樂(lè)

［關(guān)鍵詞］ArcPy腳本；小流域劃分；屬性賦值；溝道；微流域；河南省

［摘要］小流域劃分過(guò)程中河流溝道和微流域提取、面積控制、屬性批量賦值計(jì)算及質(zhì)量檢查都需要GIS工具的深度參與，現(xiàn)有GIS軟件無(wú)法直接滿足小流域劃分過(guò)程中批量化、流程化的需求?；贏rcPy設(shè)計(jì)并開(kāi)發(fā)了河南省小流域劃分輔助工具集，有針對(duì)性地實(shí)現(xiàn)流程化的溝道及微流域提取、小流域?qū)傩再x值計(jì)算、質(zhì)量檢查等功能。該工具集算法合理、操作便捷，在河南省小流域劃分工作中發(fā)揮了提質(zhì)增效的作用。

［中圖分類號(hào)］S157.1［文獻(xiàn)標(biāo)識(shí)碼］B［文章編號(hào)］1000－0941（2023）04－0019－05

以小流域?yàn)榛締卧乃亮魇ЬC合治理是我國(guó)在長(zhǎng)期水土保持工作中總結(jié)出來(lái)的一條寶貴經(jīng)驗(yàn)，是生態(tài)建設(shè)的重要手段。多年來(lái)河南省始終重視水土流失綜合治理工作，而獲取小流域信息是開(kāi)展流域治理工作的基礎(chǔ)[1]。以河南省小流域劃分為契機(jī)，逐步摸清全省小流域現(xiàn)狀，可為全省水土保持規(guī)劃制定、開(kāi)展綜合治理和信息化建設(shè)提供基礎(chǔ)數(shù)據(jù)支撐。小流域劃分的主要工作內(nèi)容包括溝道及微流域提取、小流域邊界劃分、屬性賦值計(jì)算、質(zhì)量檢查及數(shù)據(jù)入庫(kù)等[2]。在河南省16.7萬(wàn)km2范圍內(nèi)開(kāi)展小流域劃分工作，僅依靠ArcGIS自帶的地理處理工具存在重復(fù)作業(yè)、自動(dòng)化程度低等問(wèn)題[3]，且容易產(chǎn)生人為誤差，無(wú)法直接滿足實(shí)際工作需求。在ArcPy技術(shù)支持下開(kāi)發(fā)一個(gè)輔助河南省小流域劃分工作的流程化運(yùn)行、參數(shù)化選擇，集水文分析、屬性賦值計(jì)算、質(zhì)量檢查為一體的工具集尤為必要。

1ArcPy站點(diǎn)包

Python是一種功能強(qiáng)大的跨平臺(tái)語(yǔ)言，在代碼簡(jiǎn)潔性、通用性上有著顯著優(yōu)勢(shì)，通過(guò)Python可以輕松調(diào)用ArcGIS中封裝的算法來(lái)完成復(fù)雜的數(shù)據(jù)處理和分析任務(wù)[4]。ArcPy是以ArcGISscripting模塊為基礎(chǔ)并繼承了其功能而構(gòu)建的站點(diǎn)包，它繼承了基于Python語(yǔ)言編譯的數(shù)據(jù)轉(zhuǎn)換、數(shù)據(jù)管理、數(shù)據(jù)分析、自動(dòng)化制圖等功能，開(kāi)發(fā)人員可以應(yīng)用該模塊便捷地訪問(wèn)GP工具、類和模塊，創(chuàng)建自定義的空間數(shù)據(jù)操作工作流[5]；同時(shí)ArcPy還提供了多種地理處理函數(shù)，能夠更好地支持地理數(shù)據(jù)的分析計(jì)算。基于ArcPy的腳本工具開(kāi)發(fā)具有復(fù)雜度低、效率高、可移植性強(qiáng)等特點(diǎn)，使用該技術(shù)開(kāi)發(fā)河南省小流域劃分工具更加便捷高效。

2河南省小流域劃分主要內(nèi)容

河南省小流域劃分工作的內(nèi)容主要包括三部分：一是基于水文分析模型提取河南省范圍內(nèi)的溝道和微流域數(shù)據(jù)，為小流域精細(xì)劃分提供參考；二是通過(guò)人機(jī)交互方式進(jìn)行小流域劃分，并對(duì)每條小流域進(jìn)行屬性賦值和計(jì)算；三是對(duì)河南省小流域劃分成果進(jìn)行拓?fù)浜蛯傩詸z查。

2.1溝道及微流域提取

流域河網(wǎng)水系是描述區(qū)域地形地貌和水文特征的基本地理參數(shù)，是小流域劃分的重要依據(jù)[6]。在河南省小流域劃分工作中需要基于水文分析進(jìn)行溝道、微流域提取，主要流程是基于DEM數(shù)據(jù)提取作業(yè)單元的匯水關(guān)系，結(jié)合實(shí)際需求合理設(shè)置閾值提取溝道及相應(yīng)微流域。已有DEM數(shù)據(jù)精度低于預(yù)期，造成河網(wǎng)提取精度無(wú)法達(dá)到實(shí)際標(biāo)準(zhǔn)[7]，如何校正DEM精度從而提升水文分析成果質(zhì)量值得關(guān)注。閾值設(shè)定也應(yīng)該充分考慮河南省地形、地貌等自然地理狀況，并結(jié)合微流域面積最終確定。

2.2小流域劃分及屬性賦值

小流域劃分工作主要基于GIS軟件，參考吳鵬等[8]提出的河南省水土保持小流域劃分原則及流程，通過(guò)人機(jī)交互方式進(jìn)行劃分。單個(gè)小流域?qū)傩苑倍?，依?jù)《小流域劃分及編碼規(guī)范》（SL653—2013），小流域?qū)傩灾饕ㄐ×饔虼a、小流域名稱、小流域類型等24個(gè)指標(biāo)，具體屬性及相關(guān)要求見(jiàn)表1。屬性賦值過(guò)程中需要借助其他數(shù)據(jù)圖層和屬性信息參與計(jì)算（如一級(jí)流域邊界、縣域邊界、DEM數(shù)據(jù)、地形類別數(shù)據(jù)、村莊數(shù)據(jù)等），對(duì)數(shù)據(jù)質(zhì)量要求較高。個(gè)別屬性值計(jì)算過(guò)程復(fù)雜，涉及空間分析、邏輯判斷等多種處理流程，具有較強(qiáng)的專業(yè)性，各字段之間相互關(guān)聯(lián)，在計(jì)算過(guò)程中也應(yīng)有先后之分。

2.3質(zhì)量檢查

河南省小流域質(zhì)量檢查主要包括拓?fù)錂z查和屬性檢查兩部分。拓?fù)錂z查主要分析小流域矢量數(shù)據(jù)之間以及小流域與縣級(jí)行政區(qū)邊界之間的拓?fù)溴e(cuò)誤，在拓?fù)錂z查過(guò)程中檢查工具應(yīng)定義好所有類型的拓?fù)溴e(cuò)誤，同時(shí)將檢查結(jié)果直觀展示，并支持對(duì)作業(yè)單元內(nèi)部、作業(yè)單元之間以及全省小流域間的拓?fù)錂z查。屬性檢查對(duì)象包括幾何屬性以及非幾何屬性，其中幾何屬性主要包括面積和坐標(biāo)系；通過(guò)對(duì)幾何屬性檢查可以將作業(yè)單元內(nèi)的小流域面積總和與作業(yè)單元實(shí)際面積差值控制在一定容差范圍內(nèi)，這也是驗(yàn)證小流域拓?fù)溴e(cuò)誤的輔助手段，該容差需要根據(jù)實(shí)際情況設(shè)置；非幾何屬性則是面積及坐標(biāo)系之外的其他屬性，主要檢查內(nèi)容是字段值是否在指定值域內(nèi)或必填字段值是否為空（NONE）。

3基于ArcPy的功能設(shè)計(jì)與實(shí)現(xiàn)

3.1功能設(shè)計(jì)

河南省小流域劃分輔助工具集采用ArcPy腳本進(jìn)行開(kāi)發(fā)封裝。根據(jù)實(shí)際工作需要，輔助工具分模塊獨(dú)立進(jìn)行設(shè)計(jì)，主要包括水文分析模塊、屬性賦值計(jì)算模塊和數(shù)據(jù)質(zhì)檢模塊。參照水利部發(fā)布的《小流域劃分及編碼規(guī)范》（SL653—2013），以及河南省小流域綜合管理業(yè)務(wù)需要，河南省小流域劃分過(guò)程中采用縣一級(jí)的行政區(qū)劃作為作業(yè)單元，對(duì)每個(gè)縣（市、區(qū)）逐一進(jìn)行小流域劃分操作?？傮w功能框架見(jiàn)圖1，各模塊間相互獨(dú)立，在作業(yè)流程上存在邏輯關(guān)系。

3.2功能實(shí)現(xiàn)

3.2.1水文分析模塊

溝道和微流域提取的前提是確定水文分析作業(yè)單元，本模塊以縣級(jí)行政區(qū)劃作為小流域劃分作業(yè)單元。為了保證作業(yè)單元內(nèi)溝道、微流域數(shù)據(jù)的完整性，在水文分析中使用的DEM數(shù)據(jù)范圍通常大于縣區(qū)邊界，水文分析工具模塊的具體流程見(jiàn)圖2。DEM數(shù)據(jù)選用ALOSDEM，其水平及垂直精度均約為12.5m。

為了盡可能地保證水文分析結(jié)果的準(zhǔn)確性，使用AGREE算法[9]對(duì)DEM數(shù)據(jù)進(jìn)行修正，從而獲取到更高精度的DEM，該算法根據(jù)矢量河流數(shù)據(jù)降低與矢量線重疊的格網(wǎng)高程值[10]。水文分析模塊的主要流程包括：填洼、流向分析、匯流累積量計(jì)算，閾值提取溝道和河流鏈接數(shù)據(jù)，最后基于流向分析結(jié)果和河流鏈接數(shù)據(jù)獲取微流域分布。集水面積閾值是河網(wǎng)提取的關(guān)鍵參數(shù)，決定了提取河網(wǎng)的精度和與實(shí)際河網(wǎng)的匹配度[11]。根據(jù)河南省小流域劃分的實(shí)際需要，將微流域最小面積設(shè)定為0.1km2，提取作業(yè)單元內(nèi)集水區(qū)域大于0.1km2的微流域分水嶺及其對(duì)應(yīng)的溝道弧段、水流出口位置節(jié)點(diǎn)?？紤]到DEM精

度為12.5m，將微流域控制面積換算為集水面積閾值，約為640個(gè)柵格。為了保持各作業(yè)單元小流域劃分尺度的一致性，全省統(tǒng)一采用該閾值。圖3為某作業(yè)單元提取的溝道和微流域數(shù)據(jù)。

3.2.2屬性賦值計(jì)算模塊

計(jì)算小流域?qū)傩再x值時(shí)，除個(gè)別屬性可以通過(guò)簡(jiǎn)單賦值函數(shù)實(shí)現(xiàn)外，多數(shù)屬性字段涉及多類數(shù)據(jù)間的復(fù)雜計(jì)算，需要GIS空間分析、邏輯判斷、幾何計(jì)算等算法深度參與。由于屬性字段是統(tǒng)一且固定不變的，因此只需要在工具集面板中進(jìn)行參數(shù)化的圖層選擇即可。河南省小流域?qū)傩再x值計(jì)算模塊主要流程結(jié)構(gòu)見(jiàn)圖4，通過(guò)人機(jī)交互手段劃分的小流域圖層為屬性賦值計(jì)算的核心數(shù)據(jù)，除此之外還需要屬性字段表、地形類型、行政區(qū)劃、一級(jí)流域、DEM、溝道及微流域、村莊點(diǎn)、主要河流、包含河流名稱及河流代碼的流域面等輔助圖層（圖中深灰色底紋為需要計(jì)算的小流域字段，其對(duì)應(yīng)箭頭的上游方向?yàn)橛?jì)算該字段所需的圖層數(shù)據(jù)，圖中為淺色底紋圖層）。下文根據(jù)字段賦值在流程上的先后順序詳細(xì)介紹各字段的計(jì)算過(guò)程及關(guān)鍵函數(shù)。

（1）小流域?qū)傩宰侄翁砑?。通過(guò)調(diào)用外部預(yù)定義屬性表將所需字段追加至小流域圖層，除字段名稱外還包括字段別名、字段類型、長(zhǎng)度、數(shù)值位數(shù)、小數(shù)位數(shù)等信息，在程序中通過(guò)逐單元格遍歷的形式將字段信息追加至小流域?qū)傩员怼?/p>

（2）劃分人、劃分時(shí)間、縣（市、區(qū)）、縣（市、區(qū)）代碼、地（市）。在同一作業(yè)單元內(nèi)的上述5個(gè)字段屬性值通常是相同的，通過(guò)字段批量計(jì)算函數(shù)arcpy.management.CalculateField（）進(jìn)行賦值。

（3）面積最東經(jīng)度、最西經(jīng)度、最南緯度、最北緯度（四至地理坐標(biāo)）。通過(guò)arcpy.CalculateGeometryAttributes_management（）方法進(jìn)行幾何空間屬性計(jì)算，字段值采用與小流域一致的CGCS2000坐標(biāo)，考慮ArcPy計(jì)算的四至地理坐標(biāo)結(jié)果不支持度分秒格式，需要定義函數(shù)進(jìn)行相應(yīng)格式轉(zhuǎn)換。

（4）一級(jí)流域、地形類型、河流名稱、河流代碼。參照屬性賦值計(jì)算邏輯，首先使用arcpy.analysis.Intersect（）函數(shù)獲取相交圖層信息，若單個(gè)小流域被分割為多塊，則使用arcpy.analysis.Statistics（Intersect，None，["Shape_Area"，"MAX"]，LinkField）函數(shù)得到其中面積最大的要素字段信息，并通過(guò)字段鏈接形式進(jìn)行賦值。

（5）小流域類型、亞單元。不同類型的小流域和亞單元在小流域邊界與河流溝道的空間關(guān)系特征上各有差異，其中完整型小流域與溝道只有一個(gè)交點(diǎn)，區(qū)間型小流域存在兩個(gè)或多個(gè)交點(diǎn)，地形類型為丘陵或山地的坡面型小流域部分邊界通常與河流共線。亞單元小流域往往在和縣（市、區(qū)）的公共邊界上同河流溝道存在交點(diǎn)?；谶@些特征設(shè)計(jì)小流域類型及亞單元屬性賦值功能。亞單元賦值的主要Python語(yǔ)句如下：

XLY_XJ_RV_insect=arcpy.analysis.Intersect（[[Layer，1]，[XJ_Layer，1]，[GD_Layer，1]]，out_insect_in，"ONLY_FID"，""，"POINT"）#小流域、縣界、溝道相交獲取交點(diǎn)

insect_xly_id_list.append（FID）#將所有存在交點(diǎn)的小流域FID值存入類別

witharcpy.da.UpdateCursor（Layer，["OID@"，XLY]）ascursor：#遍歷目標(biāo)要素并賦值

row[1]=1

（6）小流域代碼批量賦值。小流域代碼需要滿足全省唯一性的條件，除通過(guò)縣（市、區(qū)）代碼和河流代碼進(jìn)行控制外，為避免在作業(yè)單元內(nèi)出現(xiàn)重復(fù)值，在同一縣區(qū)內(nèi)為小流域賦予唯一標(biāo)識(shí)，即ID值。小流域代碼根據(jù)規(guī)范要求采用英文大寫字母（舍棄I、O、Z）和數(shù)字的混合碼，共17位，分別由全國(guó)河流代碼、小流域順序碼、亞單元標(biāo)識(shí)碼和縣級(jí)行政區(qū)區(qū)劃代碼構(gòu)成，主要Python語(yǔ)句如下：

arcpy.management.CalculateField（Layer，"XLYDM"，"str（！HLDM！）[1：7]+str（！ID?。?zfill（4）+str（！YDY?。?str（！XQDM?。?，"PYTHON3"，""，"TEXT"）

（7）平均海拔、平均坡度?；谧鳂I(yè)單元DEM數(shù)據(jù)計(jì)算小流域平均海拔、平均坡度，調(diào)用arcpy.ddd.Slope（）函數(shù)，通過(guò)三維分析提取坡度，使用arcpy.sa.ZonalStatisticsAsTable（）函數(shù)計(jì)算每個(gè)小流域的平均海拔和坡度值，最后通過(guò)字段鏈接的方式賦值至相應(yīng)字段。

（8）溝壑密度、最長(zhǎng)匯流路徑。溝壑密度和最長(zhǎng)匯流路徑以小流域?yàn)閱挝贿M(jìn)行統(tǒng)計(jì)，利用區(qū)域統(tǒng)計(jì)分析工具實(shí)現(xiàn)，其中小流域溝壑密度為小流域內(nèi)溝道長(zhǎng)度值之和與小流域面積之比，最長(zhǎng)匯流路徑長(zhǎng)度為小流域內(nèi)最長(zhǎng)的河流長(zhǎng)度，對(duì)于無(wú)河流數(shù)據(jù)的小流域，自動(dòng)遍歷并將其最長(zhǎng)匯流路徑賦值為“0”。溝壑密度計(jì)算過(guò)程：使用小流域圖層對(duì)溝道數(shù)據(jù)進(jìn)行裁剪，通過(guò)arcpy.analysis.Identity（"Goudao_Clip"，Layer，"Gouhe_Identify"，"ALL"，None，"NO_RELATIONSHIPS"）函數(shù)進(jìn)行標(biāo)識(shí)以建立與小流域的空間關(guān)聯(lián)，利用arcpy.analysis.Statistics（"Gouhe_Identify"，"Zonal_GHCD_SUM_Statis_Tb"，"GHCDSUM"，"XLYDM"）函數(shù)計(jì)算小流域內(nèi)的溝道長(zhǎng)度總和，并將其鏈接至小流域圖層，最后基于arcpy.management.CalculateField（Layer，"GHMD"，"！SUM_GHCD！/！MJ！"，"PYTHON3"，"TEXT"）函數(shù)批量獲取作業(yè)單元小流域溝壑密度屬性。

（9）小流域名稱。根據(jù)小流域命名原則，小流域名稱應(yīng)簡(jiǎn)明確切、易于辨識(shí)，可采用當(dāng)?shù)販系?、村莊等標(biāo)志性地物命名，本研究以小流域幾何中心附近村莊名稱命名。通過(guò)arcpy.FeatureToPoint_management（）函數(shù)將小流域面要素轉(zhuǎn)換為包含于面內(nèi)的點(diǎn)要素，使用arcpy.Near_analysis（）函數(shù)捕捉鄰近村莊點(diǎn)作為小流域名稱，在捕捉的過(guò)程中使用arcpy.management.DeleteIdentical（）函數(shù)移除村名重復(fù)點(diǎn)，以避免在同一縣區(qū)范圍內(nèi)存在同名小流域。

為了提高字段賦值計(jì)算工具的操作便捷性，將開(kāi)發(fā)的功能模塊打包為工具箱，實(shí)現(xiàn)在同一個(gè)地理處理工具面板中統(tǒng)一批處理，即選擇所有參與計(jì)算的圖層數(shù)據(jù)，按照相應(yīng)計(jì)算邏輯和字段賦值順序一次性完成一個(gè)作業(yè)單元內(nèi)小流域的所有屬性批量賦值。

3.2.3數(shù)據(jù)質(zhì)檢模塊

質(zhì)量檢查是保障小流域劃分成果準(zhǔn)確性、可靠性的重要手段。在程序中檢測(cè)到錯(cuò)誤后以彈窗形式進(jìn)行提示，便于小流域劃分人員快速定位錯(cuò)誤進(jìn)行修改。

河南省小流域?qū)傩詸z查主要包括以下內(nèi)容：

（1）屬性完整性檢查。將小流域圖層中的所有屬性字段存入列表，通過(guò)遍歷與外部預(yù)定義的屬性表信息逐一對(duì)照，對(duì)缺少的字段進(jìn)行提示，并刪除多余字段。

（2）坐標(biāo)系正確性檢查。獲取作業(yè)單元小流域的元數(shù)據(jù)信息，判斷坐標(biāo)系是否為CGCS2000國(guó)家大地坐標(biāo)系，否則進(jìn)行提示。

（3）字段完整性檢查。檢查所有小流域要素的必填字段屬性值是否為空或“NONE”值，若為空或“NONE”'值則警告提示，并列出相應(yīng)小流域編碼及字段名。

（4）面積誤差控制。匯總統(tǒng)計(jì)同一縣（市、區(qū)）內(nèi)的小流域總和面積并與該縣（市、區(qū)）總面積進(jìn)行對(duì)比，絕對(duì)誤差值應(yīng)控制在0.01km2內(nèi)，否則進(jìn)行提示。

在河南省小流域劃分空間拓?fù)鋽?shù)據(jù)質(zhì)量檢查時(shí)，調(diào)用arcpy.management.AddRuleToTopology（）函數(shù)為拓?fù)錂z查定義規(guī)則。針對(duì)以下3種常見(jiàn)的拓?fù)溴e(cuò)誤類型，分別設(shè)置“MustNotOverlap（Area）”“MustNotHaveGaps（Area）”“BoundaryMustBeCoveredByBoundaryOf（Area-Area）”為對(duì)應(yīng)的拓?fù)湟?guī)則參數(shù)。

面拓?fù)溴e(cuò)誤為小流域邊界相互重疊，如圖REF_Ref85896294h*MERGEFORMAT5中小流域A與小流域B間的區(qū)域所示。線拓?fù)溴e(cuò)誤一為小流域邊界之間存在空隙，如圖５中小流域B和小流域C間的區(qū)域所示；二為縣（市、區(qū)）行政邊界沒(méi)有被小流域邊界覆蓋，如圖５中小流域A與縣區(qū)行政邊界間的區(qū)域所示。

4總結(jié)與討論

目前關(guān)于小流域劃分技術(shù)的研究較多，但多是基于水文分析模型的小流域自動(dòng)劃分技術(shù)探討，很少有對(duì)后續(xù)屬性賦值流程及結(jié)果的關(guān)注，但這也是小流域劃分成果的重要組成部分。同時(shí)，現(xiàn)有研究對(duì)小流域劃分效率與質(zhì)量方面也鮮有提及。本研究重點(diǎn)在于：一是將縣（市、區(qū)）作為小流域劃分的最小作業(yè)單元，進(jìn)一步提升了小流域劃分成果的精細(xì)化程度；二是更加關(guān)注小流域劃分成果中的屬性賦值計(jì)算和質(zhì)量檢查過(guò)程，將小流域的屬性成果放在了與空間劃分成果同等重要的位置；三是更加重視小流域劃分的效率和質(zhì)量問(wèn)題，并在梳理每個(gè)屬性賦值計(jì)算的流程基礎(chǔ)上定制開(kāi)發(fā)了對(duì)應(yīng)的自動(dòng)批處理工具，在降低小流域劃分后續(xù)工作復(fù)雜度的基礎(chǔ)上又極大地提升了效率和質(zhì)量。

基于ArcPy站點(diǎn)包構(gòu)建的小流域劃分輔助工具集在河南省小流域劃分項(xiàng)目中實(shí)現(xiàn)了流程化、專業(yè)化、高效化的批處理，通過(guò)水文分析提取了河南省溝道、微流域數(shù)據(jù)，分別計(jì)算了河南省5576個(gè)小流域的24個(gè)屬性字段，并對(duì)小流域成果進(jìn)行了質(zhì)量檢查等。該工具集具有以下特點(diǎn)：

（1）運(yùn)行效率高。結(jié)合實(shí)際工作需求對(duì)算法邏輯進(jìn)行優(yōu)化，減少了大量重復(fù)操作，能夠?qū)崿F(xiàn)數(shù)據(jù)的批量自動(dòng)化處理，與傳統(tǒng)方式相比作業(yè)效率明顯提高。

（2）運(yùn)行穩(wěn)定可靠。字段通過(guò)鏈接形式添加在屬性中，處理完成后自動(dòng)取消字段關(guān)聯(lián)，不會(huì)產(chǎn)生冗余數(shù)據(jù)；字段賦值過(guò)程嚴(yán)格按照規(guī)范流程執(zhí)行，設(shè)置異常捕捉機(jī)制，產(chǎn)生錯(cuò)誤后立即停止運(yùn)行避免產(chǎn)生錯(cuò)誤數(shù)據(jù)；空間分析的結(jié)果為新圖層，不會(huì)更改輸入數(shù)據(jù)的空間屬性。

（3）操作簡(jiǎn)單。工具集算法合理，操作便捷，只需輸入相關(guān)圖層參數(shù)即可由后臺(tái)完全自動(dòng)化處理。河南省小流域劃分輔助工具集的研發(fā)切實(shí)解決了生產(chǎn)實(shí)踐中的問(wèn)題，提高了工作效率，是小流域劃分工作按時(shí)高質(zhì)量完成的保障。利用ArcPy工具進(jìn)行高效數(shù)據(jù)分析處理是一種科學(xué)有效的方法，這種開(kāi)發(fā)思路在眾多數(shù)據(jù)生產(chǎn)項(xiàng)目中均具有借鑒意義，本研究成果具備一定的推廣價(jià)值和可移植性。

[參考文獻(xiàn)]

[1]俞春娜，馮存均.面積約束的小流域提取算法研究[J].測(cè)繪科學(xué)，2022，47（12）：216-222.

[2]張景揚(yáng)，盧遠(yuǎn)，劉斌濤.基于ArcHydroTools的廣西小流域劃分研究[J].水利水電技術(shù)，2017，48（9）：65-71.

[3]余詠勝，彭艷麗，尹言軍，等.基于Arcpy的影像地圖自動(dòng)處理技術(shù)研究[J].測(cè)繪通報(bào)，2015（3）：82-85.

[4]呂德昊，宋偉東，董山，等.ArcGIS中應(yīng)用Python腳本實(shí)現(xiàn)地理數(shù)據(jù)的批處理[J].測(cè)繪與空間地理信息，2021，44（2）：111-112，117.

[5]馮里濤，陳艷艷，李兆豐，等.基于Geoprocessing的基礎(chǔ)地理信息數(shù)據(jù)庫(kù)質(zhì)量工具集的開(kāi)發(fā)[J].測(cè)繪與空間地理信息，2013，36（9）：102-105.

[6]王奎，林孝松.基于GIS的小流域特征信息提取[J].重慶理工大學(xué)學(xué)報(bào)（自然科學(xué)），2015，29（9）：105-110.

[7]趙遠(yuǎn)洋，丁永生，孫丹，等.基于Arcpy編程修正的數(shù)字高程模型小流域水文分析[J].大連海事大學(xué)學(xué)報(bào)，2016，42（4）：55-60.

[8]吳鵬，韓錦琳.河南省水土保持小流域劃分原則與流程[J].中國(guó)水土保持，2022（5）：29-31.

[9]FERDIHellweger.AGREE-DEMSurfaceReconditioningSystem[D].Austin：UniversityofTexasatAustin，1997：17.

[10]鄭倩，史海濱，李仙岳，等.基于AGREE算法與BURN-IN算法的平原灌區(qū)DEM河網(wǎng)提取問(wèn)題剖析[J].灌溉排水學(xué)報(bào)，2019，38（5）：122-128.

[11]李照會(huì)，郭良，劉榮華，等.基于DEM數(shù)字河網(wǎng)提取時(shí)集水面積閾值與河源密度關(guān)系的研究[J].地球信息科學(xué)學(xué)報(bào)，2018，20（9）：1244-1251.

［作者簡(jiǎn)介］王松倫（1981—），男，河南扶溝人，高級(jí)工程師，碩士，主要研究方向?yàn)樗畔⒒煌ㄐ抛髡吲俗榆帲?996—），男，河南輝縣人，助理工程師，碩士，主要研究方向?yàn)樗b感監(jiān)管。［收稿日期］2022-12-01

（責(zé)任編輯楊傲秋）