張紅敏，靳國(guó)旺，2，徐 青，秦志遠(yuǎn)，孫 偉

1.信息工程大學(xué)測(cè)繪學(xué)院，河南 鄭州450052；2.中國(guó)科學(xué)院 電子學(xué)研究所 微波成像技術(shù)國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室，北京100190

中國(guó)余數(shù)定理在雙基線(xiàn)InSAR相位解纏中的應(yīng)用

張紅敏1，靳國(guó)旺1，2，徐 青1，秦志遠(yuǎn)1，孫 偉1

1.信息工程大學(xué)測(cè)繪學(xué)院，河南 鄭州450052；2.中國(guó)科學(xué)院 電子學(xué)研究所 微波成像技術(shù)國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室，北京100190

介紹利用中國(guó)余數(shù)定理求解同余方程組的基本公式，并將中國(guó)余數(shù)定理引入到雙基線(xiàn)InSAR相位解纏中。構(gòu)造關(guān)于干涉相位模糊數(shù)的同余方程組，設(shè)計(jì)基于中國(guó)余數(shù)定理的雙基線(xiàn)InSAR相位解纏方法，擴(kuò)展相位解纏的非模糊區(qū)間，為解決干涉相位欠采樣區(qū)域的相位解纏難題奠定基礎(chǔ)。利用不同地區(qū)的多套DEM仿真的雙基線(xiàn)InSAR干涉圖進(jìn)行相位解纏試驗(yàn)，驗(yàn)證該方法在干涉相位欠采樣等區(qū)域的解纏能力。

中國(guó)余數(shù)定理；合成孔徑雷達(dá)干涉測(cè)量；多基線(xiàn)；雙基線(xiàn)；相位解纏

1 引 言

中國(guó)余數(shù)定理（CRT，Chinese remainder theory）可用于解決理想兩兩互素情況下線(xiàn)性同余方程組解的計(jì)算問(wèn)題［1］。該定理在信號(hào)處理、密碼系統(tǒng)等領(lǐng)域已得到實(shí)際應(yīng)用。文獻(xiàn)［2］將中國(guó)余數(shù)定理應(yīng)用于多基線(xiàn)相位干涉儀的相位差解模糊；文獻(xiàn)［3］研究了基于中國(guó)余數(shù)定理的同步圖像壓縮與加密方法；文獻(xiàn)［4—5］應(yīng)用中國(guó)余數(shù)定理解決SAR動(dòng)目標(biāo)成像時(shí)的相位模糊問(wèn)題；文獻(xiàn)［6］設(shè)計(jì)了基于中國(guó)余數(shù)定理作為陷門(mén)的公鑰密碼系統(tǒng)；文獻(xiàn)［7］分析了中國(guó)余數(shù)定理用于分布式星載SAR－ATI系統(tǒng)解速度模糊的方法。

多基線(xiàn)InSAR技術(shù)［8－9］具有高可靠性與高精度優(yōu)勢(shì)，已經(jīng)成為InSAR技術(shù)的重要發(fā)展方向之一。多基線(xiàn)InSAR獲取高精度DEM技術(shù)的核心問(wèn)題之一是相位解纏?？紤]到多基線(xiàn)InSAR相位解纏與同余問(wèn)題的相通性，為了解決多基線(xiàn)InSAR相位解纏難題，筆者從最簡(jiǎn)單的雙基線(xiàn)InSAR基本原理入手，將中國(guó)余數(shù)定理應(yīng)用于雙基線(xiàn)InSAR相位解纏中。構(gòu)造了關(guān)于干涉相位模糊數(shù)的同余方程組，設(shè)計(jì)了基于中國(guó)余數(shù)定理的雙基線(xiàn)InSAR相位解纏方案。利用多套不同地區(qū)DEM仿真的雙基線(xiàn)InSAR干涉圖進(jìn)行了相位解纏試驗(yàn)，得到了滿(mǎn)意的解纏結(jié)果，驗(yàn)證了所設(shè)計(jì)解纏方案的有效性。

2 中國(guó)余數(shù)定理

中國(guó)余數(shù)定理又稱(chēng)中國(guó)剩余定理、孫子定理，最早見(jiàn)諸南北朝時(shí)期的數(shù)學(xué)名著《孫子算經(jīng)》［11］，用于解決理想兩兩互素情況下線(xiàn)性同余方程組求解問(wèn)題。

《孫子算經(jīng)》中記載的“物不知其數(shù)”［1］問(wèn)題是中國(guó)余數(shù)定理的典型應(yīng)用：“今有物不知其數(shù)，三三數(shù)之剩二，五五數(shù)之剩三，七七數(shù)之剩二，問(wèn)物幾何？”用同余式表示“物不知其數(shù)”問(wèn)題，就轉(zhuǎn)化為如式（1）所示的線(xiàn)性同余方程組求解問(wèn)題。

式中，x為同余方程組的解，x≡a（mod m）表示x和a模m同余，整數(shù)m稱(chēng)為同余的模。

《孫子算經(jīng)》中給出了“物不知其數(shù)”問(wèn)題的解法，程大位進(jìn)一步在《算法統(tǒng)宗》中給出了計(jì)算口訣［11］：“三人同行七十稀，五樹(shù)梅花廿一枝，七子團(tuán)圓正月半，除百零五便得知?！睌?shù)學(xué)語(yǔ)言表述，即

式中，a1、a2、a3為余數(shù)，滿(mǎn)足式（1）的最小正整數(shù)解為23。

“物不知其數(shù)”問(wèn)題描述了同余問(wèn)題的一個(gè)典型特例，其解法口訣實(shí)際上是中國(guó)余數(shù)定理的一個(gè)特定用法。中國(guó)余數(shù)定理的表述有多種，依據(jù)本文需要，表述如下［1，12］

若m1，m2，…，mn兩兩互素，a1，a2，…，an是任意整數(shù)，則同余方程組

在0≤x＜m＝m1m2…mn內(nèi)有唯一的解，并且解一定可以寫(xiě)為

式中

3 雙基線(xiàn)InSAR獲取DEM基本原理

如圖1所示，假定主天線(xiàn)相位中心S1和從天線(xiàn)相位中心S2、S3位于同一直線(xiàn)上，三天線(xiàn)相位中心分別形成兩副穩(wěn)定基線(xiàn)，構(gòu)成單發(fā)三收式雙基線(xiàn)InSAR系統(tǒng)。記S1的高程為H，S1和S2形成的基線(xiàn)為B，基線(xiàn)水平角為α。沿視線(xiàn)向的基線(xiàn)分量，即平行基線(xiàn)為B‖，垂直于視線(xiàn)向的基線(xiàn)分量，即垂直基線(xiàn)為B⊥。對(duì)目標(biāo)點(diǎn)P的側(cè)視角為θ，干涉圖中相鄰像元對(duì)應(yīng)的地面點(diǎn)為P′，P′對(duì)P的相對(duì)高程為dZ，P到天線(xiàn)相位中心S1和S2的斜距分別為R和R－ΔR，ΔR為斜距差。

圖1 雙基線(xiàn)InSAR幾何模型Fig.1 Geometry of dual－baseline InSAR

由圖1所示幾何關(guān)系可知，地面點(diǎn)P的高程Z為

式中，側(cè)視角θ與基線(xiàn)水平角α及角β的關(guān)系為

斜距R的計(jì)算公式為

式中，R0為SAR系統(tǒng)近距延遲，Mslant為SAR圖像斜距向像元大小，y為地面點(diǎn)P對(duì)應(yīng)像點(diǎn)的斜距向坐標(biāo)。

在△S1S2P中，根據(jù)余弦定理有

則

由式（6）、式（7）和式（10），得

記λ為雷達(dá)波波長(zhǎng)。對(duì)于單發(fā)多收模式，斜距差ΔR與真實(shí)干涉相位Φa的關(guān)系為

由于干涉圖中記錄的是干涉相位主值φ，其值域?yàn)椋?，2π），在InSAR處理中需要經(jīng)過(guò)相位解纏和干涉相位偏置參數(shù)定標(biāo)才能由φ得到Φa。

式（11）對(duì)ΔR求偏導(dǎo)，可得

ΔR相對(duì)于R是一個(gè)可以忽略的小量，上式可近似為

聯(lián)立式（12）和式（14），得

整理上式，可得

式中，

Z＊表示引起一個(gè)2π干涉相位變化所對(duì)應(yīng)的高程變化，稱(chēng)為高程模糊數(shù)。由式（17）不難看出，高程模糊數(shù)與垂直基線(xiàn)成反比，垂直基線(xiàn)越長(zhǎng)，高程模糊數(shù)越小，高程反演精度越高；反之，垂直基線(xiàn)越短，高程模糊數(shù)越大，高程反演精度越低。對(duì)于如圖1所示的幾何構(gòu)型，垂直基線(xiàn)之比等于基線(xiàn)長(zhǎng)度之比。為簡(jiǎn)化算法，本文以基線(xiàn)長(zhǎng)度之比代替垂直基線(xiàn)之比。

4 雙基線(xiàn)InSAR相位解纏原理

為了將中國(guó)余數(shù)定理應(yīng)用于雙基線(xiàn)InSAR相位解纏中，需要建立同余方程組形式的雙基線(xiàn)InSAR相位解纏模型。

假設(shè)在基線(xiàn)分別為B1和B2的情況下獲取的纏繞干涉相位分別為φ1和φ2，對(duì)應(yīng)的真實(shí)干涉相位分別為Φa1和Φa2，解纏后的干涉相位分別為Φ1和Φ2。根據(jù)式（15）可知，同一相對(duì)高程dZ與各干涉相位微分dφi（i＝1，2）之間的關(guān)系為

式中，ki為dφi的模糊數(shù)。

任意有限位數(shù)的實(shí)數(shù)均可以乘以10q（其中，q為該實(shí)數(shù)的小數(shù)位數(shù)），保留全部有效位數(shù)，轉(zhuǎn)換為無(wú)精度損失的整數(shù)。因此，可假定基線(xiàn)長(zhǎng)度B1、B2均為正整數(shù)。兩基線(xiàn)長(zhǎng)度的最小公倍數(shù)為B0＝［B1，B2］，記模mi（i＝1，2）為mi＝B0／Bi。則由式（18）可得

其中，ai取整數(shù)，且－mi＜ai＜mi。需要指出的是，余數(shù)a1、a2是實(shí)際干涉圖中纏繞干涉相位微分的函數(shù)，是一組實(shí)數(shù)。并且由于實(shí)際中相位噪聲的存在，式（20）可能并不成立。然而，在一定的噪聲水平下，通過(guò)余數(shù)取整，仍可將雙基線(xiàn)InSAR相位解纏問(wèn)題視為整數(shù)域內(nèi)的同余問(wèn)題，表示成線(xiàn)性同余方程組為

因?yàn)閙i＝B0／Bi，以此方式構(gòu)造的模m1、m2互素，從而滿(mǎn)足了中國(guó)余數(shù)定理的要求。根據(jù)中國(guó)余數(shù)定理，同余方程組（21）在0≤x＜m內(nèi)有唯一解，并且解可以寫(xiě)成

式中

式中，M－1i是Mi的乘法逆元，即

求解M－1i是應(yīng)用中國(guó)余數(shù)定理的關(guān)鍵之一，可利用擴(kuò)展歐幾里得算法［12］遞歸求解M－1i。

由中國(guó)余數(shù)定理求得式（21）在［0，m）的最小正整數(shù)解xmin后，可由下式計(jì)算兩幅干涉圖中對(duì)應(yīng)像元處的模糊數(shù)ki（i＝1，2），為

式中，余數(shù)ai需進(jìn)行取整，也可寫(xiě)作Int（ai）。那么，兩幅干涉圖中該像元處解纏后的干涉相位微分dΦi分別為

在兩幅干涉圖中，均以像素（x0，y0）為相位解纏種子，分別對(duì)dΦi沿某一路徑積分，就可以同時(shí)解纏兩幅干涉圖，即

式中，（x，y）為干涉圖像點(diǎn)的方位向和斜距向坐標(biāo)。

由上述分析可知，將雙基線(xiàn)InSAR相位解纏問(wèn)題轉(zhuǎn)化為線(xiàn)性同余方程組求解問(wèn)題，突破了傳統(tǒng)相位解纏方法對(duì)相鄰干涉相位差不超過(guò)半周期的限制，擴(kuò)展了正確解纏的非模糊區(qū)間，即將非模糊區(qū)間從［－π，π）擴(kuò)展到［－mπ，mπ），從而可以解決干涉相位欠采樣和頻譜混疊區(qū)域的相位解纏難題。

5 雙基線(xiàn)InSAR相位解纏流程

基于中國(guó)余數(shù)定理的雙基線(xiàn)InSAR相位解纏包括依據(jù)基線(xiàn)情況構(gòu)造互素模、建立關(guān)于纏繞相位微分模糊數(shù)的同余方程組、利用中國(guó)余數(shù)定理求解同余方程組、計(jì)算模糊數(shù)、求得解纏后的相位微分值、相位積分等關(guān)鍵步驟。其基本流程如圖2所示。

圖2 基于CRT的雙基線(xiàn)InSAR相位解纏流程Fig.2 Phase unwrapping flowchart based on CRT for dual－baseline InSAR

（1）根據(jù)配準(zhǔn)結(jié)果，獲取分別對(duì)應(yīng)于基線(xiàn)B1、B2的雙基線(xiàn)InSAR干涉圖（即纏繞干涉相位）φ1、φ2。構(gòu)造互素模m1、m2

式中，B0＝［B1，B2］。當(dāng)B1、B2為有限位數(shù)的小數(shù)時(shí)，首先將其同時(shí)擴(kuò)大整數(shù)倍，轉(zhuǎn)化為正整數(shù)。

（2）計(jì)算相鄰像素（x1，y1）和（x2，y2）之間的纏繞干涉相位微分dφ1、dφ2

（3）構(gòu)造a1、a2。由dφ1、dφ2構(gòu)造同余方程組的余數(shù)a1、a2，并進(jìn)行取整，將同余方程組納入整數(shù)環(huán)

（4）建立線(xiàn)性同余方程組。建立關(guān)于纏繞相位微分的模糊數(shù)k1、k2的同余方程組

（5）計(jì)算p。根據(jù)中國(guó)余數(shù)定理，在0≤p＜m1m2內(nèi)，解同余方程組，得到p。其中，乘法逆元利用擴(kuò)展歐幾里得算法迭代求解。

（6）計(jì)算k1、k2。模糊數(shù)k1、k2為

（7）計(jì)算dΦ1、dΦ2。解纏后的干涉相位微分dΦ1、dΦ2為

（8）進(jìn)行相位積分，計(jì)算干涉圖任意像元的解纏干涉相位Φ1、Φ2。分別以給定的相位解纏種子φ1（x0，y0）、φ2（x0，y0），對(duì)dΦ1、dΦ2沿某一路徑積分，得到兩幅干涉圖的解纏結(jié)果

6 試 驗(yàn)

為了驗(yàn)證基于中國(guó)余數(shù)定理的雙基線(xiàn)InSAR相位解纏方法的有效性，本文利用不同地區(qū)的DEM仿真了多套雙基線(xiàn)InSAR干涉圖，進(jìn)行了相位解纏試驗(yàn)。采用的DEM如圖3所示，其中（a）為由SIR－C／X－SAR干涉數(shù)據(jù)獲取的天山地區(qū)10m格網(wǎng)間距的DEM，（b）為SRTM獲取的海南地區(qū)90m格網(wǎng)間距的DEM。為了充分驗(yàn)證該方法對(duì)含任意有限位小數(shù)基線(xiàn)情況的適用性，試驗(yàn)中，選擇的基線(xiàn)長(zhǎng)度既有整數(shù)也有一般實(shí)數(shù)。所用仿真參數(shù)如表1所示，仿真的去平地效應(yīng)后的雙基線(xiàn)干涉圖如圖4所示。

圖3 DEMFig.3 DEM

表1 干涉圖仿真參數(shù)Tab.1 Parameters of interferogram simulation

圖4 仿真的雙基線(xiàn)干涉圖Fig.4 Simulated dual－baseline InSAR interferograms

由圖3和圖4不難看出，天山及海南地區(qū)的仿真干涉圖中均存在欠采樣和頻譜混疊現(xiàn)象。為了突出單基線(xiàn)InSAR技術(shù)在欠采樣和頻譜混疊區(qū)域的相位解纏問(wèn)題，在進(jìn)行相位解纏時(shí)，兩類(lèi)解纏方法均采用沿距離向直接相位積分策略。其中，單基線(xiàn)解纏結(jié)果如圖5所示；基于中國(guó)余數(shù)定理的雙基線(xiàn)相位解纏結(jié)果如圖6所示。

圖5 單基線(xiàn)相位解纏結(jié)果Fig.5 Unwrapped interferograms with single baseline

圖6 基于中國(guó)余數(shù)定理的雙基線(xiàn)相位解纏結(jié)果Fig.6 Unwrapped interferograms with dualbaseline based on CRT

由圖5不難看出，由于仿真干涉圖中存在頻譜混疊現(xiàn)象，采用單基線(xiàn)相位解纏方法，難以得到良好的相位解纏結(jié)果。并且，由于積分路徑?jīng)]有繞開(kāi)殘差點(diǎn)，誤差傳遞導(dǎo)致了整行解纏錯(cuò)誤的點(diǎn)，從而更加鮮明地反映了單基線(xiàn)InSAR技術(shù)在頻譜混疊區(qū)域的相位解纏困境。采用基于中國(guó)余數(shù)定理的雙基線(xiàn)InSAR相位解纏方法時(shí)仍采用相同的積分路徑，由圖6易知，天山和海南地區(qū)的仿真雙基線(xiàn)干涉圖均得到了良好的解纏效果，解纏后的干涉圖所顯示的地形起伏狀況分別與圖3 DEM1和DEM2所顯示的地形相吻合，較長(zhǎng)基線(xiàn)對(duì)應(yīng)的解纏后干涉圖包含更多的細(xì)部信息。為直觀比較兩種方法的解纏差異，分別計(jì)算了兩種方法對(duì)天山地區(qū)仿真干涉圖的解纏結(jié)果（圖5（a）（b）和圖6（a）（b））與理論結(jié)果的較差圖，分別如圖7和圖8所示。

對(duì)比圖7和圖8可知，基于中國(guó)余數(shù)定理的相位解纏方法通過(guò)將雙基線(xiàn)相位解纏問(wèn)題轉(zhuǎn)化為線(xiàn)性同余方程組的求解問(wèn)題，擴(kuò)展了正確解纏的非模糊區(qū)間，從而解決了頻譜混疊區(qū)域的相位解纏難題。以海南地區(qū)仿真干涉圖的相位解纏為例，采用單基線(xiàn)InSAR相位解纏方法時(shí)，正確解纏的非模糊區(qū)間為［－π，π）；而采用本文所設(shè)計(jì)的雙基線(xiàn)InSAR相位解纏方法時(shí)，正確解纏的非模糊區(qū)間擴(kuò)展為［－mπ，mπ），從而能夠正確解纏頻譜混疊區(qū)域的纏繞相位。其中，m＝35，計(jì)算過(guò)程如下：

圖7 單基線(xiàn)相位解纏結(jié)果與理論解纏結(jié)果較差圖Fig.7 Interferogram difference between unwrapped interferograms with single－baseline and unwrapped interferograms in theory

圖8 雙基線(xiàn)相位解纏結(jié)果與理論解纏結(jié)果較差圖Fig.8 Interferogram difference between unwrapped interferograms with dual－baseline and unwrapped interferograms in theory

記基線(xiàn)B1和B2擴(kuò)大取整后對(duì)應(yīng)值分別為B′1和B′2

B′1和B′2的最小公倍數(shù)B′為

構(gòu)造模m1和m2分別為

m是模m1和m2的最小公倍數(shù)，為

由圖8易知，在無(wú)噪聲情況下，采用所設(shè)計(jì)方法，不同地區(qū)仿真干涉圖均得到了良好的解纏結(jié)果，驗(yàn)證了所設(shè)計(jì)方法的正確性和有效性。

為了分析噪聲對(duì)基于中國(guó)余數(shù)定理的雙基線(xiàn)InSAR相位解纏方法的影響，又采用仿真含噪聲干涉圖進(jìn)行了試驗(yàn)。在利用天山地區(qū)DEM仿真干涉圖時(shí)，加入了均值為0，方差為0.000 2rad2的隨機(jī)相位噪聲，干涉圖仿真結(jié)果如圖9所示。

圖9 仿真的含噪聲雙基線(xiàn)干涉圖Fig.9 Simulated dual－baseline interferograms with noise

對(duì)圖9所示含噪聲干涉圖，采用直接相位積分的單基線(xiàn)InSAR相位解纏方法進(jìn)行相位解纏的結(jié)果如圖10所示；采用基于中國(guó)余數(shù)的雙基線(xiàn)InSAR相位解纏方法進(jìn)行相位解纏，結(jié)果如圖11所示。

由圖10和圖11可以看出，兩種方法的解纏結(jié)果均受到了噪聲的影響，出現(xiàn)了解纏錯(cuò)誤的點(diǎn)并導(dǎo)致了沿積分路徑的誤差積累。其中，基于中國(guó)余數(shù)定理的解纏結(jié)果中部分解纏錯(cuò)誤區(qū)域如圖11中方框所示。為了避免噪聲影響引起的誤差積累，下一步研究中將選擇避開(kāi)殘差點(diǎn)的積分路徑進(jìn)行基于中國(guó)余數(shù)定理的雙基線(xiàn)InSAR相位解纏。關(guān)于噪聲量值對(duì)本文算法的定量影響及其解決方案將在后續(xù)工作中進(jìn)行深入研究。

圖10 單基線(xiàn)相位解纏結(jié)果Fig.10 Unwrapped interferograms with single baseline

圖11 基于中國(guó)余數(shù)定理的雙基線(xiàn)相位解纏結(jié)果Fig.11 Unwrapped interferograms with dualbaseline based on CRT

7 總 結(jié)

中國(guó)余數(shù)定理是解決同余問(wèn)題的有效方法，而多基線(xiàn)InSAR相位解纏的實(shí)質(zhì)可認(rèn)為是同余問(wèn)題。為了給多基線(xiàn)InSAR相位解纏提供良好技術(shù)方案，本文著眼于最簡(jiǎn)單的雙基線(xiàn)InSAR技術(shù)，將中國(guó)余數(shù)定理應(yīng)用于雙基線(xiàn)InSAR相位解纏中，設(shè)計(jì)了基于中國(guó)余數(shù)定理的雙基線(xiàn)InSAR相位解纏方法。采用SIR－C／X－SAR獲取的天山地區(qū)DEM和SRTM海南地區(qū)DEM仿真的雙基線(xiàn)InSAR干涉圖，進(jìn)行了相位解纏試驗(yàn)，得到了理想的解纏結(jié)果。通過(guò)進(jìn)一步對(duì)仿真含噪聲干涉圖的解纏試驗(yàn)，驗(yàn)證了一定噪聲水平下，該方案的正確性和有效性。與單基線(xiàn)相位解纏結(jié)果對(duì)比，驗(yàn)證了多基線(xiàn)InSAR技術(shù)高可解纏性的優(yōu)勢(shì)。

論文成果為進(jìn)一步研究三基線(xiàn)以上的InSAR相位解纏方案提供了參考。在后續(xù)研究中，將把該方法擴(kuò)展至三基線(xiàn)以上的相位解纏中；研究基線(xiàn)組合方案對(duì)相位解纏的影響；研究噪聲對(duì)相位解纏的影響，及采用質(zhì)量圖區(qū)域生長(zhǎng)、選擇避開(kāi)殘差點(diǎn)的積分路徑以及利用優(yōu)化思想改化同余方程組等方法，減弱噪聲對(duì)解纏結(jié)果的影響。

［1］ PEI Dingyi，ZHU Yuefei.Algorithmic Number Theory［M］.Beijing：Science Press，2002.（裴定一，祝躍飛.算法數(shù)論［M］.北京：科學(xué)出版社，2002.）

［2］ ZHOU Yaqiang，CHEN Zhu，HUANG Pukan，et al.Algorithm of Solving Multi－baseline Interferometer Phase Difference Ambiguity in Noisy Circumstance［J］.Journal of Electronics and Information Technology，2005，27（2）：259－261.（周亞強(qiáng)，陳翥，皇甫堪，孫仲康.噪擾條件下多基線(xiàn)相位干涉儀解模糊算法［J］.電子與信息學(xué)報(bào)，2005，27（2）：259－261.）

［3］ JAGANNATHAN V，MAHADEVAN A，HARIHARAN R，et al.Number Theory Based Image Compression Encryption and Application to Image Multiplexing［C］∥Proceedings of International Conference on Signal Processing，Communications and Networking.Chennai：［s.n.］，2007：59－64.

［4］ XIA Xianggen，WANG Genyuan.Phase Unwrapping and a Robust Chinese Remainder Theorem［J］.IEEE Signal Processing Letters，2007，14（4）：247－250.

［5］ LI Xiaowei，XIA Xianggen.A Fast Robust Chinese Remainder Theorem Based Phase Unwrapping Algorithm［J］.IEEE Signal Processing Letters，2008，15：665－668.

［6］ YASUYUKI M，KIYOKO K，MASAO K.A New Class of Cryptosystems Based on Chinese Remainder Theorem［C］∥Proceedings of International Symposium on Information Theorem and Its Applications.Auckland：［s.n.］，2008：1－6.

［7］ QI Weikong，DANG Yawen，YU Weidong.Deblurring Velocity Ambiguity of Distributed Space－borne SAR Based on Chinese Remainder Theorem［J］.Journal of Electronics and Information Technology，2009，31（10）：2493－2497.（齊維孔，黨雅文，禹衛(wèi)東.基于中國(guó)剩余定理解分布式星載SAR－ATI測(cè)速模糊［J］.電子與信息學(xué)報(bào)，2009，31（10）：2493－2497.）

［8］ JIN Guowang.Research on Key Processing Techniques for Accurate DEM Deriving from InSAR［D］.Zhengzhou：Institute of Surveying and Mapping，Information Engineering University，2007.（靳國(guó)旺.InSAR獲取高精度DEM關(guān)鍵處理技術(shù)研究［D］.鄭州：信息工程大學(xué)測(cè)繪學(xué)院，2007.）

［9］ ZHANG Hongmin.Research on Data Simulation and Phase Unwrapping for Multi－baseline InSAR［D］.Zhengzhou：Institute of Surveying and Mapping，Information Engineering University，2010.（張紅敏.多基線(xiàn)InSAR數(shù)據(jù)仿真及其相位解纏技術(shù)研究［D］.鄭州：信息工程大學(xué)測(cè)繪學(xué)院，2010.）

［10］ JIN Guowang，XU Qing，ZHANG Yan，et al.The Zero Intermediate Frequency Vector Filtering for Interferograms［J］.Acta Geodaetica et Cartographica Sinica，2006，35（1）：24－29.（靳國(guó)旺，徐青，張燕，等.InSAR干涉圖的零中頻矢量濾波算法［J］.測(cè)繪學(xué)報(bào)，2006，35（1）：24－29.）

［11］ SHEN Hua，LIU Heguo.Application of Chinese Remainder Theorem on Some Problems［J］.High School Mathematics，2003，12：43－44.（沈華，劉合國(guó).中國(guó)剩余定理對(duì)幾道賽題的應(yīng)用［J］.中學(xué)數(shù)學(xué)，2003，12：43－44.）

［12］ KENNETH H R.Elementary Number Theory and Its Application［M］.XIA Honggang，trans.Beijing：China Machine Press，2009：104.（KENNETH H R.初等數(shù)論及其應(yīng)用［M］.夏鴻剛，譯.北京：機(jī)械工業(yè)出版社，2009：104.）

Application of Chinese Remainder Theorem in Phase Unwrapping for Dualbaseline InSAR

ZHANG Hongmin1，JIN Guowang1，2，XU Qing1，QIN Zhiyuan1，SUN Wei1
1.Institute of Surveying and Mapping，Information Engineering University，Zhengzhou 450052，China；2.National Key Laboratory of Science and Technology on Microwave Imaging，Institute of Electronics，Chinese Academy of Sciences，Beijing 100190，China

The basic formulation of solving congruence equations with Chinese Remainder Theorem（CRT）is introduced，and CRT is applied to phase unwrapping for dual－baseline InSAR.The congruence equations refer to interferometric phase fuzzy numbers are designed，and the phase unwrapping method for dual－baseline InSAR based on Chinese Remainder Theorem is proposed.In this algorithm，the un－fuzzy phase interval of phase unwrapping is extended，which helps to solve the difficult problem of phase unwrapping in sub－sampling area.The phase unwrapping experiments are done with simulated interferograms from different areas DEMs.The satisfying phase unwrapping results are gotten，as validated the phase unwrapping ability for the proposed algorithm in sub－sampling area.

Chinese remainder theorem；interferometric synthetic aperture radar；multi－baseline；dual－baseline；phase－inwrapping

ZHANG Hongmin（1984—），female，PhD candidate，major in InSAR，photogrammetry and remote sensing.

1001－1595（2011）06－0770－08

TP701

A

國(guó)家自然科學(xué)基金（40771142；40871213；41071296）；中國(guó)博士后科學(xué)基金（200801111）；國(guó)家西部1∶50 000地形圖空白區(qū)測(cè)圖工程項(xiàng)目；測(cè)繪學(xué)院學(xué)位論文創(chuàng)新與創(chuàng)優(yōu)基金

雷秀麗）

2010－09－25

2010－12－31

張紅敏（1984—），女，博士生，主要從事合成孔徑雷達(dá)干涉測(cè)量、攝影測(cè)量與遙感技術(shù)研究。

E－mail：zhmin1206＠163.com