王露 盧駿

摘 要：【目的】梳理現(xiàn)有土地利用現(xiàn)狀及規(guī)劃分類(lèi)統(tǒng)計(jì)相關(guān)研究成果，從而優(yōu)化和完美傳統(tǒng)人工分類(lèi)方式?！痉椒ā拷Y(jié)合土地利用現(xiàn)狀及規(guī)劃分類(lèi)要求，以python通用編程語(yǔ)言為基礎(chǔ)，基于arcpy站點(diǎn)包及openpyxl模塊設(shè)計(jì)一套腳本工具，引入“四舍六入五看齊、奇進(jìn)偶不進(jìn)”的進(jìn)位規(guī)則，確保成果精度?！窘Y(jié)果】解決excel數(shù)據(jù)統(tǒng)計(jì)結(jié)果不準(zhǔn)確問(wèn)題，實(shí)現(xiàn)建設(shè)項(xiàng)目土地利用現(xiàn)狀及規(guī)劃分類(lèi)自動(dòng)統(tǒng)計(jì)，通過(guò)試驗(yàn)對(duì)比驗(yàn)證工具的可靠性和穩(wěn)定性?！窘Y(jié)論】提高用地預(yù)審與規(guī)劃選址、土地利用規(guī)劃修改、用地報(bào)批等的工作效率。

關(guān)鍵詞：土地利用；腳本工具；奇進(jìn)偶舍；用地預(yù)審

中圖分類(lèi)號(hào)：P963? ? ?文獻(xiàn)標(biāo)志碼：A? ? ? 文章編號(hào)：1003-5168（2023）09-0023-04

Abstract：[Purposes] To sort out the existing problems in the statistical research on land use status and planning classification， and improve the cumbersome and error-prone status of traditional manual classification methods. [Methods] Based on the current situation of land use and the requirements of planning and classification， a set of scripting tools was designed based on the Python general programming language， the Arcpy site package， and the openpyxl module. The rounding rule of "round to six， keep to five， and odd to even" was introduced to ensure the accuracy of the results. [Findings] The problem of unbalanced statistical results from Excel data was solved， and automatic statistics of land use status and planning classification for construction projects were realized. The reliability and stability of the tool were verified through experimental comparison. [Conclusions] The tool effectively promote the efficiency of land use preliminary review and planning site selection， land use planning modification， and land use approval.

Keywords： land use; script tools; banker's rounding; preliminary review of land use

0 引言

用地預(yù)審與規(guī)劃選址、土地利用規(guī)劃修改、用地報(bào)批等前期工作是國(guó)土空間用途管制十分重要的環(huán)節(jié)，是落實(shí)國(guó)土空間規(guī)劃的主要調(diào)控措施，也是建設(shè)項(xiàng)目立項(xiàng)階段的重要工作之一［1］。土地規(guī)劃是國(guó)土資源管理的重要依據(jù)，對(duì)發(fā)揮國(guó)土宏觀調(diào)控作用和國(guó)土資源合理利用等具有重要意義，而土地利用現(xiàn)狀和規(guī)劃分類(lèi)統(tǒng)計(jì)數(shù)據(jù)又是國(guó)土資源管理的基礎(chǔ)資料。

目前，對(duì)土地利用現(xiàn)狀及規(guī)劃分類(lèi)自動(dòng)化統(tǒng)計(jì)的研究相對(duì)較少。傳統(tǒng)的人工統(tǒng)計(jì)法是先在ArcGIS中對(duì)數(shù)據(jù)進(jìn)行處理，再將處理結(jié)果圖層導(dǎo)入到excel［2］或ACCESS表格［3］中，利用excel公式或SQL函數(shù)完成相關(guān)統(tǒng)計(jì)工作。數(shù)據(jù)量較小時(shí)，該方法具有一定優(yōu)勢(shì)，但隨著項(xiàng)目地塊面積的增加，地類(lèi)數(shù)據(jù)量和復(fù)雜度增加，該方法會(huì)出現(xiàn)效率低、易出錯(cuò)等弊端。此外，excel統(tǒng)計(jì)存在一定的局限性，容易出現(xiàn)分項(xiàng)之和不等于總數(shù)的數(shù)據(jù)不準(zhǔn)確問(wèn)題［4］。有學(xué)者利用信息化技術(shù)來(lái)實(shí)現(xiàn)土地分類(lèi)的自動(dòng)統(tǒng)計(jì)，但主要面向勘測(cè)定界階段，如鄭浩［5］、高玉榮等［6］、于國(guó)強(qiáng)［7］等通過(guò)信息化技術(shù)，實(shí)現(xiàn)對(duì)勘測(cè)定界階段土地現(xiàn)狀分類(lèi)的快速統(tǒng)計(jì)。鄭浩［5］、高玉榮等［6］基于ArcGIS Engine進(jìn)行系統(tǒng)開(kāi)發(fā)，于國(guó)強(qiáng)［7］則基于AutoCAD平臺(tái)和VBA語(yǔ)言對(duì)系統(tǒng)進(jìn)行開(kāi)發(fā)，分別從手工統(tǒng)計(jì)對(duì)比和實(shí)際應(yīng)用層來(lái)驗(yàn)證系統(tǒng)的高效性、正確性。以上研究均未考慮預(yù)審、規(guī)劃選址、土地利用規(guī)劃修改階段的土地分類(lèi)統(tǒng)計(jì)要求，無(wú)法實(shí)現(xiàn)對(duì)土地規(guī)劃分類(lèi)的自動(dòng)統(tǒng)計(jì)。其中，鄭浩［5］設(shè)計(jì)的系統(tǒng)還存在軟件環(huán)境要求高、可移植性差等問(wèn)題。需要著重從用地預(yù)審與規(guī)劃選址、土地利用規(guī)劃修改等環(huán)節(jié)的土地地類(lèi)統(tǒng)計(jì)需求出發(fā)，開(kāi)發(fā)一套ArcGIS通用腳本工具，同步實(shí)現(xiàn)土地利用現(xiàn)狀及規(guī)劃分類(lèi)統(tǒng)計(jì)結(jié)果的自動(dòng)生成。該工具可移植性強(qiáng)、操作簡(jiǎn)便，生成的成果精度較高，可大大提高生產(chǎn)工作效率。

1 自動(dòng)化分類(lèi)統(tǒng)計(jì)工具實(shí)現(xiàn)

1.1 研究方法和技術(shù)路線

ArcGIS軟件提供一系列地理空間數(shù)據(jù)自動(dòng)化處理方法，包括簡(jiǎn)單批處理工具、模型構(gòu)建器、arcpy腳本語(yǔ)言［8］。研究選用arcpy腳本語(yǔ)言來(lái)構(gòu)建具備可視化對(duì)話(huà)框的定制腳本工具，對(duì)土地利用現(xiàn)狀和規(guī)劃分類(lèi)統(tǒng)計(jì)過(guò)程實(shí)現(xiàn)自動(dòng)化出表。

分類(lèi)統(tǒng)計(jì)工具的實(shí)現(xiàn)步驟如下。①按照當(dāng)前國(guó)土資源管理工作中土地利用現(xiàn)狀及規(guī)劃地類(lèi)統(tǒng)計(jì)要求，建立現(xiàn)狀地類(lèi)與統(tǒng)計(jì)地類(lèi)、規(guī)劃地類(lèi)與統(tǒng)計(jì)地類(lèi)對(duì)應(yīng)關(guān)系，完成地類(lèi)轉(zhuǎn)換；②將項(xiàng)目范圍線與行政區(qū)、現(xiàn)狀地類(lèi)、規(guī)劃地類(lèi)、永久基本農(nóng)田等數(shù)據(jù)進(jìn)行相交分析，得到項(xiàng)目涉及到行政區(qū)、現(xiàn)狀地類(lèi)、規(guī)劃地類(lèi)、永久基本農(nóng)田等信息；③對(duì)上一步驟獲得的相交成果數(shù)據(jù)，利用cursor游標(biāo)函數(shù)分行政區(qū)、分地類(lèi)讀取字段面積，將結(jié)果存儲(chǔ)到嵌套字典變量中；④利用openpyxl模塊來(lái)完成統(tǒng)計(jì)表格的讀寫(xiě)與保存；⑤利用arcpy.GetParameterAsText（）函數(shù)構(gòu)建輸入?yún)?shù)，完成腳本工具的創(chuàng)建。

1.2 土地利用現(xiàn)狀與規(guī)劃分類(lèi)體系

年度土地利用變更調(diào)查是基于土地利用現(xiàn)狀的調(diào)查，為保持變更調(diào)查成果的現(xiàn)勢(shì)性和準(zhǔn)確性而每年進(jìn)行的成果更新過(guò)程［9］，也是反映當(dāng)前土地利用現(xiàn)狀的主要基礎(chǔ)數(shù)據(jù)。

土地利用現(xiàn)狀按《土地利用現(xiàn)狀分類(lèi)》（GB/T 21010—2017）進(jìn)行分類(lèi)。土地利用總體規(guī)劃是各級(jí)行政區(qū)對(duì)土地資源進(jìn)行統(tǒng)籌布局和規(guī)劃安排而形成的基礎(chǔ)數(shù)據(jù)，其分類(lèi)體系是根據(jù)規(guī)劃管理需要，在土地利用現(xiàn)狀分類(lèi)的基礎(chǔ)上進(jìn)行歸并或細(xì)分轉(zhuǎn)換形成的。統(tǒng)計(jì)分類(lèi)是在此基礎(chǔ)上進(jìn)一步歸并形成的。因此，在進(jìn)行其他操作前，要將土地利用現(xiàn)狀地類(lèi)和規(guī)劃地類(lèi)按要求轉(zhuǎn)換為統(tǒng)計(jì)地類(lèi)。本研究將地類(lèi)轉(zhuǎn)換對(duì)應(yīng)關(guān)系存儲(chǔ)于外置excel表格中，在使用工具時(shí)，通過(guò)openpyxl模塊調(diào)用excel表格，以字典變量形式存儲(chǔ)地類(lèi)轉(zhuǎn)換關(guān)系。該方式可有效增加工具的靈活性，降低后續(xù)工具與“三調(diào)”現(xiàn)狀及國(guó)土空間規(guī)劃分類(lèi)統(tǒng)計(jì)的銜接難度。

1.3 地類(lèi)面積統(tǒng)計(jì)與存儲(chǔ)

本研究需要解決的核心問(wèn)題是各地類(lèi)面積的統(tǒng)計(jì)與存儲(chǔ)。主要統(tǒng)計(jì)以下三類(lèi)面積數(shù)據(jù)：現(xiàn)狀地類(lèi)面積、規(guī)劃地類(lèi)面積、永久基本農(nóng)田面積等。其中，現(xiàn)狀地類(lèi)和規(guī)劃地類(lèi)要分行政區(qū)、分地類(lèi)進(jìn)行面積統(tǒng)計(jì)，永久基本農(nóng)田面積要分行政區(qū)進(jìn)行面積統(tǒng)計(jì)。因此，統(tǒng)計(jì)面積時(shí)會(huì)用到行政區(qū)信息、現(xiàn)狀地類(lèi)信息、規(guī)劃地類(lèi)信息、永久基本農(nóng)田信息等。為獲取所需的數(shù)據(jù)信息，首先，將建設(shè)項(xiàng)目范圍線與行政區(qū)數(shù)據(jù)、土地利用現(xiàn)狀數(shù)據(jù)、土地利用總體規(guī)劃數(shù)據(jù)、永久基本農(nóng)田保護(hù)圖斑數(shù)據(jù)進(jìn)行疊加分析，按照要求存放分析成果，方便后續(xù)調(diào)用。其次，分別創(chuàng)建現(xiàn)狀、規(guī)劃、永久基本農(nóng)田三個(gè)字典變量，利用cursor游標(biāo)函數(shù)與循環(huán)語(yǔ)句遍歷相應(yīng)圖層數(shù)據(jù)對(duì)象，提取面積字段值。再次，按行政區(qū)、地類(lèi)、面積的對(duì)應(yīng)關(guān)系分別存儲(chǔ)到所創(chuàng)建的字典變量中。其中，規(guī)劃和現(xiàn)狀直接讀取數(shù)據(jù)圖層自帶的“shape area”面積字段，即統(tǒng)計(jì)地類(lèi)的毛面積。永久基本農(nóng)田是利用自然資源部下發(fā)的凈面積計(jì)算方法，得到永久基本農(nóng)田的凈面積，對(duì)凈面積進(jìn)行判斷，僅在凈面積大于1時(shí)進(jìn)行累加，分行政區(qū)統(tǒng)計(jì)面積之和。最后，該過(guò)程可返回{行政區(qū)：{地類(lèi)：面積}的嵌套字典。

此外，在使用excel進(jìn)行數(shù)據(jù)統(tǒng)計(jì)時(shí)，往往會(huì)出現(xiàn)分項(xiàng)之和不等于總數(shù)的情況，這是因?yàn)楹嫌?jì)數(shù)是分項(xiàng)實(shí)際數(shù)據(jù)相加得到的，不是減少位數(shù)后的分項(xiàng)數(shù)相加所得，導(dǎo)致分項(xiàng)數(shù)與合計(jì)數(shù)之間并無(wú)直接的數(shù)據(jù)關(guān)系［4］。因此，在進(jìn)行數(shù)據(jù)計(jì)算前，最好對(duì)數(shù)據(jù)成果保留精度位數(shù)?！八纳嵛迦搿笔且环N常用的精確度計(jì)數(shù)保留法，但因其規(guī)則存在“舍的少、進(jìn)的多”，因而結(jié)果往往偏大。為了能更好地滿(mǎn)足地類(lèi)面積統(tǒng)計(jì)高精度的要求，引入測(cè)量中“四舍六入五看齊、奇進(jìn)偶不進(jìn)”的進(jìn)位規(guī)則，即當(dāng)小數(shù)點(diǎn)后第五位為1到4的數(shù)字時(shí)不進(jìn)位，為6到9的數(shù)字時(shí)進(jìn)位，為5時(shí)要看小數(shù)點(diǎn)后第四位，如果是奇數(shù)，則進(jìn)位，否則不進(jìn)位。從統(tǒng)計(jì)學(xué)角度來(lái)看，該方法可避免出現(xiàn)“四舍五入”進(jìn)位制結(jié)果偏大的現(xiàn)象，將結(jié)果受到舍入誤差的影響降到最低，提高統(tǒng)計(jì)結(jié)果精度。

1.4 統(tǒng)計(jì)表格自動(dòng)化輸出

本研究設(shè)計(jì)的工具最終要實(shí)現(xiàn)對(duì)土地利用現(xiàn)狀及規(guī)劃地類(lèi)統(tǒng)計(jì)結(jié)果表的自動(dòng)生成。openpyxl模塊可通過(guò)“工作簿workbook-工作表sheet-單元格cell”的模式對(duì).xlsx文件進(jìn)行表格的創(chuàng)建、讀取、寫(xiě)入與修改，能同步調(diào)整表格樣式。首先，要通過(guò)openpyxl模塊的create_sheet（）方法創(chuàng)建一張現(xiàn)狀及規(guī)劃表；其次，通過(guò)merge_cells（）、Font（）等方法完成對(duì)單元格的合并、表頭制作、樣式設(shè)計(jì)等；再次，分行政區(qū)和地類(lèi)讀取嵌套字典中的面積數(shù)據(jù)，填寫(xiě)到對(duì)應(yīng)的單元格中；最后，通過(guò)save（）方法來(lái)完成對(duì)表格的保存。

由于不同項(xiàng)目涉及的行政區(qū)數(shù)量有所不同，導(dǎo)致統(tǒng)計(jì)結(jié)果表中的列數(shù)也不相同，固定列數(shù)的表格無(wú)法滿(mǎn)足統(tǒng)計(jì)需求，因此對(duì)不同行政區(qū)數(shù)量要靈活設(shè)計(jì)統(tǒng)計(jì)表。經(jīng)初步分析，每增加一個(gè)行政區(qū)，就要增加一列數(shù)據(jù)，套用行政區(qū)數(shù)量和列數(shù)之間的關(guān)系，可采用循環(huán)方式，以行政區(qū)數(shù)量為循環(huán)次數(shù)，按行政區(qū)逐項(xiàng)讀取相關(guān)地類(lèi)面積，填寫(xiě)到對(duì)應(yīng)單元格中。

1.5 定制腳本工具

ArcGIS為用戶(hù)提供一套可用于空間地理數(shù)據(jù)處理分析的工具集［10］，可滿(mǎn)足絕大部分空間數(shù)據(jù)的處理需求。當(dāng)現(xiàn)有工具集無(wú)法滿(mǎn)足需求時(shí)，可定制工具，以工具箱可交互界面的形式來(lái)滿(mǎn)足用戶(hù)個(gè)性化需求。本研究使用ArcGIS腳本制作工具來(lái)完成代碼包裝，在工具箱目錄中新建工具箱，并新建腳本工具，將python代碼文件導(dǎo)入腳本工具中，設(shè)置好參數(shù)。工具界面的參數(shù)通過(guò)arcpy.GetParameterAsText（）函數(shù)設(shè)置的參數(shù)與代碼一一對(duì)應(yīng)起來(lái)［11］，完成腳本工具的包裝。

2 應(yīng)用實(shí)踐

所需的基礎(chǔ)數(shù)據(jù)包括實(shí)踐區(qū)域內(nèi)2018年土地利用現(xiàn)狀、土地利用總體規(guī)劃、永久基本農(nóng)田、行政區(qū)、建設(shè)項(xiàng)目用地范圍線。其中，建設(shè)項(xiàng)目用地范圍線為試驗(yàn)區(qū)域內(nèi)自行創(chuàng)建的模擬用地范圍［5］，用地紅線總面積為3 322.360 4 hm2。建設(shè)項(xiàng)目模擬用地范圍如圖1所示。

工具運(yùn)算時(shí)間一般在30 s內(nèi)，分析出的結(jié)果未見(jiàn)總數(shù)與分項(xiàng)數(shù)不等的問(wèn)題。同時(shí)，采用傳統(tǒng)的人工統(tǒng)計(jì)方式進(jìn)行驗(yàn)算，驗(yàn)算用時(shí)為40～60 min，得到的數(shù)據(jù)基本一致，工具統(tǒng)計(jì)結(jié)果未見(jiàn)數(shù)據(jù)不平衡問(wèn)題。人工統(tǒng)計(jì)結(jié)果不可避免地存在總數(shù)與分項(xiàng)數(shù)不等的數(shù)據(jù)不平衡問(wèn)題，如行政區(qū)1中耕地、園地、林地、草地、其他農(nóng)用地之和與總數(shù)會(huì)存在0.000 1的差別，要額外進(jìn)行數(shù)據(jù)平衡處理。由此證明工具統(tǒng)計(jì)效率較高，且統(tǒng)計(jì)結(jié)果穩(wěn)定性和可靠性更強(qiáng)。

3 結(jié)論

本研究通過(guò)arcpy站點(diǎn)包設(shè)計(jì)出一套土地利用現(xiàn)狀及規(guī)劃分類(lèi)統(tǒng)計(jì)的腳本工具，該腳本工具可完成土地利用現(xiàn)狀及規(guī)劃分類(lèi)統(tǒng)計(jì)表的自動(dòng)生成，提高工作效率。將分類(lèi)表對(duì)應(yīng)存儲(chǔ)在外置excel表格中，可按需對(duì)分類(lèi)表進(jìn)行修改，滿(mǎn)足后續(xù)“三調(diào)”及國(guó)土空間規(guī)劃分類(lèi)統(tǒng)計(jì)需求。此外，在面積統(tǒng)計(jì)時(shí)考慮到“四舍五入”精度保留法存在的缺陷及excel表格統(tǒng)計(jì)結(jié)果存在數(shù)據(jù)不準(zhǔn)確的問(wèn)題，融入“奇進(jìn)偶不進(jìn)”的“四舍六入”精度保留規(guī)則，能有效提高成果精度，解決excel數(shù)據(jù)統(tǒng)計(jì)結(jié)果不平衡的問(wèn)題。通過(guò)試驗(yàn)對(duì)比，驗(yàn)證工具的可靠性和穩(wěn)定性，為提升用地預(yù)審與規(guī)劃選址、土地利用規(guī)劃修改、用地報(bào)批等工作的效率提供技術(shù)支撐。

該工具僅統(tǒng)計(jì)永久基本農(nóng)田的凈面積，其他地類(lèi)只統(tǒng)計(jì)毛面積。在實(shí)際工作中，如用地預(yù)審與規(guī)劃選址、用地報(bào)批階段，對(duì)地類(lèi)的凈面積也有統(tǒng)計(jì)要求。因此，要進(jìn)一步納入地類(lèi)凈面積計(jì)算模型，以實(shí)現(xiàn)成果的精確統(tǒng)計(jì)。

參考文獻(xiàn)：

［1］許良一，陳晨.“多審合一”背景下的建設(shè)項(xiàng)目選址和用地預(yù)審合并審批的探索和思考［J］.特區(qū)經(jīng)略.2019（8）：82-84.

［2］顧建，王芹.基于Arc GIS的土地利用分類(lèi)面積統(tǒng)計(jì)及應(yīng)用研究［C］//江蘇省測(cè)繪學(xué)會(huì)，2013年度江蘇省測(cè)繪學(xué)會(huì)年會(huì)論文集.江蘇：《現(xiàn)代測(cè)繪》編輯部，2013：113-115.

［3］胡大國(guó).巧用SQL函數(shù)實(shí)現(xiàn)土地利用數(shù)據(jù)統(tǒng)計(jì)［J］.國(guó)土資源信息化，2007（1）：33-35.

［4］張潤(rùn)舒，汪虹希.Excel輸出表格中分項(xiàng)數(shù)與合計(jì)數(shù)不平衡的解決方案［J］.中國(guó)集體經(jīng)濟(jì)，2012（19）：183-184.

［5］鄭浩.基于ArcGIS Engine的土地利用現(xiàn)狀分類(lèi)統(tǒng)計(jì)系統(tǒng)的研究與實(shí)現(xiàn)［C］//中國(guó)測(cè)繪學(xué)會(huì)科技信息網(wǎng)分會(huì)，全國(guó)測(cè)繪科技信息網(wǎng)中南分網(wǎng)第三十次學(xué)術(shù)信息交流會(huì)論文集.南寧：廣西人民出版社，2016：207-214.

［6］高玉榮，高慶強(qiáng)，劉玉春，等.面向勘測(cè)的土地利用分類(lèi)統(tǒng)計(jì)［J］.城市勘測(cè).2015（6）：28-30.

［7］于國(guó)強(qiáng).基于VBA的土地勘測(cè)定界信息系統(tǒng)開(kāi)發(fā)［J］.價(jià)值工程.2015（9）：311-312.

［8］李詩(shī)宇，鄧吉強(qiáng).基于arcpy的地學(xué)數(shù)據(jù)自動(dòng)化處理與分析方法［J］.科技創(chuàng)新與生產(chǎn)力，2018（3）：44-46，49.

［9］郭誰(shuí)瓊，黃賢金，白曉飛，等.土地利用變更調(diào)查數(shù)據(jù)的應(yīng)用研究現(xiàn)狀與前景［J］.中國(guó)土地科學(xué)，2013（12）：18-24，97.

［10］陳棟，姚維嶺.基于ArcPy與定制ArcToolbox的礦山新增圖斑自動(dòng)編號(hào)及方法改進(jìn)［J］.國(guó)土資源遙感，2021（2）：262-269.

［11］歐丁丁，張琪，劉世好，等.基于ArcPy和VBA的長(zhǎng)順縣國(guó)家儲(chǔ)備林小班設(shè)計(jì)圖ArcGis的自動(dòng)出圖方法研究［J］.中南林業(yè)調(diào)查規(guī)劃，2021（1）：32-38，44.