[發明專利]大尺度地距星載SAR圖像鑲嵌方法

有效

• 申請號: 201810299312.5
• 申請日: 2018-03-30
• 公開（公告）號: CN108564532B
• 公開（公告）日: 2022-02-15
• 發明（設計）人: 郎文輝;余不凡;石聰聰;趙子航
• 申請（專利權）人: 合肥工業大學
• 主分類號: G06T3/40
• 分類號: G06T3/40;G06T5/00;G06T5/50;G06T7/32
• 代理公司: 安徽匯樸律師事務所 34116
• 代理人: 劉海涵

1.大尺度地距星載SAR圖像鑲嵌方法，包括多幅已編號地距星載SAR圖像，其特征在于：該方法包括以下步驟：

步驟1、圖像預處理；

對每個已編號地距星載SAR圖像進行參數修正后獲得含地理定位信息的待處理SAR圖像；

步驟2、計算每個待處理SAR圖像中每個像素的經、緯度坐標；

步驟3、判定所有待處理SAR圖像兩兩之間是否重疊；

根據步驟2所得到的每個待處理SAR圖像中每個像素的經、緯度坐標，判斷所有待處理SAR圖像兩兩之間是否重疊；

存在重疊的，記錄所有兩兩之間存在重疊的待處理SAR圖像對應的編號以及兩兩之間重疊區域范圍，并生成含有重疊區域范圍信息的編號組進行存儲；

步驟4、計算每個編號組對應重疊區域范圍之間的匹配偏移量；

對于峰均比低的區域，采用互相關法計算每個編號組對應重疊區域范圍之間的匹配偏移量；

對于峰均比高的區域，采用基于過零點的二值匹配統計法來計算每個編號組對應重疊區域范圍之間的匹配偏移量；

步驟5、創建全零像素圖像和全局結構體變量

根據步驟2中得到的待處理SAR圖像中每個像素的經、緯度坐標和步驟4得到的每個編號組對應重疊區域范圍之間的匹配偏移量獲得所有待處理SAR圖像修正后的地理坐標，得到修正后每個待處理SAR圖像對應的待拼接圖像；

創建與最終要生成的鑲嵌圖像同樣大小的全零像素圖像，同時建立一個全局結構體變量；

步驟6、采用Wallis濾波器對步驟5獲得的每個編號組對應的重疊區域范圍進行勻色；

步驟7、對勻色后的所有待拼接圖像進行拼接；

拼接方法如下：

a.在全局結構體變量中為所有待拼接圖像設置拼接勻色標記，并復位為0；

b.輸入一幅待拼接圖像，基于編號組找出所有與該幅待拼接圖像重疊且勻色標記為0的待拼接圖像；

c.假設與當前輸入待拼接圖像重疊且勻色標記為0的待拼接圖像總數為Num，以當前輸入待拼接圖像為參考圖像，依次按照步驟6對這Num幅待拼接圖像進行Wallis勻色，且每一次勻色均更新Wallis濾波器中相應的參數，將勻色標記置為1并保存勻色后的圖像；

d.Num次Wallis勻色后，將該輸入待拼接圖像寫入步驟5中創建的全零像素圖像中，并將拼接勻色標記設置為1；

e.對所有待拼接圖像重復上述b、c、d步驟，處理完所有待拼接圖像后，得到最終的鑲嵌結果。

2.根據權利要求1所述的大尺度地距星載SAR圖像鑲嵌方法，其特征在于，步驟1中所述參數修正具體為如下過程：

對每個已編號地距星載SAR圖像進行相干斑濾波，以去除相干斑噪聲，獲得無相干斑噪聲的SAR圖像；

對無相干斑噪聲的SAR圖像進行輻射定標，以消除系統輻射誤差，并將定標后的SAR圖像用后向散射系數表示，獲得用后向散射系數表示的SAR圖像；

對用后向散射系數表示的SAR圖像進行橢球體校正，并將輔助信息中的地理定位信息寫入SAR圖像，得到待處理SAR圖像。

3.根據權利要求2所述的大尺度地距星載SAR圖像鑲嵌方法，其特征在于，步驟3中每個待處理SAR圖像中每個像素的經、緯度坐標的計算公式如下：

E_geo＝A(0)+X_pix×A(1)+Y_pix×A(2) (1)

N_geo＝A(3)+X_pix×A(4)+Y_pix×A(5) (2)

其中：E_geo和N_geo分別表示所述待處理SAR圖像中像素的經度和緯度坐標；

X_pix和Y_pix分別表示所述待處理SAR圖像中像素的列坐標和行坐標；

所述地理定位信息具體由地理定位數組構成，所述地理定位數組具體由六個參數A(0),A(1),...A(5)構成，其中：

A(0)和A(3)表示分別表示所述待處理圖像左上角的經、緯度坐標；

A(1)和A(5)分別表示所述待處理圖像的橫向與縱向分辨率；

A(2)和A(4)分別表示所述待處理圖像中像素在緯、經度方向上的旋轉系數。