王 巖，王曉青，竇愛霞，孫素梅，付 聰

（1. 遼寧省地震局，遼寧 沈陽 110034；2. 中國(guó)地震局地震預(yù)測(cè)研究所，北京 100036）

0　引言

強(qiáng)震發(fā)生后，評(píng)估地震災(zāi)害的范圍程度是最急迫的需求，評(píng)估過程中，對(duì)災(zāi)害進(jìn)行準(zhǔn)確、快速的提取能夠第一時(shí)間掌握全局，為進(jìn)一步的救災(zāi)工作提供最有效的依據(jù)[1]。在震害快速提取中，應(yīng)用遙感衛(wèi)星影像已成為主要手段之一，提高災(zāi)害提取的精度與速度，是研究的關(guān)鍵問題[2-3]。隨著遙感技術(shù)的日趨成熟，獲取的各類影像分辨率越來越高，其包含的空間信息越來越豐富，地物的幾何結(jié)構(gòu)、紋理信息突出，運(yùn)算量也逐步增大，造成基于像元的分類方法越來越難以適應(yīng)需求[4]。

面向?qū)ο蟮倪b感圖像分類方法，以影像對(duì)象代替像元，是一種運(yùn)用多尺度分割實(shí)現(xiàn)圖像信息的融合再分類的方法[5-6]。面向?qū)ο蟮姆诸愔?，?yīng)用信息不僅局限于光譜信息，更側(cè)重于對(duì)象與對(duì)象之間的紋理、形狀、關(guān)系等。在精度上，由于應(yīng)用了更接近人類思維模式的模糊邏輯分類而更為精確；在速度上，采用影響對(duì)象為基礎(chǔ)計(jì)算單位，大大減少了計(jì)算量，提高了速度[7-8]。

本文以汶川地震中破壞嚴(yán)重的北川縣城作為研究對(duì)象，應(yīng)用面向?qū)ο蠓诸惙椒?，?shí)現(xiàn)震后災(zāi)害區(qū)域的快速提取。

1　研究區(qū)域與數(shù)據(jù)

研究選取2008年5月12日汶川8.0級(jí)地震后，烈度為Ⅺ的北川縣城及周邊地區(qū)作為研究區(qū)域。北川位于四川盆地西北部，境內(nèi)地質(zhì)構(gòu)造復(fù)雜，震后二次災(zāi)害嚴(yán)重，是汶川地震后受災(zāi)最為嚴(yán)重的區(qū)域之一。

研究數(shù)據(jù)：福衛(wèi)二號(hào)衛(wèi)星多光譜影像，分辨率8m，獲取時(shí)間2008年5月14日。福衛(wèi)二號(hào)衛(wèi)星由中國(guó)臺(tái)灣制造，2004年5月美國(guó)發(fā)射，多光譜波段包括：0.45～0.52微米（藍(lán)），0.52～0.60微米（綠），0.63～0.69微米（紅），0.76～0.90微米（近紅外）四個(gè)波段。

2　面向?qū)ο笥跋穹诸?/h2>

2.1　目標(biāo)類別

研究主要目的是快速提取震后的受災(zāi)區(qū)域，即地震發(fā)生后因嚴(yán)重災(zāi)害而引起變化的區(qū)域，稱為“熱區(qū)”。熱區(qū)地物類型包括建筑物嚴(yán)重?fù)p毀區(qū)域和地質(zhì)災(zāi)害區(qū)域。

根據(jù)災(zāi)害評(píng)估內(nèi)容類別的需求和北川縣城及周邊地區(qū)的地物特征，在面向?qū)ο蠓诸愔袆澐至怂膫€(gè)類別：植被、水體、熱區(qū)和陰影。其中熱區(qū)類別是震前震后變化巨大的地物范圍；水體類別是震前震后未發(fā)生變化的水體部分，不包括堰塞湖；陰影主要是各山體之間的陰影部分。四個(gè)類別相互獨(dú)立，為并列關(guān)系。

2.2　多尺度分割

面向?qū)ο蠓诸愂菑膯我幌裨_始，采用自下而上區(qū)域合并的方法形成影像對(duì)象，再針對(duì)眾多對(duì)象進(jìn)行分類的方法。多尺度影像分割技術(shù)，是面向?qū)ο蠓诸惖暮诵膬?nèi)容，是完成影像對(duì)象分割的基礎(chǔ)。多尺度影像分割包含兩個(gè)重要選擇：分割尺度和同質(zhì)性系數(shù)。分割尺度的選擇是根據(jù)地物類別的顆粒大小決定分割單元大??；同質(zhì)性系數(shù)是根據(jù)地物光譜和形狀特征的不同選擇不同的要素權(quán)重。

分割尺度的大小代表對(duì)象異質(zhì)性的閾值，用無量綱的數(shù)值表示，數(shù)值越大顆粒越大。異質(zhì)性是由對(duì)象的光譜和形狀差異確定的，光譜異質(zhì)性、質(zhì)密度異質(zhì)性共同提供形狀信息，與顏色信息相互依存，組成同質(zhì)性系數(shù)，如圖1所示。

光譜異質(zhì)性的計(jì)算公式如下：

圖1　同質(zhì)性系數(shù)的組成Fig.1　 Composition of homogeneity coefficient

光滑度異質(zhì)性的計(jì)算公式為：

質(zhì)密度異質(zhì)性計(jì)算公式為：

其中，n是組成影像對(duì)象的像元個(gè)數(shù)。

形狀因子綜合了質(zhì)密度異質(zhì)性與光滑度異質(zhì)性兩種信息，計(jì)算公式為：

同質(zhì)性系數(shù)綜合了光譜異質(zhì)性與形狀異質(zhì)性兩種信息，計(jì)算公式為：

如圖2所示，不同的分割尺度下，地物分割成不同大小的單元集合。尺度為10時(shí)，地物分割太過于細(xì)碎；尺度為20時(shí)，水體、居民區(qū)分割最為合理；尺度為30時(shí)，植被分割最為合理；尺度為40、50時(shí)，過于概括而無法有效區(qū)分。為了最大程度正確分類各種地物，研究選擇兩個(gè)不同尺度進(jìn)行交叉分類：尺度為20時(shí)分割水體、熱區(qū)，尺度為30時(shí)分割植被、陰影。

圖2　不同尺度參數(shù)下的分割Fig.2　 Segmentation result of different scales

同時(shí)，從圖2中可見，同一分割尺度下，分割后的單元并不是大小一致的，而是隨著地物特征變化的，這些變化取決于同質(zhì)性系數(shù)的選擇。由于北川縣城的災(zāi)害快速提取中對(duì)形狀特征的要求不是很突出，在分割過程中形狀因子的權(quán)重選擇較小，其中植被與陰影區(qū)受形狀因素影響程度更低，具體歸一化權(quán)重選擇如表1所示。

表1　尺度參數(shù)列表

另外，在多光譜福衛(wèi)二號(hào)衛(wèi)星影像中有四個(gè)波段數(shù)據(jù)，其中不同波段對(duì)不同地物的反映不同。例如近紅外波段對(duì)植被、水體的反映就較為敏感，在這種情況下，研究通過適當(dāng)調(diào)整不同波段在不同地物分類中的權(quán)重進(jìn)行篩選。如表2所示，為實(shí)驗(yàn)獲得的最佳權(quán)重。

表2　波段權(quán)重表

3　分類結(jié)果與效能評(píng)價(jià)

3.1　分類結(jié)果

面向?qū)ο笥跋穹指钪锌梢愿鶕?jù)每一種地物的特征，分別選擇不同的分類標(biāo)準(zhǔn)，在不同的尺度下建立分類流程。如圖3，為北川縣城及周邊區(qū)域震害快速提取的工作示意圖，其中用到了比值植被指數(shù)（RVI）、差值環(huán)境指數(shù)（DVIEVI）、歸一化植被指數(shù)（NDVI）三個(gè)光譜特征來分別分類植被、熱區(qū)和水體。

RVI是比值植被指數(shù)，其公式為：

圖3　北川縣城及周邊區(qū)域震害快速提取Fig.3　 Rapid seismic damage extraction of Beichuan County and its surrounding area

其中NIR為近紅外波段的亮度值，R為紅光波段的亮度值。

DVIEVI是差值環(huán)境指數(shù)，其公式為：

NDVI是歸一化植被指數(shù)，其公式為：

選取的分類特征是取自大量不同類別的特征，選取過程中經(jīng)過對(duì)圖像的直方圖分析調(diào)整，綜合分析得到。

按照上述分類流程，得到了分類結(jié)果，如圖4（a）所示，粉色為熱區(qū)、綠色為植被、藍(lán)色為水體、黑色為陰影。面向?qū)ο蟮姆诸惤Y(jié)果是矢量化的結(jié)果，不再是單個(gè)像元分離，更易于進(jìn)行各種變換和統(tǒng)計(jì)。將提取的熱區(qū)用矢量化圖層表示，與原始影像疊合，結(jié)果如圖4（b），可見各個(gè)類別的劃分非常明確，沒有細(xì)碎的噪聲，分類結(jié)果與目視解譯基本吻合。

圖4　分類結(jié)果圖與熱區(qū)Fig.4　 Classification result and the hot-area

3.2　效能評(píng)價(jià)

面向?qū)ο蟮倪b感分類結(jié)果提供了矢量分類數(shù)據(jù)，能夠快速給出災(zāi)害范圍的區(qū)域、面積等定量化結(jié)果，便于與基礎(chǔ)生命線數(shù)據(jù)庫中的矢量信息進(jìn)行交叉計(jì)算，第一時(shí)間提供震后的評(píng)估數(shù)據(jù)。

與人工解譯結(jié)果比對(duì)分析發(fā)現(xiàn)，除了山體滑坡信息中較為細(xì)碎的小滑坡和影像中云的影響，判讀結(jié)果基本一致。去除云的影響，受災(zāi)面積正確率達(dá)90%、水體95%、植被85%。在面向?qū)ο蠓诸愔?，同一類型遙感影像數(shù)據(jù)在大范圍宏觀災(zāi)害提取中選取的各參數(shù)信息具有可遺傳性，試用廣泛，因此可以在大區(qū)域計(jì)算中減少重復(fù)工作，大大提高災(zāi)害提取速度。

21天結(jié)束后，共99人參加打卡活動(dòng)，總打卡次數(shù)1218次。根據(jù)讀者在21天的打卡表現(xiàn)，最終評(píng)選出31名全勤打卡獎(jiǎng)及28名優(yōu)質(zhì)打卡獎(jiǎng)，在寫書評(píng)的讀者中評(píng)選出五名優(yōu)秀書評(píng)獎(jiǎng)。開學(xué)后我們邀請(qǐng)三十名寒假打卡的讀者組織了“寒假21天閱讀習(xí)慣養(yǎng)成”活動(dòng)的線下交流會(huì)暨頒獎(jiǎng)儀式?；顒?dòng)中讀者交流寒假打卡的心得體會(huì)，邀請(qǐng)兩位同學(xué)就“如何挑選你要讀的書”“如何抽出時(shí)間讀更多的書”以及“讀書給你帶來了什么”這三個(gè)問題分享了自己的心得和感悟。

4　結(jié)論

在震后快速評(píng)估中，快速、準(zhǔn)確的提取有效災(zāi)害信息是評(píng)估依據(jù)的基礎(chǔ)。本文應(yīng)用面向?qū)ο蟮姆诸惙椒?，快速提取了汶川地震發(fā)生后北川縣城及周邊的宏觀災(zāi)害區(qū)域。提取過程中綜合運(yùn)用了遙感影像數(shù)據(jù)的多光譜特征和紋理特征，模擬人類思維的模式，使得分類結(jié)果最大程度接近于人工解譯結(jié)果，并大大提高了解譯速度。同時(shí)，面向?qū)ο蟮姆诸惤Y(jié)果得到的是矢量化的數(shù)據(jù)信息，相對(duì)于傳統(tǒng)的基于像元的分類結(jié)果，更便于定量化信息的分析和后續(xù)應(yīng)用計(jì)算。

參考文獻(xiàn)：

[1] 王曉青， 丁香. 地震現(xiàn)場(chǎng)災(zāi)害損失評(píng)估地理信息系統(tǒng)[M]. 北京： 地震出版社， 2002.

[2] 王曉青， 魏成階， 苗崇剛， 等 . 震害遙感快速提取研究—以2003年2月24日巴楚－伽師6.8級(jí)地震為例 [J]. 地學(xué)前緣， 2003， 10（U08）： 285-291.

[3] 王巖， 李彤霞， 李秀麗， 等 . 基于 MODIS 的日本 9.0級(jí)地震前后溫度變化研究[J]. 防災(zāi)減災(zāi)學(xué)報(bào)， 2015，31（4）： 14-19.

[4] 竇愛霞. 震害遙感圖像的變化檢測(cè)技術(shù)研究[D]. 山東： 山東科技大學(xué)， 2003.

[6] 孫曉霞， 張繼賢， 劉正軍. 利用面向?qū)ο蟮姆诸惙椒◤腎KONOS全色影像中提取河流和道路[J]. 測(cè)繪科學(xué)， 2006， 31（1）： 62-63.

[7] 王巖， 王曉青， 竇愛霞. 面向?qū)ο筮b感分類方法在汶川地震震害提取中的應(yīng)用[J]. 地震， 2009， 29（3）：54-60.

[8] 王巖. 面向?qū)ο髨D像處理方法在震害信息提取中的應(yīng)用研究[D]. 北京： 中國(guó)地震局預(yù)測(cè)研究所， 2009.