基于地理信息系統的建設用地審查報批中地類面積計算方法淺析

李良毅LI Liang-yi

（黔西南州自然資源局，興義 562400）

0 引言

近年來，貴州省社會經濟高速發展，建設用地需求大量增加，僅2022 年，經貴州省自然資源廳復審，報省人民政府（包括省人民政府授權市、州人民政府）批準，并報自然資源部備案的項目建設用地共計770 宗、用地面積9075.3835 公頃（其中：新增建設用地8284.6738 公頃，農用地7618.7817 公頃，耕地3908.3482 公頃，未利用土地677.6121 公頃），由此可以得出，在辦理建設用地審查報批中，涉及到大量各地類面積的計算。傳統的計算統計方法繁瑣復雜且不準確，因此使用現代地理信息系統技術展現出數據庫圖層集中疊加拓撲處理、地塊采集與計算處理能力強、數據精準等各類優勢，為建設用地審查報批提供強大的技術保障支撐作用。

1 地理信息系統建設用地審查報批中的應用

地理信息系統在建設用地審查報批中的作用與地位日益凸顯，而且處于基礎性地位，在建設用地報批地類審查過程可以說是不可或缺。目前，建設用地報批地類審查基于全國第三次國土調查數據庫為基礎的最新年度變更調查數據庫，該數據庫采用國家測繪標準為：《土地利用現狀分類》（GB/T 21010-2017）、《遙感影像平面圖制作規范》（GB 15968-2008）、《全球定位系統（GPS）測量規范》（GB/T 18314-2009）。采用的數學基礎為：大地基準：采用“2000國家大地坐標系”；高程基準：采用“1985 國家高程基準”；投影分帶：采用高斯-克呂格投影3°分帶。農村調查以1:5000 比例尺標準分幅圖為主，對于調查精度較高的城鎮調查部分，采用1∶2000 比例尺標準分幅。按照國家相關規范技術要求，建立起最新的年度變更調查數據庫，在數據庫建立過程中，地理信息系統尤其是坐標系、高程基準可謂發揮了基礎性、至關重要的作用。尤其地理信息系統在圖層管理、屬性字段、遙感影像的分層疊加，地塊與各類圖層的裁剪分析，數據計算方面優勢突出。由此可見，在年底變更調查數據庫基礎上開展的地類審查工作，地理信息系統發揮了不容忽視的作用。

2 名詞解釋

本文重點論述建設用地審查報批中地類面積計算方法，對國土調查相關術語進行解釋如下。

建設用地：是指建造建筑物、構筑物的土地，包括城鄉住宅和公共設施用地，工礦用地，能源、交通、水利、通信等基礎設施用地，旅游用地，軍事用地等。亦是需要付出一定投資（土地開發建設費用），通過工程手段，為各項建設提供的土地。

農用地：直接或間接為農業生產所利用的土地。包括耕地、園地、林地、牧草地、養殖水面、坑塘水面、農田水利設施用地，以及田間道路和其他一切農業生產性建筑物占用的土地等。

未利用地：指農用地和建設用地以外的土地，主要包括荒草地、鹽堿地、沼澤地、沙地、裸土地、裸巖等。

第三次國土調查：在第二次全國土地調查成果基礎上，全面細化和完善全國土地利用基礎數據，掌握翔實準確的全國國土利用現狀和國土資源變化情況，進一步完善國土調查、監測和統計制度，實現成果信息化管理與共享，滿足生態文明建設、空間規劃編制、供給側結構性改革、宏觀調控、自然資源管理體制改革和統一確權登記、國土空間用途管制、國土空間生態修復、空間治理能力現代化和國土空間規劃體系建設等各項工作的需要。

建設用地審查報批：根據2016 年修正的《建設用地審查報批管理辦法》，第四條：“在建設項目審批、核準、備案階段，建設單位應當向建設項目批準機關的同級國土資源主管部門提出建設項目用地預審申請。受理預審申請的國土資源主管部門應當依據土地利用總體規劃、土地使用標準和國家土地供應政策，對建設項目的有關事項進行預審，出具建設項目用地預審意見。”

農用地轉用：指按照土地利用總體規劃和國家規定的批準權限獲得批準后，將農用地轉變為建設用地的行為。又稱為農用地轉為建設用地。用于非農建設有以下情形之一者，應當辦理農用地轉用審批手續：①征用農村集體經濟組織農用地的；②農村集體經濟組織使用本集體農用地的；③使用國有農用地；④需要辦理農用地轉用的其他土地。

3 建設用地追溯概述

我國為落實最嚴格的耕地保護制度，毫不動搖地堅持耕地保護紅線。嚴格限制農用地轉為建設用地，控制建設用地總量，對耕地實行特殊保護，在近年來的建設用地審查報批時對建設用地進行追溯。追溯流程如圖1 所示，將審查報批范圍內最新年度國土變更數據庫中的地類進行分類，其中的農用地和未利用地按自身地類報批，建設用地根據其來源合法性進行分類處理，有合法來源按照建設用地地類報批；無合法來源的，需追溯到上一年度變更數據庫，將追溯建設用地范圍內上一年度變更數據庫中的地類再進行分類，農用地和未利用地按自身地類報批，建設用地需重復追溯過程，直到追溯地類分類中沒有建設用地或已經追溯到最早年度的國土調查數據庫（最早可追溯至第二次土地調查數據庫），才停止追溯。

圖1 建設用地追溯流程示意圖

4 追溯過程中的問題及分析

在建設用地追溯過程中，一般會追溯到2019 年之前的數據庫，這樣計算就會同時涉及到三調數據庫和二調數據庫。三調數據庫中只涉及到面狀的地類圖斑圖層，追溯時關系明晰，使用地理信息系統的空間疊加分析就能得到追溯結果，而二調數據庫中涉及到的除了面狀的地類圖斑圖層（DLTB），還涉及到線狀的線狀地物圖層（XZDW）和點狀的零星地物圖層（LXDW），如何在追溯時統籌考慮這三個圖層之間的關系，進而追溯出新增建設用地占用其他地類的面積，成為了追溯中遇到的主要難題。

以單個年度二調地類面積的計算方法為例，闡述因追溯而引起的計算問題。計算單個年度二調地類面積的方法如下：

計算單個年度二調地類面積的方法如下：

①零星地物面積SLX為LXDW 圖層的MJ 字段值；

②線狀地物面積SXZ=長度*寬度*扣除比例；

③田坎面積STK=（圖斑面積S-SLX-SXZ）*田坎系數；

④圖斑地類面積SDL=圖斑面積S-SLX-SXZ-STK。

建設用地追溯的過程，可看做是將當前的建設用地范圍套合上一年度變更數據庫中計算的過程，而從上述的二調地類面積計算方法中可知，二調中可能包含建設用地的部分有零星地物面積SLX、線狀地物面積SXZ和圖斑地類面積SDL，但這些面積卻不一定對應明確范圍。如圖2 所示，在2.45 萬平方米的建制鎮中有500 平方米的農村道路，作為建設用地需追溯的建制鎮面積為2.4 萬平方米，就沒有一個明確的范圍，故需要找到計算方法來解決此類問題。

5 地類面積計算方法及代碼實現

針對上述問題，下面提出兩種解決方法：

方法一：圖斑面積比較法。追溯的實質是當前的建設用地與上一年度相比是否增加，如果增加了，是什么地類轉變成了建設用地，可以表示為：當前的建設用地面積SJn 減去上一年度的建設用地面積SJ（n-1），即SJn-SJ（n-1）。可以把地類圖斑、零星地物、線狀地物看做是一個整體，它們共同表示了地類圖斑范圍內由什么地類的組成，圖斑面積STBn 等于建設用地面積SJn 加上農用地面積SNn 再加上未利用地面積SWn，公式為STBn=SJn+SNn+SWn（公式1）。只要圖斑范圍類有建設用地，都將整個圖斑范圍帶入上一年度數據庫進行追溯，STBn=SJ（n-1）+SN（n-1）+SW（n-1）（公式2）。通過整理公式1 和公式2 可以得出：SJn-SJ（n-1）=（SN（n-1）-SNn）+（SW（n-1）-SWn），即增加的建設用地SJn-SJ（n-1）等于減少的農用地與減少的未利用地之和。

方法二：分層比例法。在追溯過程中，將2018 年度的地類圖斑、零星地物、線狀地物統一起來，按照單個年度二調地類面積的計算方法計算出地類分類面積，然后每個圖層的建設用地再去各自圖層分別追溯，在追溯地類圖斑圖層時遇到地類是建設用地，且地類圖斑范圍內包含零星地物或線狀地物的，追溯到上一年度時，以地類圖斑范圍追溯，建設用地面積作為控制面積，各地類的面積需要乘以建設用地的面積與圖斑面積的比例。以圖2 為例，需要用建制鎮面積2.4 萬平方米作為控制面積，追溯到上一年度各地類面積需要乘以比例2.4/2.45。這兩種方法都能完成追溯，方法一的計算步驟比較直觀，容易理解，計算難度底，貼合實際生活中的體驗，但在農用地和未利用地地類上，當上一年度的面積小于當前年度面積時，計算結果中會出現負數，這會引起對此方法不熟悉人員的疑惑；方法二中計算出的結果不會出現負數，常人更容易理解，但在計算過程中需要記錄下一年度建設用地與圖斑面積的比例，計算略微復雜。通過編寫軟件完成計算過程后，方法二的結果更易理解，所以在生產環境中推薦使用方法二。

圖2 二調中的地類示例

實現核心代碼如下（本文使用Python 代碼實現了所有功能）：

將2018 年度的地類圖斑、零星地物、線狀地物統一起來計算后得出這三個圖層中各自的建設用地。

相關變量定義如下：

dltbJsyd：地類圖斑中的建設用地；

lxdwJsyd：零星地物中的建設用地；

xzdwJsyd：線狀地物中的建設用地。

# 圖層各自的建設用地從上一年度繼續追溯，以2018

追溯現狀地物核心代碼如下：

# len: 紅線范圍內線狀地物長度kd: 線狀地物寬度jsydLength:原始線狀地物中建設用地的長度 oriLength:原始線狀地物的長度

追溯零星地物核心代碼如下：

# 零星地物為點狀，只要在紅線范圍內，必定所有面積都是同一地類

resultMap.append(mj)

6 結語

本文介紹了兩種解決建設用地追溯中問題的方法，基于結果更易理解的特性，推薦在生產環境中使用方法二。筆者基于方法二制作分類面積計算軟件，在貴州省黔西南州、遵義市等地的建設用地報批中得到使用，提高了報批的工作效率。