CN101454801B - 基于列表模式数据的局部运动补偿 - Google Patents
基于列表模式数据的局部运动补偿 Download PDFInfo
- Publication number
- CN101454801B CN101454801B CN2007800070638A CN200780007063A CN101454801B CN 101454801 B CN101454801 B CN 101454801B CN 2007800070638 A CN2007800070638 A CN 2007800070638A CN 200780007063 A CN200780007063 A CN 200780007063A CN 101454801 B CN101454801 B CN 101454801B
- Authority
- CN
- China
- Prior art keywords
- interest
- area
- barycenter
- image
- line
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/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 for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/037—Emission tomography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/412—Dynamic
Abstract
一种在由成像系统(10)的重建系统(42)利用原始数据重建的图像中局部校正运动的方法,包括:估计从所述原始数据重建的图像内部的感兴趣区域的特有特征;利用所估计的区域特有特征校正与所述感兴趣区域相关联的原始数据;以及利用校正的原始数据重建与所述感兴趣区域对应的运动校正图像。
Description
说明书
下文涉及医学成像系统。其特别适用于PET成像中的运动补偿,但也适用于诊断成像的其他模态。
在医学成像(例如正电子发射断层摄影(PET)、单光子发射计算机断层摄影(SPECT)、计算机断层摄影(CT)、磁共振成像(MR)等)中公知的一个问题是解剖体的运动。随着患者的呼吸、心动周期等,解剖体的肺部和其他部分会发生周期性运动。在数据采集期间的这种运动造成运动平均的(例如模糊的)图像,其往往具有低的图像质量和有限的定量能力(例如,对比度降低、损伤可探测性降低、标准摄取值(SUV)的定量变差等)。
一种用于减小运动伪影的技术是让患者在数据采集期间屏住呼吸。然而,在一些情况下,采集时间太长,患者无法在整个数据采集期间都屏住呼吸。
在门控采集的情况下,根据患者当前的呼吸状态将事件存储在特定的时间元中,事件是利用外部传感器测量的或直接从采集的数据提取的。然后对独立的元或门独立进行重建,以给出具有较少运动伪影的图像。然而,这些图像的计数统计的降低导致噪声增大且妨碍精确的定量。此外,很多临床应用不要求在整个图像上进行运动补偿。例如,在肿瘤学研究中,临床医生常常仅对SUV和肺损伤大小的定量感兴趣。
常规的主动运动校正方案(例如变形模型)可能非常复杂且耗时,而且容易出错。此外,它们通常需要有关患者呼吸状态随时间变化的额外输入,这需要利用外部传感器来采集。
一方面,示出了一种利用原始数据在由成像系统的重建系统重建的图像中校正运动的方法。该方法包括根据所述列表模式数据估计重建图像内部的感兴趣区域的特有特征。然后利用所估计的区域特有特征对于与所述 感兴趣区域相关联的原始数据进行运动校正。然后对与具有所校正的原始数据的所述感兴趣区域对应的运动校正图像进行重建。
一个优点包括便于在图像中的运动伪影的校正。
另一个优点在于使用用于重建图像的数据来在该图像中校正运动。
另一个优点在于不需要身体运动探测器和复杂的运动校正算法。
在阅读并理解优选实施例的详细描述后,本领域的普通技术人员将明了其他优点。
附图仅用于示出实施例,不应视为限制权利要求。
图1示出了一种医学成像系统,其具有利用用于重建图像的数据校正运动的部件;
图2示出了一种校正图像内运动的示范性技术,其中使LOR组移位相应的运动校正偏移量,然后进行重建以生成运动校正图像;
图3示出了从VOI/ROI导出的多个示范性中间图像的图示范例,其用于确定运动校正偏移量;
图4示出了便于生成中间图像的示范性技术;
图5图示了有关每幅时间部分图像(time partial image)的质心位置的信息,其中可以对其间的值进行内插以生成平滑的随时间变化的平移偏移量;
图6示出了用于校正图像内运动的另一示范性技术,其中使用矢量分析来校正LOR;
图7示出了用于确定质心估计值的示范性矢量方法;
图8示出了运动球体的模拟列表模式数据的示范性结果。
图1示出了一种医学成像系统10,其包括扫描系统12和控制台14,控制台14控制扫描系统12,呈现重建图像并提供显象和数据/图像校正工具。扫描系统12可以是任何类型的医学成像扫描仪,包括正电子发射断层摄影(PET)扫描仪、单光子发射计算机断层摄影(SPECT)扫描仪、计算机断层摄影(CT)扫描仪、磁共振成像(MR)扫描仪、PET/CT扫描仪等中的一种或多种。对于该例而言,扫描系统12为PET扫描仪。同样,扫描 系统12包括多个均设置在成像区域18周围的辐射探测器模块16(例如几百个、几千个等)。定位或旋转模块16以沿着轴向方向界定一个或多个辐射探测器环(例如两个、十个、一百个等),该探测器探测成像区域18内发生的辐射事件(例如γ射线)。扫描系统12还包括支撑机构20,用于在成像之前、期间和/或之后将受检者(例如人)定位在成像区域18中。支撑机构20可以沿着大致横过辐射探测器16的轴向方向直线移动。
控制台14为基于微处理器的系统(例如主机、工作站、台式机、笔记本等)。同样,控制台14包括至少一个或多个处理器22、存储部件24(例如易失性和非易失性、驻留和可移动存储器等)、呈现部件26(例如平板监视器、液晶显示器(LCD)、阴极射线管(CRT)监视器等)、输入装置28(例如鼠标、键盘、小键盘、滚动球、数字笔、扬声器、触摸屏、磁带驱动器、磁盘驱动器等)、输入驱动器30(例如DVD、CD、磁盘、光盘、磁带等)、有线和/或无线通信部件32(例如Ethemet、USB、串行、并行、FireWire等)以及各种计算机相关部件。
控制台14用于规划患者过程(例如选择成像协议,适当设置成像参数等),利用扫描系统12开始扫描,呈现重建图像并提供各种显象和数据/图像校正能力,以及其他功能,例如终止过程、传输数据、定位支撑机构20等。适当的显象能力范例包括但不限于从二维数据集形成三维体积,界定图元(例如感兴趣的区域和/或体积),测量各种图像质量统计数据,叠加通过不同模态获得的图像,从图像中去除感兴趣解剖结构和/或界定的区域/体积等。如下文详细讨论的,利用这些显象能力便于通过由一个或多个处理器22执行的一个或多个运动估计/校正算法从图像中的多个区域(例如一个或多个感兴趣区域或体积等)去除运动伪影。这些算法使用用于重建图像的数据(例如列表模式数据)来在重建图像中校正运动。
在准备利用扫描系统12成像时,给受检者施予放射性药物,并将受检者适当放置在成像区域18之内。放射性药物经历放射性衰变,产生正电子发射。正电子与一个或多个附近的电子交互作用并湮灭,产生两个方向相反、每个能量均约为511keV的γ射线。两束方向相反的γ射线基本在同一时间同时撞击到相对的探测器模块16上。由于距一对探测器不等距的位置经过不同的距离,因此在符合事件之间有一时间偏移量。
典型的探测器模块16包括一个或多个闪炼晶体(未示出)。每个闪炼晶体在被诸如正电子湮灭产生的γ射线的放射事件撞击时会产生闪光。每个晶体产生的光被一个或多个诸如光电倍增管的光探测器(未示出)接收。每个光探测器将光转换成代表性电信号。也可以考虑使用直接响应于接收辐射产生电信号的固态探测器和其他类型的探测系统。由一个或多个处理部件34传输并处理由一个或多个光探测器所得到的电信号。处理部件34进行的适当处理包括信号放大、滤波和/或调整。然后将处理过的电信号传送到转换器36,其对信号进行数字化并打上时间戳。可以使用缓冲器和/或其他存储介质(未示出)以便于信号的预处理和数字化。
然后将数据传送给符合探测器38,其识别对应于电子正电子湮灭事件对的基本同时的符合γ射线探测对。例如,该处理可以包括能量滤波(例如,抛弃设置在511keV附近的选定能量滤波窗口之外的辐射探测事件)和/或符合性滤波(例如,抛弃彼此暂时分隔大于选定时间滤波间隔的辐射探测事件对或公共窗口内的两个以上事件、利用延迟窗口技术的随机符合性识别等)。在识别出符合事件对时,响应线(LOR)处理器40处理每对事件的空间信息,以识别在每个事件对中连接γ射线探测的空间LOR。由于正电子电子湮灭事件发射的两束γ射线是空间上相反方向的,因此知道电子正电子湮灭事件发生在LOR上的某处。在飞行时间(TOF)PET中,处理器40分析γ射线对的探测时间之间的差异,以沿着LOR确定发生湮灭事件的段。
将空间LOR和/或符合事件对传送到与至少便于LOR的运动校正的简化的重建系统43(如这里所述,其可以是伪影校正器42的一部分或如下文图2所述在伪影校正器42外部)相关的伪影校正器42(下文详细描述)、重建系统44、存储部件46和/或显象系统48。重建系统44使用任何适当的重建算法从该数据重建一个或多个图像。适当算法的范例包括但不限于有校正或无校正的滤波反向投影和/或迭代反向投影。由于运动校正可能会给感兴趣区域之外的图像区域带来伪影,因此在一个实施例中,图像组合器50利用来自运动校正图像的ROI替换未校正图像的ROI。
所得到的重建图像和/或原始数据可以被传送到控制台14并存储在一个或多个存储部件24中和/或通过一个或多个呈现部件26呈现给用户。此 外后者替代地,可以将所得到的重建图像传送给存储部件46。在这一情况下,控制台14可以从存储部件46获取重建图像和/或原始数据。此外或者替代地,可以从控制台14、重建系统D和/或存储部件46将重建图像和/或原始数据传送到显象系统48。控制台14和显象系统48都可以将数据和/或图像发送给一个或多个外围装置52,包括胶片、打印机、绘图仪、传真机、电子邮箱地址、计算机、网络等。
如上文所简述的,控制台14提供了各种显象和数据/图像校正能力。于是,当在扫描过程期间和/或之后把图像发送到控制台14时,临床医生可以利用显象和数据/图像校正能力观察和/或校正图像伪影(例如运动伪影)。这种能力可以是存储在存储部件24中并由处理器22执行的系统软件和/或应用程序的一部分。此外或替代地,显象系统48可以便于实时和/或后期处理伪影校正。显象系统48可以是具有专用硬件(未示出)和/或便于复杂且耗时的图像处理算法的专用软件应用程序的专用处理系统,专用硬件例如为图形引擎、图形加速器、更大量的存储器、更高能力的处理器、高分辨率监视器等。在备选实施例中(未示出),控制台14和显象系统48被组合到单个单元中,以提供扫描仪控制以及高级图像处理能力。
接下来描述用于校正运动伪影的适当技术的范例。然而,要理解的是,也可以考虑采用其他技术以及这里所述的技术的导出技术。此外,可以同时重建运动伪影图像和运动校正图像,或者可以首先重建运动伪影图像,随后重建对应于一个或多个ROI/VOI的运动校正图像。第一种技术(在下面的图2中给出)使用简化的基于Siddon的算法,另一种技术(在下面的图6中给出)使用矢量空间方法。
开始参考图2,LOR处理器40将剩余的LOR传送到伪影校正器42的分组器52。例如,分组器52将LOR分成跨越大约100毫秒的组。在一种情况下,时间跨度可以是固定的,每组的LOR数量可以变化。在另一种情况下,例如可以将LOR的数量固定为大约10000个,时间跨度可以变化。在又一种情况下,在约束条件之内,数量和时间跨度二者都可以有阈值限制,或者允许二者变化。滤波器54滤除不与用户或者自动定义的感兴趣体积/感兴趣区域(VOL/ROT)56相交叉的LOR。对于TOF PET而言,还消除了与感兴趣体积相交叉但其识别湮灭事件的可能位置的段落在 VOT/ROL56之外的LOR。
特征处理器58确定每个组的LOR的特有特征,例如质心或其他中心函数。在该范例中,特征处理器58的控制器60访问简化重建系统43以将LOR的每个组重建成部分图像。图3示出了VOI/ROI56和由简化重建系统43从VOI/ROI56导出的多个示范性中间图像62的图示范例。示出了十个中间图像62;不过,在其他实施例中,可以生成更多的中间图像。
在一种简化形式中,不计算每个LOR与每个体素的确切重叠长度。相反,针对每个中间图像62,对于与LOR相交叉的每个体素将VOI/ROI56内的体素值加一。这一点在图4中示出,其中LOR64与VOI/ROI56内的体素66、68、70、72、74和76相交叉。如图所示,将体素66-76中的每个的值增加1。该简化的反向投影在计算上是高效的,因为所需的除法和乘法次数是最小的。生成足够数量的中间图像62的时间典型地为几秒左右。
回到图2,在存储部件78中存储该部分图像。质心(COM)发生器80计算存储部件78中存储的每个部分图像的质心或其他特有特征。来自与碰撞VOI/ROI56的LOR相关联的VOI/ROI56外部的符合性对中间图像62有贡献。于是,上述局部受约束的反向投影可以具有与VOI/ROI56外部的事件相比明显的背景。为了使背景的影响最小化,给数据设定阈值,从而在计算质心期间仅考虑值在最大ROI/VOI值的50%(或其他临床医生认为适当的百分比)左右以上的体素。
变异处理器82确定LOR当前组的质心或其他特有特征与存储器84中存储的预先选定的质心或其他特有特征之间的平移偏移量或其他偏差。可以将预先选定的中心设置为第一组被处理的运动伪影图像等的质心或者可以由操作员选定。将每个质心的偏差存储在矢量、阵列、串、矩阵等中。图5图示了有关每幅时间部分图像的质心位置的信息。整个时间范围上的总质心偏差被用于计算时间相关的平移矢量86。
返回到图2,LOR移位器88通过将原始数据集中的每个LOR或每个LOR组移位由偏差给出的量(例如平移矢量86)以得到运动校正的列表模式数据集来校正运动列表模式数据。与仍在每个时间帧的持续时间期间上对数据求平均或模糊化的常规门控技术不同,这是以列表模式数据集的完整时间分辨率进行的。还可以将该校正集成为实时图像重建的一部分。假定可以忽略ROI/VOI内部组织的形变,使得质心运动反映ROI/VOI内部组织的运动。对于小VOI而言这种假定是合理的。虽然基于平移做了说明,但也可以考虑其他连接,例如旋转、弹性形变等。
然后由重建系统44重建移位后或校正偏差后的LOR以及来自VOI/ROI外部的平移后或未平移的LOR。因为LOR都被移位到目标的质心,所以感兴趣区域的重建图像得到了运动校正。在一种情况下,由控制台14的呈现部件26在不同的显示窗口中呈现所得到的校正图像和未校正图像。在另一种情况下,在从扫描仪读出数据时,基本实时地对数据进行运动校正和重建。由于每个额外的LOR组变成所显示图像的一部分,因此显示图像逐渐变得清晰起来。在又一种情况下,将校正图像的一部分拼接到更大的未校正图像的对应区域中。由控制台14的呈现部件26呈现各图像、图像组合和/或其导出结果。在又一种情况下,将运动校正图像提供给存储部件46和/或显象系统48。
或者,在图6所示的实施例中,由特征处理器56利用矢量分析部件96确定LOR的质心。如上所述,LOR处理器40将LOR提供给伪影校正器42。伪影校正器42的分组器52将LOR分成一个或多个组。滤波器54滤除VOI/ROI 56之外的LOR,对于TOF PET而言,滤除与感兴趣体积相交叉但其识别湮灭事件的可能位置的段落在VOI/ROI 56之外的LOR。
滤波后的LOR被提供给矢量分析部件96。对于该范例而言,假定列表模式发射数据中的每个LOR由LOR上的点 和指向LOR方向的单位矢量 定义。图7示出了两个这种LOR,其中第一LOR 98与点 和指向LOR98的方向的单位矢量 相关,第二LOR 100与点 和指向LOR 100的方向的单位矢量 相关。
其中d为连接LOR 98和100的线段长度,被定义为:
其中n为指向连接LOR 98和100的线段的方向的单位矢量,被定义为:
且b被定义为:
其中上标x和y表示对应矢量的x和y分量。
返回到图6,利用该体系,矢量分析部件96从列表模式数据生成用于给定数量的n个连续LOR对的点的集合{C1,…Cn}。质心发生器80使用该信息确定该点云的质心,这是对指定ROI/VOI中活动分布的质心的估计值。通过去除点云中的异常值(outlier)可以改进上述技术。分别处理平行的LOR或垂直于x轴(即 )的LOR。图8示出了运动球体的模拟列表模式数据的示范性结果。在图104中以图形方式给出了结果,该图104在第一轴106上示出了位置,在第二轴108上示出了时间。对于该范例而言,针对以±2cm幅度运动的活动球体(r=1cm)模拟列表模式数据。在110处指定的事件数量用于直接从列表模式数据确定球体随时间而变化的位置。
返回到图7,变异处理器82确定当前LOR组和存储器84中存储的预先选定的质心之间的适当平移函数或其他偏差。LOR移位器88利用该偏移量或其他偏差调整LOR。然后由重建系统44重建调整后或校正偏差后的LOR以及来自VOI/ROI外部的平移后或未平移的LOR。类似地,由于LOR都被移位到目标的质心,所以感兴趣区域的重建图像得到了运动校正。然后将运动校正图像提供给控制台14、存储部件46和/或显象部件48。
在一种应用中,临床医生对定量值(例如,特定肺损伤的定量值)感兴趣。在这一情况下,临床医生通过在一个或多个呈现部件26呈现的图像上(利用一个或多个输入装置28)划出几何图元来界定感兴趣区域或体积。例如,临床医生可以画出包括了受检者损伤的基本部分的图元。根据临床医生的偏好和/或损伤的属性(例如大小、形状、相对位置等),该图元可以是正方形、矩形、圆形、椭圆形、随手画成的形状、球体、立方体等,和/或任何其他控制台14和/或显象系统46给出的用于感兴趣区域定义的二维或三维形状。
或者,可以利用算法自动定义一个或多个这种感兴趣区域或体积。例如,基于对应于损伤灰度级值的值的范围(例如由临床医生定义或预先定义的),执行软件可以为图像像素(或体素)值设定阈值并画出包括该范围 内的数据的区域。这可能会在图像中产生“热点”自动绘制或较高辐射发射区域。类似地,该图元的形状可以基于已知的几何形状或不规则形状。在使用多种成像模态时,可以在来自一种模态的数据产生的图像上定义感兴趣区域或体积,并将该图像记录到从不同模态的数据产生的图像。例如,在组合PET/CT的情况下,临床医生可以在CT图像中定义感兴趣区域并随后将CT图像记录到PET图像。
上述技术的扩展包括如下内容。如果可以得到飞行时间(TOF)信息,可以利用该信息来改善运动估计。此外,可以通过图形用户接口(GUI)调整时间步长的数值,即运动估计的时间分辨率。此外,可以经由GUI调整用于估计质心的事件数量。或者,可以由算法确定这些参数的最佳设置。此外,上述技术可用于包括肺损伤定量或肺肿瘤分期的各种应用或其他在呼吸、心脏或患者运动的情况下定量局部目标(例如,大小和绝对放射性浓度)的应用。
已经参考优选实施例描述了本发明。在阅读和理解前述详细说明的前提下,其他人可以想到各种修改和变化。只要修改和变化落入所附权利要求或其等价要件的范围内,本发明意在被视为包括所有这种修改和变化。
Claims (38)
1.一种在根据原始数据重建的图像中校正运动的方法,包括:
根据所述原始数据估计所述重建图像内部的感兴趣区域的质心;
基于所述感兴趣区域的所述质心的变化计算运动校正平移矢量(86);
利用所述运动校正平移矢量对于与所述感兴趣区域相关联的所述原始数据进行运动校正;以及
利用所校正的原始数据重建与所述感兴趣区域对应的运动校正图像。
2.根据权利要求1所述的方法,其中,所述原始数据为列表模式数据。
3.根据权利要求1所述的方法,其中,所述感兴趣区域小于视场并定位至感兴趣的解剖结构。
4.根据权利要求1所述的方法,还包括:
利用所述运动校正的感兴趣区域替换所述图像内的未校正的感兴趣区域。
5.根据权利要求1所述的方法,还包括在分开的显示中同时呈现所述运动校正的感兴趣区域和未校正的图像。
6.根据权利要求1所述的方法,还包括与估计所述质心和校正所述原始数据的步骤相关联的如下步骤:
针对多个短时间帧从所述原始数据的段生成多个中间感兴趣区域图像(62);
针对所述多个中间感兴趣区域图像(62)中的每个,调整一个或多个位于感兴趣区域(56)内部并与对应的响应线(64)相交叉的体素(66-76)的体素值;
计算每个所述中间感兴趣区域图像(62)的质心;
计算每个质心和选定质心之间的差异;以及
基于所述差异计算运动校正平移矢量(86)。
7.根据权利要求6所述的方法,还包括:
在估计所述质心之前去除体素值小于感兴趣区域最大体素值的预定百分比的体素,以减小来自所述感兴趣区域(56)外部的符合的影响。
8.根据权利要求6所述的方法,还包括:
通过将初始数据中的响应线移位所述运动校正平移矢量给出的一个量而对于所述原始数据进行运动校正。
9.根据权利要求1所述的方法,还包括与估计所述质心和校正所述原始数据的步骤相关的如下步骤:
计算表示连接所述原始数据中的一对响应线的最短线段的中点的矢量。
10.根据权利要求9所述的方法,其中,所述矢量包括给定数量的连续对的响应线的点的集合。
12.根据权利要求9所述的方法,还包括计算所述感兴趣区域中的分布的质心。
13.根据权利要求12所述的方法,还包括:
通过根据所述质心的运动移位所述初始数据中的响应线来对于所述原始数据进行运动校正。
14.根据权利要求9所述的方法,还包括通过忽略所述感兴趣区域之外的数据来改善所述结果。
15.根据权利要求1所述的方法,其中,校正所述原始数据包括:
将所述原始数据分成组;
根据每个组中的数据确定所述质心在每个组中的位置;
确定所述质心在每个组中的位置和所估计的质心之间的偏差;
根据所确定的偏差调整所述原始数据。
16.根据权利要求15所述的方法,其中,所述原始数据包括多个响应线(LOR),且所述调整步骤包括:
根据所确定的偏差平移所述响应线以对准所述组的所述质心。
17.根据权利要求1所述的方法,还包括如下操作中的至少一种以改进所述运动校正:
使用飞行时间信息以便于估计运动;
调整时间分辨率以影响所述运动估计;
调整用于所述特征估计的事件数量;以及
采用为一个或多个参数确定最佳设置的算法。
18.一种在根据原始数据重建的图像中校正运动的装置,包括:
用于根据所述原始数据估计所述重建图像内部的感兴趣区域的质心的模块;
用于基于所述感兴趣区域的质心的变化计算运动校正平移矢量(86)的模块;
用于利用所述运动校正平移矢量对于与所述感兴趣区域相关联的所述原始数据进行运动校正的模块;以及
用于利用所校正的原始数据重建与所述感兴趣区域对应的运动校正图像的模块。
19.一种利用用于重建图像的数据在重建图像中校正运动的医学成像系统(10),包括:
多个探测器(16),用于探测电子正电子湮灭事件;
符合探测器(38),用于识别所述电子正电子湮灭事件中的基本同时的符合事件对;
响应线处理器(40),用于处理所述事件并识别连接每个符合事件对内的事件的响应线;
一个或多个处理器(22,42),包括:
分组器(52),用于将所识别的响应线分成多个组;
滤波器(54),用于滤除不与感兴趣区域相交叉的响应线;
特征处理器(58),用于估计每个组的经滤波的响应线的质心;
变异处理器(82),用于确定在每个组的响应线的质心和预先选择的质心之间的偏差;
LOR移位器(88),用于将每组响应线移位由所述偏差给出的量以获得运动校正后的响应线;以及
重建系统(44),用于根据所述运动校正后的响应线重建与所述感兴趣区域对应的运动校正图像。
20.根据权利要求19所述的医学成像系统(10),其中,所述响应线为列表模式数据。
21.根据权利要求19所述的医学成像系统(10),其中,将所述感兴趣区域定位至小于所述医学成像系统(10)的视场。
22.根据权利要求19所述的医学成像系统(10),其中,所述处理部件(22,42)执行简化的基于Siddon的算法和基于矢量的算法之一,以便于估计所述感兴趣区域的质心。
23.根据权利要求19所述的医学成像系统(10),还包括:
呈现如下内容中的至少一个的呈现部件(26):
所述运动校正图像,
所述运动校正图像和未校正图像的同时显示,以及
未校正图像,利用所述运动校正的感兴趣区域替代所述未校正图像的对应感兴趣区域。
24.一种对图像进行校正的方法,包括:
针对多个短时间帧从运动模糊的感兴趣区域生成多个中间感兴趣区域图像(62);
针对所述多个中间感兴趣区域图像(62)中的每个,调整一个或多个位于感兴趣区域(56)内部并与对应的响应线(64)相交叉的体素(66-76)的体素值;
计算每个所述中间感兴趣区域图像(56)的质心;
组合各个质心以计算在采集的全部时间范围内的总计质心;
计算每个质心和所述总计质心之间的差异;
基于所述差异计算运动校正平移矢量(86);
通过将数据中的响应线移位所述运动校正平移矢量给出的一个量来校正与所述感兴趣区域对应的数据;以及
利用所述校正数据重建与所述感兴趣区域对应的运动校正图像。
25.根据权利要求24所述的方法,其中,将所述感兴趣区域定位至感兴趣的解剖结构。
26.根据权利要求24所述的方法,其中,与所述感兴趣区域对应的所述数据为列表模式数据。
27.一种对图像进行校正的方法,包括:
计算表示连接一对响应线的最短线段的中点的矢量;
根据中点在所述感兴趣区域中的分布计算质心的时域运动;
通过将所述初始数据中的响应线移位到总质心来对于所述数据进行运动校正;以及
利用所述校正数据重建与所述感兴趣区域对应的运动校正图像。
28.一种对图像进行校正的装置,包括:
用于针对多个短时间帧从运动模糊的感兴趣区域生成多个中间感兴趣区域图像(62)的模块;
用于针对所述多个中间感兴趣区域图像(62)中的每个,调整一个或多个位于感兴趣区域(56)内部并与对应的响应线(64)相交叉的体素(66-76)的体素值的模块;
用于计算每个所述中间感兴趣区域图像(56)的质心的模块;
用于组合各个质心以计算在采集的全部时间范围内的总计质心的模块;
用于计算每个质心和所述总计质心之间的差异的模块;
用于基于所述差异计算运动校正平移矢量(86)的模块;
用于通过将数据中的响应线移位所述运动校正平移矢量给出的一个量来校正与所述感兴趣区域对应的数据的模块;以及
用于利用所述校正数据重建与所述感兴趣区域对应的运动校正图像的模块。
29.根据权利要求28所述的装置,其中,将所述感兴趣区域定位至感兴趣的解剖结构。
30.根据权利要求28所述的装置,其中,与所述感兴趣区域对应的所述数据为列表模式数据。
31.一种对图像进行校正的装置,包括:
用于计算表示连接一对响应线的最短线段的中点的矢量的模块;
用于根据中点在所述感兴趣区域中的分布计算质心的时域运动的模块;
用于通过将所述初始数据中的响应线移位到总质心来对于所述数据进行运动校正的模块;以及
用于利用所述校正数据重建与所述感兴趣区域对应的运动校正图像的模块。
32.一种形成运动校正图像的方法,包括:
在未校正图像内界定感兴趣区域(56);
基于所述感兴趣区域(56)的质心的变化计算平移矢量(86);以及
利用所述平移矢量(86)针对局部运动引起的差异校正所述感兴趣区域(56)。
33.根据权利要求32所述的方法,还包括:
针对多个短时间帧从所述感兴趣区域(56)生成多个中间图像(62);
针对所述多个中间图像(62)中的每个,调整一个或多个位于所述感兴趣区域(56)内部并与对应的响应线(64)相交叉的体素(66-76)的体素值;
计算每个所述中间图像(62)的质心;以及
组合各个质心以计算在采集的全部时间范围内的所述感兴趣区域(56)的所述质心。
34.根据权利要求33所述的方法,还包括:
计算所述各个质心和所述感兴趣区域(56)的质心之间的差异;以及
基于所述差异计算所述平移矢量(86)。
35.根据权利要求32所述的方法,其中,校正所述感兴趣区域(56)的步骤包括:
从与所述感兴趣区域(56)对应的列表模式数据计算所述平移矢量(86)。
36.根据权利要求32所述的方法,其中,将所述感兴趣区域定位至视场内的解剖结构。
37.一种在图像中局部校正运动的方法,包括:
在未校正图像内定位感兴趣区域;
基于所定位的感兴趣区域的质心的变化计算运动校正平移矢量(86),
利用所述运动校正平移矢量(86)对于所定位的感兴趣区域进行运动校正;以及
通过用所定位的经校正的感兴趣区域取代所定位的未校正的感兴趣区域来形成校正图像。
38.根据权利要求37所述的方法,其中,所定位的感兴趣区域对应于列表模式数据。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US77746906P | 2006-02-28 | 2006-02-28 | |
US60/777,469 | 2006-02-28 | ||
PCT/US2007/061597 WO2007100955A2 (en) | 2006-02-28 | 2007-02-05 | Local motion compensation based on list mode data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101454801A CN101454801A (zh) | 2009-06-10 |
CN101454801B true CN101454801B (zh) | 2012-12-05 |
Family
ID=38459705
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007800070638A Expired - Fee Related CN101454801B (zh) | 2006-02-28 | 2007-02-05 | 基于列表模式数据的局部运动补偿 |
Country Status (5)
Country | Link |
---|---|
US (1) | US8144962B2 (zh) |
EP (1) | EP1991959B1 (zh) |
JP (1) | JP5254810B2 (zh) |
CN (1) | CN101454801B (zh) |
WO (1) | WO2007100955A2 (zh) |
Families Citing this family (60)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010517672A (ja) * | 2007-02-07 | 2010-05-27 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 数量的データ分析の運動補償及び治療 |
CN101681520B (zh) * | 2007-05-30 | 2013-09-25 | 皇家飞利浦电子股份有限公司 | Pet局部断层摄影 |
WO2009040690A2 (en) * | 2007-09-24 | 2009-04-02 | Koninklijke Philips Electronics N.V. | Preclinical time of flight pet imaging |
US8017915B2 (en) | 2008-03-14 | 2011-09-13 | Reflexion Medical, Inc. | Method and apparatus for emission guided radiation therapy |
US8472683B2 (en) * | 2008-05-09 | 2013-06-25 | General Electric Company | Motion correction in tomographic images |
EP2286274B1 (en) * | 2008-05-28 | 2019-07-24 | Koninklijke Philips N.V. | Geometrical transformations preserving list-mode format |
CN102098964B (zh) * | 2008-06-13 | 2015-04-08 | 皇家飞利浦电子股份有限公司 | 生理列表模式核成像中用于计数的最佳时间采样的反向数据重建 |
EP2291827B1 (en) * | 2008-06-18 | 2012-11-14 | Koninklijke Philips Electronics N.V. | Radiological imaging incorporating local motion monitoring, correction, and assessment |
EP2163201A1 (en) * | 2008-09-15 | 2010-03-17 | Westfälische Wilhelms-Universität Münster | List mode-based respiratory and cardiac gating in positron emission tomography |
JP5523791B2 (ja) | 2008-10-27 | 2014-06-18 | 株式会社東芝 | X線診断装置および画像処理装置 |
CN102282589A (zh) * | 2009-01-19 | 2011-12-14 | 皇家飞利浦电子股份有限公司 | 列表模式pet成像中的区域重建和定量评估 |
US8787643B2 (en) | 2009-02-17 | 2014-07-22 | Koninklijke Philips B.V. | Functional imaging |
US9196030B2 (en) | 2009-09-18 | 2015-11-24 | Koninklijke Philips N.V. | System and method for determining a property of blur in a blurred image |
US8466419B2 (en) * | 2009-10-01 | 2013-06-18 | Kabushiki Kaisha Toshiba | System and method for enhanced sampling via helical scanning and list-mode reconstruction in positron emission tomography |
US8824757B2 (en) | 2009-12-10 | 2014-09-02 | Koninklijke Philips N.V. | Method and apparatus for using time of flight information to detect and correct for motion in imaging scans |
US9111381B2 (en) * | 2010-01-27 | 2015-08-18 | Koninklijke Philips N.V. | Shift-varying line projection using graphics hardware |
WO2012037151A2 (en) * | 2010-09-13 | 2012-03-22 | University Of Southern California | Efficient mapping of tissue properties from unregistered data with low signal-to-noise ratio |
US8761467B2 (en) * | 2010-10-04 | 2014-06-24 | General Electric Company | Method and apparatus for assessing motion correction |
US8811695B2 (en) | 2010-12-14 | 2014-08-19 | General Electric Company | Methods, apparatus and articles of manufacture to adaptively reconstruct medical diagnostic images |
US9305377B2 (en) | 2011-01-05 | 2016-04-05 | Koninklijke Philips N.V. | Method and apparatus to detect and correct motion in list-mode PET data with a gated signal |
CN103650095B (zh) | 2011-03-31 | 2016-12-07 | 反射医疗公司 | 用于在发射引导的放射治疗中使用的系统和方法 |
CN103534730B (zh) | 2011-05-12 | 2016-12-14 | 皇家飞利浦有限公司 | 列表模式动态图像重建方法和装置 |
US8421021B2 (en) * | 2011-06-21 | 2013-04-16 | General Electric Company | Motion correction of SPECT images |
CN102324089B (zh) * | 2011-07-13 | 2013-04-03 | 南方医科大学 | 基于广义熵与mr先验的pet图像最大后验重建方法 |
US9788803B2 (en) | 2011-09-15 | 2017-10-17 | Shimadzu Corporation | Medical-data processing device and radiation tomography apparatus having the same |
CN103126701B (zh) * | 2011-11-30 | 2016-04-27 | 株式会社东芝 | 正电子发射计算机断层摄影装置和图像处理装置 |
BR112015002566A2 (pt) * | 2012-08-10 | 2017-07-04 | Koninklijke Philips Nv | sistema de tomografia por emissão de pósitrons (pet), método de tomografia por emissão depósitrons (pet), e meio de armazenamento não transitório legível por computador. |
DE102012216292B4 (de) * | 2012-09-13 | 2021-02-18 | Siemens Healthcare Gmbh | Magnetresonanzbaueinheit, eine Magnetresonanzvorrichtung mit der Magnetresonanzbaueinheit sowie ein Verfahren zu einem Bestimmen einer Bewegung eines Patienten während einer Magnetresonanzuntersuchung |
KR20140042461A (ko) | 2012-09-28 | 2014-04-07 | 삼성전자주식회사 | 영상을 생성하는 방법 및 장치 |
CA2906227C (en) * | 2013-03-15 | 2020-04-28 | Synaptive Medical (Barbados) Inc. | System and method for magnetic resonance coil arrangement |
KR101604812B1 (ko) * | 2014-01-15 | 2016-03-18 | 삼성전자주식회사 | 의료 영상 처리 장치 및 그에 따른 의료 영상 처리 방법 |
CN103942763A (zh) * | 2014-05-03 | 2014-07-23 | 南方医科大学 | 一种基于mr信息引导的体素水平pet图像部分容积校正方法 |
DE102014219376A1 (de) * | 2014-09-25 | 2016-03-31 | Siemens Aktiengesellschaft | Verfahren zu einer Erfassung eines hochaufgelösten Magnetresonanzbilddatensatzes zumindest eines begrenzten Körperbereichs mit zumindest einer anatomischen Struktur eines Patienten |
US9606245B1 (en) | 2015-03-24 | 2017-03-28 | The Research Foundation For The State University Of New York | Autonomous gamma, X-ray, and particle detector |
CN107924730B (zh) | 2015-06-10 | 2021-09-28 | 反射医疗公司 | 高带宽双态多叶式准直器设计 |
CN107133549B (zh) | 2016-02-29 | 2020-11-24 | 上海联影医疗科技有限公司 | Ect运动门控信号获取方法及ect图像重建方法 |
WO2017156316A1 (en) | 2016-03-09 | 2017-09-14 | Reflexion Medical, Inc. | Fluence map generation methods for radiotherapy |
CN106097384B (zh) * | 2016-05-30 | 2018-12-18 | 江苏赛诺格兰医疗科技有限公司 | 一种确定门控信号的方法和装置 |
CN106548473B (zh) * | 2016-11-07 | 2019-03-08 | 赛诺联合医疗科技(北京)有限公司 | 一种构建相位图像的方法及装置 |
US10695586B2 (en) | 2016-11-15 | 2020-06-30 | Reflexion Medical, Inc. | System for emission-guided high-energy photon delivery |
WO2018093849A1 (en) | 2016-11-15 | 2018-05-24 | Reflexion Medical, Inc. | Methods for radiation delivery in emission-guided radiotherapy |
CN116943051A (zh) | 2016-11-15 | 2023-10-27 | 反射医疗公司 | 放射治疗患者平台 |
WO2018099772A1 (en) * | 2016-11-29 | 2018-06-07 | Koninklijke Philips N.V. | Interactive targeted ultrafast reconstruction in emission and transmission tomography |
KR101946576B1 (ko) * | 2016-12-23 | 2019-02-11 | 삼성전자주식회사 | 의료 영상 장치 및 의료 영상 처리 방법 |
US10307116B2 (en) | 2017-02-22 | 2019-06-04 | Uih America, Inc. | System and method for detecting organ motion |
US11410349B2 (en) | 2017-03-24 | 2022-08-09 | Koninklijke Philips N.V. | Methods for data driven respiratory motion estimation |
WO2018183748A1 (en) | 2017-03-30 | 2018-10-04 | Reflexion Medical, Inc. | Radiation therapy systems and methods with tumor tracking |
EP3630286A4 (en) | 2017-05-30 | 2021-03-03 | RefleXion Medical, Inc. | PROCESS FOR IMAGE-GUIDED RADIATION THERAPY IN REAL-TIME |
EP3641884B8 (en) | 2017-06-22 | 2024-01-03 | RefleXion Medical, Inc. | Methods for calculating bounded dose-volume histograms (dvh) for evaluating a treatment plan |
CN114699655A (zh) | 2017-07-11 | 2022-07-05 | 反射医疗公司 | 用于pet检测器余辉管理的方法 |
CN111095026B (zh) | 2017-07-26 | 2023-11-17 | 反射医疗公司 | 放射治疗的图形表示 |
JP7315961B2 (ja) | 2017-08-09 | 2023-07-27 | リフレクション メディカル, インコーポレイテッド | 放出誘導放射線療法における異常検出のためのシステムおよび方法 |
US10993689B2 (en) * | 2017-08-31 | 2021-05-04 | General Electric Company | Method and system for motion assessment and correction in digital breast tomosynthesis |
US11369806B2 (en) | 2017-11-14 | 2022-06-28 | Reflexion Medical, Inc. | Systems and methods for patient monitoring for radiotherapy |
US11151726B2 (en) * | 2018-01-10 | 2021-10-19 | Canon Medical Systems Corporation | Medical image processing apparatus, X-ray diagnostic apparatus, and medical image processing method |
WO2019160958A1 (en) | 2018-02-13 | 2019-08-22 | Reflexion Medical, Inc. | Beam station treatment planning and radiation delivery methods |
CN108876730B (zh) * | 2018-05-24 | 2022-03-04 | 东软医疗系统股份有限公司 | 校正运动伪影的方法、装置及设备和存储介质 |
US11426131B2 (en) * | 2018-12-17 | 2022-08-30 | Siemens Medical Solutions Usa, Inc. | Automated motion correction in PET imaging |
WO2020150505A1 (en) | 2019-01-16 | 2020-07-23 | Reflexion Medical, Inc. | Methods for setup corrections in radiation therapy |
WO2023103975A1 (zh) * | 2021-12-06 | 2023-06-15 | 柯纪纶 | 医学影像的移动检测和校正方法、系统和计算机可读介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1997005574A1 (en) * | 1995-07-27 | 1997-02-13 | Imperial Cancer Research Technology Limited | Raw data segmentation and analysis in image tomography |
US6292683B1 (en) * | 1999-05-18 | 2001-09-18 | General Electric Company | Method and apparatus for tracking motion in MR images |
US20050123183A1 (en) * | 2003-09-02 | 2005-06-09 | Paul Schleyer | Data driven motion correction for nuclear imaging |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6980683B2 (en) * | 2000-08-28 | 2005-12-27 | Cti Pet Systems, Inc. | On-line correction of patient motion in three-dimensional positron emission tomography |
US7574249B2 (en) * | 2005-02-08 | 2009-08-11 | General Electric Company | Device-less gating of physiological movement for improved image detection |
-
2007
- 2007-02-05 JP JP2008557440A patent/JP5254810B2/ja not_active Expired - Fee Related
- 2007-02-05 CN CN2007800070638A patent/CN101454801B/zh not_active Expired - Fee Related
- 2007-02-05 US US12/279,996 patent/US8144962B2/en active Active
- 2007-02-05 EP EP07717552.9A patent/EP1991959B1/en not_active Not-in-force
- 2007-02-05 WO PCT/US2007/061597 patent/WO2007100955A2/en active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1997005574A1 (en) * | 1995-07-27 | 1997-02-13 | Imperial Cancer Research Technology Limited | Raw data segmentation and analysis in image tomography |
US6292683B1 (en) * | 1999-05-18 | 2001-09-18 | General Electric Company | Method and apparatus for tracking motion in MR images |
US20050123183A1 (en) * | 2003-09-02 | 2005-06-09 | Paul Schleyer | Data driven motion correction for nuclear imaging |
Non-Patent Citations (2)
Title |
---|
GJ Klein等.Fine-Scale Motion Detection Using Intrinsic List Mode PET Information.《IEEE Workshop on Mathematical Methods in Biomedical Image Analysis》.2001,第71-78页. * |
Philippe P. Bruyant等.Correction of the Respiratory Motion of the Heart by Tracking of the Center of Mass of Thresholded Projections: A Simulation Study Using the Dynamic MCAT Phantom.《IEEE TRANSACTIONS ON NUCLEAR SCIENCE》.2002,第49卷(第5期),第2159-2166页. * |
Also Published As
Publication number | Publication date |
---|---|
WO2007100955A8 (en) | 2009-04-02 |
EP1991959B1 (en) | 2017-08-30 |
US8144962B2 (en) | 2012-03-27 |
JP5254810B2 (ja) | 2013-08-07 |
EP1991959A2 (en) | 2008-11-19 |
CN101454801A (zh) | 2009-06-10 |
WO2007100955A3 (en) | 2008-08-07 |
WO2007100955A2 (en) | 2007-09-07 |
US20100166274A1 (en) | 2010-07-01 |
JP2009528139A (ja) | 2009-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101454801B (zh) | 基于列表模式数据的局部运动补偿 | |
US9495771B2 (en) | Systems and methods for motion correction in positron emission tomography imaging | |
Tong et al. | Image reconstruction for PET/CT scanners: past achievements and future challenges | |
US8098916B2 (en) | System and method for image-based attenuation correction of PET/SPECT images | |
US9747701B2 (en) | Systems and methods for emission tomography quantitation | |
CN1809841B (zh) | 运动补偿的重建方法、设备与系统 | |
US9031300B1 (en) | System and method reconstructing a nuclear medicine image using deformed attenuation image | |
CN101681520B (zh) | Pet局部断层摄影 | |
CN102483852B (zh) | 利用基于飞行时间信息逐个事件生成的图像内容的飞行时间正电子发射断层摄影重建 | |
US7507968B2 (en) | Systems and methods for correcting a positron emission tomography emission image | |
US8575554B2 (en) | Methods and systems for scatter estimation in positron emission tomography | |
US20060214097A1 (en) | Detector head position correction for hybrid SPECT/CT imaging apparatus | |
US7057178B1 (en) | Method and system for imaging using a filter for Time-of-Flight PET | |
CN102067176A (zh) | 结合局部运动监测、校正和评估的辐射成像 | |
CN101765865A (zh) | 核成像中的运动校正 | |
GB2450073A (en) | Regional reconstructions of spatially distributed functions | |
CN105147312A (zh) | Pet图像获取方法及系统 | |
US10304219B2 (en) | Intra reconstruction motion correction | |
CN101499173A (zh) | 一种pet成像中卡尔曼滤波图像重建方法 | |
US20160163095A1 (en) | Systems and methods for selecting imaging data for principle components analysis | |
CN104274201A (zh) | 乳腺层析成像方法和系统及成像设备和图像采集处理方法 | |
RU2493572C2 (ru) | Геометрические преобразования, сохраняющие данные в виде списка | |
US10852449B2 (en) | System and method for self-time alignment calibration for a positron emission tomography system | |
US10482634B2 (en) | Systems and methods for imaging with anisotropic voxels | |
CN104000618B (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 | ||
C14 | Grant of patent or utility model | ||
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: 20121205 Termination date: 20200205 |