高分辨率遙感影像土地利用變化檢測方法研究

王 琰，舒 寧，2，龔 龑

高分辨率遙感影像土地利用變化檢測方法研究

王 琰1，舒 寧1，2，龔 龑1

(1.武漢大學遙感信息工程學院，武漢 430079;2.武漢大學測繪遙感信息工程國家重點實驗室，武漢 430079)

提出一種利用高分辨率遙感影像進行土地利用變化檢測的方法。以土地利用圖為輔助數據，通過土地利用圖和遙感影像的配準套合，獲取影像像斑;同時，對遙感影像進行基于像素的監督分類，獲取概略的類別圖;再根據像斑內像素的類別編碼完成子像斑的劃分。以子像斑為影像分析的基本單位提取特征，以相關系數為相似性測度衡量不同時期子像斑的特征相似性，用ROC曲線(接受者操作特性曲線)代替經驗選取的方法自動獲取變化閾值，確定像斑是否發生變化。以武漢市區局部QuickBird 2002年和2005年多光譜影像、相同地區2002年1∶10 000土地利用圖為實驗數據進行了算法的實驗，結果顯示絕大部分的變化區域都可以被提取出來，實驗方法可行。

變化檢測;土地利用;面向對象;高分辨率;遙感;多源數據

0 引言

利用遙感影像的宏觀性、實時性等優勢進行國土資源的變化檢測，是遙感影像分析和應用領域的研究熱點。從技術路線層面而言，利用遙感技術進行土地利用變化檢測的方法可分為基于像素和基于像斑兩種方式。基于像素的影像分析受噪聲影響較大，能夠在影像上提取的特征信息較少;基于像斑的影像分析適用于高分辨率遙感影像的分析，為近幾年來常用的一種影像分析方法［1］。像斑是影像上的灰度同質區域，在絕大多數國內外文獻中稱之為“對象”，基于像斑的分析一般稱作“面向對象”或“基于對象”的方法。以像斑為影像分析的基本單位進行變化檢測，變化檢測結果更容易和矢量數據相聯系，從而更好地生成或更新矢量專題圖。目前，像斑往往通過各種影像分割的算法獲取，但合適的分割參數的設置以及分割尺度的確定往往較為困難。在變化檢測方面，衡量兩個時期像斑變化與否的變化閾值的設定，往往采用經驗的方式，通過多次嘗試和對結果的目視判斷獲取合適的閾值，這不利于技術路線的自動化。

本文提出一種利用高分辨率影像進行土地利用變化檢測的方法。首先，綜合基于像斑和基于像素的影像分析方式，通過土地利用圖和遙感影像的配準套合獲取影像像斑，用基于像素的監督分類方法，將影像數據分成簡單的基礎地類，生成不同時期影像的分類圖;然后分別以各分類圖為底圖，根據每個像斑內像素在分類圖上的類別標識，將通過配準套合生成的像斑進行再分割，獲得子像斑。以子像斑為后續分析的基本單位，通過特征提取、特征矢量相似性分析以及變化閾值的自動判定，獲得子像斑的變化信息;最后，將以不同時期分類底圖做參照獲取的子像斑變化檢測結果求并集，獲取最終的變化檢測結果。

1 像斑的獲取及像斑特征提取

1.1 像斑及子像斑的獲取

應用多源數據可以獲得有利于遙感影像解譯的更多信息。實驗區除遙感影像外的其他數據，可以以先驗知識的形式出現，直接參與到影像分析中去，這樣的思路也更為符合人工識別地物的流程和機理［2］;而變化檢測的結果，一般用于更新現有的矢量數據圖。因此在本文的變化檢測算法中，將現有土地利用圖與遙感影像配準套合，以土地利用圖圖斑矢量邊界為準，獲取影像像斑。由于土地利用圖上的圖斑劃分和生成須遵循其成圖規范和標準，并不是以影像上光譜信息為生成依據，僅憑配準套合得到的像斑并不能保證其光譜同質性，常會形成所謂的混合像斑，使后續的分析產生困難。因此有必要對套合得到的像斑進行再劃分，生成子像斑，從而保證各像斑內的光譜同質性。運用基于像素的監督分類方法，將地物分成水體、城區、裸地和綠地4大類，使得每一個像素都有一個類別編碼與其相對應。將分類影像與土地利用圖配準套合，這樣，對于每個像斑，根據其包含像素的基礎類別編碼進行再劃分。一個像斑內，具有同一類別編碼的像素生成一個子像斑。以不同時期的分類影像圖為參照進行像斑再分割，獲取的子像斑是不相同的。即以T1時期影像分類圖為參照獲取的T1，T2時期影像子像斑，和以T2時期影像分類圖為參照獲取的T1，T2時期影像子像斑不相同。子像斑生成的過程如圖1((b)圖上不同顏色代表再劃分后生成的子像斑)所示。

圖1 子像斑生成圖Fig.1 Subsegment generation graph

1.2 子像斑特征提取

通過計算子像斑內所有像素的統計值獲取子像斑的特征。對于實驗用QuickBird影像而言，獲取的特征為藍、綠、紅、近紅外4個波段影像的灰度均值、灰度標準差和灰度信息熵，共12個特征。

2 變化檢測及其閾值的自動獲取

2.1 變化檢測方法

像斑之間相似關系，可以看作是離散化的空間面模式之間的相似性關系［3］。這種關系的表達可以使用相關系數作為測度，用于衡量不同時期影像上相應子像斑特征矢量的相似性。假如某子像斑在基準時期的特征向量為(x1，x2，…，xn)，在待檢測時期的特征向量為(y1，y2，…，yn)，則兩個時期特征向量的相關系數ρ可以用

計算，式中n為特征向量維數。

2.2 基于ROC曲線的變化閾值自動判定方法

ROC曲線即接受者操作特性曲線，在醫學診斷學中有著較為廣泛的應用。對于二分類的問題，可以通過繪制ROC曲線，在使得誤檢率和漏檢率相同的情況下，獲得一個最合理的分類閾值。Yitzhaky和Peli曾將ROC曲線應用于影像邊緣提取中參數的選擇［5］;楊朝暉等也應用ROC曲線進行SAR影像的邊緣檢測［6］;Chen等應用ROC曲線進行基于衛星影像的地震前后城區損壞估計［7］，都獲取了較為理想的結果。本文引入上述ROC曲線，嘗試將其應用到變化檢測中，自動獲取變化閾值。

利用ROC曲線進行變化閾值獲取，其原理是設定多個不同的變化閾值，從而對應不同的變化檢測結果。對上述各結果進行分析，找到最合理的一個結果作為可能的真實變化結果，最接近可能真實結果的變化檢測結果，其對應的變化閾值則被認為是最合理的變化閾值。這種方法可以取代傳統的經驗閾值獲取方法，實現閾值的自動獲取。具體方法如下:

計算基準時期影像上各子像斑與其在待檢測時期影像上對應子像斑特征矢量的相關系數，并設定N個相關系數閾值Si(i=1，2，…，N)，相關系數大于閾值認為沒有發生變化，反之則判定子像斑發生了變化。這樣將會獲得N個影像變化檢測結果Ri(i=1，2，…，N)。每一個變化檢測結果為一個和像斑數目大小相等的數組，數組元素值為1表示像斑發生了變化，0則表示沒有發生變化。可能的真實結果CGTj(candidate ground truth，CGT)(j=1，2，…，N)，則根據Ri得出。CGTj的計算通過依次統計每個像斑在所有變化檢測結果上的變化次數獲取，對于第M個像斑，

這里 j(j=1，2，…，N)，稱為相關閾值(correspondence threshold)，因此可以獲得N個候選的真實變化結果。從上式可以看出，當某個像斑在所有N個變化檢測結果中被標記為1的次數較多，即認為發生變化的次數較多時，這個像斑發生變化的可能則較大，反之則小。CGTj仍然是一個和影像像斑數目相同的數組，數組元素值為1代表像斑發生變化，0則表示沒有變化。隨著j值的不斷增大，CGTj中數組元素值為1的越來越少，判定為發生變化的像斑數目則逐步減少。

將變化檢測視為一個二分類的問題，則4種判斷決策true positive(TP)(hit)，false positive(FP)(false alarm)，true negative(TN)(correct rejection)和false nagtive(FN)(miss)可能出現。對這4種判別決策的描述如表1所示。

表1 判別決策描述Tab.1 Describe discrimination decison

各個判別決策概率的計算公式為

根據式(3)—(6)可以計算出檢測概率(靈敏度)TPR和虛驚概率(1－特異度)FPR分別為

獲取了可能的真實變化結果之后，則可將所有發生變化的檢測結果與其相比較，根據上述算法可以看出，必然有一個變化檢測結果與已經獲取的真實變化結果相同，而這個變化檢測結果所對應的變化閾值則是需要求取的閾值。

3 實驗及結果分析

本文具體技術路線如圖2所示。

圖2 技術路線圖Fig.2 Technique flow chart

選用了武漢市局部2002年、2005年QuickBird多光譜影像，以及相同區域2002年1∶10 000土地利用圖。影像分辨率經重采樣后為2 m，實驗區大小為2 155×1 633像元。這一區域在2002—2005年間土地利用變化大，較為適合變化檢測實驗。應用監督分類方法，在ERDAS軟件環境下，分別將2002年和2005年影像分成城市、綠地、裸地和水體4大類。通過2002年土地利用圖GIS數據和兩個時期影像分別配準套合，共獲取像斑256個。

第一步，先以2002年影像監督分類圖為底圖對兩個時期像斑進行再劃分，得到兩時期影像上子像斑的數目為932個;應用相關系數對兩個時期像斑的特征矢量進行相似性評價，給定變化閾值獲取范圍，用2.2節介紹的基于ROC曲線的變化閾值自動判定方法即可求出最終的變化閾值(圖3)。最接近交點是第22個數據對，其對應的候選真實結果則被認為是可能的真實結果，對應的變化檢測閾值為 0.868。

圖3 ROC曲線結果Fig.3 ROC result

根據獲取的變化檢測閾值進行影像的變化檢測，將兩幅變化檢測實驗結果影像求并集，以獲取影像上所有可能發生變化的區域。實驗結果如圖4所示。

從圖3(a)中看出，最接近交點是第20個數據對，其對應的候選真實變化結果則被認為是可能的真實變化結果，對應的變化閾值為0.880。

第二步，以2005年影像監督分類圖為底圖對兩個時期像斑進行再劃分，得到兩個時期影像子像斑數目為801個，再重復上述的實驗步驟。獲取真實結果和變化檢測閾值的ROC曲線如圖3(b)所示，

圖4 實驗區影像及實驗結果Fig.4 Image of study area and the change result

通過人工目視解譯與實驗結果相比較，發現絕大部分發生變化的部分都能夠被檢驗并提取出來。

4 結論

1)提出了一種高分辨率遙感影像土地利用變化檢測的方法。綜合以像斑和以像素為單位進行影像分析的優勢，應用多源數據獲取影像像斑并進行再分割，在此基礎上運用相關系數衡量像斑特征的相似性，并運用ROC曲線自動獲取變化閾值，得到變化檢測結果。

2)實驗結果表明，文中提出的方法簡單、有效、可行，且具有一定的通用性。在實驗過程中發現，變化閾值數組范圍的設定會影響可能的真實變化結果的獲取，繼而影響變化閾值的判定。怎樣選取合適的變化閾值數組范圍，應作進一步探討。

［1］ Blaschke T.Object Based Image Analysis for Remote Sensing［J］.ISPRS Journal of Photogrammetry and Remote Sensing，2010，65(1):2－16.

［2］ 舒 寧.關于遙感影像處理分析的理論與方法之若干問題［J］.武漢大學學報:信息科學版，2007，32(11):1007－1010.

［3］ 毛政元，李 霖.空間模式的測度及其應用［M］.北京:科學出版社，2004:74－75.

［4］ 張劍清，潘 勵，王樹根.攝影測量學［M］.武漢:武漢大學出版社，2003:109.

［5］ Yitzhaky Y，Peli E.A Method for Objective Edge Detection Evaluation and Detector Parameter Selection［J］.IEEE Transactions on Pattern Analysis and Machine Intelligence，2003，25(8):1027 －1033.

［6］ 楊朝暉，陳 鷹.基于ROC融合準則的SAR邊緣檢測算法［J］.光電子·激光，2010，21(7):1053 －1057.

［7］ Chen Z Q，Hutchinson T C.Urban Damage Estimation Using Statistical Processing of Satellite Images［J］.Journal of Computing in Civil Engineering，2007，21(3):187 －199.

A Study of Land Use Change Detection Based on High Resolution Remote Sensing Images

WANG Yan1，SHU Ning1，2，GONG Yan1
(1.School of Remote Sensing and Information Engineering，Wuhan University，Wuhan 430079，China;2.State Key Laboratory of Information Engineering in Surveying，Mapping and Remote Sensing，Wuhan University，Wuhan 430079，China)

An approach to land use change detection by using high resolution remote sensing images is put forward in this paper.With the help of GIS land use map，image objects can be obtained by the matching of land use map and remote sensing images in the same region.Meanwhile pixel－based supervised classification is conducted for each image so that each pixel has its own class code.Then image subsegments can be obtained based on the image segment and the class code of each pixel within it.Image subsegments can be regarded as the basic units for feature extraction.Correlation coefficient is used for detecting changes between the images gotten from different time periods，and instead of the empirical selection，the change threshold is founded automatically by using ROC curve(receiver operating characteristic curve).Two multispectral Quickbird images obtained in 2002 and 2005 respectively and a 1 ∶10 000 land use map of 2002 in the same region were used in the experiment.This study area is located in Wuhan City and the result shows that most land use changes can be detected，and hence this approach is effective.

change detection;land use;object based;high resolution;remote sensing;multisource data

TP 79

A

1001－070X(2012)01－0043－05

10.6046/gtzyyg.2012.01.08

2011－05－21;

2011－07－01

湖北省自然科學基金重大項目(編號:2006ABD003)及中央高校基本科研業務費專項資金(編號:3101009)共同資助。

王 琰(1984－)，女，在讀博士，研究方向為遙感影像解譯。E－mail:wangyan.jenny@163.com。

(責任編輯:李 瑜)