CN107316277B - 一种医学锥束ct图像中运动伪影修正方法 - Google Patents

一种医学锥束ct图像中运动伪影修正方法 Download PDF

Info

Publication number
CN107316277B
CN107316277B CN201710269022.1A CN201710269022A CN107316277B CN 107316277 B CN107316277 B CN 107316277B CN 201710269022 A CN201710269022 A CN 201710269022A CN 107316277 B CN107316277 B CN 107316277B
Authority
CN
China
Prior art keywords
image
artifact
cone beam
cone
shadow
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.)
Expired - Fee Related
Application number
CN201710269022.1A
Other languages
English (en)
Other versions
CN107316277A (zh
Inventor
牛田野
江阳康
薛一
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Suzhou Dongying Information Technology Co ltd
Original Assignee
Suzhou Dongying Information Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Suzhou Dongying Information Technology Co ltd filed Critical Suzhou Dongying Information Technology Co ltd
Priority to CN201710269022.1A priority Critical patent/CN107316277B/zh
Publication of CN107316277A publication Critical patent/CN107316277A/zh
Application granted granted Critical
Publication of CN107316277B publication Critical patent/CN107316277B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/168Segmentation; Edge detection involving transform domain methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/421Filtered back projection [FBP]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明公开了一种医学锥束CT图像中运动伪影修正方法,运动伪影主要包括阴影伪影和条状伪影,本方法首先利用病人或者模体锥束CT投影数据重建初始图像,其次利用图像分割技术和图像特征提取技术获取运动伪影存在的主要区域,然后利用线性插值技术消除阴影伪影,再利用频率分割技术消除条状伪影,实现锥束CT图像的运动伪影消除。本发明使用基于信号处理的方法对锥束CT图像的运动伪影进行修正,有效地消除了锥束CT图像中的运动伪影,提高了图像的质量,增加了锥束CT图像在图像引导放射治疗中的图像引导的准确性。

Description

一种医学锥束CT图像中运动伪影修正方法
技术领域:
本发明属于医学成像技术领域,具体是涉及一种医学锥束CT图像中运动伪影修正方法。
背景技术:
在图像引导放射治疗中,锥束CT成像技术是非常重要的图像引导技术。但是由于呼吸作用和肠道蠕动的存在,使得锥束CT图像存在较为严重的运动伪影,影像图像的质量,降低图像引导放射治疗中定位的准确性。
现有的阴影伪影修正方法主要有呼吸修正、屏息以及运动模型修正的方法,这些方法需要依靠复杂的数学模型进行运动估计,而且对于不规则运动的估计有局限性。
发明内容:
为此,本发明所要解决的技术问题在于现有技术中用于锥束CT图像技术的伪影修正方法需要依靠复杂的数学模型进行运动估计,且对于不规则运动的估计有局限性,从而本申请提出基于信号处理技术的一种医学锥束CT图像中运动伪影修正方法,该方法不需要依靠先验知识、原始投影以及复杂的运动模型,能够有效的消除图像中的运动伪影,提高锥束CT图像的质量。
为达到上述目的,本发明的技术方案如下:
一种医学锥束CT图像中运动伪影修正方法,包括如下步骤:
S1:初始图像重建,利用病人或者模体锥束CT投影数据重建初始图像f0
S2:区域分割提取,利用图像分割技术和图像特征提取技术对所述初始图像f0进行处理,获取气泡区域fgas、骨性组织区域fbone、阴影伪影与条状伪影存在的主要区域fr
S3:阴影伪影消除,利用线性插值技术对所述气泡区域fgas、所述阴影伪影与条状伪影存在的主要区域fr进行处理,获取阴影伪影修正后的第一锥束CT图像;
S4:条状伪影消除,利用频率分割技术对阴影伪影修正后的第一锥束CT图像、所述气泡区域fgas、所述骨性组织区域fbone进行处理,获取条状伪影修正后的第二锥束CT图像,实现锥束CT图像的运动伪影消除。
作为上述技术方案的优选,所述步骤S1中:
采用滤波反投影技术重建所述初始图像f0
作为上述技术方案的优选,所述步骤S1中:
所述初始图像f0中的运动伪影包括阴影伪影和条状伪影。
作为上述技术方案的优选,所述步骤S2中:
所述图像分割技术采用阈值分割技术,所述图像特征提取技术采用图像腐蚀和图像膨胀技术。
作为上述技术方案的优选,所述步骤S3具体包括如下步骤:
S31:将局部Radon变换作用于所述主要区域fr得到第一虚拟投影,即将图像域转换到虚拟投影域;
S32:将Radon变换作用于所述气泡区域fgas,获取含有阴影伪影污染的虚拟投影数据;
S33:在所述第一虚拟投影域中,利用线性插值技术替换含有阴影伪影污染的虚拟投影数据;
S34:重建替换后的虚拟投影数据,对重建后的虚拟投影数据进行低通滤波后得到fmask1,将fmask1加到所述初始图像f0中,获取阴影伪影修正后的第一锥束CT图像,即实现锥束CT图像的阴影伪影修正。
作为上述技术方案的优选,所述步骤S34中:
采用高斯滤波器对所述投影数据进行低通滤波处理。
作为上述技术方案的优选,所述步骤S4具体包括如下步骤:
S41:对步骤S3中获取的第一锥束CT图像进行频率分割,获取所述第一锥束CT图像的中频成份;
S42:将Radon变换作用于所述中频成份、所述气泡区域fgas和所述骨性组织区域fbone,获取第二虚拟投影;
S43:选取中频成份投影中的气泡成份和骨性组织成份,选取到的成份为条状伪影成份;
S44:重建第二虚拟投影中选取的条状伪影成份,得到fmask2,将步骤S3中获取的第一锥束CT图像减去fmask2,获取条状伪影修正后的第二锥束CT图像,实现锥束CT图像的运动伪影消除。
作为上述技术方案的优选,所述步骤S41中:
所述第一锥束CT图像的中频成份包括条状伪影、图像边界以及图像的纹理特征。
本发明的有益效果在于:本发明使用基于信号处理的方法对锥束CT图像的运动伪影进行修正,可以有效地消除锥束CT图像中的运动伪影,提高图像的质量,并且有望能够运用于现有的临床治疗过程中,提高锥束CT图像在图像引导放射治疗中的作用。
附图说明:
以下附图仅旨在于对本发明做示意性说明和解释,并不限定本发明的范围。其中:
图1为本发明一个实施例的一种医学锥束CT图像中运动伪影修正方法流程图;
图2为本发明一个实施例的一种医学锥束CT图像中阴影伪影消除的流程图;
图3本发明一个实施例的一种医学锥束CT图像中条状伪影消除的流程图;
图4为本发明一个实施例的运动伪影修正方法在病人腹部锥束CT数据上的实施结果示意图;
图5为本发明一个实施例的运动伪影修正方法在病人腹部锥束CT数据上的实施结果的矢状位图像示意图。
具体实施方式:
如图1所示,本发明的一种医学锥束CT图像中运动伪影修正方法,包括如下步骤:
S1:使用滤波反投影技术产生初始图像,利用实际测得的病人或者模体锥束CT投影数据使用滤波反投影技术重建出初始图像f0。由于呼吸作用和肠道蠕动,所述初始图像f0中往往含有严重的运动伪影,运动伪影包括阴影伪影和条状伪影。
S2:区域分割提取,利用图像分割技术和图像特征提取技术对所述初始图像f0进行处理,获取气泡区域fgas、骨性组织区域fbone、阴影伪影与条状伪影存在的主要区域fr。通过阈值分割从初始图像f0中得到气泡区域fgas和骨性组织区域fbone,通过图像腐蚀与图像膨胀得到阴影伪影与条状伪影存在的主要区域fr
S3:阴影伪影消除,利用线性插值技术对所述气泡区域fgas、所述阴影伪影与条状伪影存在的主要区域fr进行处理,获取阴影伪影修正后的第一锥束CT图像;
如图2所示,所述步骤S3具体包括如下步骤:
S31:将局部Radon变换作用于所述主要区域fr得到第一虚拟投影,即将图像域转换到虚拟投影域;
S32:将Radon变换作用于所述气泡区域fgas,获取含有阴影伪影污染的虚拟投影数据;
S33:在所述第一虚拟投影域中,利用线性插值技术替换含有阴影伪影污染的虚拟投影数据;
S34:重建替换后的虚拟投影数据,对重建后的虚拟投影数据进行低通滤波后得到fmask1,将fmask1加到所述初始图像f0中,获取阴影伪影修正后的第一锥束CT图像,即实现锥束CT图像的阴影伪影修正。本实施例中,采用高斯滤波器对所述投影数据进行低通滤波处理。
S4:条状伪影消除,利用频率分割技术对阴影伪影修正后的第一锥束CT图像、所述气泡区域fgas、所述骨性组织区域fbone进行处理,获取条状伪影修正后的第二锥束CT图像,实现锥束CT图像的运动伪影消除。
如图3所示,所述步骤S4具体包括如下步骤:
S41:对步骤S3中获取的第一锥束CT图像进行频率分割,获取所述第一锥束CT图像的中频成份;所述第一锥束CT图像的中频成份包括条状伪影、图像边界以及图像的纹理特征。
S42:将Radon变换作用于所述中频成份、所述气泡区域fgas和所述骨性组织区域fbone,获取第二虚拟投影;
S43:选取中频成份投影中的气泡成份和骨性组织成份,选取到的成份为条状伪影成份;
S44:重建第二虚拟投影中选取的条状伪影成份,得到fmask2,将步骤S3中获取的第一锥束CT图像减去fmask2,获取条状伪影修正后的第二锥束CT图像,实现锥束CT图像的运动伪影消除。
本实施例还给出该运动伪影修正方法在病人腹部锥束CT数据上的实施结果示意图,如图4所示,其中(a):使用传统滤波反投影重建的结果,(b):对重建结果进行阴影伪影修正的结果,(c):对阴影伪影修正结果进行条状伪影修正的结果,(d):计划CT的图像,显示窗为[0.015 0.025]。
本实施例还给出该运动伪影修正方法在病人腹部锥束CT数据上的实施结果的矢状位图像示意图,如图5所示,其中(a’):使用传统滤波反投影重建的结果,(b’):对重建结果进行阴影伪影修正的结果,(c’):对阴影伪影修正结果进行条状伪影修正的结果,显示窗为[0.015 0.025]。
上述技术方案中,病人上腹部锥束CT数据,使用本发明中的运动伪影修正方法得到的结果,相比较于初始滤波反投影重建得到的结果,显著的消除了原有图像中的阴影伪影和运动伪影,CT数精确度从90HU降低到了20HU,并且较好的保存了原有图像的纹理特征和边界特征。
本实施例所述的一种医学锥束CT图像中运动伪影修正方法,运动伪影主要包括阴影伪影和条状伪影,本方法首先利用病人或者模体锥束CT投影数据重建初始图像,其次利用图像分割技术和图像特征提取技术获取运动伪影存在的主要区域,然后利用线性插值技术消除阴影伪影,再利用频率分割技术消除条状伪影,实现锥束CT图像的运动伪影消除。本发明使用基于信号处理的方法对锥束CT图像的运动伪影进行修正,有效地消除了锥束CT图像中的运动伪影,提高了图像的质量,增加了锥束CT图像在图像引导放射治疗中的图像引导的准确性。
显然,上述实施例仅仅是为清楚地说明所作的举例,而并非对实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保护范围之中。

Claims (6)

1.一种医学锥束CT图像中运动伪影修正方法,其特征在于,包括如下步骤:
S1:初始图像重建,利用病人或者模体锥束CT投影数据重建初始图像f0
S2:区域分割提取,利用图像分割技术和图像特征提取技术对所述初始图像f0进行处理,获取气泡区域fgas、骨性组织区域fbone、阴影伪影与条状伪影存在的主要区域fr
S3:阴影伪影消除,利用线性插值技术对所述气泡区域fgas、所述阴影伪影与条状伪影存在的主要区域fr进行处理,获取阴影伪影修正后的第一锥束CT图像;
所述步骤S3具体包括如下步骤:
S31:将局部Radon变换作用于所述主要区域fr得到第一虚拟投影,即将图像域转换到虚拟投影域;
S32:将Radon变换作用于所述气泡区域fgas,获取含有阴影伪影污染的虚拟投影数据;
S33:在所述第一虚拟投影域中,利用线性插值技术替换含有阴影伪影污染的虚拟投影数据;
S34:重建替换后的虚拟投影数据,对重建后的虚拟投影数据进行低通滤波后得到fmask1,将fmask1加到所述初始图像f0中,获取阴影伪影修正后的第一锥束CT图像,即实现锥束CT图像的阴影伪影修正;
S4:条状伪影消除,利用频率分割技术对阴影伪影修正后的第一锥束CT图像、所述气泡区域fgas、所述骨性组织区域fbone进行处理,获取条状伪影修正后的第二锥束CT图像,实现锥束CT图像的运动伪影消除;
所述步骤S4具体包括如下步骤:
S41:对步骤S3中获取的第一锥束CT图像进行频率分割,获取第一锥束CT图像的中频成份;
S42:将Radon变换作用于所述中频成份、所述气泡区域fgas 和所述骨性组织区域fbone ,获取第二虚拟投影;
S43:选取中频成份投影中的气泡成份和骨性组织成份,选取到的成份为条状伪影成份;
S44:重建第二虚拟投影中选取的条状伪影成份,得到fmask2,将步骤S3中获取的第一锥束CT图像减去fmask2,获取条状伪影修正后的第二锥束CT图像,实现锥束CT图像的运动伪影消除。
2.根据权利要求1所述的医学锥束CT图像中运动伪影修正方法,其特征在于,所述步骤S1中:采用滤波反投影技术重建所述初始图像f0
3.根据权利要求1所述的医学锥束CT图像中运动伪影修正方法,其特征在于,所述步骤S1中:所述初始图像f0中的运动伪影包括阴影伪影和条状伪影。
4.根据权利要求1所述的医学锥束CT图像中运动伪影修正方法,其特征在于,所述步骤S2中:所述图像分割技术采用阈值分割技术,所述图像特征提取技术采用图像腐蚀和图像膨胀技术。
5.根据权利要求1所述的医学锥束CT图像中运动伪影修正方法,其特征在于,所述步骤S34中:采用高斯滤波器对所述投影数据进行低通滤波处理。
6.根据权利要求1所述的医学锥束CT图像中运动伪影修正方法,其特征在于,所述步骤S41中:所述第一锥束CT图像的中频成份包括条状伪影、图像边界以及图像的纹理特征。
CN201710269022.1A 2017-04-24 2017-04-24 一种医学锥束ct图像中运动伪影修正方法 Expired - Fee Related CN107316277B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710269022.1A CN107316277B (zh) 2017-04-24 2017-04-24 一种医学锥束ct图像中运动伪影修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710269022.1A CN107316277B (zh) 2017-04-24 2017-04-24 一种医学锥束ct图像中运动伪影修正方法

Publications (2)

Publication Number Publication Date
CN107316277A CN107316277A (zh) 2017-11-03
CN107316277B true CN107316277B (zh) 2020-12-18

Family

ID=60184773

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710269022.1A Expired - Fee Related CN107316277B (zh) 2017-04-24 2017-04-24 一种医学锥束ct图像中运动伪影修正方法

Country Status (1)

Country Link
CN (1) CN107316277B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113902737A (zh) * 2021-11-23 2022-01-07 西南医科大学附属医院 一种基于甲状腺ct图像异常的检测方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4670840A (en) * 1983-03-09 1987-06-02 Elscint, Inc. Ring artifact correction for computerized tomography
CN101842807A (zh) * 2007-05-31 2010-09-22 伊利克塔股份有限公司 Ct扫描中运动伪影的减少
CN102567958A (zh) * 2010-12-31 2012-07-11 东软飞利浦医疗设备系统有限责任公司 一种去除ct图像中金属伪影的图像后处理方法
CN103366389A (zh) * 2013-04-27 2013-10-23 中国人民解放军北京军区总医院 Ct图像重建方法
CN105264358A (zh) * 2013-02-18 2016-01-20 赛拉诺斯股份有限公司 生物样品的图像分析及测量
CN105310726A (zh) * 2014-06-18 2016-02-10 株式会社东芝 超声波诊断装置、图像处理装置以及图像处理方法
CN105701778A (zh) * 2016-01-11 2016-06-22 赛诺威盛科技(北京)有限公司 一种从ct图像中去除金属伪影的方法
CN106485680A (zh) * 2016-10-13 2017-03-08 上海联影医疗科技有限公司 图像校正方法及装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104135931B (zh) * 2012-02-22 2017-06-20 皇家飞利浦有限公司 用于减少成像数据中的局部伪影的方法和系统
CN106551702B (zh) * 2015-09-30 2021-05-14 通用电气公司 Ct图像的射束硬化伪影校正装置及方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4670840A (en) * 1983-03-09 1987-06-02 Elscint, Inc. Ring artifact correction for computerized tomography
CN101842807A (zh) * 2007-05-31 2010-09-22 伊利克塔股份有限公司 Ct扫描中运动伪影的减少
CN102567958A (zh) * 2010-12-31 2012-07-11 东软飞利浦医疗设备系统有限责任公司 一种去除ct图像中金属伪影的图像后处理方法
CN105264358A (zh) * 2013-02-18 2016-01-20 赛拉诺斯股份有限公司 生物样品的图像分析及测量
CN103366389A (zh) * 2013-04-27 2013-10-23 中国人民解放军北京军区总医院 Ct图像重建方法
CN105310726A (zh) * 2014-06-18 2016-02-10 株式会社东芝 超声波诊断装置、图像处理装置以及图像处理方法
CN105701778A (zh) * 2016-01-11 2016-06-22 赛诺威盛科技(北京)有限公司 一种从ct图像中去除金属伪影的方法
CN106485680A (zh) * 2016-10-13 2017-03-08 上海联影医疗科技有限公司 图像校正方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CT图像环形伪影校正方法研究;黄苏红;《中国优秀硕士学位论文全文数据库(信息科技辑)》;20120115(第01期);第I138-552页 *
Cupping artifact correction and automated classification for high-resolution dedicated breast CT images;Xiaofeng Yang等;《Medical physics》;20121031;第6397-6405页 *

Also Published As

Publication number Publication date
CN107316277A (zh) 2017-11-03

Similar Documents

Publication Publication Date Title
US9934597B2 (en) Metal artifacts reduction in cone beam reconstruction
RU2605519C2 (ru) Двухпроходная коррекция металлического артефакта с компенсацией движения для изображений срезов компьютерной томографии
Meilinger et al. Metal artifact reduction in cone beam computed tomography using forward projected reconstruction information
CN107106102B (zh) 数字减影血管造影
US8855394B2 (en) Methods and apparatus for texture based filter fusion for CBCT system and cone-beam image reconstruction
CN106530236B (zh) 一种医学图像处理方法及系统
US8565502B2 (en) Method and system for reconstruction of tomographic images
CN112017131B (zh) Ct图像金属伪影去除方法、装置及计算机可读存储介质
CN110458762B (zh) 一种基于可调双因子的ct图像射束硬化伪影校正系统
CN111223156B (zh) 一种用于牙科锥束ct系统的金属伪影消除方法
CN109523458A (zh) 一种结合稀疏诱导动态引导滤波的高精度稀疏角度ct重建方法
CN113706409A (zh) 基于人工智能的cbct图像增强方法、装置和存储介质
EP3209207B1 (en) Methods and systems for normalizing contrast across multiple acquisitions
CN110599530B (zh) 基于双正则约束的mvct图像纹理增强方法
CN107316277B (zh) 一种医学锥束ct图像中运动伪影修正方法
Chen et al. Novel method for reducing high-attenuation object artifacts in CT reconstructions
CN107492132B (zh) 一种数字层析成像系统的图像伪影去除方法
JP5317612B2 (ja) 断層像処理装置、x線ct装置およびプログラム
Anhaus et al. Nonlinearly scaled prior image‐controlled frequency split for high‐frequency metal artifact reduction in computed tomography
CN113570586A (zh) 神经网络系统的创建、处理ct图像的方法及其装置
KR102248191B1 (ko) 콘빔 컴퓨터 단층촬영 영상에서 콘빔아티팩트를 보정하는 방법 및 장치, 그것을 포함하는 콘빔 컴퓨터 단층촬영 장치
Tang et al. Efficient metal artifact reduction method based on improved total variation regularization
CN106875342B (zh) 一种计算机断层图像处理方法和装置
TWI593392B (zh) Metal Detection and Artifact Removal Methods in Computerized Tomography Images
Rodríguez-Gallo et al. Inpainting-filtering for metal artifact reduction (IMIF-MAR) in computed tomography

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20201218

Termination date: 20210424