[發明專利]一種基于空間分割的矩匹配衛星影像條帶噪聲去除方法

有效

• 申請號: 201310577168.4
• 申請日: 2013-11-18
• 公開（公告）號: CN103679648A
• 公開（公告）日: 2014-03-26
• 發明（設計）人: 邢坤;岳春宇;劉爽;何紅艷;周楠;鮑云飛;齊文雯;李巖;李方琦
• 申請（專利權）人: 北京空間機電研究所
• 主分類號: G06T5/00
• 分類號: G06T5/00
• 代理公司: 中國航天科技專利中心 11009
• 代理人: 安麗

1.一種基于空間分割的矩匹配衛星影像條帶噪聲去除方法，其特征在于步驟如下：

1）對待處理衛星影像的每個像素(x,y)建立鄰域范圍k×k，其中k為正奇數，即以像素點(x,y)為中心，選擇k×k大小的區域L；獲取區域L中所有像素的灰度值的平均值f(x,y)，獲取區域L中所有像素的灰度值的中間值g(x,y)，獲取梯度G(x,y)=Gx2+Gy2,]]>其中Gx=-101-202-101*L,Gy=121000-1-2-1*L,]]>其中*為卷積符號；

以每個像素點的均值f、中值g、梯度G為該像素點的三維坐標，映射到以f、g、G分別為三軸的三維直角坐標系中，構建待處理衛星影像的三維共生矩陣模型；并設置分割次數N，其中N為大于1的正整數；

2）在步驟1）建立的三維共生矩陣模型中設置初始分割點s(f₀,g₀,G₀)；以過點s(f₀,g₀,G₀)且平行于fOg的平面、過點s(f₀,g₀,G₀)且平行于gOG的平面和過點s(f₀,g₀,G₀)且平行于fOG的平面將三維共生矩陣分為八個子空間；定義在fOg平面內的四個子空間為下層子空間，其余四個子空間為上層子空間；在上層四個子空間中，定義與G軸相交的子空間，以及與該子空間對稱的子空間為參考區域A₀和B₀；在下層四個子空間中，定義與G軸相交的子空間，以及與該子空間對稱的子空間為計算區域C₀和D₀；

3）沿G軸移動s點，在移動過程中實時獲取參考區域A₀、B₀內的熵的和，取和的最大值點為記G軸上的最佳分割點s_G(f₀,g₀,G₀')，并重新獲得參考區域A_0G和B_0G；沿f軸移動s_G點，在移動過程中實時獲取參考區域A_0G和B_0G內的熵的和，取和的最大值點記為f軸上的最佳分割點s_fG(f₀',g₀,G₀')，并重新獲得參考區域A_0fG和B_0fG；沿g軸移動s_fG點，在移動過程中實時獲取參考區域A_0fG和B_0fG的熵的和，取和的最大值點記為g軸上的最佳分割點s_fgG(f₀',g₀',G₀')，并重新獲得參考區域A_0fgG和B_0fgG；將最佳分割點s_fgG記為s'，將參考區域A_0fgG和B_0fgG，記為參考區域A₁和B₁，獲得三維共生矩陣第一級分割的八個子空間，并得到第一級分割的計算區域C₁和D₁；

4）在步驟3）中得到的計算區域C₁、D₁中的每個區域內重復步驟2）、步驟3），得到四個第二級分割的計算區域；

5）重復步驟4）直到設置的迭代次數N執行完畢，得到最終2N的計算區域；

6）計算步驟5）中獲得的每個計算區域所對應的像素灰度值的均值和方差，采用標準矩匹配方法，重新計算該區域所對應的像素的灰度值，最終去除待處理衛星影像的條帶噪聲。