[發明專利]CT圖像偽影去除方法

有效

• 申請號: 202110059221.6
• 申請日: 2021-01-14
• 公開（公告）號: CN112785520B
• 公開（公告）日: 2022-10-21
• 發明（設計）人: 徐亦飛;陳偉;潘華素
• 申請（專利權）人: 明峰醫療系統股份有限公司
• 主分類號: G06T5/00
• 分類號: G06T5/00;G06T5/50
• 代理公司: 紹興市越興專利事務所(普通合伙) 33220
• 代理人: 蔣衛東

1.一種CT圖像偽影去除方法，其特征在于包括如下步驟：

(1)輸入圖像，對輸入圖像進行濾波得到帶有邊緣保護區域的圖像紋理；

(2)獲得圖像紋理在頻率閾的相位分布和幅值分布；

(3)計算幅值分布中相同頻率成分的均值并做歸一化，得到歸一化的幅值分布，并根據設定的閾值和低頻掩膜得到峰值的分布；

(4)對峰值分布進行高斯濾波得到濾波矩陣；將濾波矩陣與紋理幅值分布做矩陣點乘，得到無混疊的紋理頻域圖像；

(5)結合步驟(2)圖像紋理在頻率閾的相位分布，對無混疊的紋理頻域圖像處理，得到無混疊的紋理空間域圖像；然后將步驟(1)得到的圖像紋理與無混疊的紋理空間域圖像相減，得到混疊空間域圖像；

(6)將步驟(1)輸入圖像減去混疊偽影圖像，得到消除混疊偽影之后的圖像；

所述步驟(1)中，首先對輸入圖像分別做中值濾波和拉普拉斯濾波，將輸入圖像與中值濾波的結果相減得到圖像紋理分布；輸入圖像通過拉普拉斯濾波可以得到圖像的邊緣分布；根據拉普拉斯濾波的結果，將圖像紋理分布中較強的邊緣設為保護區域，得到帶有保護區域的圖像紋理分布。

2.如權利要求1所述的一種CT圖像偽影去除方法，其特征在于：所述步驟(2)中，對帶有邊緣保護區域的圖像紋理分布做二維傅里葉變換，得到圖像紋理分布在頻率域的相位分布和幅值分布。

3.如權利要求1所述的一種CT圖像偽影去除方法，其特征在于：所述步驟(4)中，濾波核為51x51像素的高斯形凹陷。

4.如權利要求2所述的一種CT圖像偽影去除方法，其特征在于：或者，所述步驟(5)為：將步驟(1)得到的圖像紋理減去無混疊的紋理頻域圖像得到混疊頻域圖像；結合步驟(2)圖像紋理在頻率閾的相位分布，對混疊頻域圖像處理得到混疊空間域圖像。

5.如權利要求4所述的一種CT圖像偽影去除方法，其特征在于：所述步驟(5)中的圖像處理具體為對圖像做二維傅里葉逆變換。