CN103065279B - 圆轨道锥束计算机断层摄影(ct)中用于大幅度地减轻伪影的方法以及*** - Google Patents
圆轨道锥束计算机断层摄影(ct)中用于大幅度地减轻伪影的方法以及*** Download PDFInfo
- Publication number
- CN103065279B CN103065279B CN201210400916.7A CN201210400916A CN103065279B CN 103065279 B CN103065279 B CN 103065279B CN 201210400916 A CN201210400916 A CN 201210400916A CN 103065279 B CN103065279 B CN 103065279B
- Authority
- CN
- China
- Prior art keywords
- mentioned
- image
- circular
- artifact
- circular orbit
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000002591 computed tomography Methods 0.000 title claims abstract description 36
- 238000012937 correction Methods 0.000 claims abstract description 60
- 230000008569 process Effects 0.000 claims description 26
- 230000033228 biological regulation Effects 0.000 claims description 25
- 230000000116 mitigating effect Effects 0.000 claims description 22
- 238000001914 filtration Methods 0.000 claims description 15
- 230000000007 visual effect Effects 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 11
- 238000005259 measurement Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 claims description 2
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims 2
- 239000013589 supplement Substances 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 230000009467 reduction Effects 0.000 description 4
- 230000005855 radiation Effects 0.000 description 3
- VWDWKYIASSYTQR-UHFFFAOYSA-N sodium nitrate Chemical compound [Na+].[O-][N+]([O-])=O VWDWKYIASSYTQR-UHFFFAOYSA-N 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 230000003321 amplification Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 239000002872 contrast media Substances 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000003014 reinforcing effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005266 casting Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 210000004247 hand Anatomy 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000032696 parturition Effects 0.000 description 1
- 229910052704 radon Inorganic materials 0.000 description 1
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical group [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
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/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/416—Exact reconstruction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
本发明的目的在于在圆轨道锥束计算机断层摄影(CT)中大幅度地减轻伪影。在圆轨道计算机断层摄影中校正锥束伪影的方法具备:a)根据规定的圆形射线源轨道中的实际测量到的投影数据来重建参照图像的步骤;b)通过按照与上述圆形射线源轨道对应的理论上的完整的圆形轨道对上述参照图像进行正投影来生成投影数据的步骤;c)根据上述生成的投影数据来重建校正图像的步骤;以及d)通过使用上述参照图像和上述校正图像来生成最终图像从而大幅度地减轻上述锥束伪影的步骤。
Description
技术领域
本发明一般涉及图像处理以及***,更具体而言,涉及圆轨道锥束计算机断层摄影(CT)中伪影的大幅度减轻。
背景技术
锥束伪影是在圆轨道计算机断层摄影中被熟知的问题。由于东芝的AquilionONE等大部分的先进的CT***中的X射线源的锥角非常大,氡区域的数据发生缺损,因此,在CFK图像中产生锥束伪影。
在现有技术的一个方法中,圆轨道锥束CT中的伪影通过应用根据圆形轨道与线轨道的组合等理论上完整的轨道的准确的重建来大幅度地除去。通过追加的线扫描,得到具有用于准确的重建的圆轨道的理论上完整的轨道,但大多数情况下,追加的扫描是收集不能够使用或者非实用的之一的情况。并且,由于不能同时得到圆形数据和线数据,因此,在两次扫描间运动或者造影剂的强化效果(agent enhancement)发生变化时,在两个数据组间数据发生不一致,图像的精度受到影响。最后,通过追加扫描,增加辐射患者的放射线量。
由于上述理由,特别关注只来自圆形数据的图像体的准确的重建。在现有技术的另一方法中,为了推定线数据而使用扫描图(scanogram)。在该方法中患者的放射线量不增加,但即使在推定出的线数据有益于减轻大半的伪影的情况下,依然能够识别锥束伪影。同时,另外,运动或者造影剂的强化效果发生变化时,作为结果而得到的图像稍微变得不准确。
发明内容
目的在于在圆轨道锥束计算机断层摄影(CT)中大幅度地减轻伪 影。
在本实施方式的圆轨道计算机断层摄影中校正锥束伪影的方法具备:
a)根据规定的圆形射线源轨道中的实际测量的投影数据重建参照图像的步骤;
b)通过按照与上述圆形射线源轨道对应的理论上的完整圆形轨道对上述参照图像进行正投影从而生成投影数据的步骤;
c)根据上述生成的投影数据重建校正图像的步骤;以及
d)通过使用上述参照图像和上述校正图像生成最终图像从而大幅度地减轻上述锥束伪影的步骤。
附图说明
图1是表示用于通过本发明在根据在圆轨道上收集到的数据而重建出的图像中大幅度地减轻锥束伪影的多切片X射线CT装置或者扫描仪的一实施方式的图。
图2是表示根据本发明一实施方式来大幅度地减轻锥束(CB)伪影的理由的一侧面的图。
图3是表示与根据本发明来大幅度地减轻锥束(CB)伪影的示例性的处理相关的全部步骤的流程图。
图4是表示用于根据本发明一实施方式来大幅度地减轻锥束伪影的特定的概念性的解决策略的图的集合。
图5是表示与本发明一实施方式中的锥束伪影减轻处理相关的步骤的流程图。
图6A是表示为了评价本发明的示例性的实施方式中的大幅度的锥束伪影减轻的效果而使用的盘模型(ディスクファントム)的侧视图的图。
图6B是表示为了评价本发明的示例性的实施方式中的大幅度的锥束伪影减轻的效果而使用的盘模型的俯视图的图。
图7A是未进行锥束伪影校正而根据盘模型重建的以往的FDK图像。
图7B是未进行锥束伪影校正而根据重复了500次之后的盘模型重建的SART图像。
图7C是通过用于通过本发明来大幅度地减轻锥束伪影的1个示例性的处理来根据盘模型重建的校正后的图像。
图8A是未进行重组(rebinning)而使用了CFK的重建图像。
图8B是未进行重组但使用了进行了OSR滤波处理的CFK的重建图像。
图8C是为了进行比较而使用了HFK的重建图像。
图8D是使用了进行了重组的CFK的重建图像。
图8E是使用了通过本发明的大幅度地减轻锥束伪影的处理来进行了重组和一次正投影的线数据校正的CFK的重建图像。
图8F是使用了通过本发明的大幅度地减轻锥束伪影的处理来进行了重组和三次重复正投影的线数据校正的CFK的重建图像。
具体实施方式
在此,参照附图,相同的符号表示在所有的附图中对应的构造,但特别地参照图1,图中示出包含机架(ガントリ)100和其他的设备或者单元的、本发明的多切片X射线CT装置或者扫描仪的一实施方式。机架100由侧视图示出,还包含X射线管101、环状架102、多列或者二维排列型的X射线检测器103。X射线管101以及X射线检测器103安装在横穿环状架102上的被检体(subject)S的对角线上,环状架102能够围绕旋转轴RA旋转地被支承。旋转单元107使架102以0.4秒/旋转等的高速旋转,被检体S沿着轴RA,向相对于图示纸面的深度方向或者向前方方向(手前方向)移动。
多切片X射线CT装置还包含生成经由滑动环108应该对X射线管101施加的管电压的高电压产生装置109,以使X射线管101生成X射线。朝向被检体S照射X射线,被检体S的剖面区域由圆来表示。X射线检测器103为了检测透过了被检体S的照射X射线,隔着被检体S位于X射线管101的相反侧。
另外,参照图1,X射线CT装置或者扫描仪还包含用于处理由 X射线检测器103检测到的信号的其他的设备。数据收集电路或者数据收集***(DAS)104对于各通道将从X射线检测器103输出的信号转换成电压信号,并将其放大,再转换成数字信号。X射线检测器103以及DAS104构成为,为了最大能够使每一次旋转的规定的全部投影数(TPPR:total numberof projections per rotation)为900TPPR、从900TPPR到1800TPPR之间、以及从900TPPR到3600TPPR之间而进行处理。
上述数据通过非接触数据传送装置105,发送至机架100的外部的收容于控制台的前处理设备106。前处理设备106对于原始数据执行灵敏度校正等特定的校正。接着,存储设备112在重建处理即将开始之前的阶段存储作为其结果而生成的数据(也被称为投影数据的数据)。存储设备112经由数据/控制总线,与重建设备114、输入设备115、显示设备116、多尺度(Multiscale)处理设备117、以及扫描计划辅助装置200一起,与***控制器110连接。扫描计划辅助装置200包含用于对摄像技师制定扫描计划进行辅助的功能。
重建设备114的一实施方式还包含各种软件构成要素和硬件构成要素,对于投影数据执行规定的解析重建处理。根据本发明的一方式,CT装置的重建设备114有利的是使用规定的滤波校正反投影(FBP)法来重建图像体。
根据本发明的另一方式,CT装置的重建设备114有利的是使用迭代重建法来将总变分(TV)最小化。一般而言,本发明的一实施方式的重建设备114实施总变分迭代重建(TVIR)算法。TVIR对于投影数据,执行有序子集同时代数重建法(OS-SART)步骤等的同时代数重建法和TV最小化步骤等的调整。在一实施方式中,这两个步骤在预先设定了重复次数的主循环中连续实施。
在TV最小化步骤之前,投影数据被提供给有序子集同时代数重建法(OS-SART)。投影数据被分组为分别具有特定数的视角(ビュ一)的规定数的子集N。在一实施方式中,在正在实施有序子集同时代数重建法(OS-SART)中,也可以连续地处理各子集。在另一实施方式中,多个子集也可以利用多个中央处理装置(CPU)或者一个图 形处理装置(GPU)等的特定微处理器来并列进行处理。在总变分(TV)最小化步骤中,在重建设备114的一实施方式中,为了以当前的图像体的目标函数比之前的图像体的目标函数小的方式探索正的步长(stepsize)而使用直线探索方式(line search strategy)。
重建设备114还在实施有序子集同时代数重建法(OS-SART)中,执行两个主要的运算。即,重建设备114对于各子集N,为了生成计算用的投影数据而对图像体进行重投影,为了重建更新图像体,对测量到的投影数据和计算上的投影数据的归一化后的差分进行反投影。更具体而言,在重建设备114的一实施方式中,通过使用***矩阵的系数未被高速缓存的射线追踪法,来对图像体进行重投影。另外,在重建设备114的一实施方式中,同时对子集内的所有的射线进行重投影。该处理任意选择,并列地实施。在反投影中,在重建设备114的一实施方式中,以生成所希望的更新图像体为目的,为了对子集内的所有的归一化后的差分投影数据进行反投影而使用像素驱动(pixel-driven)法。由于重建设备114为了生成图像体而对子集内的所有的射线总数(即差分投影数据)进行反投影,因此,该运算也任意选择,并列地实施。为了结束单一的OS-SART步骤,对于所有的子集N应用这些运算。另外,AWAD任意选择进行组合。
本发明的一实施方式除了上述构成要素之外,还具备用于执行锥束伪影减轻的各种其他的软件模块和硬件构成要素。根据本发明的一方式,CT装置的锥束(CB)伪影减轻设备117有利的是在特定的状况下执行用于大幅度地减轻锥束伪影的锥束伪影减轻功能。一般而言,CB伪影具有使画质劣化的、包含阴影(shading)锥束和高对比度锥束的两个分量。在1个示例的处理中,为了提高画质,通过对重组进行滤波处理来校正阴影(shading),高对比度由线图像来校正。针对上述的伪影,参照本申请的其他的图进一步进行说明。
在本发明的另一实施方式中,CT装置的锥束(CB)伪影减轻设备117有利的是为了大幅度地减轻锥束伪影而将准确的重建的要素与SART等的迭代重复重建的要素结合。如更详细地说明的那样,重建设备114根据使用圆形射线源轨道上的锥束射线源而收集到的测量投影数据来重建圆形图像。之后,CB伪影减轻设备117对来自该圆形图像的线数据进行正投影,接着,重建设备114根据该正投影后的线数据来重建线图像。根据上述的步骤,CB伪影减轻设备117合成圆形图像和线图像。该合成图像包含对于应该在以往方法的重建图像中存在的CB伪影而言大幅度地减轻了CB伪影的CB伪影。最后,CB伪影减轻设备117输出校正后的图像。
在本发明的一实施方式中,锥束伪影减轻设备117与其他的软件模块和/或存储设备112、重建设备114、显示设备116、以及输入设备115等***构成要素经由数据/控制总线动作性地连接。关于该点,在本发明的另一实施方式中,锥束伪影减轻设备117不必独立地执行锥束伪影功能和/或与其相关的任务。并且,在本发明的代替实施方式中,锥束伪影减轻设备117任意选择,是重建设备114等其他的设备的一部分。
图2是表示通过本发明一实施方式来大幅度地减轻锥束(CB)伪影的理由的一个侧面的图。该图表示一系列的细小的对象0位于相对于存在于规定的圆轨道CT上的射线源而言相对的位置的示例性的状况。位于射线源位置S1的锥束以规定的锥束角朝向对象0进行照射。该示例的状况还示出了图像为了进行诊断而具有良好的分辨率的被放大或者在所希望的摄像视野(FOV)内被重建的情况很多。在正投影中,由于需要与使X射线束衰减的对象整体相关的信息,因此,如图2所示的那样生成两个体图像CFK_A以及CFK_B。阴影区域与非充分地收集到的实际测量到的投影数据对应。图像CFK_A具有最终结果所应该使用的所希望的FOV,图像CFK_B为了生成缺损的线数据而具有全(full)FOV。所需的范围2H沿着用于生成该线数据的线数据方向LTD而示出。
接着,参照图3,流程图表示与根据本发明来大幅度地减轻锥束(CB)伪影的示例性的处理相关的全部步骤。实际上,该流程图是通过本发明来大幅度地减轻锥束(CB)伪影的示例性的处理的概念性的方案,本发明未必限定于以下的说明。在步骤S100中,收集实际测量到的投影数据。一般而言,收集到的实际测量到的投影数据由 于锥束圆形扫描几何图形由于某些原因不足,因此,当通过以往的方法进行重建时,易于产生锥束伪影。在步骤S110中,圆形图像根据使用圆形射线源轨道上的锥束射线源而收集到的实际测量到的投影数据来重建。重建后的圆形图像任意选择,为了之后取出而存储于规定的存储区域。在步骤S120中,根据在步骤S110中重建的圆形图像,将线数据进行正投影或者重投影。针对步骤S120中的线数据生成的细节,之后参照其他的流程图进行说明。在步骤S130中,线图像根据在步骤S120中生成的被正投影后的线数据来重建。根据上述的步骤S110、S120、以及S130,在步骤S140中,在此合成圆形图像与线图像。该合成图像包含对于应该存在于以往方法的重建图像中的CB伪影而言大幅度地减轻了的CB伪影。最后,在步骤S150中,为了对校正后的图像进行显示或者解析而进行输出。在大幅度地减轻CB伪影的另一实施方式中,为了提高线数据的精度,并进一步减轻输出图像的伪影,使用SART等已知的迭代的技术反复地重复上述的步骤中的几个。
图4是表示用于通过本发明一实施方式来大幅度地减轻锥束伪影的特定的概念性的解决策略的图的集合。如关于图2已经说明了的那样,在示例的状况中,为了进行诊断,需要具有分辨率的放大或者所希望的摄像视野(FOV)。在正投影中,由于需要与使X射线束衰减的对象整体相关的信息,因此,图4A以及图4B所示的那样的两个图像CFK_A以及CFK_B如图4C所示的那样根据使用圆轨道CT上的锥束射线源CBS而收集到的对应的圆形数据使用圆形费尔德坎普(CFK:circular Feld-Kamp)法来重建。图4C还示出用于图像体CFK_A的圆形数据是用于图像CFK_B的圆形数据的放大后的一部分的情况。也就是说,图像CFK_A在图像CFK_B的全FOV或者大的FOV中具有所希望的FOV,该全FOV沿着Z轴向延伸。
根据一实施方式,图4D所示的线数据根据进行了本发明的滤波校正反投影的体图像CFK_B的正投影而生成。关于该点,在本申请中,重投影在上述线数据生成中与正投影同义地使用。如上述那样,图像CFK_B具有全FOV。由于具有大的锥角的X射线能够超过重建 图像透过Z方向的空间,因此,图像CFK_B任意选择,超过由扫描***的锥角确定的某个Z范围而延伸。在另一实施方式中,为了减轻图像CFK_B中的锥束伪影,在沿着线轨道LT的正投影之前,对图像CFK_B应用规定的自适应低通3D滤波器。或者,在另一实施方式中,在沿着线轨道LT的正投影之前,对图像CFK_B应用规定的因数分解方法。
在本发明的锥束伪影减轻处理的一实施方式中,在通过图像CFK_B的重投影从线轨道得到了线数据之后,具有如图4E所示的所希望的FOV的线图像根据图4D所示的那样的线数据来重建。最后,具有如图4E所示的所希望的FOV的线图像如由+记号表示的那样与如图4A所示的图像体CFK_A合成,生成如图4F所示的校正后的图像。该合成图像具有大幅地减轻了的锥束伪影。在一实施方式中,图像体CFK_B任意地选择,通过追加根据具有全FOV的线数据重建后的图像来更新。并且,由于根据具有伪影的体图像进行正投影,因此,该线数据是近似的数据。因此,在另一实施方式中,为了改良线数据,并且提高锥束伪影的减轻,使用重复方法。
图5是表示与本发明一实施方式中的锥束伪影减轻处理相关的步骤的流程图。在大幅度地减轻锥束伪影的示例性的处理中,假设实际测量到的投影数据使用特定的锥束角通过规定的圆形数据收集技术来收集。一般而言,在实施方式中的锥束伪影减轻处理中,生成分别具有所希望的小的摄像视野和来自对应的圆轨道锥束数据的全摄像视野的两个参照图像CFK_A和CFK_B。在实施方式中的锥束伪影减轻处理中,通过沿着规定的射线源轨道对具有全摄像视野的参照图像CFK_B进行正投影来生成投影数据,接着,根据该投影数据重建校正图像。该校正通过线轨道等规定的射线源轨道补充参照图像CFK_B的圆形射线源轨道来实现。并且,在实施方式中的锥束伪影减轻处理中,为了使步骤S360之后的校正后的图像具有通过本发明大幅度地减轻的锥束伪影,关于线图像反复地重复从步骤S230到S290等的特定步骤。在代替实施方式中,在锥束伪影减轻处理中,代替重复而以规定的方法执行线图像的操作。
另外,参照图5,针对大幅度减轻锥束伪影的处理进一步详细地进行说明。一般而言,通过沿着由Cseg+z/cos(gamma)提供的规定的滤波处理方向进行滤波处理从而根据滤波校正反投影算法来得到参照图像,在此,gamma是锥角,z是来自Cseg的垂直方向的距离,Cseg通过用检测器列数-1后除以2来规定。在步骤S210中,实际测量到的投影数据被提供给规定的Hconv步骤,在该步骤中,为了重建两个参照图像而输出卷积(convolution)数据。在一实施方式中,在Hconv步骤中,利用Ramp核心数据(カ一ネル)+Hilbert核心数据的混合核心数据。在另一实施方式中,为了生成两个参照图像,使用圆形费尔德坎普(CFK)法。并且,在用于头部成像的脑部的阴影等中,为了在某种程度上提高画质,以任意选择的方式对Hconv步骤追加重组步骤。在代替实施方式中,卷积数据的重组步骤以及逆重组(inverserebinning)分别为了在某种程度上提高画质而在Hconv步骤之前和之后以任意选择的方式执行。在步骤S300中,第1参照图像CFK A使用圆形费尔德坎普(CFK)法根据圆轨道锥束数据来重建,第1参照体图像CFK_A具有第1摄像视野(FFOV)。FFOV一般而言是所希望的反投影视野(back projection field of view)。同样地,在步骤S220中,第2参照图像CFK_B根据圆轨道锥束数据来重建,第2参照体图像CFK_B具有第2摄像视野(SFOV),SFOV比FFOV大,以任意选择的方式将CT图像***的机架作为对象。接着,双方的参照图像CFK_A以及CFK_B以任意选择的方式,为了之后取出分别在步骤S300以及S220中进行存储。并且,在步骤S220中,为了与重复相关的步骤之后的示例对重复计数Iter进行初始化处理。
另外,参照图5,针对大幅度地减轻锥束伪影的处理的重复步骤进一步详细地进行说明。在步骤S230中,为了追踪规定数的重复,重复计数Iter仅增加1。在步骤S240中,对第2参照体图像CFK_B应用OSR滤波等规定的滤波,生成图像CFK_C,但该图像CFK_C是被LLFDK校正后的CFK图像。之后,接着,在步骤S250中,图像CFK_C为了进行正投影而在扩张切片CFK_D上进行体扩张(volume expanded)。在步骤S260中,切片CFK_D正投影为线数据, 在步骤S310以及步骤S270中,相同的线数据分别重建为线图像A以及线图像B。之后,线图像B在步骤S270中每当重复就根据线图像换算系数(scaling factor)进行换算(scaling),在步骤S280中与第2参照图像CFK B合成,在步骤S290中生成了LC图像B之后,在步骤S230中使用LC图像B开始下一重复,在该步骤中重复计数Iter增加。
同样地,线图像A也在S310中根据线图像换算系数进行换算,在步骤S320中与第1参照图像CFK_A合成,在步骤S330中生成LC图像A。在步骤S340中,判断重复计数Iter的值是否比规定的重复次数Niter大。在步骤S340中判断为Iter计数的值不比最大重复值Niter大时,规定的重复次数还没有结束,大幅度地减轻锥束伪影的处理为了进一步重复而进入步骤S260。另一方面,在步骤S340中判断为Iter计数值比最大重复值Niter大时,规定的重复次数结束,大幅度地减轻锥束伪影的处理任意选择,在步骤S350中应用了LLFDK校正之后,在步骤S360中生成校正后的最终图像。根据本发明,大幅度地减轻锥束伪影的处理并不限定于上述的说明,还包含其他的实施步骤。
图6A以及图6B分别表示为了评价本发明的示例的实施方式中的大幅度的锥束伪影减轻的效果而使用的盘模型的侧视图和俯视图。该模型(ファントム)沿着Z方向具有盘状的裂口。
图7A、图7B以及图7C表示本发明的一比较例中的大幅度的锥束伪影减轻的效果。图7A表示未进行锥束伪影校正而根据盘模型重建的以往的FDK图像。能够识别出明显的锥束伪影。图7B表示根据没有进行锥束伪影校正的500次重复之后的盘模型重建的SART图像。图7C表示根据进行了用于根据本发明来大幅度地减轻锥束伪影的1个示例的处理的盘模型重建的校正后的图像。大幅度地减轻了明显的锥束伪影。
图8A至图8F表示本发明的另一比较例中的大幅度的锥束伪影减轻的效果。图8A至图8F中的示例的图像是根据临床上得到的头部数据的厚的MPR切片使用CFK重建法或者HFK重建法而生成的 图像。并且,图8A至图8F中的示例性的图像使用线校正、数据的重组、以及OSR校正的组合来生成。另外,图8A表示未进行重组而使用了CFK的重建图像。图8B表示未进行重组但使用了进行了OSR滤波处理的CFK的重建图像。图8C表示为了进行比较而使用了HFK的重建图像。图8D表示使用了进行了重组的CFK的重建图像。图8E表示使用了通过本发明的大幅度地减轻锥束伪影的处理来进行了重组和一次正投影的线数据校正的CFK的重建图像。图8F表示使用了通过本发明的大幅度地减轻锥束伪影的处理来进行了重组和三次重复正投影的线数据校正的CFK的重建图像。
但是,本发明的多个特征以及优点与本发明的构造以及功能的细节一起在上述的说明进行了叙述,但本公开只不过是个示例,另外,关于细节,特别是关于部件的形状、大小、以及结构和软件、硬件、或者其双方的组合中的安装方式,能够进行变更,但要理解为该变更包含于表现附加的权利要求的范围的用语的广泛的一般意思所最大限定的本发明的原理内。
Claims (26)
1.一种在圆轨道计算机断层摄影中校正锥束伪影的方法,具备:
a)根据规定的圆形射线源轨道中的实际测量的投影数据来重建圆形图像的步骤;
b)通过沿着线轨道对上述圆形图像进行正投影从而生成线数据的步骤;
c)根据所生成的上述线数据来重建包含伪影的线图像的步骤;以及
d)使用上述圆形图像和上述线图像来生成校正图像的步骤。
2.根据权利要求1所述的圆轨道计算机断层摄影中校正锥束伪影的方法,其中,
上述圆形图像通过根据滤波校正反投影算法沿着规定的滤波处理方向进行滤波处理来得到。
3.根据权利要求2所述的圆轨道计算机断层摄影中校正锥束伪影的方法,其中,
上述规定的滤波处理方向由Cseg+z/cos(gamma)来提供,在此,gamma是锥角,z是自Cseg起的垂直方向的距离,Cseg通过用检测器列的数量-1后除以2来规定。
4.根据权利要求2所述的圆轨道计算机断层摄影中校正锥束伪影的方法,其中,
上述圆形图像是圆形费尔德坎普图像与混合卷积重建图像中之一。
5.根据权利要求1所述的圆轨道计算机断层摄影中校正锥束伪影的方法,其中,
上述校正图像是线扫描重建图像。
6.根据权利要求1所述的圆轨道计算机断层摄影中校正锥束伪影的方法,其中,
还具备追加的步骤f),使用校正后的上述图像作为上述步骤b)中的上述圆形图像,并将上述步骤b)到d)重复规定的次数。
7.根据权利要求1所述的圆轨道计算机断层摄影中校正锥束伪影的方法,其中,
在上述步骤b)中应用低通滤波器。
8.根据权利要求1所述的圆轨道计算机断层摄影中校正锥束伪影的方法,其中,
在上述步骤a)中,根据上述测量圆形数据生成第1圆形图像和第2圆形图像,上述第1圆形图像具有全摄像视野,上述第2圆形图像是所希望的摄像视野。
9.根据权利要求8所述的圆轨道计算机断层摄影中校正锥束伪影的方法,其中,
在上述步骤b)中,对上述第1圆形图像应用规定的OSR滤波。
10.根据权利要求9所述的圆轨道计算机断层摄影中校正锥束伪影的方法,其中,
在上述步骤b)中,为了生成被正投影后的线数据而对上述第1圆形图像进行正投影。
11.根据权利要求10所述的圆轨道计算机断层摄影中校正锥束伪影的方法,其中,
在上述步骤c)中,根据上述被正投影后的线数据来重建线图像。
12.根据权利要求11所述的圆轨道计算机断层摄影中校正锥束伪影的方法,其中,
在上述步骤d)中,通过使用上述第2圆形图像和上述线图像生成校正后的上述图像从而大幅度地减轻上述锥束伪影。
13.根据权利要求12所述的圆轨道计算机断层摄影中校正锥束伪影的方法,其中,
还具备追加的步骤f),使用校正后的上述图像作为上述步骤b)中的上述圆形图像,并将上述步骤b)到d)重复规定的次数。
14.一种在圆轨道计算机断层摄影中校正锥束伪影的***,具备:
重建设备,用于根据规定的圆形射线源轨道中的实际测量到的投影数据重建圆形图像;以及
锥束伪影减轻设备,用于通过沿着线轨道对上述圆形图像进行正投影来生成线数据,并与上述重建设备连接,
上述重建设备根据所生成的上述线数据来重建包含伪影的线图像,
上述锥束伪影减轻设备使用上述圆形图像和上述线图像来生成校正图像。
15.根据权利要求14所述的圆轨道计算机断层摄影中校正锥束伪影的***,其中,
上述圆形图像通过基于滤波校正反投影算法沿着规定的滤波处理方向进行滤波处理来得到。
16.根据权利要求15所述的圆轨道计算机断层摄影中校正锥束伪影的***,其中,
上述规定的滤波处理方向由Cseg+z/cos(gamma)来提供,在此,gamma是锥角,z是自Cseg起的垂直方向的距离,Cseg是通过用检测器列的数量-1后除以2来规定。
17.根据权利要求15所述的圆轨道计算机断层摄影中校正锥束伪影的***,其中,
上述圆形图像是圆形费尔德坎普图像与混合卷积重建图像中的一个图像。
18.根据权利要求14所述的圆轨道计算机断层摄影中校正锥束伪影的***,其中,
上述校正图像是线扫描重建图像。
19.根据权利要求14所述的圆轨道计算机断层摄影中校正锥束伪影的***,其中,
上述重建设备以及上述锥束伪影减轻设备使用校正后的上述图像作为上述圆形图像,并重复规定的重复次数。
20.根据权利要求14所述的圆轨道计算机断层摄影中校正锥束伪影的***,其中,
上述重建设备利用低通滤波器。
21.根据权利要求14所述的圆轨道计算机断层摄影中校正锥束伪影的***,其中,
上述重建设备根据上述测量圆形数据生成第1圆形图像和第2圆形图像,上述第1圆形图像具有全摄像视野,上述第2圆形图像是所希望的摄像视野。
22.根据权利要求21所述的圆轨道计算机断层摄影中校正锥束伪影的***,其中,
上述锥束伪影减轻设备对上述第1圆形图像应用规定的OSR滤波。
23.根据权利要求22所述的圆轨道计算机断层摄影中校正锥束伪影的***,其中,
上述锥束伪影减轻设备为了生成被正投影后的线数据而对上述第1圆形图像进行正投影。
24.根据权利要求23所述的圆轨道计算机断层摄影中校正锥束伪影的***,其中,
上述重建设备根据上述被正投影后的线数据重建线图像。
25.根据权利要求24所述的圆轨道计算机断层摄影中校正锥束伪影的***,其中,
上述锥束伪影减轻设备通过使用上述第2圆形图像与上述线图像生成校正后的上述图像从而大幅度地减轻上述锥束伪影。
26.根据权利要求25所述的圆轨道计算机断层摄影中校正锥束伪影的***,其中,
上述重建设备以及上述锥束伪影减轻设备使用校正后的上述图像作为上述圆形图像,并重复规定的重复次数。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/276,841 US8731269B2 (en) | 2011-10-19 | 2011-10-19 | Method and system for substantially reducing artifacts in circular cone beam computer tomography (CT) |
US13/276,841 | 2011-10-19 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103065279A CN103065279A (zh) | 2013-04-24 |
CN103065279B true CN103065279B (zh) | 2017-04-12 |
Family
ID=48107896
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210400916.7A Active CN103065279B (zh) | 2011-10-19 | 2012-10-19 | 圆轨道锥束计算机断层摄影(ct)中用于大幅度地减轻伪影的方法以及*** |
Country Status (3)
Country | Link |
---|---|
US (1) | US8731269B2 (zh) |
JP (1) | JP6109524B2 (zh) |
CN (1) | CN103065279B (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6158910B2 (ja) * | 2012-03-29 | 2017-07-05 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 正則化による反復画像再構成 |
CN104182932B (zh) * | 2013-05-27 | 2017-04-12 | 株式会社日立制作所 | Ct 装置、ct 图像***及ct 图像生成方法 |
CN106255994B (zh) * | 2014-02-18 | 2019-09-17 | 皇家飞利浦有限公司 | 针对正电子发射断层摄影(pet)列表模式迭代重建的重建中滤波 |
CN104952043B (zh) * | 2014-03-27 | 2017-10-24 | 株式会社日立制作所 | 图像滤波方法及ct *** |
US9704223B2 (en) * | 2014-06-10 | 2017-07-11 | Toshiba Medical Systems Corporation | Method and system for substantially reducing cone beam artifacts based upon adaptive scaling factor in circular computer tomography (CT) |
US9558569B2 (en) * | 2014-08-12 | 2017-01-31 | Toshiba Medical Systems Corporation | Method and system for substantially reducing cone beam artifacts based upon image domain differentiation in circular computer tomography (CT) |
CN104574292B (zh) | 2014-11-26 | 2018-06-26 | 沈阳东软医疗***有限公司 | 一种ct图像的校正方法和装置 |
US10019818B2 (en) * | 2015-03-12 | 2018-07-10 | Toshiba Medical Systems Corporation | Concurrent forward- and back-projection iterative reconstruction system |
US10210633B2 (en) * | 2015-06-12 | 2019-02-19 | Hitachi, Ltd. | X-ray CT device and sequential correction parameter determination method |
DE102015007934A1 (de) * | 2015-06-19 | 2016-12-22 | Universität Stuttgart | Verfahren und Computerprogrammprodukt zum Erzeugen eines artefaktreduzierten Voxeldatensatzes |
CN105118039B (zh) * | 2015-09-17 | 2018-04-10 | 广州华端科技有限公司 | 实现锥束ct图像重建的方法及*** |
WO2017129459A1 (en) * | 2016-01-29 | 2017-08-03 | Koninklijke Philips N.V. | Cone beam computed tomography projection values providing system and method |
CN105608720B (zh) * | 2016-01-30 | 2018-01-16 | 上海联影医疗科技有限公司 | 计算机断层成像伪影校正方法及装置 |
WO2017128892A1 (zh) | 2016-01-30 | 2017-08-03 | 上海联影医疗科技有限公司 | 计算机断层成像伪影校正方法及*** |
US10977840B2 (en) * | 2018-08-21 | 2021-04-13 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for iterative reconstruction |
US11026642B2 (en) * | 2019-03-29 | 2021-06-08 | Canon Medical Systems Corporation | Apparatuses and a method for artifact reduction in medical images using a neural network |
CN110766629B (zh) * | 2019-10-17 | 2022-03-01 | 广州华端科技有限公司 | Cbct***几何校正方法、装置、计算机设备和存储介质 |
EP4012648A1 (en) * | 2020-12-10 | 2022-06-15 | Koninklijke Philips N.V. | Reduction of artefacts in medical images |
CN116580115B (zh) * | 2023-03-29 | 2023-12-05 | 赛诺威盛医疗科技(扬州)有限公司 | 螺旋ct图像迭代重建方法、装置及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1541374A (zh) * | 2001-08-10 | 2004-10-27 | �ʼҷ����ֵ��ӹɷ�����˾ | 用于根据投影图像重建三维数据集的x-射线检查装置 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5841829A (en) * | 1997-05-13 | 1998-11-24 | Analogic Corporation | Optimal channel filter for CT system with wobbling focal spot |
JP4282110B2 (ja) * | 1998-06-04 | 2009-06-17 | 株式会社東芝 | X線コンピュータ断層撮影装置 |
US6560308B1 (en) * | 2001-10-26 | 2003-05-06 | Kabushiki Kaisha Toshiba | Method and system for approximating missing data in cone beam x-ray CT reconstruction |
ATE354838T1 (de) * | 2003-02-14 | 2007-03-15 | Koninkl Philips Electronics Nv | System und verfahren für exakte bildrekonstruktion für spiral-kegelstrahl- computertomographie mit redundanten daten |
WO2006116316A2 (en) * | 2005-04-22 | 2006-11-02 | University Of Chicago | Open source trajectory method and apparatus for interior imaging |
WO2007034342A2 (en) * | 2005-09-26 | 2007-03-29 | Koninklijke Philips Electronics, N.V. | Iterative reconstruction with enhanced noise control filtering |
AU2007217369A1 (en) * | 2006-02-27 | 2007-08-30 | University Of Rochester | Method and apparatus for cone beam CT dynamic imaging |
US8582855B2 (en) * | 2007-01-04 | 2013-11-12 | Koninklijke Philips N.V. | Apparatus, method and computer program for producing a corrected image of a region of interest from acquired projection data |
JP2008237886A (ja) * | 2007-02-28 | 2008-10-09 | Toshiba Corp | X線ct装置及びその制御方法 |
JP2009006133A (ja) * | 2007-05-31 | 2009-01-15 | Toshiba Corp | X線ct装置及びその制御方法 |
JP5525706B2 (ja) * | 2007-07-25 | 2014-06-18 | 株式会社東芝 | X線コンピュータ断層撮影装置 |
WO2009115982A1 (en) * | 2008-03-21 | 2009-09-24 | Koninklijke Philips Electronics N.V. | Computed tomography scanner apparatus and method for ct-based image acquisition based on spatially distributed x-ray microsources of the cone-beam type |
US8194961B2 (en) * | 2008-04-21 | 2012-06-05 | Kabushiki Kaisha Toshiba | Method, apparatus, and computer-readable medium for pre-reconstruction decomposition and calibration in dual energy computed tomography |
US8577113B2 (en) * | 2009-10-22 | 2013-11-05 | Koninklijke Philips N.V. | Image reconstruction for truncated detector array |
US8559687B2 (en) * | 2009-10-30 | 2013-10-15 | Kabushiki Kaisha Toshiba | Method for circular scan reconstruction in computed tomography and computed tomographic device |
JP5406063B2 (ja) * | 2010-01-27 | 2014-02-05 | 株式会社日立メディコ | 再構成演算装置、再構成演算方法、及びx線ct装置 |
US8031828B1 (en) * | 2010-04-30 | 2011-10-04 | General Electric Company | Method and apparatus for computed tomography |
US8805037B2 (en) * | 2011-05-31 | 2014-08-12 | General Electric Company | Method and system for reconstruction of tomographic images |
-
2011
- 2011-10-19 US US13/276,841 patent/US8731269B2/en active Active
-
2012
- 2012-10-19 JP JP2012232359A patent/JP6109524B2/ja active Active
- 2012-10-19 CN CN201210400916.7A patent/CN103065279B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1541374A (zh) * | 2001-08-10 | 2004-10-27 | �ʼҷ����ֵ��ӹɷ�����˾ | 用于根据投影图像重建三维数据集的x-射线检查装置 |
Also Published As
Publication number | Publication date |
---|---|
JP2013085965A (ja) | 2013-05-13 |
CN103065279A (zh) | 2013-04-24 |
JP6109524B2 (ja) | 2017-04-05 |
US8731269B2 (en) | 2014-05-20 |
US20130101192A1 (en) | 2013-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103065279B (zh) | 圆轨道锥束计算机断层摄影(ct)中用于大幅度地减轻伪影的方法以及*** | |
CN105118039B (zh) | 实现锥束ct图像重建的方法及*** | |
CN102270350B (zh) | 结合四维噪声滤波器的迭代ct图像重建 | |
CN106683144A (zh) | 一种图像迭代重建方法及装置 | |
CN109949411B (zh) | 一种基于三维加权滤波反投影和统计迭代的图像重建方法 | |
EP3673457B1 (en) | A method of generating an enhanced tomographic image of an object | |
CN103054600B (zh) | X射线ct装置、医用图像处理装置及医用图像处理方法 | |
CN103180879B (zh) | 用于从投影数据对对象进行混合重建的设备和方法 | |
CN1980604A (zh) | 心脏锥面光束ct重建中条纹伪影的减少 | |
JP6505513B2 (ja) | X線コンピュータ断層撮像装置及び医用画像処理装置 | |
CN102270349A (zh) | 在没有调节项的情况下对ct 图像的迭代重建 | |
CN103156629A (zh) | 图像处理设备和图像处理方法 | |
CN103054599B (zh) | X射线ct装置及x射线ct装置的动作方法 | |
EP3847623B1 (en) | A method of generating an enhanced tomographic image of an object | |
US8798350B2 (en) | Method and system for reconstruction algorithm in cone beam CT with differentiation in one direction on detector | |
US20180286087A1 (en) | Volume image reconstruction using projection decomposition | |
JP6615531B2 (ja) | X線コンピュータ断層撮影装置及び医用画像処理装置 | |
JP4222930B2 (ja) | 3次元逆投影方法および装置並びにx線ct装置 | |
CN107251095A (zh) | 图像重建***、方法和计算机程序 | |
Jakab et al. | High quality cone-beam CT reconstruction on the GPU | |
Sun | Rigid motion correction for head CT imaging | |
Mahmoudi et al. | Recent Advances in X-Ray CT Image Reconstruction Techniques | |
WO2024008764A1 (en) | Cone beam artifact reduction | |
WO2023227511A1 (en) | Simulating x-ray from low dose ct | |
CN117830187A (zh) | 医学图像处理方法、装置以及医疗设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C41 | Transfer of patent application or patent right or utility model | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20160719 Address after: Tochigi County, Japan Applicant after: Toshiba Medical System Co., Ltd. Address before: Tokyo, Japan, Japan Applicant before: Toshiba Corp Applicant before: Toshiba Medical System Co., Ltd. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |