CN110889878A - 结合图像平滑处理的ct图像环形伪影抑制方法 - Google Patents
结合图像平滑处理的ct图像环形伪影抑制方法 Download PDFInfo
- Publication number
- CN110889878A CN110889878A CN201911087706.5A CN201911087706A CN110889878A CN 110889878 A CN110889878 A CN 110889878A CN 201911087706 A CN201911087706 A CN 201911087706A CN 110889878 A CN110889878 A CN 110889878A
- Authority
- CN
- China
- Prior art keywords
- image
- annular
- smoothing
- strip
- artifact
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
本发明提供结合图像平滑处理的CT图像环形伪影抑制方法,属于医学图像处理技术领域。该结合图像平滑处理的CT图像环形伪影抑制方法包括如下步骤:S1:将直角坐标下的第一环形图像转化为极坐标下的第一条形图像,通过图像纹理平滑处理方法将第一条形图像处理为第二条形图像;S2:将极坐标下的第二条形图像转化为直角坐标下的第二环形图像,通过图像纹理平滑处理方法将第二环形图像处理为环形伪影抑制后的第三环形图像。本发明能够在很好的保护图像细节和图像清晰度的情况下,能够对图像中的环形伪影有比较好的抑制作用。
Description
技术领域
本发明属于医学图像处理技术领域,涉及结合图像平滑处理的CT图像环形伪影抑制方法。
背景技术
图像伪影是 CT 图像在重建过程中受到的一些干扰信号产生的或者是原来物体中所不存在的影像。CT 图像的伪影常常表现为是环状伪影。图像质量的下降正是由于这些环状伪影的出现会导致,从而造成一些误判的发生。虽然环形伪影的出现只是图像中的一个干扰,但是它如果出现在一些研究方向的重要位置上就可能会产生比较严重的问题。尤其在医学上,环状伪影的出现会影响 CT 图像中的一些病灶处的影像,甚至还会掩盖一些有效区域,从而导致医生对疾病的诊断,造成一些重大的医疗事故。而且 CT 设备是一个很复杂的系统,由于使用不当、保养不及时或未及时利用标准模型校准等均可引起一些伪影。由这些因素所导致的环形伪影通常很容易被看出,但如果所造成的伪影不明显或是由于一些机器伤影而产生的时,很容易把它当成图像本身的特点而造成误判,如当环形伪影出现在脊髓中央时,会导致诊断错误等。此外,在环形伪影中有一个比较特殊表现形式,即它在图像中表现为环的中央模糊,这种环形伪影一般表现在小部分范围内的。又由于它的表现形态与病理结构很相似,从而可能导致误判。
发明内容
本发明针对现有的技术存在的上述问题,提供结合图像平滑处理的CT图像环形伪影抑制方法,本发明所要解决的技术问题是:如何提供一种可以结合图像平滑处理CT图像以对CT图像中的环形伪影进行抑制的CT图像环形伪影抑制方法。
本发明的目的可通过下列技术方案来实现:
结合图像平滑处理的CT图像环形伪影抑制方法,包括如下步骤:
S1:将直角坐标下的第一环形图像转化为极坐标下的第一条形图像,通过图像纹理平滑处理方法将第一条形图像处理为第二条形图像;
S2:将极坐标下的第二条形图像转化为直角坐标下的第二环形图像,通过图像纹理平滑处理方法将第二环形图像处理为环形伪影抑制后的第三环形图像。
优选的,步骤S1中第一环形图像转化为第一条形图像,第一环形图像中的第一环形伪影转化为第一条形伪影,通过对第一条形图像进行图像纹理平滑处理以将第一条形图像中的第一条形伪影消除。
优选的,步骤S2中第二条形图像转化第二环形图像,第二条形图像中的第二条形伪影转化为第二环形伪影,通过对第二环形图像进行图像纹理平滑处理以将第二环形图像中的第二环形伪影消除。
优选的,图像纹理平滑处理过程具体包括:将图像窗口的固有变差和图像窗口
总变差结合起来,形成结构-纹理分解正则化器,结构-纹理分解正则化器的目标函数为,其中使得输入和输出的结果图像不发生变化,正则项为,正则项
为图像窗口的相对总变差,为权重值,像素级窗口总变差分别为、,窗口固
有变化分别为、,为根据空间相关关系定义的
权重函数,,为图像窗口的固有变化,大于0。
优选的,步骤S1中直角坐标下的第一环形图像转化为极坐标下的第一条形图像时对于超出第一环形图像的采样点通过最邻近插值法确定其超出预设范围的值。
优选的,步骤S2中极坐标下的第二条形图像转化为直角坐标下的第二环形图像时对于超出第二条形图像的位置通过最邻近插值法进行填补。
本发明中首先将直角坐标下的第一环形图像转化为极坐标下的第一条形图像,通过图像纹理平滑处理方法将第一条形图像处理为第二条形图像,第一环形图像中的部分环形伪影的程度得到减弱,接着将极坐标下的第二条形图像转化为直角坐标下的第二环形图像,通过图像纹理平滑处理方法将第二环形图像处理为环形伪影抑制后的第三环形图像,未能被完全抑制的环形伪影能够有效得到抑制,能够在很好的保护图像细节和图像清晰度的情况下,能够对图像中的环形伪影有比较好的抑制作用。
附图说明
图1是本发明的流程示意图。
具体实施方式
以下是本发明的具体实施例并结合附图,对本发明的技术方案作进一步的描述,但本发明并不限于这些实施例。
请参阅图1,本实施例中的结合图像平滑处理的CT图像环形伪影抑制方法,包括如下步骤:
S1:将直角坐标下的第一环形图像转化为极坐标下的第一条形图像,通过图像纹理平滑处理方法将第一条形图像处理为第二条形图像;
S2:将极坐标下的第二条形图像转化为直角坐标下的第二环形图像,通过图像纹理平滑处理方法将第二环形图像处理为环形伪影抑制后的第三环形图像。
此处,首先将直角坐标下的第一环形图像转化为极坐标下的第一条形图像,通过图像纹理平滑处理方法将第一条形图像处理为第二条形图像,第一环形图像中的部分环形伪影的程度得到减弱,接着将极坐标下的第二条形图像转化为直角坐标下的第二环形图像,通过图像纹理平滑处理方法将第二环形图像处理为环形伪影抑制后的第三环形图像,未能被完全抑制的环形伪影能够有效得到抑制,能够在很好的保护图像细节和图像清晰度的情况下,能够对图像中的环形伪影有比较好的抑制作用。
步骤S1中第一环形图像转化为第一条形图像,第一环形图像中的第一环形伪影转化为第一条形伪影,通过对第一条形图像进行图像纹理平滑处理以将第一条形图像中的第一条形伪影消除。
步骤S2中第二条形图像转化第二环形图像,第二条形图像中的第二条形伪影转化为第二环形伪影,通过对第二环形图像进行图像纹理平滑处理以将第二环形图像中的第二环形伪影消除。
图像纹理平滑处理过程可以具体包括:将图像窗口的固有变差和图像窗口总变
差结合起来,形成结构-纹理分解正则化器,结构-纹理分解正则化器的目标函数为,其中使得输入和输出的结果图像不发生变化,正则项为,正则项为图
像窗口的相对总变差,为权重值,像素级窗口总变差分别为、,窗口固有
变化分别为、,为根据空间相关关系定义的权重
函数,,为图
像窗口的固有变化,大于0,是一个正数,为了避免被除数出现为 0 的情况。只包含纹
理的窗口的结果的通常比包含结构边缘的窗口的结果的小。
步骤S1中直角坐标下的第一环形图像转化为极坐标下的第一条形图像时对于超出第一环形图像的采样点可以通过最邻近插值法确定其超出预设范围的值。
步骤S2中极坐标下的第二条形图像转化为直角坐标下的第二环形图像时对于超出第二条形图像的位置可以通过最邻近插值法进行填补。
参数的值可以大于0.09,这样可以保护图像中的细节部分。在对数据图像处理
时,参数在大于 0.025 和参数在小于 0.09 时,在对极坐标下处理后的图像再在直
角坐标下处理中,图像会变得很模糊,甚至图像中的细节部分也消失了,反之,图像中的环
形伪影处理效果不是特别的好。对数据图像处理时, 参数在大于 0.05和参数在小
于 0.03 时,在对极坐标下处理后的图像再在直角坐标下处理中,图像会变得很模糊,甚至
图像中的细节部分也消失了,反之,图像中的环形伪影处理效果不是特别的好。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (10)
1.结合图像平滑处理的CT图像环形伪影抑制方法,其特征在于,包括如下步骤:
S1:将直角坐标下的第一环形图像转化为极坐标下的第一条形图像,通过图像纹理平滑处理方法将第一条形图像处理为第二条形图像;
S2:将极坐标下的第二条形图像转化为直角坐标下的第二环形图像,通过图像纹理平滑处理方法将第二环形图像处理为环形伪影抑制后的第三环形图像。
2.如权利要求1所述的结合图像平滑处理的CT图像环形伪影抑制方法,其特征在于:步骤S1中第一环形图像转化为第一条形图像,第一环形图像中的第一环形伪影转化为第一条形伪影,通过对第一条形图像进行图像纹理平滑处理以将第一条形图像中的第一条形伪影消除。
3.如权利要求1或2所述的结合图像平滑处理的CT图像环形伪影抑制方法,其特征在于:步骤S2中第二条形图像转化第二环形图像,第二条形图像中的第二条形伪影转化为第二环形伪影,通过对第二环形图像进行图像纹理平滑处理以将第二环形图像中的第二环形伪影消除。
6.如权利要求1或2所述的结合图像平滑处理的CT图像环形伪影抑制方法,其特征在于:步骤S1中直角坐标下的第一环形图像转化为极坐标下的第一条形图像时对于超出第一环形图像的采样点通过最邻近插值法确定其超出预设范围的值。
7.如权利要求1或2所述的结合图像平滑处理的CT图像环形伪影抑制方法,其特征在于:步骤S2中极坐标下的第二条形图像转化为直角坐标下的第二环形图像时对于超出第二条形图像的位置通过最邻近插值法进行填补。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911087706.5A CN110889878B (zh) | 2019-11-08 | 2019-11-08 | 结合图像平滑处理的ct图像环形伪影抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911087706.5A CN110889878B (zh) | 2019-11-08 | 2019-11-08 | 结合图像平滑处理的ct图像环形伪影抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110889878A true CN110889878A (zh) | 2020-03-17 |
CN110889878B CN110889878B (zh) | 2023-06-09 |
Family
ID=69747126
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911087706.5A Active CN110889878B (zh) | 2019-11-08 | 2019-11-08 | 结合图像平滑处理的ct图像环形伪影抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110889878B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112233027A (zh) * | 2020-09-30 | 2021-01-15 | 西北工业大学 | 一种ct图像环形伪影的迭代后处理去除方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080019607A1 (en) * | 2006-07-21 | 2008-01-24 | Josh Star-Lack | System and method for correcting for ring artifacts in an image |
US20100284599A1 (en) * | 2008-01-11 | 2010-11-11 | Akinori Fujita | Image processing method, an apparatus therefor and a tomographic apparatus |
CN102274040A (zh) * | 2010-06-10 | 2011-12-14 | 东软飞利浦医疗设备系统有限责任公司 | 一种去除ct图像中环状伪影的后处理方法 |
CN102521801A (zh) * | 2011-11-23 | 2012-06-27 | 中国科学院深圳先进技术研究院 | Ct图像的环形伪影及弧形伪影校正方法 |
CN102800073A (zh) * | 2012-06-28 | 2012-11-28 | 西北工业大学 | 一种锥束ct环形伪影的自动判别与校正方法 |
CN105321155A (zh) * | 2015-10-29 | 2016-02-10 | 北京理工大学 | 一种cbct图像环形伪影消除方法 |
CN106296763A (zh) * | 2016-07-20 | 2017-01-04 | 中国兵器科学研究院宁波分院 | 一种金属材料工业ct图像质量快速校正方法 |
CN106651981A (zh) * | 2016-12-07 | 2017-05-10 | 深圳先进技术研究院 | 一种环形伪影修正的方法及装置 |
WO2018103015A1 (zh) * | 2016-12-07 | 2018-06-14 | 深圳先进技术研究院 | 一种环形伪影修正的方法及装置 |
CN109636872A (zh) * | 2018-12-10 | 2019-04-16 | 合肥中科离子医学技术装备有限公司 | 基于滑动窗口差值和条带噪声检测的cbct环形伪影消除方法 |
CN110060316A (zh) * | 2019-04-29 | 2019-07-26 | 重庆邮电大学 | 一种用于ct重建中多区域分割的环伪影校正方法 |
-
2019
- 2019-11-08 CN CN201911087706.5A patent/CN110889878B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080019607A1 (en) * | 2006-07-21 | 2008-01-24 | Josh Star-Lack | System and method for correcting for ring artifacts in an image |
US20100284599A1 (en) * | 2008-01-11 | 2010-11-11 | Akinori Fujita | Image processing method, an apparatus therefor and a tomographic apparatus |
CN102274040A (zh) * | 2010-06-10 | 2011-12-14 | 东软飞利浦医疗设备系统有限责任公司 | 一种去除ct图像中环状伪影的后处理方法 |
CN102521801A (zh) * | 2011-11-23 | 2012-06-27 | 中国科学院深圳先进技术研究院 | Ct图像的环形伪影及弧形伪影校正方法 |
CN102800073A (zh) * | 2012-06-28 | 2012-11-28 | 西北工业大学 | 一种锥束ct环形伪影的自动判别与校正方法 |
CN105321155A (zh) * | 2015-10-29 | 2016-02-10 | 北京理工大学 | 一种cbct图像环形伪影消除方法 |
CN106296763A (zh) * | 2016-07-20 | 2017-01-04 | 中国兵器科学研究院宁波分院 | 一种金属材料工业ct图像质量快速校正方法 |
CN106651981A (zh) * | 2016-12-07 | 2017-05-10 | 深圳先进技术研究院 | 一种环形伪影修正的方法及装置 |
WO2018103015A1 (zh) * | 2016-12-07 | 2018-06-14 | 深圳先进技术研究院 | 一种环形伪影修正的方法及装置 |
CN109636872A (zh) * | 2018-12-10 | 2019-04-16 | 合肥中科离子医学技术装备有限公司 | 基于滑动窗口差值和条带噪声检测的cbct环形伪影消除方法 |
CN110060316A (zh) * | 2019-04-29 | 2019-07-26 | 重庆邮电大学 | 一种用于ct重建中多区域分割的环伪影校正方法 |
Non-Patent Citations (2)
Title |
---|
LI XU ET AL: "Structure Extraction from Texture via Relative Total Variation", 《ACM TRANSACTIONS ON GRAPHICS》 * |
LI XU ET AL: "Structure Extraction from Texture via Relative Total Variation", 《ACM TRANSACTIONS ON GRAPHICS》, 1 November 2012 (2012-11-01), pages 1 - 7 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112233027A (zh) * | 2020-09-30 | 2021-01-15 | 西北工业大学 | 一种ct图像环形伪影的迭代后处理去除方法 |
CN112233027B (zh) * | 2020-09-30 | 2022-12-09 | 西北工业大学 | 一种ct图像环形伪影的迭代后处理去除方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110889878B (zh) | 2023-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3995854B2 (ja) | 画像処理方法および装置並びに記録媒体 | |
CN109801343B (zh) | 基于重建前后图像的环形伪影校正方法、ct控制系统 | |
JP2011028588A (ja) | 情報処理装置、ライン状ノイズ低減処理方法、及びプログラム | |
US20130121611A1 (en) | Image processing device and method and image display device | |
CN102521801B (zh) | Ct图像的环形伪影及弧形伪影校正方法 | |
JP6071444B2 (ja) | 画像処理装置及びその作動方法、プログラム | |
CN103413280A (zh) | 一种低剂量x射线ct图像重建方法 | |
WO2022253093A1 (zh) | 肠道内窥镜观察视频中图像的处理方法、装置及存储介质 | |
JP2012208553A (ja) | 画像処理装置、および画像処理方法、並びにプログラム | |
Zhang et al. | Image denoising and zooming under the linear minimum mean square-error estimation framework | |
CN110889878A (zh) | 结合图像平滑处理的ct图像环形伪影抑制方法 | |
CN111553960A (zh) | 一种基于投影均值图像的环状伪影快速校正方法 | |
CN105809633A (zh) | 去除颜色噪声的方法及装置 | |
CN103971335B (zh) | 一种基于置信度核回归的图像超分辨率重建方法 | |
CN106651981B (zh) | 一种环形伪影修正的方法及装置 | |
GB2549592B (en) | Method and system for reducing grid line artifacts in x-ray image | |
Li et al. | Application of adaptive histogram equalization to x-ray chest images | |
CN111091517A (zh) | 一种残差加权成像方法和装置 | |
Goyal et al. | SGS-SRAD filter for denoising and edge preservation of ultrasound images | |
KR20100097858A (ko) | 예제 기반 신경회로망을 이용한 고해상도 영상 확대 | |
Ziou et al. | Reducing aliasing in images: a PDE-based diffusion revisited | |
CN115588060A (zh) | 基于投影平均图像的快速ct环形伪影校正方法及系统 | |
JP2013250893A (ja) | 情報処理装置、情報処理方法およびプログラム | |
Ahmed | Image enhancement and noise removal by using new spatial filters | |
Nasonov et al. | Edge quality metrics for image enhancement |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |