段 偉，呂孝雷

(中國(guó)科學(xué)院電子學(xué)研究所 中國(guó)科學(xué)院空間信息處理與應(yīng)用系統(tǒng)技術(shù)重點(diǎn)實(shí)驗(yàn)室， 北京 100190； 中國(guó)科學(xué)院大學(xué)， 北京 100049)

時(shí)間序列InSAR技術(shù)[1-2]，可以對(duì)區(qū)域和人工建筑進(jìn)行長(zhǎng)時(shí)間觀測(cè)，以較低的成本，獲取區(qū)域和建筑大面積、高精度的微小形變信息，在地理測(cè)量領(lǐng)域有著廣泛應(yīng)用。而且目前由于高分辨率SAR影像的使用，時(shí)間序列InSAR技術(shù)甚至可以對(duì)單個(gè)具體建筑進(jìn)行形變分析。相比于傳統(tǒng)測(cè)量手段，這種空間遙感技術(shù)有許多優(yōu)點(diǎn)，在橋梁形變檢測(cè)領(lǐng)域有著廣闊的應(yīng)用前景，對(duì)橋梁安全監(jiān)測(cè)有著重大意義。但目前在行業(yè)內(nèi)，大型斜拉橋的時(shí)間序列InSAR形變檢測(cè)仍然面臨較大困難，現(xiàn)有處理算法無(wú)法獲取正確的大型橋梁形變結(jié)果，極大地限制了時(shí)間序列InSAR技術(shù)在橋梁領(lǐng)域的應(yīng)用。

本文以洞庭湖大橋?yàn)槔钊敕治瞿壳皶r(shí)間序列InSAR處理方法在大型斜拉橋上遇到的問(wèn)題，證明該問(wèn)題主要產(chǎn)生于相位解纏過(guò)程。由于大橋所在區(qū)域殘差點(diǎn)分布比較復(fù)雜，目前常規(guī)的稀疏網(wǎng)絡(luò)MCF相位解纏方法[3-7]在橋塔附近的殘差點(diǎn)會(huì)產(chǎn)生錯(cuò)誤的連接，導(dǎo)致該區(qū)域產(chǎn)生的枝切線會(huì)穿過(guò)橋面，帶來(lái)較大的解纏誤差。因此，本文提出一種新的相位解纏方法，首先構(gòu)建Delaunay約束三角網(wǎng)絡(luò)，減少殘差點(diǎn)數(shù)目；然后為避免出現(xiàn)較長(zhǎng)的枝切線，同時(shí)阻止枝切線穿過(guò)橋面，對(duì)殘差點(diǎn)按距離進(jìn)行聚類，從而將區(qū)域分成若干個(gè)殘差點(diǎn)聚集區(qū)；最后對(duì)每塊中的正負(fù)殘差點(diǎn)進(jìn)行連接，為平衡孤立殘差點(diǎn)，引入一部分接地點(diǎn)作為虛擬殘差點(diǎn)，采用二分圖最優(yōu)權(quán)匹配方法，對(duì)正負(fù)殘差點(diǎn)以及接地點(diǎn)進(jìn)行最優(yōu)匹配連接，極大地降低了枝切線的長(zhǎng)度；最后利用路徑積分方法得到正確的解纏相位。通過(guò)在洞庭湖大橋?qū)嶋H星載SAR數(shù)據(jù)上進(jìn)行實(shí)驗(yàn)驗(yàn)證，并與MCF算法解纏結(jié)果進(jìn)行對(duì)比，證明該算法有效改善了解纏結(jié)果，減少了解纏誤差，適合于橋梁線狀區(qū)域的相位解纏。

1 稀疏網(wǎng)絡(luò)MCF相位解纏

基于永久散射體的干涉(PSI)技術(shù)，是目前時(shí)間序列InSAR技術(shù)最為常用的處理方法，其通過(guò)對(duì)干涉圖中的高相干性點(diǎn)(PS)進(jìn)行分析，經(jīng)過(guò)干涉處理、相位解纏和形變解算操作，可以獲取這些PS點(diǎn)的形變結(jié)果，從而得到整個(gè)區(qū)域的形變信息。

由于PSI技術(shù)只對(duì)PS點(diǎn)進(jìn)行處理，并且PS點(diǎn)在干涉圖中是稀疏分布，而MCF相位解纏方法是基于網(wǎng)絡(luò)的方法，必須先對(duì)PS點(diǎn)構(gòu)建網(wǎng)絡(luò)，一般使用Delaunay三角化方法來(lái)構(gòu)建稀疏網(wǎng)絡(luò)。因此基于Delaunay稀疏三角網(wǎng)絡(luò)的MCF算法是PSI技術(shù)中使用最多的相位解纏方法，該算法成熟可靠，計(jì)算效率高，在各類商用軟件中也廣泛使用，下面簡(jiǎn)單介紹該算法原理與步驟。

1.1 基于Delaunay稀疏三角網(wǎng)絡(luò)的MCF算法

定義纏繞相位梯度φij，解纏相位梯度φij，需要補(bǔ)償?shù)南辔徽麛?shù)模糊數(shù)kij,即有

φij=φij+2πkij.

(1)

相位解纏過(guò)程也就是求kij，在MCF相位解纏算法中，kij表示每條邊上的網(wǎng)絡(luò)流大小，如圖1所示。

圖1 網(wǎng)絡(luò)流算法模型Fig.1 Network flow algorithm model

MCF算法通過(guò)把正負(fù)殘差點(diǎn)視作供給點(diǎn)和需求點(diǎn)，將殘差點(diǎn)的連接問(wèn)題，轉(zhuǎn)換成網(wǎng)絡(luò)流問(wèn)題，在構(gòu)建的Delaunay稀疏三角網(wǎng)絡(luò)中，采用L1范數(shù)優(yōu)化模型，使得總的補(bǔ)償數(shù)最小。

(2)

若采用有向圖模型，可將公式(2)轉(zhuǎn)換成線性約束形式，本質(zhì)上來(lái)說(shuō)，該約束等效于要求在一定加權(quán)條件下使得總的枝切線長(zhǎng)度最小，而MCF方法就利用最小費(fèi)用流算法求得最優(yōu)加權(quán)枝切線。

該算法的具體步驟: 首先對(duì)PS點(diǎn)構(gòu)建Delaunay三角網(wǎng)絡(luò)；然后對(duì)每個(gè)三角形，根據(jù)閉合路徑相位梯度積分是否為0，識(shí)別正負(fù)殘差點(diǎn)；再采用最小費(fèi)用流算法，得到正負(fù)殘差點(diǎn)對(duì)的最優(yōu)連接方式，確定枝切線位置；最后采用flood-fill方法，選擇合適的參考點(diǎn)，對(duì)纏繞相位梯度沿路徑積分，當(dāng)積分路徑經(jīng)過(guò)枝切線時(shí)，要根據(jù)枝切線上流的大小，對(duì)積分相位進(jìn)行補(bǔ)償2π的整數(shù)倍操作，最終獲得解纏相位圖。

1.2 實(shí)驗(yàn)結(jié)果

本文實(shí)驗(yàn)數(shù)據(jù)采用2016年上半年共17幅洞庭湖大橋區(qū)域的TerraSAR聚束模式影像，從17幅影像的兩兩干涉對(duì)中挑選出56幅相干性比較好的干涉對(duì)，根據(jù)PSI處理方法進(jìn)行形變分析，采用幅度離差法和譜相關(guān)系數(shù)法選出的橋上PS點(diǎn)，并對(duì)PS點(diǎn)相位進(jìn)行干涉處理，得到橋面PS點(diǎn)的差分干涉相位圖，取某幅干涉圖如圖2所示(顏色變化表示干涉相位值變化)。

雖然在PSI處理流程中，可以先對(duì)差分干涉相位進(jìn)行形變解算，然后再進(jìn)行相位解纏，但由于斜拉橋受到較大的熱形變和負(fù)載影響，帶來(lái)極大的橋面PS點(diǎn)形變相位變化，此時(shí)，若直接采用形變解算算法會(huì)遇到較大困難，甚至無(wú)法求解形變結(jié)果。因此，為正確進(jìn)行形變分析，必須先對(duì)橋面相位進(jìn)行相位解纏，然后對(duì)解纏相位進(jìn)行形變解算。根據(jù)1.1節(jié)中的MCF相位解纏算法步驟，對(duì)洞庭湖大橋區(qū)域差分干涉相位進(jìn)行相位解纏，得到如圖3所示的結(jié)果，其中圖3(d)為了分析解纏相位誤差，沒(méi)有給出解纏相位干涉圖，而是給出解纏相位在空間分布的情況，可以更為直觀地發(fā)現(xiàn)誤差問(wèn)題。

圖2 洞庭湖大橋某差分干涉相位圖Fig.2 The difference interferometry phase image of the Dongtinghu Bridge

從圖3(a)可看到，利用Delaunay三角化方法構(gòu)建的三角網(wǎng)絡(luò)存在一些很長(zhǎng)的邊，而且這些邊的位置都在水面上，顯然這些邊是不合理的，很可能會(huì)出現(xiàn)相位梯度變化太大的問(wèn)題，導(dǎo)致相位解纏出錯(cuò)；圖3(b)為某干涉對(duì)相位利用閉合相位梯度積分是否為0，識(shí)別出的殘差點(diǎn)，正負(fù)殘差點(diǎn)分別用“+”和“-”表示，可發(fā)現(xiàn)三角網(wǎng)絡(luò)中較長(zhǎng)的邊會(huì)帶來(lái)許多殘差點(diǎn)；圖3(c)為利用MCF得到的非零網(wǎng)絡(luò)流(枝切線)結(jié)果，用紅色標(biāo)注；圖3(d)為利用路徑積分方法，選擇橋面第一個(gè)點(diǎn)為參考點(diǎn)，對(duì)纏繞相位梯度積分的解纏相位。由于從干涉圖上難以直觀地發(fā)現(xiàn)解纏誤差，因此給出的是橋面PS點(diǎn)解纏相位在空間上分布結(jié)果，可看到解纏相位中紅圈所示區(qū)域出現(xiàn)明顯的跳變不連續(xù)情況，反映的正是MCF相位解纏過(guò)程中存在解纏誤差。在56對(duì)解纏相位結(jié)果中，僅有4對(duì)相位變化較小的干涉圖解纏結(jié)果是正確的，其余解纏相位均存在大量跳變誤差。這樣，在后續(xù)對(duì)解纏相位進(jìn)行形變解算時(shí)，這些相位跳變點(diǎn)給形變解算過(guò)程帶來(lái)較大困難，導(dǎo)致出現(xiàn)錯(cuò)誤的形變結(jié)果。

對(duì)正負(fù)殘差點(diǎn)的連接方式進(jìn)行分析，可發(fā)現(xiàn)一旦有連接正負(fù)殘差點(diǎn)的枝切線穿過(guò)橋面，就會(huì)導(dǎo)致解纏相位出現(xiàn)誤差，枝切線和解纏誤差對(duì)應(yīng)關(guān)系如圖4所示。出現(xiàn)這種情況的主要原因在于，斜拉橋橋塔區(qū)域?qū)蛎娈a(chǎn)生遮擋，在橋塔附近，橋面的點(diǎn)與橋塔點(diǎn)混在一起，而橋塔點(diǎn)本身存在較大的高程相位變化，在該區(qū)域極易出現(xiàn)孤立殘差點(diǎn)，導(dǎo)致正負(fù)殘差點(diǎn)出現(xiàn)不合理的連接，帶來(lái)解纏誤差；另外，橋梁區(qū)域?yàn)榫€狀區(qū)域，一旦某點(diǎn)出現(xiàn)解纏誤差，會(huì)導(dǎo)致誤差傳播到整個(gè)橋上，因此只要橋面存在不合理的枝切線，就會(huì)帶來(lái)較大的解纏誤差，并導(dǎo)致橋面的形變結(jié)果無(wú)法正確獲取，這也是MCF算法在橋梁區(qū)域解纏結(jié)果不如城區(qū)的原因。

從洞庭湖大橋的處理結(jié)果來(lái)看，MCF算法在進(jìn)行正負(fù)殘差點(diǎn)連接時(shí)，僅僅以加權(quán)連接線長(zhǎng)短為目的，沒(méi)有考慮區(qū)域具體特征，在斜拉橋橋塔附近的殘差點(diǎn)容易產(chǎn)生錯(cuò)誤的連接。同理，目前其他的算法[8-18]，如枝切法、最小生成樹算法等這類路徑跟蹤算法，還有最優(yōu)化方法，在處理大型橋梁時(shí)都會(huì)遇到較大的困難。本文提出一種基于二分圖的分塊殘差點(diǎn)最優(yōu)匹配算法，可在一定程度上約束殘差點(diǎn)的連接，獲取正確的解纏結(jié)果。

2 基于二分圖的分塊殘差點(diǎn)最優(yōu)匹配算法

2.1 算法流程

基于二分圖的分塊殘差點(diǎn)最優(yōu)匹配算法主要有3個(gè)步驟，構(gòu)建約束三角網(wǎng)絡(luò)、殘差點(diǎn)分塊和殘差點(diǎn)二分圖最優(yōu)匹配，下面對(duì)這3點(diǎn)分別進(jìn)行詳細(xì)介紹。

1)構(gòu)建約束三角網(wǎng)絡(luò)

從橋面PS點(diǎn)構(gòu)成的Delaunay三角網(wǎng)絡(luò)可看到，由于Delaunay構(gòu)網(wǎng)只考慮拓?fù)浣Y(jié)構(gòu)，保證構(gòu)建的網(wǎng)絡(luò)結(jié)構(gòu)滿足凸包性質(zhì)，在橋面上一些凸出部分周圍構(gòu)建的三角形會(huì)出現(xiàn)不合理的長(zhǎng)邊，而這類區(qū)域極易產(chǎn)生殘差點(diǎn)。因此需要對(duì)構(gòu)建的Delaunay三角網(wǎng)絡(luò)進(jìn)行約束，在保證網(wǎng)絡(luò)連通性的同時(shí)，設(shè)置一個(gè)合理的邊長(zhǎng)閾值，構(gòu)建Delaunay約束三角網(wǎng)絡(luò)。

2)殘差點(diǎn)分塊

構(gòu)建約束三角網(wǎng)絡(luò)，會(huì)在原有的Delaunay三角網(wǎng)絡(luò)上刪除很多三角邊，在殘差點(diǎn)數(shù)目會(huì)降低同時(shí)，也給正負(fù)殘差點(diǎn)的匹配帶來(lái)一些問(wèn)題，由于橋面的殘差點(diǎn)分布變得更加稀疏，很可能會(huì)導(dǎo)致許多過(guò)長(zhǎng)的連接線出現(xiàn)。使用分塊方法可以解決該問(wèn)題，傳統(tǒng)上分塊操作，主要是按影像數(shù)據(jù)塊大小進(jìn)行處理的，但由于橋梁所在區(qū)域存在一些凸出部分，會(huì)給數(shù)據(jù)分塊時(shí)帶來(lái)一些問(wèn)題。本文采用基于殘差點(diǎn)距離的分塊方法，對(duì)殘差點(diǎn)進(jìn)行聚類，具體操作如下：

①計(jì)算每個(gè)殘差點(diǎn)到其他殘差點(diǎn)的距離，包括歐式距離和連通性距離(即連接兩殘差點(diǎn)需要經(jīng)過(guò)的最小三角形邊數(shù))，構(gòu)建殘差點(diǎn)距離矩陣。

②根據(jù)歐氏距離閾值和連通性距離閾值，判斷兩殘差點(diǎn)是否分在同一塊，若在，則為1，否則為0，從而構(gòu)建分塊殘差點(diǎn)0～1矩陣。

③遍歷殘差點(diǎn)0～1矩陣的各連通分量，其中每個(gè)連通分量中的殘差點(diǎn)即分在同一塊，得到最終殘差點(diǎn)的分塊結(jié)果如圖5所示。

圖5 殘差點(diǎn)分塊示意圖Fig.5 Residual points partition

3)殘差點(diǎn)二分圖最優(yōu)匹配

殘差點(diǎn)的最優(yōu)連接屬于組合優(yōu)化問(wèn)題，選擇合適的正負(fù)殘差點(diǎn)的連接方式，能使得總得枝切線長(zhǎng)度最小，枝切線的總長(zhǎng)度一般定義如下：

(3)

公式(3)只考慮正負(fù)殘差點(diǎn)完全匹配的情況。實(shí)際殘差點(diǎn)分塊后，為保證每塊中正負(fù)殘差點(diǎn)數(shù)目保持平衡，需要適當(dāng)引入一部分接地點(diǎn)來(lái)平衡殘差點(diǎn)的極性。若將接地點(diǎn)視作虛擬的正負(fù)殘差點(diǎn)，與相應(yīng)的殘差點(diǎn)進(jìn)行匹配，枝切線長(zhǎng)度就仍能保持公式(3)中的形式，即N個(gè)正殘差點(diǎn)與N個(gè)負(fù)殘差點(diǎn)匹配，這樣正負(fù)殘差點(diǎn)的連接可轉(zhuǎn)換成二分圖最優(yōu)匹配模型，如圖6所示。

利用KM算法來(lái)獲取最優(yōu)匹配結(jié)果，假定正殘差點(diǎn)個(gè)數(shù)大于負(fù)殘差點(diǎn)，算法具體步驟如下：

①將正殘差點(diǎn)構(gòu)成點(diǎn)集{X}，負(fù)殘差點(diǎn)構(gòu)成點(diǎn)集{Y}；

②計(jì)算每個(gè)正殘差點(diǎn)到每個(gè)負(fù)殘差點(diǎn)的距離，由于正殘差點(diǎn)多，要加入一部分接地點(diǎn)作為虛擬的負(fù)殘差點(diǎn)，因此需要計(jì)算每個(gè)正殘差點(diǎn)到最近接地點(diǎn)的距離。取所有的這些距離的負(fù)值作為正負(fù)殘差點(diǎn)二分圖G(X,Y)的匹配權(quán)值矩陣；

③定義頂標(biāo)，正殘差點(diǎn)頂標(biāo){Cx}，負(fù)殘差點(diǎn)頂標(biāo){Cy}；

④初始化頂標(biāo)：

cx(i)=max{G(i,j)}，j=1,2…N，

cy(i)=0；

(4)

⑤使用匈牙利匹配算法求最優(yōu)匹配；

⑥對(duì)匹配結(jié)果中任意一條匹配邊，判斷其兩端頂點(diǎn)對(duì)應(yīng)的頂標(biāo)是否滿足

cx(i)+cy(j)==G(i,j),

(5)

若相等，則該匹配結(jié)果即為最優(yōu)權(quán)匹配，停止迭代，否則繼續(xù)下一步；

⑦對(duì)不滿足⑥中條件的匹配節(jié)點(diǎn)更新頂標(biāo):

cx(i)=cx(i)-min(cx(i)+cy(j)-G(i,j)),

cy(j)=cy(j)+min(cx(i)+cy(j)-G(i,j)).

(6)

⑧跳轉(zhuǎn)到步驟⑤，直到算法結(jié)束。

該算法時(shí)間復(fù)雜度為O(N4)，但由于采用分塊方法對(duì)殘差點(diǎn)進(jìn)行了細(xì)分，每塊中要匹配的殘差點(diǎn)數(shù)目比較小，因此算法計(jì)算效率很高。使用匹配算法得到最優(yōu)殘差點(diǎn)匹配連接后，就可以將相應(yīng)的殘差點(diǎn)連接起來(lái)，得到最優(yōu)枝切線，后續(xù)對(duì)纏繞相位梯度進(jìn)行路徑積分就能得到解纏相位。

圖6 殘差點(diǎn)連接的二分圖匹配模型Fig.6 The bipartite graph matching model for residual points

2.2 實(shí)驗(yàn)驗(yàn)證

1)約束三角網(wǎng)絡(luò)

對(duì)洞庭湖大橋PS點(diǎn)構(gòu)建約束三角網(wǎng)絡(luò)，采用的邊長(zhǎng)閾值為100，得到的三角網(wǎng)絡(luò)如圖7所示。

從圖7(b)可看到，較長(zhǎng)的邊都被排除了，剩下的三角形都是比較合理的，都集中在橋上，而不再有邊出現(xiàn)在水面上。另外，采用約束三角網(wǎng)絡(luò)去除不合理邊的同時(shí)，殘差點(diǎn)也極大地減少，平均而言，殘差點(diǎn)數(shù)目減少約1/5，殘差點(diǎn)數(shù)目隨三角網(wǎng)絡(luò)變化如表1所示。

2)殘差點(diǎn)分塊

在分塊過(guò)程中，為避免有正負(fù)殘差點(diǎn)的連線穿過(guò)橋面，需要設(shè)置較小的距離閾值，但該距離閾值同時(shí)也要足夠大，以保證塊內(nèi)的正負(fù)殘差點(diǎn)能夠正常匹配。取洞庭湖大橋某干涉圖，對(duì)其構(gòu)建約束三角網(wǎng)絡(luò)后，對(duì)每個(gè)三角網(wǎng)，利用閉合相位梯度是否為0，來(lái)識(shí)別出其中的正負(fù)殘差點(diǎn)，然后選擇好合適的閾值，根據(jù)2.1節(jié)中的分塊算法，對(duì)殘差點(diǎn)進(jìn)行分塊，得如圖8所示結(jié)果。

圖7 (a)原始三角網(wǎng)絡(luò)和(b)約束三角網(wǎng)絡(luò)Fig.7 Delaunay triangulation network (a) and the constrained triangulation network (b)

于涉圖序號(hào)殘差點(diǎn)數(shù)目Delaunay三角網(wǎng)絡(luò)約束三角網(wǎng)絡(luò)168482776037256464425534065843

從圖8可看到，采用分塊算法將殘差點(diǎn)分成了很多小的聚集區(qū)域，但這不會(huì)影響塊內(nèi)正負(fù)殘差點(diǎn)極性的平衡，因?yàn)闃蛄簠^(qū)域?yàn)榫€狀分布的特點(diǎn)，會(huì)存在很多邊界，有豐富的接地點(diǎn)，通過(guò)加入接地點(diǎn)可以平衡多余極性的殘差點(diǎn)。對(duì)殘差點(diǎn)進(jìn)行分塊，可以避免較長(zhǎng)的連接線，同時(shí)減小進(jìn)行匹配的正負(fù)殘差點(diǎn)數(shù)目，降低后續(xù)匹配優(yōu)化的難度。

圖8 殘差點(diǎn)分塊結(jié)果示例圖Fig.8 Residual points partition results

3)殘差點(diǎn)二分圖最優(yōu)匹配

根據(jù)前面的KM算法流程，對(duì)洞庭湖大橋的某干涉圖分塊殘差點(diǎn)進(jìn)行匹配，得到枝切線結(jié)果和MCF算法對(duì)比如圖9所示。

圖9 MCF算法得到的枝切線分布圖(a)和本文算法得到的枝切線分布圖(b)Fig.9 The branch cuts by MCF method (a) and by the proposed method (b)

從圖9可看到，本文提出的算法產(chǎn)生枝切線長(zhǎng)度比較小，明顯優(yōu)于MCF算法結(jié)果。通過(guò)采用殘差點(diǎn)分塊方式，正負(fù)殘差點(diǎn)在保持MCF算法中較短枝切線的同時(shí)，對(duì)長(zhǎng)枝切線處正負(fù)殘差點(diǎn)進(jìn)行重新連接，從而獲得整體比較短的枝切線結(jié)果。

完成殘差點(diǎn)的連接后，就可以選擇積分參考點(diǎn)，采用路徑積分方法對(duì)纏繞相位梯度進(jìn)行積分，得到解纏相位。按照該方法對(duì)56對(duì)PS點(diǎn)差分干涉圖進(jìn)行相位解纏，得到的解纏結(jié)果相對(duì)于MCF算法而言，解纏誤差得到了很大改善，有20對(duì)解纏結(jié)果是完全正確的，如圖10所示。可看到本文得到的解纏相位比MCF結(jié)果更加平滑，原來(lái)的跳變點(diǎn)處的相位也變得連續(xù)了，證明本算法能有效解決橋梁區(qū)域相位解纏問(wèn)題。

圖10 MCF算法解纏結(jié)果(a-c)和對(duì)應(yīng)的本文算法解纏結(jié)果(d-f)Fig.10 The unwrapped results by MCF method (a-c) and the proposed method (d-f)

除上述正確的干涉圖外，剩下的干涉圖只解決了部分誤差問(wèn)題，仍存在一些較小的誤差，主要由于橋面殘差點(diǎn)分布比較復(fù)雜，很難完全避免枝切線穿過(guò)橋面。對(duì)橋梁上的相位而言，必須保證完全沒(méi)有誤差，才能正常獲取形變信息。為獲得更多完全正確的解纏相位圖，后續(xù)還需要采用更為有效的約束條件，來(lái)限制殘差點(diǎn)的連接方式，從而解決橋梁形變測(cè)量的難題。

3 總結(jié)

本文針對(duì)現(xiàn)有相位解纏方法在大型橋梁區(qū)域存在較大解纏誤差的問(wèn)題，提出一種改進(jìn)的解纏算法，通過(guò)構(gòu)建約束三角網(wǎng)絡(luò)，降低殘差點(diǎn)數(shù)量，隨后又利用殘差點(diǎn)分塊方法，避免殘差點(diǎn)穿過(guò)橋面，保持較短的枝切線，同時(shí)降低匹配優(yōu)化難度，最后采用二分圖最優(yōu)權(quán)匹配方法，得到殘差點(diǎn)的最優(yōu)連接。該算法尤其適合在橋梁這類線狀區(qū)域進(jìn)行相位解纏，能有效降低解纏誤差，提高解纏相位精度。

由于橋面殘差點(diǎn)分布復(fù)雜，即便采用本文提出的約束方法，仍不能確保連接線不穿過(guò)橋面，需要后續(xù)對(duì)該問(wèn)題進(jìn)行更深入的分析，尋找更好的約束條件，來(lái)解決橋面PS點(diǎn)相位解纏誤差問(wèn)題，從而更好地獲取大型橋梁形變信息。