CN102648857A - 在多源ct中用于散射束校正的方法和计算机*** - Google Patents

在多源ct中用于散射束校正的方法和计算机*** Download PDF

Info

Publication number
CN102648857A
CN102648857A CN2012100282475A CN201210028247A CN102648857A CN 102648857 A CN102648857 A CN 102648857A CN 2012100282475 A CN2012100282475 A CN 2012100282475A CN 201210028247 A CN201210028247 A CN 201210028247A CN 102648857 A CN102648857 A CN 102648857A
Authority
CN
China
Prior art keywords
scattering
detector
projection data
require
data set
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
Application number
CN2012100282475A
Other languages
English (en)
Other versions
CN102648857B (zh
Inventor
T.弗洛尔
M.彼得希尔卡
K.斯蒂尔斯托佛
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.)
Siemens Healthineers AG
Original Assignee
Siemens AG
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 Siemens AG filed Critical Siemens AG
Publication of CN102648857A publication Critical patent/CN102648857A/zh
Application granted granted Critical
Publication of CN102648857B publication Critical patent/CN102648857B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G01N23/046Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/40Arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4007Arrangements for generating radiation specially adapted for radiation diagnosis characterised by using a plurality of source units
    • A61B6/4014Arrangements for generating radiation specially adapted for radiation diagnosis characterised by using a plurality of source units arranged in multiple source-detector units
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • A61B6/5282Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to scatter
    • 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
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/06Diaphragms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/42Arrangements for detecting radiation specially adapted for radiation diagnosis
    • A61B6/4291Arrangements for detecting radiation specially adapted for radiation diagnosis the detector being combined with a grid or grating
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/40Imaging
    • G01N2223/419Imaging computed tomograph
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/60Specific applications or type of materials
    • G01N2223/612Specific applications or type of materials biological material
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/424Iterative

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Optics & Photonics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Pulmonology (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明涉及用于在多源CT中对对象进行CT检查时的散射束校正的一种方法和一种计算机***(10),具有如下的方法步骤:-产生原始投影数据组,-利用至少一个探测器(3,5)的原始投影数据组重建对象(P),-确定由每个辐射器(2,4)所产生的、仅在该至少一个其它辐射器的朝向其相对布置的探测器(3,5)的原始射线(O)的方向上的散射,-通过从所述原始投影数据组中去除所计算出的散射,产生校正后的投影数据组,-利用所述校正后的投影数据组重建对象(P),并且-如果符合至少一个预先给定的中断标准,则在确定散射的情况下执行该方法的进一步迭代或者输出重建结果。

Description

在多源CT中用于散射束校正的方法和计算机***
技术领域
本发明涉及用于在多源CT中对对象进行CT检查时的散射束校正的一种方法和一种计算机***,其中散射束校正在投影数据上迭代地执行。
背景技术
在双源CT中通过同时使用两个辐射器-探测器***,X射线量子在待测量的对象处从第一辐射器-探测器***的X射线辐射器散射到在同一机架上角度错开地布置的第二辐射器-探测器***的探测器,并且反之亦然。为了避免伪影必须通过合适的方法校正这些所谓的横向散射。为此,需要识别在每个投影中的散射束强度的轮廓(Profil)。
为此以下三种方法基本上是公知的:
(i)根据关于校正时刻时间上在前的数据,从正弦图数据(Sinugrammdaten)中基于模型地估计散射。
(ii)借助专用传感器测量散射。
(iii)交替消隐相应的原始辐射并且在消隐阶段直接地测量横向散射。
所有三个方案的共同点在于,在以原始数据计算CT图像之前对原始数据执行一次性的散射束校正。基于模型地估计以及利用专用传感器测量的方法成功地应用在可商业化使用的双源***中。
基于模型地估计的方法(i)的严格前提一方面在于以顺序模式进行扫描或者在沿着进给方向(=z方向)螺旋形地扫描时不改变或仅不太多地改变对象。但随着探测器的z覆盖的增加,在同样高的值的情况下,对于在螺旋形地扫描时的节距来说,对散射束轮廓的估计逐渐变得不正确。另一方面,例如在文献DE 10 2007 014 829 B3中公开的、用于分类散射的方法中,引入的对象切线在确定散射平面时出现多义性,其导致不足以估计在对像上形成的散射束轮廓。
利用专用传感器测量的方法(ii)一方面需要附加的硬件,该硬件可以构成双源设备的制造开销的不小部分。另一方面,在该方法中,散射束轮廓不直接在随后在其上执行校正的z位置处测量。也就是,不采集或者仅空间上强烈欠扫描地采集散射的z轮廓的可能的变化。这点在探测器的z覆盖增加的情况下导致测量出的散射束轮廓与必要的数据校正逐渐缺少一致。
交替消隐相应的原始辐射并且在消隐阶段直接测量横向散射的第三方法(iii)虽然在该处在z方向上的数据可能不缺少一致,但在此的数据不能任意细微地按照角方向进行扫描,这点导致在散射数据中以及在起始数据中对应的混淆差错。此外,通过消隐原始辐射增加了图像噪声。同样,剂量效率变得更低,因为在实际上辐射器的开阶段和关阶段不可能任意短。特别地,快速消隐X射线辐射对CT设备的X射线辐射器和高压发电机提出了高要求。
发明内容
因此,本发明要解决的技术问题是,找到用于在多源CT中对对象进行CT检查时的散射束校正的一种改善的方法以及一种改善的计算机***。
发明人认识到,在迭代的图像重建的范围内可以根据重建的图像数据确定在双源扫描中出现的横向散射强度的轮廓,并且可以迭代地用于校正原始数据。也就是,通过每个迭代步骤得出更好地近似于在原始数据校正中使用的散射束强度。也就是成立:
fk+1=fk-Q(A·fk+S(fk)-t)-R(fk)
在此,fk是在迭代步骤k中的图像,A是前向投影的算子,Q表示例如根据WFBP方法(WFBP=Weighted Filtered Back-Projektion,加权滤波反投影)的反投影,t表示带有散射的原始数据以及R表示合适的正则化项。
按照本发明该迭代方法可以通过下面描述的散射束估计S(fk)的不同方案来执行。
方案1:根据对从探测器元件到初始图像的射线的追溯计算横向散射。
在此,假定具有准直器的探测器装置,其中单个的探测器元件准直到相对布置的管的焦点,例如为此可以使用二维散射束栅格。通过该准直器对于每个探测器元件定义原始射线,沿着这些原始射线在对象或在重建图像中识别散射中心,并且必要时关于其散射性能进行分类。由此,对于每个散射中心i对于特定的管位置,散射角θi是已知的,同样,散射射线的入射强度和能量分布也是已知的,从而可以由此计算各个微分的散射横截面。在第二步骤中,每个该计算出的散射幅度经受单独的衰减,该衰减由对从散射中心位置出发直至相应的探测器元件的对象衰减执行线性积分得出。由此,计算一阶散射分布,也就是每量子的一次散射事件。考虑限制散射中心的数量作为进一步简化,其也可以取决于对探测器进行线性积分的值。在最简单的情况下可以仅考虑在对象表面上的散射中心。
更高阶,也就是每量子的多次散射事件,可以通过所描述的措施的级联包含在内,但计算时间是开销非常大的。
方案2:从初始重建的图像中确定散射平面。
该方案2基于如下的认识,即图像有效的横向散射主要在散射对象的表面附近形成。也就是对对象表面的了解越精确,则对对象的横向散射分布的确定就会越精确。在迭代重建的范围内,在此可以考虑如下的子方案:
子方案a:根据查表(Tabellenzugriff)确定散射束轮廓。该方法表示基于模型的横向散射校正的扩展。在此,对每个投影的特定的、列表的散射束分布的选择通过分置三个参数来进行,这三个参数例如是在对象与横向散射的探测器之间的气隙的宽度、在对象与引起横向散射的X射线源之间的气隙的宽度以及在散射表面的法线与入射射线之间的角度。原则上也可以考虑其它参数,例如还可以考虑在散射平面的平面法线与散射射线之间的角度。但是按照本发明,感兴越的对象边缘直接从初始或最后可用的重建图像中得出。此外,根据重建的层厚也可以考虑横向散射的z变化,因为对于每个重建的层给出不同的参数组。
子方案b:表面散射的模型。作为方案1的简化可以根据函数
Figure BDA0000134696510000031
确定向探测器元件发射的原始射线中的每一个的散射。在此,
Figure BDA0000134696510000032
是探测器元件的原始射线与在弯曲的对象表面上的参考点(Aufpunkt)的平面法线之间的角度。角度
Figure BDA0000134696510000033
是在该平面法线与到X射线源的连接线之间的角度。函数表征了发生的横向散射过程。例如,如果横向散射的基本过程简单地是反射,(入射角=出射角)则函数与狄拉克δ函数成比例:
Figure BDA0000134696510000036
但X射线量子与对象表面的相互作用是更复杂的。由光子效应和康普顿散射组成的、在对象表面上发生的X射线量子的散射过程例如可以通过在函数
Figure BDA0000134696510000037
中解析的模型来采集。第二可能性是通过蒙特卡罗法(Monte-Carlo-Verfahren)对情况进行仿真。第三可能性在于,直接测量不同角度
Figure BDA0000134696510000038
的函数并且将其存储在二维表格中。在此,函数可以作为分量
Figure BDA00001346965100000310
的加权和来理解,从而成立:
Figure BDA00001346965100000311
在此,每个分量n表示在平面上特定阶的散射:例如具有不同半径的平面、球体等。在迭代重建期间,表面的每个点的相应的权重an可以根据该点的局部环境来确定。在最简单的情况下,
Figure BDA0000134696510000041
也可以仅近似于平面上的散射,具有:
Figure BDA0000134696510000042
方案3:利用蒙特卡罗法计算横向散射。
也可以通过蒙特卡罗法对横向散射进行预测,其中对从源出发的、穿过对象并到达探测器的大量光子的路径进行仿真。在此,蒙特卡罗法可以直接将对象的没有散射校正的初始重建的三维衰减分布作为散射对象使用。这允许得出对散射束分布非常精确地预测。但该方法非常麻烦,因为一方面穿过按像素的对象跟踪射线要求大量计算操作,并且另一方面仅有很少的散射光子到达探测器,而此必须对极多的原始光子进行仿真。补充地,为了确定物质分布也可以以公知的方式对初始或最后重建的取决于能量的图像数据进行多成分地分解。为了加速并简化可以采用如下的子方案:
子方案a:为了简化和加速仿真,对象可以通过由例如水的单一物质组成的对象替代。
子方案b:实际的对象可以通过类似的、例如椭圆圆柱体的几何对象或者通过几个几何对象的组合替代。
子方案c:在子方案b中附加地,对象可以由例如水和骨骼的几种物质组成。
子方案d:光子仅能跟踪至给定的散射阶,优选仅跟踪至第一散射阶。
子方案e:为了加速仿真以及为了降低噪声,基于到达探测器的光子的少的数量,按照公知的有效横截面和概率根据散射光子的预计值关闭每个相互作用点的所有探测器信道。
在所有的方法中,也能单独外加地处理在患者卧榻处的散射。针对测量场中的每个卧榻位置,例如可以通过直接测量预先确定患者卧榻的横向散射。虽然患者卧榻的衰减通常可以忽略,但是在卧榻的下侧和侧棱上的横向散射可以明显地为整个横向散射强度提供份额。例如如果在图像中这样选择阈值,使得仅考虑将在患者表面上和在患者内部的点作为散射中心,则该过程提供了散射束模型的改善。
此处提供的方法的主要优点在于,在CT设备的辐射器或发生器上无需附加的传感器或者没有附加的要求就已足够。与基于模型的方法相比优点在于,由通过少量的相对于当前测量时间上在前的切线的正弦图数据形成的对象表面的近似特征可以通过由当前重建的CT图像形成的对象表面的更稳健的特征代替。由此可以避免多义性,并且能回避在z方向上近似均匀的对象的前提条件。另一个优点是对横向散射的更灵活的建模。
相应于上面提到的基本构思,发明人建议一种用于在多源CT(特别是双源CT)中对对象(特别是患者)进行CT检查时的散射束校正的方法,具有如下的方法步骤:
-通过至少两个在平面上角度错开地布置的辐射器-探测器***对对象进行扫描,其中每个辐射器将射线束分别发射到相对布置的探测器并且在每个探测器上采用了散射格栅,并且产生原始投影数据组,
-利用至少一个探测器的原始投影数据组重建对象,
-确定由每个辐射器所产生的、仅在该至少一个其它辐射器的朝向其相对布置的探测器的原始射线的方向上的散射,
-通过从原始投影数据组中去除所计算出的散射,产生校正后的投影数据组,
-利用校正后的投影数据组重建对象,
-检查至少一个中断标准,并且
-在确定散射的情况下执行该方法的进一步迭代或者输出重建结果,如果对至少一个中断标准的检查的结果得到这点的话。
按照本发明可以将一个或多个如下标准的组合用作中断标准:迭代的最大次数、使图像结果与预定值相适应和/或使确定的散射结果与预定值相适应。
在此具有优势地,在该方法的替换中建议仅考虑在对象的表面上形成的散射。
在此,为了确定在对象表面上的散射可以考虑对象的表面的局部曲率、入射射线的入射角以及在至少一个其它探测器的方向上的至少一个反射的散射的至少一个反射角。在此,还可以采用三维表格,其中录入了来自预先试验的基于入射角、反射角以及散射的对象的表面曲率的散射值。此外,可以借助之前重建的图像数据确定散射的对象的表面曲率。
按照该方法的另一种替换,可以在整个对象中在产生散射的射线束的原始射线与对于各个其它射线路径来说至少一个其它射线束的其它射线路径的至少一部分交点上,在考虑在对象中衰减的情况下,确定原始射线的强度和仅在至少一个其它射线路径的方向上形成的散射及其在探测器上的强度。
在此,为了确定在观察的交点上的散射,还可以直接通过考虑到入射角、散射角以及由最后的重建数据所确定的散射物质的函数来计算该散射。
在一种特别简单的实施中,为了确定在形成散射的位置上的原始射线的强度和/或能量,可以假设均匀的物质分布。此外,为了确定在形成散射的位置上的原始射线的强度,可以假设预定的与身体类似的物质分布或者假设至少近似于从对象的上述(特别是最后)重建得出的物质分布。
作为对直接按照公式计算的替换,为了确定在观察的交点上的散射,也可以通过蒙特卡罗仿真来确定其。在此,同样可以在对象中假设均匀的物质分布,更好地假设预定的与身体类似的物质分布或者还更好地假设至少近似于从对象的上述(优选最后)重建得出的物质分布。
为了避免过长的计算时间,还有利的是,在蒙特卡罗仿真中光子仅跟踪至预定的散射阶,特别是跟踪至第一散射阶。
此外,在所有上面提到的方法中,为了确定关于对象的散射分布,还可以加上预先已知的患者卧榻的散射分布。
如果在CT检查的范围内实施上面描述的方法,其中或者通过采用不同的X射线谱和/或通过采用能量分辨的探测器开辟对已知物质的分解的可能性,则为了确定散射可以在CT图像中考虑这种此前计算出的物质分布。
除了上述方法之外,本发明还涉及一种具有计算机程序存储器的多源或双源CT***的计算机***,其中存储器也包含如下的计算机程序,其在运行时执行如下的方法步骤:
-通过至少两个在平面上角度错开地布置的辐射器-探测器***对对象进行扫描,其中每个辐射器将射线束分别发射到相对布置的探测器并且在每个探测器上采用了散射格栅,并且产生原始投影数据组,
-利用至少一个探测器的原始投影数据组重建对象,
-确定由每个辐射器所产生的、仅在该至少一个其它辐射器的朝向其相对布置的探测器的原始射线的方向上的散射,
-通过从原始投影数据组中去除所计算出的散射,产生校正后的投影数据组,
-利用校正后的投影数据组重建对象,
-检查至少一个中断标准,并且
-在确定散射的情况下执行该方法的进一步迭代或者输出重建结果,如果对至少一个中断标准的检查的结果得到这点的话。
本发明还涉及一种不必直接与CT***相连的计算机***,并且被构造为具有计算机程序存储器以用于图像处理,其中存储器也包含如下的计算机程序,其在运行时执行如下的方法步骤:
-接收由对对象扫描得到的探测器数据或投影数据,
-利用根据探测器数据计算出的或接收的至少一个探测器的原始投影数据组重建对象,
-确定由每个辐射器所产生的、仅在该至少一个其它辐射器的朝向其相对布置的探测器的原始射线的方向上的散射,
-通过从原始投影数据组中去除所计算出的散射,产生校正后的投影数据组,
-利用校正后的投影数据组重建对象,
-检查至少一个中断标准,并且
-在确定散射的情况下执行该方法的进一步迭代或者输出重建结果,如果对至少一个中断标准的检查的结果得到这点的话。
当然,在这种计算机***中存储器中的计算机程序也可以这样构造,使得在运行时执行其它上面详细描述的方法步骤。
附图说明
下面结合优选的实施例借助附图对本发明作进一步描述,其中仅示出为了理解本发明所必要的特征。采用了如下的附图标记:1:CT***;2:第一X射线管;3:第一探测器;3.1:散射格栅(Streustrahlungsgitter);4:第二X射线管;5:第二探测器;5.1:散射格栅;6:机架壳体;7:造影剂施放器;8:患者卧榻;9:***轴;10:计算机***;Ba:第一射线束;Bb:第二射线束;i:散射中心;K:机架的运转圆(Laufkreis);N:面法线;O:原始射线;P:患者/对象;Prg1-Prgn:计算机程序;S:散射束;S1:对对象进行扫描;S2:产生原始投影数据组;S3:利用原始投影数据组重建对象;S4:确定由每个辐射器产生的、仅在分别向其它探测器发射的方向上的散射;S5:产生校正后的投影数据组;S6:利用校正后的投影数据组重建对象;S7对至少一个中断标准进行检查;S8:输出、存储或进一步处理重建结果;Z:产生散射的射线;θi:在散射中心i处的散射角;
Figure BDA0000134696510000081
入射角;
Figure BDA0000134696510000082
散射角。
附图中:
图1示出了双源CT***,
图2示出了按照本发明的方法的流程图,
图3示出了计算对象的相对于表面的横向散射的示意图,以及
图4示出了计算对象的相对于体积的横向散射的示意图。
具体实施方式
图1示出了具有机架壳体6的双源CT***(=具有两个辐射器-探测器***的CT***)1的示例性的附图,其中在未详细示出的机架上固定了两个角度错开地布置的辐射器-探测器***。辐射器-探测器***一方面由第一X射线管2与对应于第一X射线管相对设置的探测器3的组成,而在另一方面由第二X射线管4与对应于第二X射线管相对设置的探测器5组成。两个辐射器-探测器***都扫掠在中央圆形开口中设置的测量场。患者P可以借助患者卧榻8沿着***轴9移动穿过该测量场。原则上因此既执行螺旋扫描也执行顺序扫描。为了改善血管或其它构造的成像也可以通过造影剂施放器7向患者注射造影剂。
CT***1的控制以及对患者P的扫描的分析通过与此相连接的计算机***10来执行,其中该计算机***具有至少一个在其中存储了计算机程序Prg1-Prgn的存储器。按照本发明,在其中也包含或存储这样构造的程序,即其在***的运行时执行按照本发明的方法的不同实施方式。
为了描述按照本发明的方法在图2中示出了流程图,该流程图示出了基本的方法步骤S1至S8。在步骤S1中通过至少两个在平面上角度错开地布置的辐射器-探测器***对对象进行扫描。在此假定布置类似于图1的辐射器-探测器***,其中每个辐射器将原始射线束分别发射到相对布置的探测器。同时该原始射线束在对象上产生散射,但是其中在一个探测器上仅探测由于在每个探测器上所使用的散射格栅在相应的原始射线的方向上延伸的散射束。
在第二步骤S2中,从所测量的探测器数据中产生原始投影数据组。
在步骤S3中,利用至少一个探测器的原始投影数据组来重建被扫描对象的CT图像。
然后在步骤S4中,基于通过前面重建所提供的关于对象的实施方式的信息,即特别是其几何数据和/或物质分布,确定由每个辐射器所产生的、仅在该至少一个其它辐射器的朝向其相对布置的探测器的原始射线的方向上的散射。
在下一个步骤S5中,将此前计算出的散射的影响从原始投影数据组中去除。
此时在步骤S6中,利用散射校正后的投影数据组进行对象的重建。
然后在步骤S7中,对至少一个迭代的中断标准进行检查,并且判断应当在步骤S4中利用改变了的CT图像数据再次重新确定更好的散射数据还是应当中断迭代。
最后在步骤S8中,输出、存储或进一步处理重建结果,也就是摆脱了散射影响的CT图像数据。
关于对在各个不同的辐射器-探测器***的探测器中可探测的散射进行确定,在图3和图4中分别示出了在仅考虑表面的情况下或在考虑整个对象的情况下进行确定的方案。
在图3中可以看出机架的运转圆K,在该机架上可以看出角度错开地布置的第一和第二辐射器-探测器***。属于第一辐射器-探测器***的是X射线管2和相对布置的探测器3,在其之间示出了具有边界射线的射线束Ba。属于第二辐射器-探测器***的是X射线管4和相对布置的探测器5,在其之间示出了具有其边界射线的射线束Bb。此外,可以看出一条从第一管2向探测器3的所选择的原始射线O,该原始射线通过由第二管4发出的射线的从散射中心i处起始的散射束S来叠加。同样画出了在原始射线O与第二射线Z之间的散射角θi
也就是说,图3示出了在根据对从探测器元件到初始图像的原始射线的追溯计算横向散射时的几何关系。沿着该原始射线O存在多个散射中心i。对于每个该散射中心已知散射角θi,在该散射角之下散射将(不期望的)份额分别提供到其它辐射器-探测器***。对于关于该原始射线的一阶散射份额来说,所有的点i的份额可以按照沿着原始射线穿过对象的衰减的加权地来计算并合计。如果这点对于所有散射中心或对于至少代表性的多个散射中心进行,则在考虑整个对象的情况下得出良好地近似于实际出现的散射的强度分布。
在图4中示出了原则上类似的情况,其中示出了仅在对象P的表面上确定散射。在此,辐射器-探测器***的结构与图3一致,但其中仅考虑在对象P的表面上的散射中心,并且相应地需要计算更少的散射中心及其影响。此处,对于所有在表面上的散射中心来说,相应于原始射线O和产生横向散射的射线Z的方向,相对于面法线N的入射角
Figure BDA0000134696510000101
和散射角
Figure BDA0000134696510000102
是已知的。相应地,也能在各个对探测器元件有效的散射中心处计算并且由此从重建的CT图像数据组中确定那里所产生的散射强度。
总之,本发明因此介绍了用于在多源CT中对对象进行CT检查时的散射束校正的一种方法和一种计算机***,具有如下的主要方法步骤:
-产生原始投影数据组,
-利用至少一个探测器的原始投影图像数据组重建对象,
-确定由每个辐射器所产生的、仅在该至少一个其它辐射器的朝向其相对布置的探测器的原始射线的方向上的散射,
-通过从原始投影数据组中去除所计算出的散射,产生校正后的投影数据组,
-利用校正后的投影数据组重建对象,并且
-如果符合至少一个预先给定的中断标准,则在确定散射的情况下执行该方法的进一步迭代或者输出重建结果。
可以理解的是,上述提到的本发明的特征不仅按照分别给出的组合,而且按照其它组合或在单独设置时也是适用的,而不会脱离本发明的范围。

Claims (20)

1.一种用于在多源CT、特别是双源CT中对对象进行CT检查时的散射束校正的方法,具有如下的方法步骤:
1.1.通过至少两个在平面上角度错开地布置的辐射器-探测器***(2,3;4,5)对对象(P)进行扫描,其中每个辐射器(2,4)将射线束(Ba,Bb)分别发射到相对布置的探测器(3,5)并且在每个探测器(3,5)上采用了散射格栅(3.1,5.1),并且产生原始投影数据组,
1.2.利用至少一个探测器(3,5)的原始投影数据组重建对象(P),
1.3.确定由每个辐射器(2,4)所产生的、仅在该至少一个其它辐射器的朝向其相对布置的探测器(3,5)的原始射线(O)的方向上的散射,
1.4.通过从所述原始投影数据组中去除所计算出的散射,产生校正后的投影数据组,
1.5.利用所述校正后的投影数据组重建对象(P),
1.6.检查至少一个中断标准,并且
1.7.从特征1.3开始执行所述方法的进一步迭代,或者根据由特征1.6得出的结果输出重建结果。
2.根据上述权利要求1所述的方法,其特征在于,仅考虑在对象(P)的表面上形成的散射。
3.根据上述权利要求2所述的方法,其特征在于,为了确定在对象(P)的表面上的散射,可以考虑对象表面的曲率、入射射线(Z)的入射角
Figure FDA0000134696500000011
以及在至少一个其它探测器(3)的方向上的至少一个反射的散射(S)的至少一个反射角
4.根据上述权利要求3所述的方法,其特征在于,为了确定散射采用三维表格,其中录入了来自预先试验的基于入射角
Figure FDA0000134696500000013
反射角
Figure FDA0000134696500000014
以及散射的对象的表面曲率的散射值。
5.根据上述权利要求4所述的方法,其特征在于,所述散射的对象(P)的表面曲率借助此前重建的图像数据来确定。
6.根据上述权利要求1所述的方法,其特征在于,在整个对象(P)中在产生散射的射线束的原始射线(O)与对于各个其它射线路径来说至少一个其它射线束的其它射线路径的至少一部分交点上,在考虑在对象(P)中衰减的情况下确定所述原始射线(O)的强度和仅在该至少一个其它射线路径的方向上形成的散射及其在探测器上的强度。
7.根据上述权利要求6所述的方法,其特征在于,为了确定在观察的交点上的散射,直接通过考虑到入射角散射角
Figure FDA0000134696500000022
以及由最后的重建数据所确定的散射物质的函数来计算所述散射。
8.根据上述权利要求7所述的方法,其特征在于,为了确定在形成散射的位置上的原始射线(O)的强度和/或能量,假设均匀的物质分布。
9.根据上述权利要求7所述的方法,其特征在于,为了确定在形成散射的位置上的原始射线(O)的强度,假设预定的与身体类似的物质分布。
10.根据上述权利要求7所述的方法,其特征在于,为了确定在形成散射的位置上的原始射线(O)的强度,假设至少近似于从对象(P)的最后重建得出的物质分布。
11.根据上述权利要求6所述的方法,其特征在于,为了确定在所观察的交点上的散射,通过蒙特卡罗仿真来确定所述散射。
12.根据上述权利要求11所述的方法,其特征在于,为了计算蒙特卡罗仿真,在对象(P)中假设均匀的物质分布。
13.根据上述权利要求11所述的方法,其特征在于,为了计算蒙特卡罗仿真,在对象(P)中假设预定的与身体类似的物质分布。
14.根据上述权利要求11所述的方法,其特征在于,为了计算蒙特卡罗仿真,假设至少近似于从对象(P)的最后重建中得出的物质分布。
15.根据上述权利要求11至14中任一项所述的方法,其特征在于,在蒙特卡罗仿真中光子仅跟踪至预定的散射阶,特别是跟踪至第一散射阶。
16.根据上述权利要求1至15中任一项所述的方法,其特征在于,为了确定关于对象的散射分布,加上预先已知的患者卧榻(8)的散射分布。
17.根据上述权利要求1至16中任一项所述的方法,其特征在于,为了确定散射,通过基于取决于能量的衰减值进行物质成分分解,来考虑物质分布。
18.一种多源或双源CT***的计算机***(10),具有计算机程序存储器,其特征在于,所述存储器也包含如下的计算机程序(Prg1-Prgn),其在运行时执行如下的方法步骤:
18.1.通过至少两个在平面上角度错开地布置的辐射器-探测器***(2,3;4,5)对对象(P)进行扫描,其中每个辐射器(2,4)将射线束(Ba,Bb)分别发射到相对布置的探测器(3,5)并且在每个探测器(3,5)上采用了散射格栅(3.1,5.1),并且产生原始投影数据组,
18.2.利用至少一个探测器(3,5)的原始投影数据组重建对象(P),
18.3.确定由每个辐射器(2,4)产生的、仅在该至少一个其它辐射器的朝向其相对布置的探测器(3,5)的原始射线(O)的方向上的散射,
18.4.通过从所述原始投影数据组中去除计算出的散射,产生校正后的投影数据组,
18.5.利用所述校正后的投影数据组重建对象(P),
18.6.检查至少一个中断标准,并且
18.7.从特征18.3开始执行所述方法的进一步迭代,或者根据由特征18.6得出的结果输出重建结果。
19.一种计算机***(10),被构造为用于图像处理,具有计算机程序存储器,其特征在于,所述存储器也包含如下的计算机程序(Prg1-Prgn),其在运行时执行如下的方法步骤:
19.1.接收由对对象(P)扫描得到的探测器数据或原始投影数据,
19.2.利用从探测器数据中计算出的或接收的至少一个探测器(3,5)的原始投影数据组重建对象(P),
19.3.确定由每个辐射器(2,4)产生的、仅在该至少一个其它辐射器的朝向其相对布置的探测器(3,5)的原始射线(O)的方向上的散射,
19.4.通过从所述原始投影数据组中去除计算出的散射,产生校正后的投影数据组,
19.5.利用所述校正后的投影数据组重建对象(P),
19.6.检查至少一个中断标准,并且
19.7.从特征19.3开始执行所述方法的进一步迭代,或者根据由特征19.6得出的结果输出重建结果。
20.根据上述权利要求18或19所述的计算机***(10),其特征在于,所述存储器中的计算机程序(Prg1-Prgn)被这样构造,使得其在运行时执行上述权利要求2至17中任一项所述的方法步骤。
CN201210028247.5A 2011-02-23 2012-02-09 在多源ct中用于散射束校正的方法和计算机*** Expired - Fee Related CN102648857B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102011004598.8 2011-02-23
DE102011004598.8A DE102011004598B4 (de) 2011-02-23 2011-02-23 Verfahren und Computersystem zur Streustrahlkorrektur in einem Multi-Source-CT

Publications (2)

Publication Number Publication Date
CN102648857A true CN102648857A (zh) 2012-08-29
CN102648857B CN102648857B (zh) 2015-10-28

Family

ID=46604805

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210028247.5A Expired - Fee Related CN102648857B (zh) 2011-02-23 2012-02-09 在多源ct中用于散射束校正的方法和计算机***

Country Status (3)

Country Link
US (1) US8744161B2 (zh)
CN (1) CN102648857B (zh)
DE (1) DE102011004598B4 (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106709913A (zh) * 2016-12-30 2017-05-24 上海联影医疗科技有限公司 锥束ct入射强度计算方法和装置
CN107837090A (zh) * 2016-09-20 2018-03-27 西门子医疗有限公司 计算机断层成像中基于正弦图的散射射线校正
CN108351429A (zh) * 2015-10-28 2018-07-31 皇家飞利浦有限公司 Ct***和ct方法
WO2020082280A1 (en) * 2018-10-25 2020-04-30 Shanghai United Imaging Healthcare Co., Ltd. System and method for scatter correction
CN111096761A (zh) * 2018-10-29 2020-05-05 上海西门子医疗器械有限公司 修正楔形滤波器散射的方法、装置和相关设备
CN111657979A (zh) * 2019-03-08 2020-09-15 江苏一影医疗设备有限公司 Ct成像***及其成像方法
CN114113173A (zh) * 2021-11-18 2022-03-01 上海联影医疗科技股份有限公司 一种x射线设备、应用于x射线设备中的散射校正方法

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8437525B2 (en) * 2011-03-28 2013-05-07 Siemens Medical Solutions Usa, Inc. Method and system for using a modified ordered subsets scheme for attenuation weighted reconstruction
DE102012204980B4 (de) 2012-03-28 2021-09-30 Siemens Healthcare Gmbh Verfahren zur Rekonstruktion von CT-Bildern mit Streustrahlenkorrektur, insbesondere für Dual-Source CT-Geräte
DE102013104720A1 (de) * 2012-05-07 2013-11-07 Werth Messtechnik Gmbh Computertomografie-Verfahren und Anordnung zur Bestimmung von Merkmalen an einem Messobjekt
US9724056B2 (en) * 2013-11-28 2017-08-08 Toshiba Medical Systems Corporation Method and system for spectral computed tomography (CT) with inner ring geometry
JP6783702B2 (ja) * 2017-05-25 2020-11-11 東芝Itコントロールシステム株式会社 X線断層撮影装置
JP6925294B2 (ja) * 2018-03-23 2021-08-25 富士フイルム株式会社 画像処理装置、放射線画像撮影システム、画像処理方法、及び画像処理プログラム
US11166690B2 (en) * 2020-03-19 2021-11-09 Accuray, Inc. Noise and artifact reduction for image scatter correction
CN114549360A (zh) * 2022-02-25 2022-05-27 清华大学 多能量ct去散射降噪物质分解方法、装置、设备及介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020048339A1 (en) * 2000-09-27 2002-04-25 Stefan Schneider Device and method for forming a computed X-ray tomogram with scatter correction
CN101028195A (zh) * 2005-10-10 2007-09-05 西门子公司 对计算机断层造影***进行辐射校正的方法
US20080159469A1 (en) * 2006-09-28 2008-07-03 Ernst-Peter Ruhrnschopf Method for combined bone hardening and scattered radiation correction in X-ray computed tomography
US20080292170A1 (en) * 2007-05-15 2008-11-27 Herbert Bruder Method for scattered radiation correction in an X-ray computed tomography system, and method for generating a tomographic display corrected for scattered radiation, and/or an X-ray computed tomography system
CN101810488A (zh) * 2009-01-29 2010-08-25 皇家飞利浦电子股份有限公司 探测值校正设备

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6421412B1 (en) * 1998-12-31 2002-07-16 General Electric Company Dual cardiac CT scanner
US20030128801A1 (en) * 2002-01-07 2003-07-10 Multi-Dimensional Imaging, Inc. Multi-modality apparatus for dynamic anatomical, physiological and molecular imaging
CN1481756A (zh) * 2002-08-14 2004-03-17 ��ʽ���綫֥ 集中照射型放疗设备
JP4314008B2 (ja) * 2002-10-01 2009-08-12 株式会社東芝 X線ctスキャナ
DE10302567A1 (de) * 2003-01-22 2004-08-12 Siemens Ag Bildgebendes Tomographiegerät mit wenigstens zwei Strahler-Detektor-Systemen und Verfahren zum Betrieb eines solchen Tomographiegeräts
US20060083351A1 (en) * 2004-10-18 2006-04-20 Ge Medical Systems Global Technology Company, Llc Method and System for Scatter Correction During Bi-Plane Imaging with Simultaneous Exposure
EP1743579A1 (en) * 2005-07-15 2007-01-17 Kabushi Kaisha Toshiba X-Ray CT apparatus
DE102006019920B4 (de) * 2006-04-28 2008-04-10 Siemens Ag Verfahren zur Streustrahlungskorrektur eines CT-Systems mit mindestens zwei winkelversetzt angeordneten Fokus-Detektor-Systemen und Röntgen-CT-System
DE102006019923A1 (de) * 2006-04-28 2007-11-15 Siemens Ag Verfahren zur Streustrahlungskorrektur bei einem Röntgen-CT und Röntgen-CT zur Anwendung dieses Verfahrens
EP2034899B1 (en) 2006-06-22 2016-03-16 Koninklijke Philips N.V. Multi-tube imaging system scatter correction
DE102007014829B3 (de) 2007-03-28 2008-09-11 Siemens Ag Verfahren zur Streustrahlungskorrektur in bildgebenden Röntgengeräten sowie Röntgenbildgebungssystem
JP5214942B2 (ja) * 2007-10-02 2013-06-19 東芝メディカルシステムズ株式会社 多管球x線ctにおける散乱線強度分布のスケーリング方法および多管球x線ct装置
US7912180B2 (en) * 2009-02-19 2011-03-22 Kabushiki Kaisha Toshiba Scattered radiation correction method and scattered radiation correction apparatus

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020048339A1 (en) * 2000-09-27 2002-04-25 Stefan Schneider Device and method for forming a computed X-ray tomogram with scatter correction
CN101028195A (zh) * 2005-10-10 2007-09-05 西门子公司 对计算机断层造影***进行辐射校正的方法
US20080159469A1 (en) * 2006-09-28 2008-07-03 Ernst-Peter Ruhrnschopf Method for combined bone hardening and scattered radiation correction in X-ray computed tomography
US20080292170A1 (en) * 2007-05-15 2008-11-27 Herbert Bruder Method for scattered radiation correction in an X-ray computed tomography system, and method for generating a tomographic display corrected for scattered radiation, and/or an X-ray computed tomography system
CN101810488A (zh) * 2009-01-29 2010-08-25 皇家飞利浦电子股份有限公司 探测值校正设备

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108351429A (zh) * 2015-10-28 2018-07-31 皇家飞利浦有限公司 Ct***和ct方法
CN108351429B (zh) * 2015-10-28 2022-04-15 皇家飞利浦有限公司 Ct***和ct方法
CN107837090A (zh) * 2016-09-20 2018-03-27 西门子医疗有限公司 计算机断层成像中基于正弦图的散射射线校正
CN106709913A (zh) * 2016-12-30 2017-05-24 上海联影医疗科技有限公司 锥束ct入射强度计算方法和装置
CN106709913B (zh) * 2016-12-30 2019-10-08 上海联影医疗科技有限公司 锥束ct入射强度计算方法和装置
WO2020082280A1 (en) * 2018-10-25 2020-04-30 Shanghai United Imaging Healthcare Co., Ltd. System and method for scatter correction
CN111096761A (zh) * 2018-10-29 2020-05-05 上海西门子医疗器械有限公司 修正楔形滤波器散射的方法、装置和相关设备
CN111096761B (zh) * 2018-10-29 2024-03-08 上海西门子医疗器械有限公司 修正楔形滤波器散射的方法、装置和相关设备
CN111657979A (zh) * 2019-03-08 2020-09-15 江苏一影医疗设备有限公司 Ct成像***及其成像方法
CN114113173A (zh) * 2021-11-18 2022-03-01 上海联影医疗科技股份有限公司 一种x射线设备、应用于x射线设备中的散射校正方法

Also Published As

Publication number Publication date
CN102648857B (zh) 2015-10-28
US8744161B2 (en) 2014-06-03
US20120213424A1 (en) 2012-08-23
DE102011004598B4 (de) 2019-07-11
DE102011004598A1 (de) 2012-08-23

Similar Documents

Publication Publication Date Title
CN102648857A (zh) 在多源ct中用于散射束校正的方法和计算机***
US8989342B2 (en) Methods and systems for volumetric reconstruction using radiography
US7308072B2 (en) Device and method for x-ray scatter correction in computer tomography
JP3476381B2 (ja) 二元エネルギーx線撮像のための装置および方法
CN101416073B (zh) 用于重建图像的双能量衰减数据的信噪比的动态优化
CN101489486B (zh) 基于散射辐射分数的x-射线探测器增益校准
EP3400876A1 (en) Automatic estimating and reducing scattering in computed tomography scans
US20100208964A1 (en) Method for eliminating scatter artefacts
CN101473348A (zh) 用于误差补偿的方法和***
CN101011253A (zh) 产生投影或断层造影的相位对比照片的焦点-检测器装置
CN101111758A (zh) 用于x射线投影的校正或扩展的设备和方法
CN106233335A (zh) X射线光谱成像
US9125286B2 (en) X-ray dose estimation technique
WO2019083014A1 (ja) 光子計数型のx線検出データを処理する方法及び装置、並びに、x線装置
US20210307714A1 (en) X-ray imaging system
CN110389138A (zh) 使用Alvarez-Macovski衰减模型在断层摄影重建中进行X射线束硬化校正
CA3008893A1 (en) Three-dimensional image reconstruction using transmission and scatter radiography methods
US9615807B2 (en) Systems and methods for improving image quality in cone beam computed tomography
US20170086757A1 (en) Multiple Emission Energies in Single Photon Emission Computed Tomography
CN110197515A (zh) 基于康普顿背散射扫描的医学ct机的图像重建方法和***
US11464470B2 (en) X-ray backscatter systems and methods for performing imaging tomosynthesis
Ordoñez et al. Fast in situ image reconstruction for proton radiography
US7385200B2 (en) Re-binning method for nuclear medicine imaging devices
CN107212898A (zh) 图像重建方法
Cruvinel et al. Compton scattering tomography for agricultural measurements

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220127

Address after: Erlangen

Patentee after: Siemens Healthineers AG

Address before: Munich, Germany

Patentee before: SIEMENS AG

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: 20151028

Termination date: 20220209