[發(fā)明專(zhuān)利]一種InSAR時(shí)序相位的優(yōu)化方法及裝置

在審

• 申請(qǐng)?zhí)枺?/td>: 202210097356.6
• 申請(qǐng)日: 2022-01-26
• 公開(kāi)（公告）號(hào): CN114488150A
• 公開(kāi)（公告）日: 2022-05-13
• 發(fā)明（設(shè)計(jì)）人: 蔣彌;程曉;梁琦
• 申請(qǐng)（專(zhuān)利權(quán)）人: 中山大學(xué)
• 主分類(lèi)號(hào): G01S13/90
• 分類(lèi)號(hào): G01S13/90
• 代理公司: 廣州三環(huán)專(zhuān)利商標(biāo)代理有限公司 44202
• 代理人: 呂金金

1.一種InSAR時(shí)序相位的優(yōu)化方法，其特征在于，包括：

獲取時(shí)序SAR數(shù)據(jù)集，并對(duì)所述時(shí)序SAR數(shù)據(jù)集進(jìn)行L視處理，分別獲得L視強(qiáng)度數(shù)據(jù)集和干涉數(shù)據(jù)集；

以所述L視強(qiáng)度數(shù)據(jù)集為參考，獲取并分別對(duì)預(yù)設(shè)數(shù)字高程模型和預(yù)設(shè)土地覆蓋圖進(jìn)行配準(zhǔn)與地理編碼，分別得到SAR圖像坐標(biāo)系下的數(shù)字高程和SAR圖像坐標(biāo)系下的土地覆蓋圖；

根據(jù)所述SAR圖像坐標(biāo)系下的數(shù)字高程，對(duì)所述干涉數(shù)據(jù)集進(jìn)行差分操作，得到差分干涉數(shù)據(jù)集；

根據(jù)所述SAR圖像坐標(biāo)系下的土地覆蓋圖和所述差分干涉數(shù)據(jù)集，估算所述差分干涉數(shù)據(jù)集的協(xié)方差矩陣，并根據(jù)預(yù)設(shè)的時(shí)序相位極大似然估計(jì)公式，估計(jì)并獲得優(yōu)化后的時(shí)序相位。

2.如權(quán)利要求1所述的一種InSAR時(shí)序相位的優(yōu)化方法，其特征在于，所述獲取時(shí)序SAR數(shù)據(jù)集，并對(duì)所述時(shí)序SAR數(shù)據(jù)集進(jìn)行L視處理，分別獲得L視強(qiáng)度數(shù)據(jù)集和干涉數(shù)據(jù)集，具體為：

獲取具有n幅時(shí)序SAR單視影像的時(shí)序SAR數(shù)據(jù)集，并對(duì)所述n幅時(shí)序SAR單視影像進(jìn)行配準(zhǔn)，對(duì)配準(zhǔn)后的n幅時(shí)序SAR單視影像進(jìn)行L視處理，分別獲得n幅SAR強(qiáng)度序列和n(n－1)/2個(gè)干涉數(shù)據(jù)集；其中，n為大于1的正整數(shù)。

3.如權(quán)利要求1所述的一種InSAR時(shí)序相位的優(yōu)化方法，其特征在于，所述以所述L視強(qiáng)度數(shù)據(jù)集為參考，獲取并分別對(duì)預(yù)設(shè)數(shù)字高程模型和預(yù)設(shè)土地覆蓋圖進(jìn)行配準(zhǔn)與地理編碼，分別得到SAR圖像坐標(biāo)系下的數(shù)字高程和SAR圖像坐標(biāo)系下的土地覆蓋圖，具體為：

以所述L視強(qiáng)度數(shù)據(jù)集為參考，采用預(yù)設(shè)的互相關(guān)最大化算法，獲取并對(duì)預(yù)設(shè)數(shù)字高程模型與所述L視強(qiáng)度數(shù)據(jù)集進(jìn)行配準(zhǔn)和地理編碼，得到SAR圖像坐標(biāo)系下的數(shù)字高程；獲取并對(duì)預(yù)設(shè)的衛(wèi)星軌道數(shù)據(jù)信息與預(yù)設(shè)土地覆蓋圖進(jìn)行配準(zhǔn)和地理編碼，得到SAR圖像坐標(biāo)系下的土地覆蓋圖。

4.如權(quán)利要求2所述的一種InSAR時(shí)序相位的優(yōu)化方法，其特征在于，所述根據(jù)所述SAR圖像坐標(biāo)系下的數(shù)字高程，對(duì)所述干涉數(shù)據(jù)集進(jìn)行差分操作，得到差分干涉數(shù)據(jù)集，具體為：

根據(jù)所述SAR圖像坐標(biāo)系下的數(shù)字高程，對(duì)所述n(n－1)/2個(gè)干涉數(shù)據(jù)集進(jìn)行差分操作，得到n(n－1)/2幅差分干涉數(shù)據(jù)集；其中，所述差分干涉數(shù)據(jù)集與所述SAR圖像坐標(biāo)系下的土地覆蓋圖的像素位置一一對(duì)應(yīng)。

5.如權(quán)利要求4所述的一種InSAR時(shí)序相位的優(yōu)化方法，其特征在于，所述根據(jù)所述SAR圖像坐標(biāo)系下的土地覆蓋圖和所述差分干涉數(shù)據(jù)集，估算所述差分干涉數(shù)據(jù)集的協(xié)方差矩陣，具體為：

依次對(duì)所述n(n－1)/2幅差分干涉數(shù)據(jù)集建立滑動(dòng)窗口，根據(jù)所述SAR圖像坐標(biāo)系下的土地覆蓋圖，依次標(biāo)記每幅差分干涉數(shù)據(jù)集對(duì)應(yīng)的滑動(dòng)窗口內(nèi)與預(yù)設(shè)參考像素相同屬性的像素，從而依次統(tǒng)計(jì)所述滑動(dòng)窗口內(nèi)與中心像素具有相同類(lèi)別標(biāo)記的像素，計(jì)算協(xié)方差矩陣的元素，從而估算出所述差分干涉數(shù)據(jù)集的協(xié)方差矩陣。

6.如權(quán)利要求1所述的一種InSAR時(shí)序相位的優(yōu)化方法，其特征在于，所述根據(jù)預(yù)設(shè)的時(shí)序相位極大似然估計(jì)公式，估計(jì)并獲得優(yōu)化后的時(shí)序相位，具體為：

根據(jù)預(yù)設(shè)的時(shí)序相位極大似然估計(jì)公式：從而估計(jì)并獲得優(yōu)化后的時(shí)序相位；

其中，為優(yōu)化后的時(shí)序相位，θ為時(shí)序相位的正值，表示取實(shí)部算子，∑為所述差分干涉數(shù)據(jù)集的協(xié)方差矩陣。