[發明專利]一種圖像處理方法及裝置

有效

• 申請號: 202010526199.7
• 申請日: 2020-06-09
• 公開（公告）號: CN111784617B
• 公開（公告）日: 2023-08-15
• 發明（設計）人: 孫凌;張群燕;張鵬
• 申請（專利權）人: 國家衛星氣象中心（國家空間天氣監測預警中心）
• 主分類號: G06T5/50
• 分類號: G06T5/50;G06T5/00
• 代理公司: 北京路浩知識產權代理有限公司 11002
• 代理人: 馬瑞

1.一種圖像處理方法，其特征在于，包括：

獲取遙感圖像中的待去條紋域，其中所述待去條紋域中包含待去條紋的像素；

對所述待去條紋域中的原始圖像進行迭代分解，得到迭代后的存在條紋的第一目標分量圖像和不存在條紋的第二目標分量圖像；

所述對所述待去條紋域中的原始圖像進行迭代分解，得到迭代后的存在條紋的第一目標分量圖像和不存在條紋的第二目標分量圖像，包括：

通過迭代計算，從所述待去條紋域中的原始圖像中逐步提取出不存在條紋的分量圖像，并將逐步提取得到的不存在條紋的分量圖像進行合并得到所述第二目標分量圖像；

所述通過迭代計算，從所述待去條紋域中的原始圖像中逐步提取出不存在條紋的分量圖像，包括：

將存在條紋的分量圖像初始化為所述待去條紋域中的原始圖像，并將不存在條紋的分量圖像初始化為0；

通過下述第一公式，對所述存在條紋的分量圖像計算受限梯度的拉普拉斯算子：

L_i(x,y)＝[v_i-1(x-1,y)-2v_i-1(x,y)+v_i-1(x+1,y)]+[v_i-1(x,y+1)

-v_i-1(x,y)]M(x,y)+[v_i-1(x,y-1)-v_i-1(x,y)]M(x,y-1)；

通過下述第二公式，對所述不存在條紋的分量圖像計算得到離散泊松方程：

L_i(x,y)＝u_i(x-1,y)+u_i(x+1,y)+u_i(x,y-1)+u_i(x,y+1)-4u_i(x,y)；

對所述第一公式進行離散傅里葉變換，得到頻域方程：

L_i(k_x,k_y)＝[2cos(πk_x/N_x)+2cos(πk_y/N_y)-4]u_i(k_x,k_y)；

基于所述頻域方程和所述第二公式，得到第i次迭代后的不存在條紋的分量圖像u_i(x,y)；

基于下述第三公式，從存在條紋的分量圖像中剔除不存在條紋的分量圖像，得到第i次迭代后的存在條紋的分量圖像：

v_i(x，y)＝v_i-1(x，y)-u_i(x，y)；

其中，L_i(x，y)表示第i次迭代時所對應的拉普拉斯算子，v_i-1(x，y)表示第(i-1)次迭代后所對應的存在條紋的分量圖像，M(x，y)的取值為1，u_i(x，y)表示第i次迭代后所對應的不存在條紋的分量圖像，v_i(x，y)表示第i次迭代后所對應的存在條紋的分量圖像，L_i(k_x，k_y)表示頻域方程；

將所述待去條紋域中的原始圖像中除所述第二目標分量圖像之外的圖像確定為所述第一目標分量圖像；

對所述第一目標分量圖像進行處理，得到低頻圖像；

所述對所述第一目標分量圖像進行處理，得到低頻圖像，包括：

采用非線性濾波器，通過下述公式，對所述第一目標分量圖像進行處理，得到所述低頻圖像：

C1＝e^C2；

C2＝-[V_N(x，z)-V_N(x，y)]²/(2δ²)；

其中，V_N(x，y)表示所述低頻圖像，x表示掃描方向上的像素，y表示跨掃描方向上的像素，H表示跨掃描方向上的濾波域范圍；

將所述低頻圖像和所述第二目標分量圖像進行合并，得到目標圖像。