[發明專利]一種基于方向金字塔濾波的圖像處理方法及裝置

無效

• 申請號: 201210468306.0
• 申請日: 2012-11-20
• 公開（公告）號: CN102945548A
• 公開（公告）日: 2013-02-27
• 發明（設計）人: 劉力輝
• 申請（專利權）人: 成都晶石石油科技有限公司
• 主分類號: G06T5/00
• 分類號: G06T5/00
• 代理公司: 成都行之專利代理事務所(普通合伙) 51220
• 代理人: 梁田

1.一種基于方向金字塔濾波的圖像處理方法，包括如下步驟，

步驟1.對輸入圖像信號XIN進行傅里葉變換后輸出圖像信號A₁

步驟2.對圖像信號A₁，分別經過低通濾波器組（L₀…L_M-1）的前兩個低通濾波器L₀和L₁，得到濾波后的圖像信號P₀，P₁，所述低通濾波器組（L₀…L_M-1）中各低通濾波器的通帶的帶寬是等比遞減的；

步驟3.用P₁減去P₀得到圖像信號VB₁,將圖像信號VB₁分別通過方向濾波器組(B₀…B_N)，后再經過反傅里葉變換，得到N+1個帶通濾波圖像信號Q₀₀~Q_0N；

步驟4.包括M-2個如下的循環步驟：

從m=2開始，將上一步的圖像信號P_m-1通過低通濾波器組的第m+1個低通濾波器L_m,得到圖像信號P_m；

用圖像信號P_m減去P_m-1，得到圖像信號VB_m,將圖像信號VB_m分別通過方向濾波器組(B₀…B_N)后再經過一次反傅里葉變換，得到(N+1)個帶通濾波圖像信號Q_(m-1)0~Q_(m-1)N；

步驟4中M為預先設定的尺度分解參數,是大于2的整數，m為大于1但不大于M的整數，在各個循環步驟中，m依次等于2,3，…,（M-1）；

步驟5.步驟4中最后一次循環得到的圖像信號P_M-1直接進行反傅里葉變換，得到粗尺度層圖像信號Q；

步驟6.將步驟3和步驟4中得到的全部（N+1）*（M-1）個帶通濾波圖像信號與步驟5中得到的粗尺度層圖像信號Q相加得到最終圖像信號；

其中步驟3和4中所述方向濾波器組(B₀…B_N)符合下式：

??-----②

,為方向濾波器的輸入信號頻率在x，y方向的頻率分量，②式中j為不大于N的非負整數；

上述各步驟中N為預先設定的方向分解參數，N為正整數。

2.如權利要求1所述的基于方向金字塔濾波的圖像處理方法，其特征在于：所述步驟2中的低通濾波器組（L₀…L_M-1）的各個濾波器符合：

??------①

其中，k為輸入信號的頻率，分別為低通濾波器組（L₀…L_M-1）的各低通濾波器錐形區域的起始點和截止點,對第j個低通濾波器Lj，，KA,KB為輸入圖像信號的頻率范圍端點，且KA>KB,β為預先設定的截止波數比；

j為不大于(M-1)的非負整數。

3.如權利要求2所述的基于方向金字塔濾波的圖像處理方法，其特征在于：所述步驟2中的截止波數比β=2。

4.如權利要求1所述的基于方向金字塔濾波的圖像處理方法，其特征在于：所述方向分解參數N=2。

5.如權利要求1所述的基于方向金字塔濾波的圖像處理方法，其特征在于：所述尺度分解參數M=3。

6.如權利要求1所述的基于方向金字塔濾波的圖像處理方法，其特征在于：所述步驟6在對各個圖像信號合成相加之前，還包括對各個圖像信號的去噪處理，去噪處理采用如下的步驟：

其中式中分別為去噪后的信號和去噪前的信號；

σn為噪聲的方差，σx?為無噪信號的方差；

????????，???

其中median()函數是中值估計函數，M為以為中心點，?2N_k為邊長的正方形區域內的數據點個數。N_k為預先設定的誤差評估區域參數。