[發明專利]一種識別濾線柵安裝狀態并進行柵影抑制的方法

有效

• 申請號: 201811275720.3
• 申請日: 2018-10-30
• 公開（公告）號: CN111110259B
• 公開（公告）日: 2023-06-30
• 發明（設計）人: 申果;金榮飛;周傳麗
• 申請（專利權）人: 遼寧開普醫療系統有限公司
• 主分類號: A61B6/00
• 分類號: A61B6/00
• 代理公司: 沈陽科苑專利商標代理有限公司 21002
• 代理人: 王倩

本發明公開了一種識別濾線柵安裝狀態并進行柵影抑制的方法，所述方法包括：獲取DR原始圖像；從原始圖像中抽取關鍵數據；根據關鍵數據識別濾線柵安裝狀態；對已有濾線柵圖像計算濾線柵柵影頻率；利用濾線柵柵影頻率創建濾波器；根據濾波器進行柵影抑制；輸出處理后圖像。本發明采用的方法可以有效識別出濾線柵安裝狀態，便于反饋濾線柵是否正確安裝；該發明同時對已進行安裝的濾線柵裝置所產生的柵影進行抑制，從而提高圖像質量。

技術領域

本發明涉及X射線數字圖像處理技術領域，特別是涉及一種識別濾線柵安裝狀態并進行柵影抑制的方法。

背景技術

在X射線成像過程中，入射X線光子與人體發生相互作用，并產生散射X線光子。散射光子會給圖像增加一層隨機的灰度信息，降低影像中組織之間的對比度，增加閱片的難度。

濾線柵裝置可以濾除散射線，從而提高X光圖像的對比度。然而濾線柵本身也會在X光圖像上產生相應的影像，這種影像被稱作濾線柵影或濾線柵偽影。濾線柵柵影的存在會影響組織的正常顯示，增加醫生閱片的難度，影響醫生對病癥的診斷。

目前普遍使用的濾線柵柵影處理方法大多是在已知濾線柵已正確安裝的情況下，且濾線柵的頻率和方向已知且固定不可改變的情況下進行的。有些方法在計算柵影頻率時，僅考慮了高頻部分的柵頻，沒有考慮柵影頻率會出現在2個或2個以上柵頻段中。在使用頻域處理時，大多數方法采用直接對柵影頻率截斷、或是矩形窗卷積，該方法會導致變換中帶進了高頻干擾和泄漏，甚至出現負譜現象。有些本發明是完全不知道濾線柵安裝狀態的情況下，首先進行濾線柵安裝狀態的識別，如未安裝濾線柵、濾線柵安裝成反面，或是濾線柵安裝為縱向或是橫向。然后對確認濾線柵已正確安裝的情況，自動計算濾線柵頻率，采用濾波器進行濾線柵抑制。

發明內容

本發明的目的是提供一種未知濾線柵安裝狀態，且未知濾線柵型號和平板探測器像素大小的濾線柵狀態識別和柵影抑制方法，該方法能夠有效的識別出當前設備安裝濾線柵的狀態，且對已安裝濾線柵的圖像進行柵影頻段確定，并進行柵影抑制。

本發明提供了一種識別濾線柵安裝狀態的方法，包括：識別濾線柵是否已安裝、濾線柵正反面是否顛倒、濾線柵的安裝方向為縱向還是橫向的方法。本發明同時提供了一種對已正確安裝濾線柵的情況下，進行柵影抑制的方法。

本發明采用如下技術方案：一種識別濾線柵安裝狀態并進行柵影抑制的方法，包括以下步驟：

A.獲取原始圖像；

B.從原始數據中抽取關鍵行和關鍵列；

C.計算關鍵行和關鍵列的頻域曲線；

D.根據頻域曲線識別濾線柵安裝狀態；

E.對安裝濾線柵的圖像，計算濾線柵柵影頻率；

F.利用濾線柵柵影頻率創建濾波器；

H.利用濾波器進行柵影抑制。

所述步驟B包括：

設定初始需要計算濾線柵頻率所占原始圖像行數或列數的比率a、b、c，abc；設定計算圖像中行方向濾線柵頻率和列方向濾線柵頻率的個數；

對于行方向和列方向，分別根據設定的個數，在比率[a,b]和[b,c]中插值，使總的行方向比率個數、列方向比率個數等于設定的個數；

用圖像總行數乘以根據插值得到的行方向的各個比率，即得到參與濾線柵頻率計算的關鍵行；

用圖像總列數乘以根據插值得到的列方向的各個比率，即得到參與濾線柵頻率計算的關鍵列。

所述步驟C包括：

計算行方向柵頻率的長度L_{Row_FFT}，計算公式如下：