摘" 要：掌握城市生態(tài)用地結(jié)構(gòu)變化對(duì)保護(hù)生態(tài)用地、建設(shè)新型城市具有重要意義。該文選取畢節(jié)七星關(guān)城區(qū)為研究區(qū)，以2016年、2021年的GF-2遙感影像為數(shù)據(jù)源，通過(guò)監(jiān)督分類(lèi)、非監(jiān)督分類(lèi)以及面向?qū)ο蠓诸?lèi)方法進(jìn)行生態(tài)用地遙感監(jiān)測(cè)及變化分析。研究結(jié)果表明，面向?qū)ο蟮姆诸?lèi)方法能夠更好地實(shí)現(xiàn)對(duì)高分影像分類(lèi)，其中，監(jiān)督分類(lèi)中最大似然法分類(lèi)結(jié)果的總體精度為80.47%，非監(jiān)督分類(lèi)中ISODATA分類(lèi)結(jié)果的總體精度為65.56%，而面向?qū)ο蟮姆诸?lèi)方法總體精度可達(dá)到91.82%；面向?qū)ο蟮姆诸?lèi)方法結(jié)合了光譜特征、形狀特征、紋理特征以及對(duì)象之間的拓?fù)湫畔?，?dāng)分割尺度為20時(shí)可以將影像對(duì)象很好地分割，解譯結(jié)果更精準(zhǔn)，分類(lèi)效果更好；2016—2021年間，七星關(guān)區(qū)生態(tài)用地范圍持續(xù)減小，共減少面積19.32 km2。其中生態(tài)用地轉(zhuǎn)化為建設(shè)用地的面積為17.78 km2，占比92.03%；生態(tài)用地轉(zhuǎn)化為道路的面積為1.54 km2，占比7.97%。

關(guān)鍵詞：高分二號(hào)；城市生態(tài)用地；面向?qū)ο蠓诸?lèi)方法；監(jiān)督分類(lèi)；非監(jiān)督分類(lèi)

中圖分類(lèi)號(hào)：P237 文獻(xiàn)標(biāo)志碼：A" " " " " 文章編號(hào)：2095-2945（2023）16-0079-05

Abstract： It is of great significance to master the structural change of urban ecological land for protecting ecological land and building a new city. This paper selects Qixingguan urban area of Bijie City as the study area， takes the GF-2 remote sensing images of 2016 and 2021 as the data source， and carries on the remote sensing monitoring and change analysis of ecological land through supervised classification， unsupervised classification and object-oriented classification. The results show that， the object-oriented classification method can better achieve the classification of high-score images， in which the overall accuracy of the maximum likelihood classification results in supervised classification is 80.47%， the overall accuracy of ISODATA classification results in unsupervised classification is 65.56%， and the overall accuracy of object-oriented classification method can reach 91.82%. The object-oriented classification method combines spectral features， shape features， texture features and topological information between objects. When the segmentation scale is 20， the image object can be segmented well， the interpretation result is more accurate， and the classification effect is better. From 2016 to 2021， the scope of ecological land in Qixingguan District continued to decrease， with a total reduction area of 19.32 km2. Among them， the area of ecological land converted into construction land is 17.78 km2， accounting for 92.03%; and the area of ecological land converted into roads is 1.54 km2， accounting for 7.97%.

Keywords： GF-2; urban ecological land; object-oriented classification; supervised classification; unsupervised classification

城市生態(tài)用地作為一種土地利用方式是城市復(fù)合生態(tài)系統(tǒng)的重要組成部分[1]，在維持生態(tài)平衡、保障生態(tài)安全、應(yīng)對(duì)全球氣候變化中具有特殊地位[2]。然而，隨著城市化進(jìn)程的加速，城市生態(tài)用地大量被占用，對(duì)城市的生態(tài)環(huán)境造成了一定的損害。

隨著遙感技術(shù)的發(fā)展，基于遙感的城市用地及生態(tài)環(huán)境研究逐漸成為熱點(diǎn)。在遙感影像分類(lèi)研究中，Baatz等[3]針對(duì)高空間分辨率遙感影像提出了面向?qū)ο蟮倪b感信息分類(lèi)方法；Shackelford等[4]指出面向?qū)ο蠓诸?lèi)方法比較適合于城市或城郊區(qū)域分類(lèi)；王建如等[5]提出了一種高空間分辨率遙感數(shù)據(jù)最優(yōu)分割結(jié)果確定方法；黃邵東等[6]使用面向?qū)ο蟮姆诸?lèi)方法對(duì)研究區(qū)茶園進(jìn)行提?。环冬摿盏萚7]提出了一種基于多特征的面向?qū)ο蠓诸?lèi)方法；劉凌露等[8]在遙感影像分類(lèi)中明確了分割精度及分類(lèi)規(guī)則集，便于進(jìn)行城市生態(tài)用地監(jiān)測(cè)；孟晉杰等[9]建立相應(yīng)的規(guī)則集以提取遙感影像中生態(tài)用地信息；劉蕓等[10]根據(jù)NDVI特征、光譜特征、紋理特征及形狀特征提取信息，認(rèn)為面向?qū)ο蠓诸?lèi)方法分類(lèi)質(zhì)量達(dá)到了很好的水平；馬飛虎等[11]在不同層次進(jìn)行最優(yōu)尺度地物分割，認(rèn)為多尺度多層次分割分類(lèi)的總體分類(lèi)精度更高；胡金梅等[12]對(duì)影像進(jìn)行分割，綜合利用SVM和隨機(jī)森林分類(lèi)算法以實(shí)現(xiàn)土地利用分類(lèi)。

本研究基于前人研究成果，擬探索適用于高分二號(hào)（GF-2）影像數(shù)據(jù)的生態(tài)用地分類(lèi)方法，以期為國(guó)產(chǎn)GF-2衛(wèi)星監(jiān)測(cè)城市生態(tài)用地提供參考。

1" 研究區(qū)概況及衛(wèi)星影像預(yù)處理

1.1" 研究區(qū)概況

本文選取位于典型喀斯特山區(qū)的畢節(jié)市七星關(guān)區(qū)為研究區(qū)域。七星關(guān)區(qū)是畢節(jié)市的政治、經(jīng)濟(jì)、文化中心，是交通運(yùn)輸物流中心，位于貴州省西北部。西接赫章縣，東部與金沙縣相連，南鄰納雍縣。東西長(zhǎng)105 km，南北寬80 km，總面積3 412.2 km2。該區(qū)地形西高東低，地勢(shì)由西南向東北傾斜，境內(nèi)喀斯特地貌發(fā)育良好。

1.2" 數(shù)據(jù)源

本文數(shù)據(jù)源主要有：①遙感影像數(shù)據(jù)，本研究遙感影像數(shù)據(jù)分別為2016年8月24日、2021年3月31日成像的GF-2影像，包括紅、綠、藍(lán)和近紅外4個(gè)波段；②矢量數(shù)據(jù)，主要包括畢節(jié)市中心城區(qū)。

GF-2衛(wèi)星2014年8月19日發(fā)射成功，其星下點(diǎn)全色波段空間分辨率可達(dá)0.8 m，標(biāo)志著我國(guó)遙感衛(wèi)星進(jìn)入了亞米級(jí)“高分時(shí)代”[13]。GF-2配備了2個(gè)1 m全色、4 m多波段高分辨率攝像機(jī)，具有空間分辨率高、定位精度高、快速姿態(tài)機(jī)動(dòng)能力等特點(diǎn)。

1.3" 數(shù)據(jù)預(yù)處理

1.3.1" 正射校正

正射校正可以消除由于衛(wèi)星相機(jī)傾斜、地形起伏及地物所引起的圖像畸變[14]。正射校正通常是在圖像上選擇一定的地面控制點(diǎn)，然后利用DEM資料，對(duì)圖像進(jìn)行傾斜修正和投影誤差修正，使圖像重新取樣為正射圖像。本文選擇PIE-Basic 6.3平臺(tái)上“正射校正”工具，基于GF-2自帶的RPC文件進(jìn)行正射校正。

1.3.2" 圖像配準(zhǔn)

圖像配準(zhǔn)是圖像融合的基礎(chǔ)，也是圖像預(yù)處理中的基本環(huán)節(jié)之一。本文使用PIE-Basic 6.3中“影像配準(zhǔn)”工具，以全色數(shù)據(jù)作為基準(zhǔn)影像對(duì)多光譜數(shù)據(jù)進(jìn)行配準(zhǔn)，為下一步圖像融合奠定基礎(chǔ)。

1.3.3" 圖像融合

利用空間分辨率較低的多譜影像或高光譜資料與高空間分辨率全色波段影像進(jìn)行再取樣，生成高分辨率的多光譜影像，稱(chēng)為影像融合。本文使用PIE-Basic 6.3平臺(tái)中的Pansharp融合方法對(duì)GF-2影像進(jìn)行融合處理。

1.3.4" 圖像裁剪

由于本研究區(qū)同時(shí)覆蓋兩幅影像，所以在裁剪之前，需要將兩幅影像鑲嵌，將鑲嵌后的影像進(jìn)行裁剪。圖像裁剪是為了減少后期數(shù)據(jù)的處理量，裁剪分為規(guī)則裁剪與不規(guī)則裁剪，本文借助PIE-Basic 6.3平臺(tái)“圖像裁剪”工具，使用畢節(jié)市中心城區(qū)矢量圖進(jìn)行裁剪，得到研究區(qū)影像。

2" 基于像元的城市生態(tài)用地遙感分類(lèi)

2.1" 監(jiān)督分類(lèi)

監(jiān)督分類(lèi)又稱(chēng)訓(xùn)練場(chǎng)法，它是一種以統(tǒng)計(jì)為基礎(chǔ)的遙感分類(lèi)法[14]。蘇志強(qiáng)[15]發(fā)現(xiàn)最大似然法能夠最大化地發(fā)揮GF-2影像在土地利用分類(lèi)中的潛力，且其分類(lèi)精度較高，因此，本文選用最大似然法將預(yù)處理好的影像進(jìn)行監(jiān)督分類(lèi)。

最大似然分類(lèi)（Maximum Likelihood Classification）的原理是求出每一個(gè)像元對(duì)于各類(lèi)別的歸屬概率，把該像元?jiǎng)澐值綒w屬概率最大的類(lèi)別中去[14]。七星關(guān)區(qū)土地用地類(lèi)型復(fù)雜，包括綠地、水體、裸地、建設(shè)用地、道路，本研究選取GF-2影像光譜特征參數(shù)見(jiàn)表1。

本文使用PIE-Basic 6.3平臺(tái)中的“最大似然分類(lèi)”模塊進(jìn)行監(jiān)督分類(lèi)，最終得到了畢節(jié)市七星關(guān)城區(qū)土地利用現(xiàn)狀的監(jiān)督分類(lèi)圖，如圖1所示。

2.2" 非監(jiān)督分類(lèi)

非監(jiān)督分類(lèi)通常使用迭代自組織數(shù)據(jù)分析（ISODATA）和K均值聚類(lèi)算法[14]。胡文敏等[16]發(fā)現(xiàn)ISODATA分類(lèi)方法精度較高，故本文使用ISODATA分類(lèi)方法進(jìn)行非監(jiān)督分類(lèi)。為了保證分類(lèi)精度，增加了迭代次數(shù)為10，并將高分影像分為15類(lèi)，得到七星關(guān)區(qū)非監(jiān)督分類(lèi)結(jié)果如圖2所示。

2.3" 面向?qū)ο蠓诸?lèi)

面向?qū)ο蟮姆诸?lèi)方法不同于傳統(tǒng)的分類(lèi)方法，其分類(lèi)流程主要包括影像分割和分類(lèi)兩步。影像分割的目的是為了將具有相同性質(zhì)的影像劃分在一起，從而提取目標(biāo)地物的光譜、紋理、幾何等特征，以適應(yīng)解譯精度的要求。

2.3.1" 影像分割

影像分割是面向?qū)ο蟮幕A(chǔ)和關(guān)鍵，影像分割方法分為兩大類(lèi)，一種是自上而下的知識(shí)驅(qū)動(dòng)方法，一種是自下而上的數(shù)據(jù)驅(qū)動(dòng)方法[17]。在多次實(shí)驗(yàn)的基礎(chǔ)上，結(jié)合影像對(duì)象的特點(diǎn)，選擇適合于本地區(qū)圖像的分割尺度。結(jié)果表明，分割尺度為20對(duì)不同類(lèi)型的地物邊界具有良好的識(shí)別效果，分割時(shí)影像各波段權(quán)重比例設(shè)為（1，1，1，3），緊實(shí)度為0.5，形狀參數(shù)設(shè)為0.2。

2.3.2" 分類(lèi)方法規(guī)則集

經(jīng)過(guò)多次實(shí)驗(yàn)，本文提出的規(guī)則集主要是光譜均值和圖像亮度，形狀的規(guī)則即長(zhǎng)寬比（length/width）；地形的規(guī)則包括DEM等。根據(jù)上述的特征值，在結(jié)合規(guī)則集進(jìn)行分類(lèi)時(shí)，選用通過(guò)多次試驗(yàn)分類(lèi)精度較高的最鄰近分類(lèi)法。

在進(jìn)行最鄰近分類(lèi)法分類(lèi)后，發(fā)現(xiàn)明顯錯(cuò)誤和由于兩種地物光譜特征極為相似所帶來(lái)的錯(cuò)分較多，本研究利用易康（eCognition）的人工編輯模塊進(jìn)行手動(dòng)修改。修改完成后，為了提高影像分類(lèi)的完整性，還需要使用eCognition中的Merge選項(xiàng)對(duì)被分割的影像進(jìn)行合并。影像合并完成后發(fā)現(xiàn)小圖斑較多，因此進(jìn)行小圖斑處理，最后將分類(lèi)處理好的影像以raster格式輸出，即輸出柵格文件，得到如圖3所示分類(lèi)結(jié)果。

2.4" 分類(lèi)結(jié)果及精度評(píng)價(jià)

在分類(lèi)結(jié)果精度評(píng)價(jià)中，選擇156個(gè)樣本，計(jì)算分類(lèi)結(jié)果的總體精度和Kappa系數(shù)（表2）。其中，監(jiān)督分類(lèi)的最大似然法總體精度為80.47%，Kappa指數(shù)分別為0.714 7。經(jīng)分析可知，裸露土地與建設(shè)用地、道路等都很難分辨，其主要原因在于建筑用地的影子與道路的光譜特性十分相近，導(dǎo)致道路與建設(shè)用地之間的差異很大，監(jiān)督分類(lèi)對(duì)綠地和水體的分類(lèi)效果較好。非監(jiān)督分類(lèi)ISODATA分類(lèi)算法總體精度為65.56%，Kappa系數(shù)為0.506 9。經(jīng)分析可知，非監(jiān)督分類(lèi)中建設(shè)用地、裸地以及道路難以區(qū)分，裸地與道路的分類(lèi)精度不超過(guò)50%，因此，非監(jiān)督分類(lèi)不適合高分影像。面向?qū)ο蠓诸?lèi)的總體分類(lèi)精度和Kappa系數(shù)分別為91.82%和0.890 0。綜上所述，面向?qū)ο蟮姆诸?lèi)方法分類(lèi)精度最高，因此，本文選擇面向?qū)ο蠓诸?lèi)方法對(duì)研究區(qū)兩期影像進(jìn)行解譯，以進(jìn)行城市生態(tài)用地變化監(jiān)測(cè)。

3" 基于GF-2的城市生態(tài)用地變化分析

3.1" 遙感影像解譯結(jié)果

本研究對(duì)2016年研究區(qū)GF-2影像進(jìn)行面向?qū)ο蟮纳鷳B(tài)用地監(jiān)測(cè)，分割尺度設(shè)置為20，各波段權(quán)重比率（1，1，1，1）、緊密性0.5，形狀參數(shù)0.2，根據(jù)上述規(guī)則，對(duì)圖像進(jìn)行嚴(yán)格的分類(lèi)，確保了一定的分類(lèi)準(zhǔn)確性，分類(lèi)結(jié)果如圖4所示。

通過(guò)分析圖3、圖4可以看出，城區(qū)植被分布較多，呈現(xiàn)較大的片狀區(qū)域分布，與各種地類(lèi)交叉分布，建設(shè)用地呈大面積條帶狀分布于城市中心，水體邊界較為明顯，主要呈條狀穿插在建設(shè)用地中，道路多呈網(wǎng)格狀分布，裸地較少，呈點(diǎn)狀分布于建設(shè)用地邊緣。在分類(lèi)過(guò)程中，發(fā)現(xiàn)建筑陰影與水體的光譜特征極其類(lèi)似，對(duì)分類(lèi)結(jié)果存在一定的影響，但影響不大，并且經(jīng)過(guò)人工校正后錯(cuò)分現(xiàn)象得到極大改善。

3.2" 城市生態(tài)用地變化分析

結(jié)合2016年和2021年七星關(guān)區(qū)土地利用分類(lèi)圖進(jìn)行分析，統(tǒng)計(jì)各種地類(lèi)的用地面積得到七星關(guān)地區(qū)土地面積的變化情況，如表3和圖5所示。

經(jīng)分析可知：①研究區(qū)植被明顯減少，七星關(guān)區(qū)的植被用地從76.33 km2降為58.91 km2。其中，16.10 km2變化為建設(shè)用地，占總變化量的92.42%，0.90 km2變化為道路，占比5.17%。②建設(shè)用地面積增加，建設(shè)用地由17.89 km2上升為36.25 km2，增加幅度為18.32%。③水體少量減少，水體由2016年的1.66 km2下降為1.63 km2，下降幅度為0.03%。其中，0.02 km2變化為建設(shè)用地，占比66.67%，0.01 km2變化為道路，占比33.33%。④裸地略有減少，2016年裸地面積為2.85 km2，而到2021年下降為0.98 km2，下降幅度為1.86%。其中，1.66 km2變化為建設(shè)用地，占比88.77%。0.19 km2變化為道路，占比10.16%。⑤道路略微增加，2016—2021年來(lái)由于交通主干道的修建，道路面積從原來(lái)的1.52 km2增加為2.49 km2，增加幅度為0.96%。

4" 結(jié)論與討論

本文以位于典型喀斯特山區(qū)的畢節(jié)市七星關(guān)城區(qū)為研究區(qū)，以GF-2影像為基礎(chǔ)數(shù)據(jù)，采用監(jiān)督分類(lèi)中最大似然分類(lèi)法、非監(jiān)督分類(lèi)中ISODATA分類(lèi)法和面向?qū)ο蠓诸?lèi)對(duì)其土地利用進(jìn)行分類(lèi)，提出了適合于城市生態(tài)用地的分類(lèi)準(zhǔn)則，并將其應(yīng)用于城市生態(tài)用地監(jiān)測(cè)。通過(guò)研究，主要獲得以下結(jié)論。

1）面向?qū)ο蟮姆诸?lèi)方法對(duì)影像分類(lèi)結(jié)果的精度評(píng)價(jià)最高。傳統(tǒng)的分類(lèi)方法具有一定的局限性，其中監(jiān)督分類(lèi)中最大似然法分類(lèi)結(jié)果的總體精度為80.47%，非監(jiān)督分類(lèi)中ISODATA法分類(lèi)結(jié)果的總體精度為65.56%，而面向?qū)ο蟮淖钹徑诸?lèi)方法總體精度可以達(dá)到91.82%。

2）在進(jìn)行面向?qū)ο蟮亩喑叨确指顣r(shí)，分別嘗試了尺度為80、50、20，最終發(fā)現(xiàn)，尺度為80時(shí)，可以明顯區(qū)分大面積湖泊，而河道則不能區(qū)分；尺度為50時(shí)，可以明顯區(qū)分較大面積的地物類(lèi)型；只有在尺度為20時(shí)，才可以將各種面積較小的地物進(jìn)行區(qū)分。因此，在進(jìn)行面向?qū)ο蠓治鰰r(shí)，分割尺度為20時(shí)分類(lèi)效果最好。

3）通過(guò)對(duì)七星關(guān)區(qū)的生態(tài)用地監(jiān)測(cè)，生態(tài)用地范圍持續(xù)減小。其中生態(tài)用地轉(zhuǎn)化為建設(shè)用地的面積為17.78 km2，占比92.03%；生態(tài)用地轉(zhuǎn)化為道路的面積為1.54 km2，占比7.97%。

通過(guò)面向?qū)ο蟮姆诸?lèi)方法對(duì)七星關(guān)區(qū)城市生態(tài)用地進(jìn)行監(jiān)測(cè)，發(fā)現(xiàn)在面向?qū)ο蠓诸?lèi)過(guò)程中，由于建筑用地陰影的光譜特征與水體的光譜特征相似，在分類(lèi)過(guò)程中出現(xiàn)錯(cuò)分現(xiàn)象；由于兩期影像獲取的月份不同，在植被的分類(lèi)中有所偏差。下一步將通過(guò)去除陰影的方法對(duì)城市建筑用地的陰影進(jìn)行去除，在進(jìn)行面向?qū)ο蠓诸?lèi)時(shí)能夠更加準(zhǔn)確。

參考文獻(xiàn)：

[1] 彭建，汪安，劉焱序，等.城市生態(tài)用地需求測(cè)算研究進(jìn)展與展望[J].地理學(xué)報(bào)，2015，70（2）：333-346.

[2] 王芳，楊武年，鄧曉宇，等.高分二號(hào)數(shù)據(jù)的城市生態(tài)用地分類(lèi)方法探[J].測(cè)繪科學(xué)，2018，43（3）：71-76.

[3] BAATZ M， SCHAPE A. Multiresolution segmentation： An optimization approach for high quality mult-i scale image segmentation[J]. Angewandte Geographische Information sverarbeitung， 2000（12）：12-23.

[4] SHACKELFORD" A K， Davis C H A. Combined" Fuzzy" Pixel-Based" and" Object-Based Approach for Classification of High-Resolution Multispectral Data over Urban Areas[C]//IEEE Transactions on Geoscience and Remote Sensing，2003，41（10）：2354-2363.

[5] 王建如，程結(jié)海.一種高空間分辨率遙感數(shù)據(jù)最優(yōu)分割結(jié)果自動(dòng)化確定方法[J].測(cè)繪地理信息，2021，46（3）：83-87.

[6] 黃邵東，徐偉恒，熊源，等.結(jié)合紋理和空間特征的多光譜影像面向?qū)ο蟛鑸@提取[J].光譜學(xué)與光譜分析，2021，41（8）：2565-2571.

[7] 范瑩琳，婁德波，張長(zhǎng)青，等.基于面向?qū)ο蟮蔫F尾礦信息提取技術(shù)研究——以遷西地區(qū)北京二號(hào)遙感影像為例[J].自然資源遙感，2021，33（4）：153-161.

[8] 劉凌露.基于高分一號(hào)影像的城市生態(tài)用地監(jiān)測(cè)研究[D].北京：中國(guó)地質(zhì)大學(xué)（北京），2015.

[9] 孟晉杰，王建華.高分一號(hào)影像數(shù)據(jù)城市生態(tài)用地監(jiān)測(cè)[J].測(cè)繪科學(xué)，2016，41（9）：33-37.

[10] 劉蕓，李雪，廖瑤，等.基于面向?qū)ο蠓诸?lèi)法的貴州烤煙種植區(qū)域提取——以烤煙成熟期為例[J].貴州科學(xué)，2021，39（2）：82-86.

[11] 馬飛虎，徐發(fā)東，孫翠羽.面向?qū)ο蟮臒o(wú)人機(jī)影像地物分類(lèi)[J].應(yīng)用科學(xué)學(xué)報(bào)，2021，39（2）：312-320.

[12] 胡金梅，董張玉，楊學(xué)志.一種面向?qū)ο蟮母叻直媛蔬b感影像信息提取方法[J].地理空間信息，2021，19（9）：10-13，18.

[13] 潘騰.高分二號(hào)衛(wèi)星的技術(shù)特點(diǎn)[J].中國(guó)航天，2015（1）：3-9.

[14] 鄧書(shū)斌.ENVI遙感圖像處理方法[M].北京：高等教育出版社，2014：188.

[15] 蘇志強(qiáng).基于GF-2遙感影像的博湖縣監(jiān)督分類(lèi)方法對(duì)比研究[J].西部大開(kāi)發(fā)（土地開(kāi)發(fā)工程研究），2020，5（1）：1-5.

[16] 胡文敏，殷梓強(qiáng)，楊睿瀚，等.高分二號(hào)衛(wèi)星在湖南雷家山滑坡災(zāi)害調(diào)查中的應(yīng)用[J].衛(wèi)星應(yīng)用，2022（1）：50-54.

[17] 關(guān)元秀.eCognition基于對(duì)象影像分析教程[M].北京：科學(xué)出版社，2019：6.