TW201624416A - 去除dr圖像靜態濾線柵柵影的方法 - Google Patents
去除dr圖像靜態濾線柵柵影的方法 Download PDFInfo
- Publication number
- TW201624416A TW201624416A TW103145350A TW103145350A TW201624416A TW 201624416 A TW201624416 A TW 201624416A TW 103145350 A TW103145350 A TW 103145350A TW 103145350 A TW103145350 A TW 103145350A TW 201624416 A TW201624416 A TW 201624416A
- Authority
- TW
- Taiwan
- Prior art keywords
- image
- grid
- frequency domain
- filter
- shadow
- Prior art date
Links
Landscapes
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本發明公開一種去除DR圖像靜態濾線柵柵影的方法，包括如下步驟：製作空域卷積核心模板，將其存儲在DR控制系統中；使用空域卷積核心模板與需處理圖像進行卷積，得到處理結果圖像。本發明的技術方案主要是爲了解决現有技術中濾線柵柵影使得醫學影像質量下降甚至影響醫生診斷以及現有解决方案算法複雜、可靠性不高的問題，採用本發明的技術方案，能夠有效的去除濾線柵柵影，提高醫學影像質量，更好的輔助醫生診斷治療。
Description
本發明涉及一種圖像處理技術，特別涉及一種去除DR圖像靜態濾線柵柵影的方法。
在X射線設備成像過程中散射線的存在使圖像整體對比度降低，影響圖像質量。而濾線柵的作用是將散射體發出的散射線在到達照片之前過濾掉。但部分濾線柵因爲材質和分辨率的原因會給最終成像結果帶來條紋狀僞影，稱之爲柵影，濾線柵的柵影在圖像中會影響正常組織顯示，進而影響醫生對病症的診斷。當前DR設備針對濾線柵的柵影問題一般使用數字圖像處理技術進行改善，通常分爲空域處理和頻域處理方法。空域處理採用圖像插值方法去除柵影，由於插值對圖像的像素進行修改，在組織區域勢必影響圖像的質量，進而影響醫生對病症的診斷。而現有頻域去柵影技術由於濾波器設置複雜，且僞影在頻域的能量分布並不能與通用的濾波器形態很好的匹配，因此也不能將僞影在頻域內的能量很好地衰减。
目前能夠檢索到的近似相關專利有：《一種濾線柵僞影檢測及抑制方法及其裝置》，申請人：深圳邁瑞，申請號：201010154719。該發明實用的技術方案通過在傅立葉頻譜特徵區域根據頻譜區域周圍鄰域內的頻譜信息對所述特徵頻譜區域進行插值處理；對插值後的特徵頻譜圖進行反傅立葉變換的方法一直濾線柵的柵影。但是其方法使用的傅立葉頻譜插值算法複雜度較高，插值不當很可能造成反變換後出現更多的人造僞影。
爲了解决上述的技術問題，本發明提出一種去除DR圖像靜態濾線柵柵影的方法，本發明採用如下技術方案來解决上述技術問題：
一種去除DR圖像靜態濾線柵柵影的方法，其包括如下步驟：
a. 製作空域卷積核心模板，將其存儲在DR控制系統中；
b. 使用空域卷積核心模板與需處理圖像進行卷積，得到處理結果圖像。
進一步，在上述的一種去除DR圖像靜態濾線柵柵影的方法中，該步驟a包括如下步驟：
a1. 採集空曝光圖像：採集一張沒有放置任何物體的空曝光DR圖像；
a2. 判別柵影頻譜區域：在得到的空曝光圖像中進行傅立葉變換，將圖像從時空域轉換到頻率域；以頻譜原點爲中心，向水平和垂直兩個方向查找濾線柵的柵影在頻譜中的位置，並測量柵影在頻譜圖像中所佔面積；
a3. 生成頻域濾波器：得到濾線柵柵影在頻譜圖像中的位置和面積後，在柵影頻譜位置的中心使用頻域低通濾波器，而濾波器的半徑大於柵影頻譜位置的中心值邊緣的最大距離；在非柵影頻譜位置區域使用1進行填充。使得在柵影頻譜位置爲低通濾波，在其他位置爲1的頻域濾波器；
頻域低通濾波器公式如下：
其中u,v爲頻率域坐標，，M、N爲原圖像長寬，爲頻域下柵影特徵的直徑。
a4. 將頻域濾波器轉換爲空域卷積核心模板：對上一步産生的頻域濾波器進行傅立葉反變換，得到其對應的空域卷積核心模板。
與現有技術相比，本發明具有如下有益效果：
使用本發明的技術方案，對應DR濾線柵的柵影一般呈水平或者垂直的特點，首先通過一張空曝光圖像對圖像柵影在頻譜中的位置進行自動檢測。然後生成與柵影頻譜對應的濾波器，接著對濾波器進行反變換得到其柵影對應的空域卷積核心。最後每次實際拍攝病人時將此低通濾波器應用於圖像中達到完美去除柵影且不影響圖像質量的目的。
下面結合附圖對本發明作進一步詳細說明，但是本發明不僅僅侷限於以下實施例。如第1圖所述，本發明所述的一種去除DR圖像靜態濾線柵柵影的方法的優選實施例具體包括以下步驟：
1.採集空曝光圖像S1
首先採集一張沒有放置任何物體的空曝光DR圖像。
2.柵影頻譜區域判別
在得到的空曝光圖像中進行傅立葉變換S2，將圖像從時空域轉換到頻率域。以頻譜原點爲中心，向水平和垂直兩個方向查找濾線柵的柵影在頻譜中的位置，並測量柵影在頻譜圖像中所佔面積。
3.頻域濾波器生成
得到濾線柵柵影在頻譜圖像中的位置和面積後S3，在柵影頻譜位置的中心使用頻域低通濾波器，而濾波器的半徑大於柵影頻譜位置的中心值邊緣的最大距離。在非柵影頻譜位置區域使用1進行填充。使得在柵影頻譜位置爲低通濾波，在其他位置爲1的頻域濾波器。頻域低通濾波器公式如下：
其中u,v爲頻率域坐標，，M、N爲原圖像長寬，爲頻域下柵影特徵的直徑。
4.頻域濾波器轉換爲空域卷積核心
對上一步産生的頻域濾波器S4進行傅立葉反變換S5，得到其對應的空域卷積核心模板。
5.空域卷積核心模板使用S6
得到空域卷積核心目標之後，將其存儲在DR控制系統中，拍攝實際病人圖像S7時，在採集到圖像後首先將採集圖像與空域卷積核心模板進行卷積就可以消除濾線柵柵影對圖像造成的不良影響S8。
使用本發明的技術方案，對應DR濾線柵的柵影一般呈水平或者垂直的特點，首先通過一張空曝光圖像對圖像柵影在頻譜中的位置進行自動檢測。然後生成與柵影頻譜對應的濾波器，接著對濾波器進行反變換得到其柵影對應的空域卷積核心。最後每次實際拍攝病人時將此低通濾波器應用於圖像中達到完美去除柵影且不影響圖像質量的目的。
以上所述僅爲本發明的較佳實施例而已，並不用以限制本發明，凡在本發明的精神和原則之內所作的任何修改、等同替換和改進等，均應包含在本發明的保護範圍之內。
S1~S8‧‧‧步驟
第1圖爲本發明流程方框圖；
第2圖爲使用本發明方法處理前的DR圖像；
第3圖爲使用本發明方法處理後的DR圖像；
第4圖爲使用本發明方法處理前的DR圖像細節；
第5圖爲使用本發明方法處理後的DR圖像細節。
S1~S8‧‧‧步驟
Claims (2)
- 一種去除DR圖像靜態濾線柵柵影的方法，其包括如下步驟： a. 製作一空域卷積核心模板，將其存儲在一DR控制系統中；以及 b. 使用該空域卷積核心模板與需處理圖像進行卷積，得到處理結果圖像。
- 如申請專利範圍第1項所述之去除DR圖像靜態濾線柵柵影的方法，其中該步驟a包括如下步驟： a1. 採集空曝光圖像：採集一張沒有放置任何物體的空曝光DR圖像； a2. 判別柵影頻譜區域：在得到的空曝光圖像中進行傅立葉變換，將圖像從時空域轉換到頻率域；以頻譜原點爲中心，向水平和垂直兩個方向查找該濾線柵的該柵影在頻譜中的位置，並測量該柵影在頻譜圖像中所佔面積；以及 a3. 生成頻域濾波器：得到該濾線柵之該柵影在頻譜圖像中的位置和面積後，在柵影頻譜位置的中心使用頻域低通濾波器，而濾波器的半徑大於柵影頻譜位置的中心值邊緣的最大距離；在非柵影頻譜位置區域使用1進行填充；使得在柵影頻譜位置爲低通濾波，在其他位置爲1的頻域濾波器； 頻域低通濾波器公式如下：其中u,v爲頻率域坐標，，M、N爲原圖像長寬，爲頻域下該柵影特徵的直徑； a4. 將頻域濾波器轉換爲該空域卷積核心模板：對上一步産生的頻域濾波器進行傅立葉反變換，得到其對應的該空域卷積核心模板。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| TW103145350A TW201624416A (zh) | 2014-12-24 | 2014-12-24 | 去除dr圖像靜態濾線柵柵影的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| TW103145350A TW201624416A (zh) | 2014-12-24 | 2014-12-24 | 去除dr圖像靜態濾線柵柵影的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| TW201624416A true TW201624416A (zh) | 2016-07-01 |
Family
ID=56984793
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| TW103145350A TW201624416A (zh) | 2014-12-24 | 2014-12-24 | 去除dr圖像靜態濾線柵柵影的方法 |
Country Status (1)
| Country | Link |
|---|---|
| TW (1) | TW201624416A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111110259A (zh) * | 2018-10-30 | 2020-05-08 | 辽宁开普医疗系统有限公司 | 一种识别滤线栅安装状态并进行栅影抑制的方法 |
-
2014
- 2014-12-24 TW TW103145350A patent/TW201624416A/zh unknown
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111110259A (zh) * | 2018-10-30 | 2020-05-08 | 辽宁开普医疗系统有限公司 | 一种识别滤线栅安装状态并进行栅影抑制的方法 |
| CN111110259B (zh) * | 2018-10-30 | 2023-06-30 | 辽宁开普医疗系统有限公司 | 一种识别滤线栅安装状态并进行栅影抑制的方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102737392B (zh) | 一种低剂量x线ct图像的非局部正则化先验重建方法 | |
| KR101389841B1 (ko) | 파노라마 영상 데이터 제공 방법 및 장치 | |
| JP2019051315A (ja) | 医用画像処理装置、医用画像撮像装置及び医用画像処理プログラム | |
| US8855394B2 (en) | Methods and apparatus for texture based filter fusion for CBCT system and cone-beam image reconstruction | |
| CN103810733B (zh) | 一种稀疏角度x射线ct图像的统计迭代重建方法 | |
| JP2014236810A (ja) | 画像処理装置、断層撮影装置、画像処理方法およびプログラム | |
| US20180211421A1 (en) | Method, x-ray unit and computer program product for determining a three-dimensional image data set | |
| CN107274350B (zh) | 用于减少x射线图像中的振铃效应的方法及系统 | |
| EP2191441A2 (en) | Projection-based removal of high-contrast objects | |
| CN104574411A (zh) | 一种基于离体牙cbct的龋齿图像处理方法 | |
| US9867586B2 (en) | Stereo X-ray tube based suppression of outside body high contrast objects | |
| GB2549592B (en) | Method and system for reducing grid line artifacts in x-ray image | |
| CN111223156A (zh) | 一种用于牙科锥束ct系统的金属伪影消除方法 | |
| CN105787887A (zh) | 一种去除dr图像静态滤线栅栅影的方法 | |
| TW201624416A (zh) | 去除dr圖像靜態濾線柵柵影的方法 | |
| KR20110020969A (ko) | 단층촬영 시스템과 토모신세시스 시스템에서 평판형 엑스선 감지기의 불량화소로 인한 영상 아티팩트 제거 방법 및 장치 | |
| JP6478774B2 (ja) | 撮像装置、撮像方法及びプログラム | |
| JP2014176565A (ja) | 画像処理装置、放射線撮影装置、画像処理方法、コンピュータプログラム及び記憶媒体 | |
| CN114299081B (zh) | 一种上颌窦cbct图像分割方法、装置、存储介质及电子设备 | |
| KR101034270B1 (ko) | 고대조도 단층영상을 재구성하기 위한 장치 및 그 방법 | |
| KR102003197B1 (ko) | 단일 방사선 영상의 산란선 교정 시스템 및 방법 | |
| CN111904445A (zh) | 一种乳腺对比增强图像的获取及处理方法、装置 | |
| TWI593392B (zh) | Metal Detection and Artifact Removal Methods in Computerized Tomography Images | |
| Zhang et al. | Metal artifacts reduction for tomosynthesis | |
| US20130051514A1 (en) | Radiographic device |