熊文豪　伊雨花

摘 要：地表覆蓋是地理國情監(jiān)測中重要的部分，其主要依賴人工解譯和自動人工結(jié)合的兩種方式。ArcGIS具備分類以及矢量處理等強大的功能，本文嘗試用ArcGIS自動人工相結(jié)合的方式進行地表覆蓋的解譯工作。詳述ArcGIS中ISO非監(jiān)督法分類在地理國情監(jiān)測地表覆蓋中的應(yīng)用，分析不同參數(shù)帶來的影響，并根據(jù)地理國情監(jiān)測的要求進行具體的實施流程確定。

關(guān)鍵詞：地理國情監(jiān)測；地表覆蓋；ISO非監(jiān)督法分類

中圖分類號：P237 文獻標(biāo)識碼：A 文章編號：1671-2064（2017）16-0249-02

在遙感影像分類中，根據(jù)是否已知訓(xùn)練樣本的數(shù)據(jù)，分類方法分為監(jiān)督分類和非監(jiān)督分類。非監(jiān)督分類的結(jié)果，只是對不同類別達到了區(qū)分，但是不確定類別的屬性，其類屬是通過事后對各類進行分析后確定的。在地理國情監(jiān)測地表覆蓋的解譯中就采用了非監(jiān)督分類，然后人工進行再分類。而ArcGIS的分類工具就具備ISO聚類非監(jiān)督分類。本文重點介紹了使用ArcGIS軟件按照地理國情的要求進行影像的分類和處理，以烏梁素海地理國情監(jiān)測為例，確定其流程和方法。

1 數(shù)據(jù)的制作

烏梁素海主要以挺水植物、沉水植物、浮游植物、水面為主，大部分屬于植物類，通過影像以及植被特性分析，重點進行植物和水面的劃分，最后在植物里區(qū)分挺水植物、沉水植物、浮游植物。這里采用一幅1萬的圖幅數(shù)據(jù)25平方公里數(shù)據(jù)的分類，并說明參數(shù)不同對分類的影像[1]。主要流程如圖1所示。

（1）首先在ArcGIS工具欄中使用影像分類工具，進行ISO聚類非監(jiān)督分類。ISO聚類非監(jiān)督分類共包含3個參數(shù)類數(shù)目、最小類大小、采樣間隔。其中類數(shù)目指要將像元劃分成的類數(shù)目，最小類的大小指一個有效類中的最小像元數(shù)，采樣間隔是指采樣所采用的間隔。因地理國情監(jiān)測數(shù)據(jù)要求矢量數(shù)據(jù)與影像的精度在5個像元內(nèi)，且此次分類最小面積1600平米，為達到較高的精度要求，最小類大小為300，采樣間隔默認10，類數(shù)目分別進行10、20、30的分類。經(jīng)過該分類后獲得的柵格數(shù)據(jù)分別為ISO10、ISO20、ISO30。

（2）數(shù)據(jù)的矢量化操作。地理國情監(jiān)測數(shù)據(jù)為矢量數(shù)據(jù)LCA層，并具備CC及TAG兩個屬性字段，因此需要轉(zhuǎn)換為矢量數(shù)據(jù)，方便后期操作。采用ArcGIS功能“轉(zhuǎn)換工具-由柵格轉(zhuǎn)出-柵格轉(zhuǎn)面”，分別獲得矢量數(shù)據(jù)。地表覆蓋圖斑數(shù)分別如表1所示。

其中可以發(fā)現(xiàn)隨著分類數(shù)的增加，圖斑數(shù)也聚增，數(shù)據(jù)量會變的很大，操作時間會變長，對計算機內(nèi)存等也有著一定的要求。分類數(shù)的增加，針對影像的分類也會更加精細，后期的處理也變得復(fù)雜。

（3）人工進行屬性分類。LCA層中的CC屬性字段代表著地物的分類屬性編碼。因整個圖幅基本可以劃分為水面、蘆葦，人為的選擇一級分類進行合并劃分為水面1001，其他的為蘆葦0431。在各個屬性字段中添加CC字符串字段，并根據(jù)grid_code進行計算。ISO10選擇1、2、5類，ISO20選擇1、2、3、5、8類，ISO30選擇1、2、3、4、8、13類，該類分別賦為水面1001。因需要使用邏輯處理，故使用代碼實現(xiàn)，如下：

If [grid_code]= 1 or [grid_code]= 2 or [grid_code]= 5 Then

[CC] = “1001”

else [CC] =”0431”

end if

分類數(shù)目增加，分類更加詳細，水面的紋理的區(qū)分更加細致，需要人為判定哪些類為水面。采用ArcGIS工具“數(shù)據(jù)管理工具-制圖綜合-融合”，針對同CC的屬性融合。

（4）數(shù)據(jù)的綜合。自動解譯過程中會產(chǎn)生大量的極小面，這是不符合地理國情監(jiān)測數(shù)據(jù)指標(biāo)的，因此需要進行小面的合并操作，分別合并至大面中。蘆葦0431和水面1001的面積指標(biāo)為1600平方米，在地理國情中面積指標(biāo)可以縮小30%，為1120平米，再次以1120為面積的限定條件。建立CC1字符串字段，將小面積的0431改為1001，小面積的1001改為0431，計算代碼如下：

If （[CC] = “1001” And [Shape_Area] <1120） Then

[CC] = “0431”

elseif （[CC] = “0431” And [Shape_Area] <1120） Then

[CC] = “1001”

end if

該操作極大的減少了小面積的噪聲，一般進行兩次操作后基本不存在極小面積等問題，最終三種分類方式的圖斑數(shù)為340，圖斑數(shù)從幾十萬到幾百，而且基本將比較明顯的類別區(qū)分出來。經(jīng)過劃分后比較三種方式，發(fā)現(xiàn)第2步人工確定分類，極大的影像了分類的精度，此處選擇錯誤將導(dǎo)致很多邊界不準(zhǔn)確或?qū)傩藻e誤，而隨著分類數(shù)目的增加，確定類的屬性變得困難，但是分類效果會更理想，尤其對于比較離散的紋理，ISO30分類效果較好，ISO20、ISO10存在屬性劃分錯誤[2]。

（5）去除齒狀節(jié)點。自動解譯是以像元為基礎(chǔ)進行分割處理，經(jīng)過此操作后，邊界齒狀明顯，在地理國情中要求圖斑節(jié)點數(shù)不能超過10萬個節(jié)點，為了減少節(jié)點數(shù)，進行簡化面操作。在ArcGIS工具中“制圖工具-制圖綜合-簡化面”，簡化面是在不改變基本幾何形狀的情況下，通過移除多余的彎曲簡化面要素。采用POINT_REMOVE簡化算法，包含兩個簡化容差，最大允許偏移量及最小面積，其中最大允許偏移量指轉(zhuǎn)換前后節(jié)點的偏移大小。偏移量設(shè)置的越大形狀改變越大，與影像的套和程度越低。因影像分辨率為2米，這里參數(shù)設(shè)置為2至10之間確保精度，當(dāng)偏移量為5時，改變了較多的形狀，邊界存在部分不準(zhǔn)確。經(jīng)測試3能達到較好的效果，節(jié)點減少較多且與影像套和效果好。在ArcGIS10.0版本以上簡化面工具新增了處理拓撲錯誤工具，該功能解決了以往簡化面中存在的大量重疊面以及簡化面操作。

（6）最后在此基礎(chǔ)上劃分浮游藻類以及沉水植物，并進行數(shù)據(jù)的整理和檢查工作，解決同普異物等問題，確保地理國情數(shù)據(jù)的精確以及可靠。

2 結(jié)語

通過ISO聚類非監(jiān)督分類極大的減小了工作量，不用人工進行面分割處理，提高了效率。在分類過程中發(fā)現(xiàn)不同的參數(shù)對自動解譯精度有較大的影響，然而分類數(shù)增加電腦的處理時間會更長，直接影響到后續(xù)的操作。因此在進行某一地區(qū)項目時，因根據(jù)影像情況制定參數(shù)方案，最后進行整個項目區(qū)的處理。

參考文獻

[1]蔣好忱，任宏權(quán)，秦先鋒，等.eCognition影像自動解譯及精度評定[J].測繪通報，2015，（10）：81-84.

[2]覃宇，張迅.地表覆蓋自動解譯方案在地理國情普查項目中的應(yīng)用[J].測繪與空間地理信息，2015，（06）：164-166.endprint