荒漠化和沙化監測數據與國土三調數據的銜接方法
——以北票市為例

2021-09-14 02:28:38姜黎黎

內蒙古林業調查設計 2021年4期

姜黎黎

（遼寧省林業調查規劃監測院，沈陽 110112）

1 數據預處理

銜接以縣為單位，對兩套數據進行預處理。

1.1 坐標系統

荒漠化和沙化監測數據與國土三調數據的坐標系統均采用CGCS2000坐標系統。文章使用三度帶投影坐標系，遼寧省共涉及3 個三度帶：CGCS2000-3-Degree-GK-Zone-40，CGCS2000-3-Degree-GK-Zone-41，CGCS2000-3-Degree-GKZone-42。北票市所在的區域為40 度帶，通過投影變換將荒漠化和沙化監測數據與國土三調數據統一為CGCS2000-3-Degree-GK-Zone-40坐標系。

1.2 國土三調數據（數據A）添加字段

添加標識ID、國土圖形面積2個字段，并計算字段填入相應的值。

1.3 荒漠化和沙化監測數據（數據B）添加字段

添加荒漠化和沙化圖形面積字段并計算字段填入相應的值。

2 兩套數據處理

2.1 Identity工具處理

ArcGIS 的Identity 工具計算輸入要素和標識要素的幾何交集，輸入要素中與標識要素不相交的部分不變，與標識要素相交的部分獲取標識要素的屬性值[1]。數據A 為輸入要素，數據B 為標識要素,處理后生成數據C。

2.2 數據C處理

圖形吻合度，以量化的面積比例來描述，即數據C 的圖形面積/國土圖形面積×100。對數據C 求圖形面積、圖形吻合度。

2.3 數據C國土地類歸組到荒漠化和沙化地類

按林沙綜字[2020]5 號文件中的附件:第三次全國國土調查工作分類與荒漠化和沙化監測土地利用對照表，建立國土地類與荒漠化和沙化地類對應關系，運用python 語句構建模型工具進行國土地類歸組，如圖1。

圖1 國土地類歸組工具

3 數據分析及處理

3.1 選定圖形吻合度閾值

設定100，95，90，85，80，70共6個圖形吻合度閾值，對數據C 進行吻合度計算，并通過標識ID 對數據A 進行分類，數據A 共212 168 個圖斑，分為大于等于閾值、小于閾值兩類，數據A圖斑個數分布情況如表2，在綜合考慮數據精度及工作量的情況下，選定80作為圖形吻合度閾值。

表2 圖斑吻合度閥值

3.2 制定數據銜接方案

數據銜接時，將數據A分為四類進行銜接：圖形吻合度大于等于閾值的圖斑166 195個；圖形吻合度小于閾值，荒漠化和沙化地類為非林草地類的圖斑11 538個；圖形吻合度小于閾值，荒漠化和沙化地類為林草地類但國土地類編碼是500以上的圖斑（建設、工礦、交通、水利用地等無植被覆蓋地類）12 009個；圖形吻合度小于閾值，荒漠化和沙化地類為林草地類但國土地類編碼為500 以下的圖斑（濕地、耕地、園地、林地、草地等有植被覆蓋地類）22 426個。

圖形吻合度大于等于閾值的圖斑，采用直接掛接的方法，即圖形采用數據A的圖形，屬性采用數據C 的荒漠化和沙化監測數據屬性，使用ArcGIS 的JOIN 功能，通過標識ID 將數據A 和數據C 進行連接，生成銜接后的新數據，下文論述的掛接也是這種方法。

圖形吻合度小于閾值，荒漠化和沙化地類為非林草地類的圖斑先掛接地類相同的圖斑（前文對國土地類進行了歸組），如有未掛接上的，再掛接圖形吻合度最大的那個圖斑的屬性值。

圖形吻合度小于閾值，荒漠化和沙化地類為林草地類，但國土地類編碼是500 以上的圖斑掛接圖形吻合度最大的那個圖斑的屬性值。

圖形吻合度小于閾值，荒漠化和沙化地類為林草地類，但國土地類編碼為500以下的圖斑，按如下方法處理。

（1）分解多部件。在Identity處理時會產生多部件要素，使用 ArcGIS 的多部件至單部件MultipartToSinglepart工具對荒漠化和沙化地類是林草地類，國土地類500以下的數據，分解多部件。

（2）消除Eliminate 工具處理。ARCMAP 的Eliminate工具通過將面與具有最大面積或最長公用邊界的鄰近面合并來消除面。通過篩選選中面積小于0.066 7 hm2的圖斑，使用Eliminate工具處理，排除圖層為國土三調數據轉換的線圖層，這樣可以保證消除合并不突破國土三調數據圖斑界線。如圖2。