CN107072627A - 用于针对心脏和肺图像的光子计数ct的心脏重建 - Google Patents
用于针对心脏和肺图像的光子计数ct的心脏重建 Download PDFInfo
- Publication number
- CN107072627A CN107072627A CN201580056917.6A CN201580056917A CN107072627A CN 107072627 A CN107072627 A CN 107072627A CN 201580056917 A CN201580056917 A CN 201580056917A CN 107072627 A CN107072627 A CN 107072627A
- Authority
- CN
- China
- Prior art keywords
- weight
- data
- cardiac motion
- projection
- flux
- 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
- 230000000747 cardiac effect Effects 0.000 title claims abstract description 96
- 210000004072 lung Anatomy 0.000 title description 18
- 230000004907 flux Effects 0.000 claims abstract description 123
- 238000000034 method Methods 0.000 claims abstract description 21
- 230000005855 radiation Effects 0.000 claims description 42
- 238000003384 imaging method Methods 0.000 claims description 33
- 238000007689 inspection Methods 0.000 claims description 14
- 238000001514 detection method Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 230000001932 seasonal effect Effects 0.000 claims description 2
- 230000001678 irradiating effect Effects 0.000 claims 2
- 238000002591 computed tomography Methods 0.000 description 14
- 230000029058 respiratory gaseous exchange Effects 0.000 description 7
- 238000001228 spectrum Methods 0.000 description 7
- 230000002861 ventricular Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 238000012986 modification Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000000718 qrs complex Methods 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 230000008602 contraction Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000000241 respiratory effect Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 229910004613 CdTe Inorganic materials 0.000 description 1
- 229910004611 CdZnTe Inorganic materials 0.000 description 1
- 229910001218 Gallium arsenide Inorganic materials 0.000 description 1
- 230000001746 atrial effect Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- 229910052732 germanium Inorganic materials 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- QIQXTHQIDYTFRH-UHFFFAOYSA-N octadecanoic acid Chemical compound CCCCCCCCCCCCCCCCCC(O)=O QIQXTHQIDYTFRH-UHFFFAOYSA-N 0.000 description 1
- 210000003540 papillary muscle Anatomy 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- 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/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/405—Source units specially adapted to modify characteristics of the beam during the data acquisition process
-
- 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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4241—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using energy resolving detectors, e.g. photon counting
-
- 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/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
-
- 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/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/503—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
-
- 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/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- 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/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
-
- 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/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5288—Devices using data or image processing specially adapted for radiation diagnosis involving retrospective matching to a physiological signal
-
- 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/54—Control of apparatus or devices for radiation diagnosis
- A61B6/541—Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
一种方法包括,与心动周期信号协作地在第一较低通量水平与第二较高通量水平之间调制发射辐射的通量,使得所述通量在具有第一较高心脏运动的第一心脏运动阶段期间处于第一较低通量水平处,并且在具有第二较低心脏运动的第二心脏运动阶段期间处于第二较高通量水平处。所述方法还包括利用第一重建窗口重建所述投影数据,所述第一重建窗口将第一权重应用于与所述第一心脏运动阶段和较低的第一通量水平相对应的所述投影数据的第一子集,并且将不同的第二权重应用于与所述第二心脏运动阶段和较高的第二通量水平相对应的所述投影数据的第二子集,以生成第一体积图像数据。
Description
技术领域
以下总体涉及利用光子计数探测器的心脏成像,并且被描述为具体应用于计算机断层摄影(CT)。然而,以下还涉及其他成像应用。
背景技术
计算机断层摄影(CT)扫描器包括发射X射线辐射的X射线管。辐射穿过位于视场中的对象或物体,并且由此被衰减。探测器阵列探测穿过视场的辐射,并且产生表示其的信号。重建器重建信号以产生体积图像数据。
被配置用于谱成像的计算机断层摄影扫描器已经包括光子计数探测器像素。在这样像素的情况下,生成的信号包括均具有指示探测到的光子的峰值幅度或峰值高度的电流或电压。该信号是被求解并且随后被处理的能量。重建器重建谱和非谱体积图像数据。
心脏CT需要高通量以在如在扫描期间视觉地冻结心脏运动所需的快速机架旋转期间得到足够的信噪比。在高通量心脏CT扫描的情况下,一些探测器将会接收高通量。例如,接收穿过患者的周边和/或患者的肺的辐射的探测器将会接收高通量,其中,患者仅具有微小的衰减。
遗憾地,光子计数探测器像素具有有限的计数速率。因此,被配置为具有光子计数探测器像素的CT扫描器不会非常适合于高通量心脏CT研究。
发明内容
在一个方面中,一种方法包括与心动周期信号协作地在第一通量水平与不同的第二通量水平之间调制发射辐射的通量,使得所述通量在具有第一心脏运动的第一心脏运动阶段期间处于所述第一通量水平处,并且在具有第二心脏运动的第二心脏运动阶段期间处于所述第二通量水平处,其中,所述第一通量水平小于所述第二通量水平,并且其中,所述第一心脏运动大于所述第二心脏运动。所述方法还包括探测经调制的发射辐射并且产生指示探测到的经调制的发射辐射的投影数据。所述方法还包括利用第一重建窗口重建所述投影数据,所述第一重建窗口将第一权重应用于对应于所述第一心脏运动阶段和较低的第一通量水平的所述投影数据的第一子集,并且将不同的第二权重应用于对应于所述第二心脏运动阶段和较高的第二通量水平的所述投影数据的第二子集,以生成第一体积图像数据,其中,所述第一权重大于所述第二权重。
在另一方面中,一种成像系统,包括:辐射源,其被配置为关于检查区域旋转并且发射穿过所述检查区域的辐射。辐射源控制器被配置为控制所述辐射源与心动周期信号协作地在第一通量水平与不同的第二通量水平之间调制发射辐射的通量,使得所述通量在具有第一心脏运动的第一心脏运动阶段期间处于所述第一通量水平处,并且在具有第二心脏运动的第二心脏运动阶段期间处于所述第二通量水平处,其中,所述第一通量水平小于所述第二通量水平,并且其中,所述第一心脏运动大于所述第二心脏运动。辐射敏感像素的阵列被配置为探测穿过所述检查区域的辐射并且生成指示探测到的辐射的信号。重建器被配置为利用第一重建窗口重建所述投影数据,所述第一重建窗口将第一权重应用于对应于所述第一心脏运动阶段和较低的第一通量水平的所述投影数据的第一子集,并且将不同的第二权重应用于对应于所述第二心脏运动阶段和较高的第二通量水平的所述投影数据的第二子集,以生成第一体积图像数据,其中,所述第一权重大于所述第二权重。
在又一方面中,一种被编码有计算机可读指令的计算机可读存储介质,所述计算机可读指令在由处理器运行时令所述处理器:与心动周期信号协作地在第一通量水平与不同的第二通量水平之间调制发射辐射的通量,使得所述通量在具有第一心脏运动的第一心脏运动阶段期间处于所述第一通量水平处,并且在具有第二心脏运动的第二心脏运动阶段期间处于所述第二通量水平处,其中,所述第一通量水平小于所述第二通量水平,并且其中,所述第一心脏运动大于所述第二心脏运动;接收由探测经调制的发射辐射的探测器生成的投影数据;并且利用第一重建窗口重建所述投影数据,所述第一重建窗口将第一权重应用于对应于所述第一心脏运动阶段和较低的第一通量水平的所述投影数据的第一子集,并且将不同的第二权重应用于对应于所述第二心脏运动阶段和较高的第二通量水平的所述投影数据的第二子集,以生成第一体积图像数据,其中,所述第一权重大于所述第二权重。
附图说明
本发明可以采取各种部件和部件的布置,以及各种步骤和步骤的安排的形式。附图仅出于图示优选实施例的目的,并且不得被解释为对本发明的限制。
图1示意地图示了具有光子计数探测器像素的阵列的范例成像系统。
图2图示了ECG信号的范例。
图3图示了结合范例ECG信号的范例通量调制模式。
图4图示了结合范例ECG信号的范例心脏重建窗口。
图5图示了结合范例ECG信号的范例肺重建窗口。
图6图示了结合范例ECG信号的另一范例肺重建窗口。
图7图示了结合范例ECG信号的另一范例肺重建窗口。
图8图示了根据本文中的实施例的方法。
具体实施方式
图1示意地图示了范例成像系统100,诸如计算机断层摄影(CT)扫描器。成像系统100包括旋转机架102和固定机架104。旋转机架102由固定机架104可旋转地支撑。旋转机架102被配置为围绕纵轴或z-轴在检查区域106周围旋转。成像系统100还包括在扫描之前、期间和/或之后支撑检查区域106中的对象或物体的对象支撑物107。对象支撑物107也能够用于将对象或物体装入检查区域106和/或从检查区域106取出。
在图示的范例中,心动周期监测器108和特征识别器110与范例成像系统100一起用于心脏扫描。在一个实例中,心动周期监测器108和/或特征识别器110是成像系统100的一部分。在另一实例中,心动周期监测器108和/或心动周期特征识别器110与成像系统100分离。心动周期监测器108从被设置在对象上和/或附近的传感器接收电信号,并且生成指示心脏的电活动的心动周期信号,诸如心电图(ECG或EKG)。
图2图示了具有两个心动周期202和204的ECG信号200的子部分的范例。在图2中,第一或y-轴206表示以毫伏(mV)为单位的电压,而第二或x-轴208表示时间。心动周期202和204中的每一个包括心室收缩期210和心室舒张期212。通常,心室收缩期210在心动周期202和204的QRS波群216的R峰214之后不久开始,并且在在T波218的末端处结束。在收缩期210期间,心室收缩(QRS波群216)并且然然后随后再极化(T波218)。
心室舒张期包括U波220(乳头肌再极化)和P波222(心房收缩)。在心室舒张期212期间,心脏松弛,并且心室再填充循环的血液。R-R间隔224表示相对于心动周期202和204的R峰214的心动周期的持续时间。区域226对应于其中心脏的运动相对于心动周期的其他区域最小的区域,并且已经称为“安静阶段”。通常,这样的区域已经位于进入R-R间隔224的大约40%、70%处。
参考图1和2,特征识别器110被配置为识别由心动周期监测器108生成的心动周期信号的预定特征。例如,特征识别器110被配置为识别具体波(例如,QRS波群216、T波218、U波220、P波222、和/或其他波)的峰。在一个非限制性实例中,特征识别器110被配置为识别R峰214。这能够电压水平阈值化和/或其他方法来实现。一旦被识别,特征识别器110生成指示心动周期的预定特征被识别的触发信号。
参考图1,成像系统100还包括由旋转机架102可旋转地支撑的辐射源112,诸如X射线管。辐射源112与旋转机架102一起在检查区域106周围旋转,并且被配置为生成并发射穿过检查区域106的辐射。成像系统100还包括辐射源控制器114。辐射源控制器114被配置为调制生成的辐射的通量。例如,辐射控制器114能够选择性地改变辐射源112的阴极加热电流,应用电荷以抑制辐射源112的电流,对发射的辐射等进行滤波等以调制通量。在图示的范例中,辐射源控制器114基于预定调制模式来调制通量。
图3示出了结合ECG信号200的范例调制模式300。在图3中,第一或y-轴302表示以毫安(mA)为单位的电流。图示的调制模式300与ECG信号200共享x-(时间)轴208。在该范例中,调制模式300在较高通量水平304与较低(非零)通量水平306之间切换生成的辐射的通量。此外,调制模式300针对对应于安静阶段226的时间窗口308将生成的辐射的通量切换到较高通量水平304,并且在安静阶段226之外切换到较低通量水平306。在其他实施例中,调制模式300能够针对(一个或多个)不同的和/或额外的阶段在多于两个的不同通量水平之间以离散和/或连续方式等进行切换。在变型中,调制模式还可以包括其中管电流并且因此通量为零的一个或多个区域。这能够利用能够被定位在射束路径中或外的网格切换管和/或物理滤波器来实现。
返回到图1,成像系统100还包括辐射敏感探测器像素116的一或二维阵列115。像素116跨检查区域106位于辐射源112对面,探测穿过检查区域106的辐射,并且生成指示其的电信号(投影数据)。在图示的范例中,像素116包括直接转换光子计数探测器像素。利用这样的像素,生成的信号包括具有指示探测到的光子的能量的峰值幅度或峰值高度的电流或电压。直接转换光子计数探测器像素可以包括任何合适的直接转换材料,诸如CdTe、CdZnTe、Si、Ge、GaAs或其他直接转换材料。
成像系统100还包括脉冲整形器118,所述脉冲整形器处理由探测器像素116输出的电信号并且生成指示探测到的光子的能量的脉冲(诸如电压或其他脉冲)。在变型中,由探测器像素116输出的电信号被预处理(例如,被放大等),并且脉冲整形器118处理经预处理的信号。成像系统100还可以包括能够分辨脉冲的能量分辨器120。在图示的范例中,能量分辨器120包括至少一个比较器122,所述比较器比较脉冲的幅度与对应于感兴趣能量的至少一个能量阈值。比较器122产生指示探测到的光子的能量是在阈值之上还是在阈值之下的输出信号。
成像系统100还包括针对每个阈值增大(或减小)计数值的计数器124。例如,当针对具体阈值的比较器122的输出指示脉冲的幅度超过对应的阈值时,针对该阈值的计数值被增大。成像系统100还包括将计数的脉冲分配到对应于不同能量范围的能量分箱的分箱器126。例如,分箱可以被定义用于两个阈值之间的能量范围。在该范例的情况下,分箱器126将导致针对较低阈值而非针对较高阈值的计数的光子分配到被定义用于两个阈值之间的能量范围的分箱。
成像系统100还包括重建器128,所述重建器利用谱和/或非谱算法重建经分箱的数据并生成谱和/或非谱体积图像数据。重建器128能够采用经滤波的反向投影、迭代、和/或其他重建算法。这包括加权和非加权重建算法。如在下面更详细地描述的,范例加权重建算法针对当通量处于较高通量水平304时采集的投影数据在更大程度上对投影数据加权。在另一范例中,加权重建算法针对当通量处于较低通量水平306处时采集的投影数据在更大程度上对投影数据加权。
这样的算法非常适合于利用图3所示的调制模式300和/或控制通量使得其在心动周期的安静阶段期间处于高水平的其他通量调制模式的心脏研究。例如,针对当通量在安静阶段226内处于较高通量水平304时采集的投影数据在更大程度上对投影数据加权的加权重建算法能够用于生成心脏的体积图像数据,并且针对当通量在其他阶段内处于较低通量水平306时采集的投影数据在更大程度上对投影数据加权的加权重建算法能够用于生成肺(并且任选地对象的周边)的体积图像数据。
这能够在下文中如所描述的那样来实现。在扫描之前,管电流调制模式被确定。例如,模式可以在最后的R峰之后使电流坡升300毫秒(300msec),并且随后再次使电流降低另一300msec。该模式假定某一长度的近似恒定的心跳。对于重建,图像在心动周期的70%处被重建。较高通量水平304的持续时间能够长于所需的用于重建的选通窗口宽度。如果这样的话,当通量处于较高通量水平304处时权重不必总是高的。然而,当通量处于较高通量水平304时,权重能够总是高的。
以此方式,心脏的体积图像数据具有用于诊断目的合适的信噪比。此外,肺的体积图像数据将不会包括由辐射敏感像素116的阵列的光子计数探测器像素的有限计数速率引起的伪影。此外,对象计量相对于完全在较高水平的通量304处的扫描被减少。
心脏的体积图像数据和肺的体积图像数据能够被同时视觉地显示在显示屏的不同、交叠、或相同的显示窗口中。此外,心脏的体积图像数据和肺的体积图像数据能够被备选地视觉地显示。图像处理器130能够用于组合心脏的体积图像数据和肺的体积图像数据以形成的单个体积图像数据集。为此,在一个范例中,心脏或肺能够被分割,并且与肺的体积图像数据或心脏的体积图像数据组合。
成像系统100还包括计算系统,其充当操作者控制台132并且包括输出设备(诸如显示器)和输入设备(诸如键盘、鼠标等)。驻留在控制台132上的软件允许用户控制系统100的操作。例如,控制台132能够用于识别成像协议(例如,具有通量调制和加权重建的心脏成像),初始化扫描,并且基于识别的成像协议控制辐射源控制器114和重建器128以调制通量并采用具体的重建算法。
应认识到,心脏监测器108、特征识别器110、脉冲整形器118、分辨器120、比较器122、计数器124、分箱器126、重建器128、图像处理器130和控制台132能够经由电路和/或软件/固件来实施。例如,控制台132包括一个或多个处理器(例如,中央处理单元、微处理器等)和被编码有计算机可执行指令的计算机可读介质(其不包括暂态存储器,而包括物理存储器),当所述计算机可执行指令被一个或多个处理器执行时,所述计算机可执行指令令一个或多个处理器执行本文中描述的功能中的一个或多个。
如在上面简单描述的,重建器128能够采用加权重建算法。以下描述了加权重建算法的范例。通常,每个重建窗口的宽度通过重建的完整性情况来确定,使得每一个体素由源辐射至少一个180度(180°)加上扇形角(或在针对平行锥形射束重分箱的扇形射束之后的180°)。
图4示出了具有结合ECG 200和调制模式300的重建窗口400的范例重建算法。在该范例中,重建窗口400是用于心脏的重建窗口。在该范例中,心脏重建窗口400对应于安静阶段226,并且位于其中辐射源电流在较高通量水平304处被调制的时间段内。在其他实施例中,心脏重建窗口400能够更窄或更宽,包括覆盖当辐射源电流在较低通量水平306处被调制时的时间。
图5、6和7分别示出了具有结合ECG200和调制模式300的重建窗口500、600和700的范例重建算法。在这些图中,重建窗口500、600和700是用于肺的呼吸重建窗口。
在图5中,呼吸重建窗口500覆盖了辐射源电流在较低通量水平306处被调制的整个时间。在图6中,呼吸重建窗口600覆盖辐射源电流在较低通量水平306处被调制的时间的仅一部分。在图7中,呼吸重建窗口700接近心脏重建窗口400(图4),但是只在辐射源电流在较低通量水平306处被调制之处。
因此,在图7的呼吸重建窗口700的情况下,心脏和肺图像的运动状态将会具有极小差异。呼吸重建窗口700可以被定位在心脏重建窗口400之前或之后(如图7所示)。在变型中,呼吸重建窗口700也可以在心脏重建窗口400的部分之前或在心脏重建窗口400的部分之后被分裂。
在Koken的“Aperture weighted cardiac reconstruction for cone-beam CT”(Phys.Med.Biol.51(2006)3433–3448)中描述了非迭代重建的非限制性范例。在Nielsen的“Cardiac cone-beam CT volume reconstruction using ART”(Med.Phys.32(4)(2005)851–860)中描述了迭代重建的非限制性范例。根据Koken,合适的重建能够如等式1中所示被公式化:
等式1:
其中表示目标点,表示目标函数,l和表示指数,表示加权函数,表示预处理的投影数据,φ表示被投影到xy-平面上的平行线的角度,并且和表示目标点的投影的探测器坐标。
对于心脏重建,其中,ωc表示心脏权重,ωap表示孔权重,并且对于肺重建,其中,对于所有{φ|I(φ)=I低},ωl(φ)=1,并且否者为0。
在另一实施例中,对于其中到探测器的X射线通量不超过指示探测器的计数速率限制的预定阈值的所有投影,ωl(φ)=1,并且对于其中到探测器的X射线通量超过指示探测器的计数速率限制的预定阈值的所有投影,ωl(φ)=0。该实施例非常适合于其中过多通量水平能够通过硬件、软件、和/或以另外的方式在探测器本身内被探测到的构造。
在其他实施例中,ωl(φ)包括边界上的cos2平滑、到心脏相位点的最小时间距离、或确保180度加上数据的扇形角度但并非更多的时间宽度。在另一实例中,其中,对所有{(φ,u,h)|Idet(φ,u,h)<Isat},为1,并且否者为0。在该实例中,Isat是其中探测器饱和的强度,并且ωl(φ)是本文中描述的加权函数。
图8图示了根据本文描述的实施例的方法。
应认识到,下面动作的顺序是用于解释性目的,而不是进行限制。这样,其他顺序在本文中也被考虑。此外,动作中的一个或多个可以被省略,和/或一个或多个其他动作可以被包括。
在802处,患者被加载到具有光子计数探测器像素的成像系统的对象支撑物上并且关于成像区域被定位以便进行扫描。
在804处,辐射发射通量在至少第一通量水平与第二通量水平之间被调制,其中,第一通量水平高于第二通量水平,并且所述调制与患者的心动周期的安静阶段协调,使得安静阶段期间的辐射发射通量为较高的通量。
在806处,扫描被执行,生成投影数据。应认识到,扫描能够是用于圆形、螺旋、分步和拍摄(例如,轴向)的选通重建、和/或其他扫描。
在808处,投影数据利用心脏加权重建算法来重建,所述心脏加权重建算法对对应于安静阶段的投影数据的第一部分比对应于安静阶段之外的一个或多个阶段的投影数据的第二部分更高地加权,从而生成心脏体积图像数据。
在810处,投影数据利用肺加权重建算法来重建,所述肺加权重建算法对对应于安静阶段之外的一个或多个阶段的投影数据的第二部分比对应于安静阶段的投影数据的第一部分更高地加权,从而生成呼吸体积图像数据。
在812处,心脏体积图像数据和呼吸体积图像数据如在本文中讨论的那样和/或以另外的方式被显示。
本文中讨论的方法中的至少一部分可以通过被编码或被嵌入在计算机可读介质(其不包括暂态介质)上的计算机可读指令来实施,当所述计算机可读指令被(一个或多个)计算机处理器执行时,令(一个或多个)处理器执行所描述的动作。额外地或备选地,计算机可读指令中的至少一个由信号、载波或其他暂态介质承载。
已经参考优选实施例描述了本发明。他人在阅读和理解以上具体实施方式的情况下能够想到修改或替代。本发明旨在被解释为包括所有这种修改和替代,只要它们落入权利要求书及其等价方案的范围之内。
Claims (21)
1.一种方法,包括:
与心动周期信号协作地在第一通量水平与不同的第二通量水平之间调制发射辐射的通量,使得所述通量在具有第一心脏运动的第一心脏运动阶段期间处于所述第一通量水平处,并且在具有第二心脏运动的第二心脏运动阶段期间处于所述第二通量水平处,其中,所述第一通量水平小于所述第二通量水平,并且其中,所述第一心脏运动大于所述第二心脏运动;
探测经调制的发射辐射;
产生指示探测到的经调制的发射辐射的投影数据;并且
利用第一重建窗口重建所述投影数据,所述第一重建窗口将第一权重应用于与所述第一心脏运动阶段和较低的第一通量水平相对应的所述投影数据的第一子集,并且将不同的第二权重应用于与所述第二心脏运动阶段和较高的第二通量水平相对应的所述投影数据的第二子集,以生成第一体积图像数据,其中,所述第一权重大于所述第二权重。
2.根据权利要求1所述的方法,其中,所述第一权重具有非零值,并且所述第二权重具有零值。
3.根据权利要求1至2中的任一项所述的方法,还包括:
将所述第一权重应用于与所述第一心脏运动阶段相对应的投影数据的整个集合。
4.根据权利要求1至2中的任一项所述的方法,还包括:
将所述第一权重应用于与所述第一心脏运动阶段相对应的所述投影数据的第一子集。
5.根据权利要求3至4中的任一项所述的方法,还包括:
利用第二重建窗口重建所述投影数据,所述第二重建窗口将第三权重应用于与所述第二心脏运动阶段和所述较高的第二通量水平相对应的所述投影数据的所述第二子集,并且将第四权重应用于与所述第一心脏运动阶段和所述较低的第一通量水平相对应的所述投影数据的所述第一子集,其中,所述第三权重大于所述第四权重,以生成第二体积图像数据。
6.根据权利要求1所述的方法,还包括:
响应于辐照所述探测器的X射线通量不超过指示所述探测器的计数速率限制的预定阈值,而将等于一的第一权重应用于所述投影数据;并且
响应于辐照所述探测器的所述X射线通量超过指示所述探测器的所述计数速率限制的所述预定阈值,而将等于零的第二权重应用于所述投影数据。
7.根据权利要求2所述的方法,还包括:
利用第二重建窗口重建所述投影数据,所述第二重建窗口将第三权重应用于与所述第二心脏运动阶段和所述较高的第二通量水平相对应的所述投影数据的所述第二子集,并且将第四权重应用于与所述第一心脏运动阶段和所述较低的第一通量水平相对应的所述投影数据的所述第一子集,其中,所述第三权重大于所述第四权重,以生成第二体积图像数据。
8.根据权利要求7所述的方法,还包括:
将所述第一权重应用于与邻近于所述第二心脏运动阶段的一个或多个心脏运动阶段相对应的投影数据。
9.根据权利要求8所述的方法,其中,所述一个或多个心脏运动阶段在所述第二心脏运动阶段之前发生。
10.根据权利要求8所述的方法,其中,所述一个或多个心脏运动阶段在所述第二心脏运动阶段之后发生。
11.根据权利要求8所述的方法,其中,所述一个或多个心脏运动阶段的第一子集在所述第二心脏运动阶段之前发生,并且所述一个或多个心脏运动阶段的第二子集在所述第二心脏运动阶段之后发生。
12.一种成像系统(100),包括:
辐射源(112),其被配置为关于检查区域(106)旋转并且发射穿过所述检查区域的辐射;
辐射源控制器(114),其被配置为控制所述辐射源以与心动周期信号协作地在第一通量水平与不同的第二通量水平之间调制发射辐射的通量,使得所述通量在具有第一心脏运动的第一心脏运动阶段期间处于所述第一通量水平处,并且在具有第二心脏运动的第二心脏运动阶段期间处于所述第二通量水平处,其中,所述第一通量水平小于所述第二通量水平,并且其中,所述第一心脏运动大于所述第二心脏运动;以及
辐射敏感像素(116)的阵列(115),其被配置为探测穿过所述检查区域的辐射并且生成指示探测到的辐射的信号;以及
重建器(128),其被配置为利用第一重建窗口重建所述投影数据,所述第一重建窗口将第一权重应用于与所述第一心脏运动阶段和较低的第一通量水平相对应的所述投影数据的第一子集,并且将不同的第二权重应用于与所述第二心脏运动阶段和较高的第二通量水平相对应的所述投影数据的第二子集,以生成第一体积图像数据,其中,所述第一权重大于所述第二权重。
13.根据权利要求12所述的成像系统,其中,所述第一权重具有非零值,并且所述第二权重为零。
14.根据权利要求12至13中的任一项所述的成像系统,其中,所述第一权重被应用于与所述第一心脏运动阶段相对应的投影数据的整个集合。
15.根据权利要求12至13中的任一项所述的成像系统,其中,所述第一权重被应用于与所述第一心脏运动阶段相对应的所述投影数据的第一子集。
16.根据权利要求13至14中的任一项所述的成像系统,其中,所述重建器还被配置为利用第二重建窗口重建所述投影数据,所述第二重建窗口将第三权重应用于与所述第二心脏运动阶段和所述较高的第二通量水平相对应的所述投影数据的所述第二子集,并且将第四权重应用于与所述第一心脏运动阶段和所述较低的第一通量水平相对应的所述投影数据的所述第一子集,其中,所述第三权重大于所述第四权重,以生成第二体积图像数据。
17.根据权利要求13所述的成像系统,其中,所述重建器还被配置为利用第二重建窗口重建所述投影数据,所述第二重建窗口将第三权重应用于与所述第二心脏运动阶段相对应的所述投影数据的所述第二子集,并且将第四权重应用于与所述第一心脏运动阶段相对应的所述投影数据的所述第一子集,其中,所述第三权重大于所述第四权重,以生成第二体积图像数据。
18.根据权利要求17所述的成像系统,其中,所述第一权重被应用于与邻近于所述第二心脏运动阶段的一个或多个心脏运动阶段相对应的投影数据。
19.根据权利要求18所述的成像系统,其中,所述一个或多个心脏运动阶段在所述第二心脏运动阶段之前或之后中的至少一个处发生。
20.根据权利要求18所述的成像系统,其中,所述一个或多个心脏运动阶段的第一子集在所述第一心脏运动阶段之前发生,并且所述一个或多个心脏运动阶段的第二子集在所述第一心脏运动阶段之后发生。
21.一种被编码有计算机可读指令的计算机可读存储介质,所述计算机可读指令在由处理器运行时令所述处理器:
与心动周期信号协作地在第一通量水平与不同的第二通量水平之间调制发射辐射的通量,使得所述通量在具有第一心脏运动的第一心脏运动阶段期间处于所述第一通量水平处,并且在具有第二心脏运动的第二心脏运动阶段期间处于所述第二通量水平处,其中,所述第一通量水平小于所述第二通量水平,并且其中,所述第一心脏运动大于所述第二心脏运动;
接收由探测经调制的发射辐射的探测器生成的投影数据;并且
利用第一重建窗口重建所述投影数据,所述第一重建窗口将第一权重应用于与所述第一心脏运动阶段和较低的第一通量水平相对应的所述投影数据的第一子集,并且将不同的第二权重应用于与所述第二心脏运动阶段和较高的第二通量水平相对应的所述投影数据的第二子集,以生成第一体积图像数据,其中,所述第一权重大于所述第二权重。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201462065939P | 2014-10-20 | 2014-10-20 | |
US62/065,939 | 2014-10-20 | ||
PCT/IB2015/057791 WO2016063171A1 (en) | 2014-10-20 | 2015-10-12 | Cardiac reconstruction for photon counting ct for heart and lung images |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107072627A true CN107072627A (zh) | 2017-08-18 |
CN107072627B CN107072627B (zh) | 2020-10-16 |
Family
ID=54396923
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201580056917.6A Active CN107072627B (zh) | 2014-10-20 | 2015-10-12 | 用于针对心脏和肺图像的光子计数ct的心脏重建 |
Country Status (5)
Country | Link |
---|---|
US (1) | US9775580B2 (zh) |
EP (1) | EP3209214B1 (zh) |
JP (1) | JP6243580B2 (zh) |
CN (1) | CN107072627B (zh) |
WO (1) | WO2016063171A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108065952A (zh) * | 2018-01-03 | 2018-05-25 | 沈阳东软医疗系统有限公司 | 一种剂量调制方法及装置 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11200709B2 (en) * | 2016-12-27 | 2021-12-14 | Canon Medical Systems Corporation | Radiation image diagnostic apparatus and medical image processing apparatus |
JP7066359B2 (ja) * | 2017-09-19 | 2022-05-13 | キヤノンメディカルシステムズ株式会社 | X線ct装置 |
WO2019149762A1 (en) | 2018-01-31 | 2019-08-08 | Koninklijke Philips N.V. | Non-spectral computed tomography (ct) scanner configured to generate spectral volumetric image data |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101495886A (zh) * | 2006-08-01 | 2009-07-29 | 皇家飞利浦电子股份有限公司 | 立体管计算机断层摄影 |
WO2010070583A1 (en) * | 2008-12-17 | 2010-06-24 | Koninklijke Philips Electronics N.V. | X-ray examination device and method |
US20100166277A1 (en) * | 2008-12-30 | 2010-07-01 | Siemens Aktiengesellschaft | Method for noise reduction of CT images and image processing system |
JP2011056167A (ja) * | 2009-09-14 | 2011-03-24 | Toshiba Corp | 放射線画像撮影装置及び放射線画像撮影方法 |
EP2452626A1 (en) * | 2005-09-07 | 2012-05-16 | Kabushiki Kaisha Toshiba | X-Ray computed tomography apparatus |
WO2012104751A1 (en) * | 2011-01-31 | 2012-08-09 | Koninklijke Philips Electronics N.V. | Detection values processing apparatus |
JP2013522657A (ja) * | 2010-03-11 | 2013-06-13 | フラウンホーファーゲゼルシャフト ツール フォルデルング デル アンゲヴァンテン フォルシユング エー.フアー. | 信号処理器、窓提供部、符号化されたメディア信号、信号を処理するための方法および窓を提供するための方法 |
CN103366389A (zh) * | 2013-04-27 | 2013-10-23 | 中国人民解放军北京军区总医院 | Ct图像重建方法 |
CN103959329A (zh) * | 2011-11-23 | 2014-07-30 | 皇家飞利浦有限公司 | 图像域去噪 |
US20140236488A1 (en) * | 2013-02-19 | 2014-08-21 | Mindways Software, Inc. | Computed Tomography Calibration Systems and Methods |
CN104105445A (zh) * | 2012-08-30 | 2014-10-15 | 株式会社东芝 | X射线ct装置、图像处理装置以及图像处理方法 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6931094B2 (en) * | 2003-10-10 | 2005-08-16 | Ge Medical Systems Global Technology Company, Llc | Methods and systems for smoothing |
JP2006101926A (ja) * | 2004-09-30 | 2006-04-20 | M & C:Kk | 放射線検出装置、放射線画像診断装置、及び放射線画像の生成方法 |
WO2006067671A2 (en) * | 2004-12-22 | 2006-06-29 | Philips Intellectual Property & Standards Gmbh | Method and apparatus for cardiac computed tomography |
JP4945203B2 (ja) * | 2005-09-07 | 2012-06-06 | 株式会社東芝 | X線コンピュータ断層撮影装置 |
US7532702B2 (en) | 2005-11-23 | 2009-05-12 | General Electric Company | Method and system for performing CT image reconstruction with motion artifact correction |
US20070147579A1 (en) * | 2005-12-23 | 2007-06-28 | De Man Bruno K B | Method and system for radiographic imaging with organ-based radiation profile prescription |
JP2011067333A (ja) * | 2009-09-25 | 2011-04-07 | Fujifilm Corp | 放射線画像撮影装置及び撮影制御装置 |
JP2015516832A (ja) | 2012-03-27 | 2015-06-18 | コーニンクレッカ フィリップス エヌ ヴェ | 光子計数検出器を備える撮像システムによる従来型のイメージング |
US9111337B2 (en) * | 2012-10-12 | 2015-08-18 | Mayo Foundation For Medical Education And Research | Low dose cardiac CT imaging with time-adaptive filtration |
JP6329490B2 (ja) * | 2013-02-05 | 2018-05-23 | 株式会社日立製作所 | X線ct装置及び画像再構成方法 |
KR102205353B1 (ko) * | 2013-09-26 | 2021-01-20 | 삼성전자주식회사 | 엑스선 촬영 장치 및 그 제어 방법 |
WO2015108097A1 (ja) * | 2014-01-20 | 2015-07-23 | 株式会社 日立メディコ | X線ct装置、画像処理装置、及び画像再構成方法 |
-
2015
- 2015-10-12 JP JP2017519490A patent/JP6243580B2/ja active Active
- 2015-10-12 EP EP15790279.2A patent/EP3209214B1/en active Active
- 2015-10-12 CN CN201580056917.6A patent/CN107072627B/zh active Active
- 2015-10-12 US US15/519,176 patent/US9775580B2/en active Active
- 2015-10-12 WO PCT/IB2015/057791 patent/WO2016063171A1/en active Application Filing
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2452626A1 (en) * | 2005-09-07 | 2012-05-16 | Kabushiki Kaisha Toshiba | X-Ray computed tomography apparatus |
CN101495886A (zh) * | 2006-08-01 | 2009-07-29 | 皇家飞利浦电子股份有限公司 | 立体管计算机断层摄影 |
WO2010070583A1 (en) * | 2008-12-17 | 2010-06-24 | Koninklijke Philips Electronics N.V. | X-ray examination device and method |
US20100166277A1 (en) * | 2008-12-30 | 2010-07-01 | Siemens Aktiengesellschaft | Method for noise reduction of CT images and image processing system |
JP2011056167A (ja) * | 2009-09-14 | 2011-03-24 | Toshiba Corp | 放射線画像撮影装置及び放射線画像撮影方法 |
JP2013522657A (ja) * | 2010-03-11 | 2013-06-13 | フラウンホーファーゲゼルシャフト ツール フォルデルング デル アンゲヴァンテン フォルシユング エー.フアー. | 信号処理器、窓提供部、符号化されたメディア信号、信号を処理するための方法および窓を提供するための方法 |
WO2012104751A1 (en) * | 2011-01-31 | 2012-08-09 | Koninklijke Philips Electronics N.V. | Detection values processing apparatus |
CN103959329A (zh) * | 2011-11-23 | 2014-07-30 | 皇家飞利浦有限公司 | 图像域去噪 |
CN104105445A (zh) * | 2012-08-30 | 2014-10-15 | 株式会社东芝 | X射线ct装置、图像处理装置以及图像处理方法 |
US20140236488A1 (en) * | 2013-02-19 | 2014-08-21 | Mindways Software, Inc. | Computed Tomography Calibration Systems and Methods |
CN103366389A (zh) * | 2013-04-27 | 2013-10-23 | 中国人民解放军北京军区总医院 | Ct图像重建方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108065952A (zh) * | 2018-01-03 | 2018-05-25 | 沈阳东软医疗系统有限公司 | 一种剂量调制方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
EP3209214A1 (en) | 2017-08-30 |
JP2017529978A (ja) | 2017-10-12 |
EP3209214B1 (en) | 2019-06-05 |
JP6243580B2 (ja) | 2017-12-06 |
US9775580B2 (en) | 2017-10-03 |
CN107072627B (zh) | 2020-10-16 |
WO2016063171A1 (en) | 2016-04-28 |
US20170231590A1 (en) | 2017-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9848843B2 (en) | Dynamic perfusion imaging | |
EP1522045B1 (en) | Motion artifact correction of tomographical images | |
JP5551436B2 (ja) | ステレオ管コンピュータ断層撮像 | |
JP6021311B2 (ja) | X線コンピュータ断層撮影装置 | |
US10420519B2 (en) | Computed tomography system | |
CN107072627A (zh) | 用于针对心脏和肺图像的光子计数ct的心脏重建 | |
CN110536640B (zh) | 从pet列表数据中的呼吸运动信号的噪声鲁棒的实时提取 | |
EP2512345A1 (en) | Computed tomography apparatus | |
KR102205353B1 (ko) | 엑스선 촬영 장치 및 그 제어 방법 | |
KR20170105876A (ko) | 단층 촬영 장치 및 그에 따른 단층 영상 재구성 방법 | |
US10984564B2 (en) | Image noise estimation using alternating negation | |
US20210386392A1 (en) | Systems and methods for four-dimensional ct scan | |
WO2018087049A1 (en) | Dose reduction in dynamic radiography | |
US20220180576A1 (en) | Cone beam artifact correction for gated imaging | |
CN102258383A (zh) | X射线计算机断层摄影装置以及图像处理方法 | |
US20230320687A1 (en) | Controller for a ct scanner | |
US20160203621A1 (en) | Medical diagnostic imaging apparatus, image processing apparatus and image processing method | |
Koken et al. | Aperture weighted cardiac cone-beam reconstruction using retrospective ECG gating | |
Zamyatin et al. | Implementation of a circle and helix reconstruction algorithm for 256-slice CT | |
Manzke et al. | Image reconstruction with multi-slice CT |
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 |