[發(fā)明專利]多基線/多頻段干涉相位解纏頻域快速算法

有效

• 申請?zhí)枺?/td>: 201110312439.4
• 申請日: 2011-10-14
• 公開（公告）號: CN102621549A
• 公開（公告）日: 2012-08-01
• 發(fā)明（設(shè)計）人: 黃海風(fēng);張永勝;王青松;何志華;孫造宇;金光虎;董臻;余安喜;何峰;杜湘瑜;陳祺
• 申請（專利權(quán)）人: 中國人民解放軍國防科學(xué)技術(shù)大學(xué)
• 主分類號: G01S13/90
• 分類號: G01S13/90
• 代理公司: 國防科技大學(xué)專利服務(wù)中心 43202
• 代理人: 王文惠

技術(shù)領(lǐng)域

本發(fā)明屬于遙感和信號處理的交叉技術(shù)領(lǐng)域，特別涉及一種利用干涉合成孔徑雷達(dá)進(jìn)行多基線/多頻段條件下的干涉相位解纏的頻域快速方法。?

背景技術(shù)

一切將相位由主值或相位差值恢復(fù)為真實值的過程統(tǒng)稱為相位解纏，除了干涉合成孔徑雷達(dá)應(yīng)用外，相位解纏在合成孔徑聲納、自適應(yīng)光學(xué)、核磁共振、地震處理等方面都有重要應(yīng)用。利用多部干涉合成孔徑雷達(dá)形成多基線或多頻段干涉可提高相位解纏的性能。?

本發(fā)明以干涉合成孔徑雷達(dá)應(yīng)用為例。傳統(tǒng)的單基線或單頻段干涉合成孔徑雷達(dá)系統(tǒng)受干涉相位模糊和高程疊掩影響，在復(fù)雜地形區(qū)域相位解纏難度較大，極大地限制了單基線或單頻段干涉合成孔徑雷達(dá)系統(tǒng)的高精度全球測繪能力。多基線或多頻段干涉合成孔徑雷達(dá)系統(tǒng)的提出與實現(xiàn)則有效地提高了干涉合成孔徑雷達(dá)對復(fù)雜地形的測量精度和測量覆蓋能力。多基線或多頻段干涉合成孔徑雷達(dá)系統(tǒng)的最大優(yōu)點就是可以充分利用其長短基線或高低頻段獲取疏密不同的干涉相位條紋來提高相位解纏的性能，短基線或低頻段可以保證相位解纏的可靠性，長基線或高頻段可以提高測量精度。因此多基線/多頻段干涉合成孔徑雷達(dá)系統(tǒng)更具吸引力，是未來發(fā)展的趨勢。?

目前多基線/多頻段相位解纏方法主要有：中國余數(shù)定律法、投影法以及線性組合法、迭代法、時域最小二乘法、Kalman濾波法、最大似然法、最大后驗法、空-像域聯(lián)合子空間正交投影法和網(wǎng)絡(luò)流法等，其中多基線/多頻段時域最小二乘法的基本思想是使相位梯度估計值與多個相位梯度值加權(quán)和之差的平方和最小，等效于求解具有牛曼邊界的泊松方程，這種方法本質(zhì)是對誤差進(jìn)行平均，?特點是十分穩(wěn)健，但效率不高。?

發(fā)明內(nèi)容

本發(fā)明提出一種計算量小的多基線/多頻段相位解纏的頻域算法，該算法的基本思想是在頻域范圍內(nèi)使相位梯度估計值與多個相位梯度值加權(quán)和之差的平方和最小。?

本發(fā)明技術(shù)方案的思路是：首先，針對各基線/各頻段的干涉相位圖，利用干涉相位圖的干涉相位觀測值計算在水平方向及垂直方向上的相位梯度值，并根據(jù)時域無旋條件計算相位梯度圖邊界上的相位梯度值。然后，對兩個方向上的相位梯度值進(jìn)行傅立葉變換得到相位梯度值的頻域函數(shù)，在相位梯度估計值的頻域函數(shù)與各基線/各頻段相位梯度值的頻域函數(shù)加權(quán)和之差的平方和最小約束條件下計算出滿足無旋條件的傅立葉系數(shù)的近似值。其后，對計算得到的滿足無旋條件的傅立葉系數(shù)近似值，進(jìn)行傅立葉反變換，得到滿足相位無旋條件的相位梯度估計值。最后，對得到的滿足無旋條件的相位梯度估計值進(jìn)行沿著任意路徑的積分，從而得到相位解纏值。?

本發(fā)明技術(shù)方案是：?

假定通過多基線/多頻段觀測得到K幅干涉相位圖，其中第k幅干涉相位圖對應(yīng)第k條基線/第k個頻段，第k條基線/第k個頻段對應(yīng)的纏繞相位函數(shù)為：?k=1,2,糑,;m=0,1,2,糓,-1;n=0,1,2,糔,-1，M和N分別表示干涉相位圖的方位向和距離向點數(shù)。記第k幅干涉相位圖的垂直有效基線為b_k（或記第k幅干涉相位圖的波長為l_k），記第k幅干涉相位圖與第一幅干涉相位圖的垂直有效基線之比為a_k=b₁/b_k（或記第k幅干涉相位圖與第一幅干涉相位圖的波長之比為a_k=l_k/l₁）。利用上述的觀測信息，完成以下步驟：?

第一步：干涉相位圖的相位梯度值計算。?

本步驟對多基線/多頻段干涉相位進(jìn)行相位梯度值計算，同時根據(jù)相位的時域無旋條件計算相位梯度圖邊界上的相位梯度值。對每一幅干涉相位圖進(jìn)行下述計算：?

第（1）步，計算非邊界上的相位梯度值。?

對第k幅干涉相位圖，用下式計算在水平方向即x方向的相位梯度?在垂直方向即y方向的相位梯度?

上式中W{×}表示取相位主值運算。?

第（2）步：計算邊界上的相位梯度值。?

第二步，傅立葉系數(shù)近似值計算。?

計算第k幅干涉相位圖x方向相位梯度值的傅立葉變換系數(shù)?和y方向相位梯度值的傅立葉變換系數(shù)?計算公式為：?

其中，p=0,1,K,M-1;q=0,1，K,N-1。?

計算x方向相位梯度的傅立葉變換系數(shù)近似值?和y方向相位梯度的傅立葉變換系數(shù)近似值?計算公式為：?