[發明專利]X光機、X光圖像中靜態濾線柵偽影的去除方法及裝置

有效

• 申請號: 201410275968.5
• 申請日: 2014-06-19
• 公開（公告）號: CN105335936B
• 公開（公告）日: 2018-07-31
• 發明（設計）人: 唐慧;仝丹;李莉
• 申請（專利權）人: 上海西門子醫療器械有限公司
• 主分類號: G06T5/00
• 分類號: G06T5/00
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種X光圖像中靜態濾線柵偽影的去除方法，包括：

A、根據靜態濾線柵的柵線間隔和采樣頻率f_s確定圖像中柵線的預估頻率f_estimate；

B、對含柵線偽影的圖像I進行二維小波分解，得到四個小波分量：水平方向和豎直方向均是低頻的LL分量A_i、豎直方向高頻和水平方向低頻的HL分量H_i、水平方向高頻和豎直方向低頻的LH分量V_i、和水平方向和豎直方向均是高頻的HH分量D_i；其中，i＝1，為小波分解的層數；所述水平方向對應柵線偽影的方向；

C、根據所述柵線的預估頻率f_estimate，判斷柵線偽影信號是位于LL分量A_i上還是位于HL分量H_i上，并將LL分量A_i和HL分量H_i中含有柵線偽影信號的分量記為含柵線分量G_i；

D、計算在所述含柵線分量G_i中，柵線偽影信號是否為主要信號，如果是，則執行步驟F；否則，執行步驟E；

E、對所述含柵線分量G_i進行二維小波分解，并令i＝i+1，得到四個小波分量：水平方向和豎直方向均是低頻的LL分量A_i、豎直方向高頻和水平方向低頻的HL分量H_i、水平方向高頻和豎直方向低頻的LH分量V_i和水平方向和豎直方向均是高頻的HH分量D_i；并返回執行步驟C；

F、設計一濾波器對所述含柵線分量G_i進行濾波，得到去除柵線的分量G_i’；所述濾波器的頻譜函數二階連續，有且僅有一個拐點，且所述濾波器對應的時域函數同樣二階連續，有且僅有一個拐點；

G、將所述去除柵線的分量G_i’和其它保留的小波分量逐層進行小波反變換，得到恢復的去除柵線偽影的圖像I’。

2.根據權利要求1所述的方法，其特征在于，所述步驟A包括：

根據靜態濾線柵的柵線間隔確定柵線的理論頻率f_g，根據圖像的像素分辨率計算得到采樣頻率f_s；

判斷所述采樣頻率f_s是否滿足奈奎斯特采樣率，使得f_s≥2f_g，如果是，則確定柵線的預估頻率f_estimate＝f_g；否則，根據如下表達式確定柵線的預估頻率：

其中，k₁和k₂為整數，且k₁、k₂分別滿足下列不等式：

3.根據權利要求1所述的方法，其特征在于，所述步驟C包括：

根據所述采樣頻率f_s計算頻率分辨率f′＝f_s/N，則在頻域中所述柵線的預估頻率f_estimate對應的位置為其中N為采樣點數；

判斷所述柵線的預估頻率f_estimate對應的位置N_p與N/2的大小，如果N_p＜N/2，則確定柵線偽影信號位于LL分量A_i上，G_i＝A_i；如果N_p≥N/2，則確定柵線偽影信號位于HL分量H_i上，G_i＝H_i。