城市建筑物永久散射體識(shí)別策略研究

楊魁，劉俊衛(wèi)（天津市測(cè)繪院，天津 300381）

城市建筑物永久散射體識(shí)別策略研究

楊魁?，劉俊衛(wèi)
（天津市測(cè)繪院，天津 300381）

針對(duì)城市建筑物永久散射體識(shí)別的應(yīng)用需求，在有效分析建筑物永久散射體原理基礎(chǔ)上，根據(jù)幅離差閾值法、子視相關(guān)法、相干系數(shù)法應(yīng)用于建筑物永久散射體識(shí)別的優(yōu)缺點(diǎn)，采用多種方法相結(jié)合的永久散射體識(shí)別策略來實(shí)現(xiàn)建筑物的監(jiān)測(cè)。并以天津市不同類型建筑物為例，從密度、點(diǎn)位分布、質(zhì)量等方面分析和對(duì)比永久散射體識(shí)別策略相對(duì)于其他單一方法的優(yōu)勢(shì)，驗(yàn)證了識(shí)別策略的有效性。

建筑物；永久散射體識(shí)別；InSAR

1 引 言

隨著天津社會(huì)經(jīng)濟(jì)的快速發(fā)展，城市建設(shè)日益增多。大量施工建設(shè)會(huì)對(duì)已有的建筑物的安全性產(chǎn)生影響，且新建的建筑物在前幾年內(nèi)也會(huì)存在一定的沉降，因此為實(shí)現(xiàn)對(duì)建筑物的安全鑒定，需要采用高精度的方法對(duì)建筑物進(jìn)行監(jiān)測(cè)。目前常用的方法主要為水準(zhǔn)測(cè)量，其具有高精度特性，但需要耗費(fèi)大量人力物力。隨著近些年來衛(wèi)星遙感技術(shù)的發(fā)展，InSAR開始作為一種重要的沉降監(jiān)測(cè)手段正逐步走向應(yīng)用，但是其針對(duì)建筑物監(jiān)測(cè)上的應(yīng)用研究較少，專門針對(duì)建筑物永久散射體的識(shí)別也很少見諸于報(bào)告。

目前對(duì)永久散射體的識(shí)別，A.Ferretti、F.Rocca等針對(duì)永久散射體目標(biāo)的識(shí)別探討了相關(guān)的探測(cè)與識(shí)別方法，主要包括幅度閾值法、幅度離差閾值法、后向散射強(qiáng)度離差閾值法等［1，2］；A.Hooper提出相位穩(wěn)定性分析法［3］，N.Adam 提出信噪比閾值法［3］，Urs Wegmüller等提出了點(diǎn)目標(biāo)的識(shí)別方法［4，5］，Mora O等提出相干系數(shù)閾值法［3］，這些方法各有優(yōu)缺點(diǎn)。本文在分析建筑物永久散射體原理基礎(chǔ)上，采用幅度離差閾值法、點(diǎn)目標(biāo)檢測(cè)法、相干系數(shù)法等多種方法相結(jié)合的永久散射體識(shí)別策略來實(shí)現(xiàn)建筑物的高密度監(jiān)測(cè)，并以天津市市內(nèi)六區(qū)的多景TerraSAR數(shù)據(jù)為例分析該識(shí)別策略的有效性。

2 主要研究?jī)?nèi)容

2.1 建筑物永久散射體原理分析

在SAR影像上建筑物一般具有較強(qiáng)回波信號(hào)。以一小區(qū)的建筑物群為例分析，在升軌方式下，建筑物的西面和南面表現(xiàn)為強(qiáng)亮度，屋頂有部分亮點(diǎn)，建筑物的東面和北面則表現(xiàn)為暗色調(diào)（如圖1所示），不同位置亮度產(chǎn)生差異的原因主要是角反射器效應(yīng)。當(dāng)?shù)匚锬繕?biāo)具有兩個(gè)互相垂直的光滑表面（房屋墻面與地表面）或有是三個(gè)互相垂直的光滑表面（建筑物的凹不與地表面）時(shí)就是所謂角反射器，它有二面角反射器和三面角角反射器之分。當(dāng)雷達(dá)波束遇到這種目標(biāo)時(shí)，由于角反射器每個(gè)表面的鏡面反射，使波束最后反轉(zhuǎn)180°方向向來波方向傳播，這樣就產(chǎn)生各條射線在反射回去的時(shí)候方向相同，相位也相同，故而信號(hào)互相增強(qiáng)的現(xiàn)象，致使回波信號(hào)極強(qiáng)，從而可以視為理想的永久散射體被有效識(shí)別出來。

圖1　建筑物解譯標(biāo)志

2.2 建筑物永久散射體識(shí)別策略

針對(duì)建筑物自身特點(diǎn)來充分研究已有永久散射體識(shí)別方法的基礎(chǔ)上，采用高可信度、高密度的建筑物永久散射體識(shí)別策略來實(shí)現(xiàn)建筑物的有效監(jiān)測(cè)。

（1）幅度離差閾值法

作為永久散射體識(shí)別技術(shù)中的經(jīng)典方法，其主要思想是基于大數(shù)據(jù)量的SAR影像下幅度離差和相位偏差間的統(tǒng)計(jì)關(guān)系來選擇高可信度的永久散射體。對(duì)經(jīng)過配準(zhǔn)處理后的時(shí)序SAR數(shù)據(jù)，計(jì)算同一點(diǎn)在多幅影像上的均差、方差，將兩者之比作為永久散射體識(shí)別方法的測(cè)度，選擇大于一定閾值的點(diǎn)為建筑物永久散射體。該方法能有效反映建筑物目標(biāo)在時(shí)間序列上的穩(wěn)定性特征，但是會(huì)錯(cuò)誤地引入建筑物群間存在的陰影。

? 收稿日期:2015—12—03

作者簡(jiǎn)介:楊魁（1987—），男，碩士，工程師，主要從事影像解譯、SAR應(yīng)用工作。

基金項(xiàng)目：國(guó)土資源部公益性行業(yè)科研專項(xiàng)（201311045）

（2）點(diǎn)目標(biāo)檢測(cè)法

通過上述對(duì)建筑物永久散射體原理分析可知，其由于角反射器效應(yīng)具有高亮度，因此可以將之歸為點(diǎn)目標(biāo)，它在獲取時(shí)間內(nèi)的回波信號(hào)保持恒定，相位特征也很穩(wěn)定，可以被識(shí)別為有效的永久散射體。在數(shù)據(jù)處理過程中，通過將SAR影像進(jìn)行子孔徑處理來獲取光譜相關(guān)性，從而可以利用子視相關(guān)理論來識(shí)別出屬于建筑物的永久散射體［6，7］。該方法物理意義明確，但是由于噪聲影響，其在時(shí)間序列上的穩(wěn)定性不夠。

（3）相干系數(shù)法

經(jīng)過統(tǒng)計(jì)分析可認(rèn)為具有SAR影像上高空間相干性的點(diǎn)對(duì)應(yīng)著地面上的建筑物等高相干區(qū)域。因此選擇合適的窗口采用相干系數(shù)來評(píng)估同一目標(biāo)在不同SAR影像上的空間類似程度以生成相干圖序列，并進(jìn)一步計(jì)算平均相干系數(shù)以選擇相干系數(shù)較高的建筑物點(diǎn)目標(biāo)。該方法可以有效剔除建筑物的陰影，提高可信度，但受估算窗口的影響容易引入低相干目標(biāo)。

（4）永久散射體識(shí)別策略

通過上述分析可知，單一的永久散射體識(shí)別方法用于永久散射體的識(shí)別都有其局限性。因此本文針對(duì)InSAR在建筑物監(jiān)測(cè)的應(yīng)用需要，在上述已有的永久散射體目標(biāo)識(shí)別理論基礎(chǔ)上采用了永久散射體的識(shí)別策略，使之適應(yīng)于建筑物永久散射體目標(biāo)的識(shí)別，識(shí)別流程如圖2所示。首先基于幅度離差閾值法來選擇時(shí)間序列上的變化小的永久散射體，以減少時(shí)間失相干的影響；其次基于點(diǎn)目標(biāo)檢測(cè)方法選擇高亮度的點(diǎn)目標(biāo)，減少空間失相干的影響；最后用相干系數(shù)閾值法來快速剔除影像中那些失相關(guān)非常嚴(yán)重的非永久散射體目標(biāo)，確保最終獲取高密度、高質(zhì)量的永久散射體。

本文所采用的建筑物永久散射體識(shí)別策略通過將振幅離差法、點(diǎn)目標(biāo)檢測(cè)方法可以選擇出受時(shí)間失相干或空間失相干影響較小的永久散射體候選點(diǎn)，大大提高單一監(jiān)測(cè)方法的永久散射體密度，實(shí)現(xiàn)建筑物的高密度監(jiān)測(cè)。在此基礎(chǔ)上采用相干系數(shù)閾值法則有效剔除建筑物陰影導(dǎo)致的錯(cuò)誤點(diǎn)，提高上述識(shí)別方法的可靠性，實(shí)現(xiàn)建筑物的高質(zhì)量監(jiān)測(cè)。因此本文所采用的多種方法相結(jié)合的識(shí)別策略為建筑物的高精度監(jiān)測(cè)提供有效的基礎(chǔ)，下面以天津市的SAR影像集和分析結(jié)果為例進(jìn)行翔實(shí)的驗(yàn)證分析。

圖2　永久散射體識(shí)別流程圖

3 實(shí)驗(yàn)及結(jié)果分析

為驗(yàn)證本項(xiàng)目永久散射體識(shí)別策略的有效性，利用11景X波段的TerraSAR-X的單視復(fù)影像，以振幅離差閾值法、子視相關(guān)法、相干系數(shù)法等三種方法為基礎(chǔ)進(jìn)行永久散射體的識(shí)別對(duì)比實(shí)驗(yàn)。

振幅離差閾值法識(shí)別的永久散射體如圖3所示，識(shí)別出來的永久散射體共174.6萬。從圖中可以看出，永久散射體分布幾乎都是獨(dú)立存在的，而且多數(shù)均分布在建筑物密集的中心城區(qū)，沿人工建筑物、構(gòu)筑物以及道路兩側(cè)分布，表現(xiàn)為高亮度的點(diǎn)或線。

子視相關(guān)法取的永久散射體如圖4所示，識(shí)別出來的永久散射體共280.4萬。從圖中可以看出，永久散射體分布比較均勻，點(diǎn)密度較高，尤其是建筑物的監(jiān)測(cè)點(diǎn)密度較高，但是由于其算法的限制，其在建筑物陰影等不應(yīng)該存在永久散射體的區(qū)域識(shí)別出一定數(shù)量的永久散射體。

圖3　振幅離差法識(shí)別永久散射體圖

圖4　子視相關(guān)法識(shí)別永久散射體圖

對(duì)比振幅離差閾值法和子視相關(guān)法識(shí)別的永久散射體分布圖，兩者在整體上表現(xiàn)出分布的一致性，但子視相關(guān)法獲取的永久散射體分布相對(duì)較均勻。局部進(jìn)行分析，兩者識(shí)別的點(diǎn)的數(shù)目和分布存在一定的差異，圖5和圖6所示的天津市標(biāo)志性建筑之一的津塔，振幅離差法識(shí)別的永久散射體在其投影面上高密度、均勻分布，而子視相關(guān)法識(shí)別的永久散射體則數(shù)量?jī)H有振幅離差法的1/3，且在投影上主要分布在頂部和底部，中間區(qū)域點(diǎn)較少，這說明振幅離差法在高層建筑物永久散射體識(shí)別上具有優(yōu)勢(shì)。但是對(duì)其周圍的眾多低中層建筑物的永久散射體分布進(jìn)行分析，振幅離差法獲取的永久散射體數(shù)量明顯少于子視相關(guān)法，且在很多建筑物上沒有永久散射體，表明中低層建筑物永久散射體識(shí)別上子視相關(guān)法更有效。

圖5　振幅離差法識(shí)別永久散射體局部圖

圖6　子視相關(guān)法識(shí)別永久散射體局部圖

因此為了實(shí)現(xiàn)對(duì)實(shí)驗(yàn)區(qū)域內(nèi)不同類型建筑物的全面檢測(cè)，需要將兩種方法進(jìn)行組合處理來實(shí)現(xiàn)優(yōu)勢(shì)互補(bǔ)。最終共識(shí)別出永久散射體375.3萬，點(diǎn)密度增大將近一倍，分布也更加均勻。

但是在上述識(shí)別出來的永久散射體中，存在部分錯(cuò)誤點(diǎn)。如陰影處的一些點(diǎn)，由于平均幅度較低，其統(tǒng)計(jì)特性較為穩(wěn)定，容易被歸為永久散射體。因此本項(xiàng)目在此基礎(chǔ)上利用相關(guān)系數(shù)法來剔除明顯的錯(cuò)誤點(diǎn)（圖7所示）39.7萬個(gè)左右，最終選擇候選點(diǎn)共有335.6萬左右，分布如圖8所示。

從點(diǎn)密度、分布、質(zhì)量（假設(shè)最終識(shí)別成果為正確值）等方面對(duì)上述分析過程進(jìn)一步詳細(xì)分析，如表1所示。采用本文所采用的永久散射體識(shí)別策略獲取的點(diǎn)密度明顯增加，與經(jīng)典的振幅離差方法對(duì)比，密度將近增加1倍；從方法適用性分析，振幅離差適用于高層建筑物，子視相關(guān)法在中、低層建筑物上優(yōu)勢(shì)明顯，永久散射體識(shí)別策略則充分利用兩者的優(yōu)勢(shì)，適用于絕大部分建筑物；從質(zhì)量進(jìn)行分析，永久散射體識(shí)別策略采用相關(guān)系數(shù)剔除建筑物陰影的影響，可靠性最高。

圖7　錯(cuò)誤永久散射體示意圖

圖8　最終識(shí)別永久散射體分布圖

四種永久散射體目標(biāo)識(shí)別方法的性能比較 表1

4 結(jié) 論

本文主要針對(duì)InSAR在城市建筑物沉降監(jiān)測(cè)中的應(yīng)用需求，重點(diǎn)對(duì)建筑物永久散射體的識(shí)別策略進(jìn)行研究。根據(jù)建筑物永久散射體原理，選擇適用于建筑物永久散射體識(shí)別的幅離差閾值法、子視相關(guān)法、相干系數(shù)法，在對(duì)這些方法原理、優(yōu)缺點(diǎn)分析的基礎(chǔ)上，采用這三種方法相結(jié)合的永久散射體識(shí)別策略對(duì)建筑物監(jiān)測(cè)，實(shí)現(xiàn)優(yōu)勢(shì)互補(bǔ)。并以天津市市內(nèi)六區(qū)的不同類型建筑物為例，分析和對(duì)比不同永久散射體識(shí)別方法與本文識(shí)別策略下永久散射體的密度、分布、質(zhì)量等特征，驗(yàn)證了永久散射體識(shí)別策略能夠?qū)崿F(xiàn)高密度、高可信度、多類型建筑物永久散射體識(shí)別。

［1］A.Ferretti，C.Prati，and F.Rocca.Permanent scatterers in SAR interferometry［J］.IEEE Trans.Geosci.Remote Sensing，2000，38：2202～2212.

［2］Ferretti A，Prati C，Rocca F.Nonlinear Subsidence Rate Esitimation Using Permanent Scatterers in Differential SAR Interferometry［J］.2001，38（5）：2202～2212.

［3］范洪冬.InSAR若干關(guān)鍵算法及其在地表沉降監(jiān)測(cè)中的應(yīng)用研究［D］.徐州：中國(guó)礦業(yè)大學(xué)，2010.

［4］盧麗君.基于時(shí)序SAR影像的地表形變監(jiān)測(cè)方法研究與應(yīng)用［D］.武漢：武漢大學(xué)，2008.

［5］U.Wegmüller，D.Walter，V.Spreckels，and Charles L.Werner.Nonuniform Ground Motion Monitoring With TerraSARX Persistent Scatterer Interferometry［J］.IEEE Trans. Geosci.Remote Sensing，2010，48（2）：895～904.

［6］Yan Wang，Daqing，Ge.Surface Subsidence Monitoring With Coherent Point Target Sar Ineterferometry［C］，IGARSS 2008.

［7］丁尚起，楊魁，陳楚.PSC選擇策略及在京津高鐵沉降監(jiān)測(cè)中的應(yīng)用研究［J］.城市勘測(cè)，2013，6：79～83.

The Study of Permanent Scatter Identification Strategy in Urban Building

Yang Kui，Liu Junwei
（Tianjin Institute of Surveying and Mapping，Tianjin 300381，China）

In order to apply InSAR to monitor urban building，a Permanent Scatter identification strategy is proposed for the building.First，the PS principle in the building is effectively analyzed.And the theory of amplitude deviation threshold method，sub-look coherent method，coherent coefficient method are explained to get corresponding advantage and disadvantage in the building.Then a new PS identification strategy utilizing these methods is advanced to solve the building monitoring problem.Different types of building in Tianjin are used to verify the validity of this identification strategy.Different PS identification method are compared from density，distribution and quality to get a conclusion that the PS identification strategy is very effectively.

building；permanent scatter identification；InSAR

1672-8262（2016）02-84-04中圖分類號(hào)：P237，TP753

A