用于运动校正的医学图像的方法和器件的制作方法

专利名称：用于运动校正的医学图像的方法和器件的制作方法
技术领域：
本文公开的主题通常涉及成像系统，并且更特别地涉及用于运动校正的医学图像的器件和方法。
背景技术：
存在使用不同模态(modality)扫描的多模态成像系统，例如，计算机断层照相(Computed Tomography, CT)、磁共振成像(Magnetic Resonance Imaging, MRI)、正电子发射断层照相(Positron Emission Tomography, PET)、以及单光子发射计算机断层照相(Single Photon Emission Computed Tomography, SPECT)。在操作期间，常规的成像系统的图像质量可能受正在成像的对象的运动影响。特别地，成像对象的运动可能降低图像质量。更具体地，在图像采集期间由对象的运动产生图像伪影(image artifact)。呼吸运动·是在医学成像系统中遇到的哺乳动物(例如，人和动物)中的非自愿的运动的公共来源。在图像检查期间(例如当内科医生确定损伤的尺寸、确定损伤的位置、或量化损伤时)，呼吸运动可能带来误差。此外，在多模态系统(例如集成PET/CT系统)中，PET和CT图像应该彼此配准。然而，由于CT图像典型地是在短时间段期间采集的，由CT图像生成的衰减图表示在没有呼吸运动的完全呼气期间患者的衰减特性。与此相反，PET图像典型地是在相对长的时间段之上采集的，该时间段中由于长的采集时间而允许患者自由呼吸。两个数据采集模式之间的失配可能导致衰减校正的PET图像中的图像伪影。一种用于减少图像伪影的已知方法是将多个呼吸阶段的CT图像(或最大强度CT图像)平均来模仿在多个呼吸循环上收集的PET采集的效果。另一种用于减少图像伪影的已知方法是使用呼吸门控CT采集来生成衰减校正图，其更好地匹配呼吸门控PET采集的呼吸特性。还有一种方法可以包含请求患者在扫描期间屏住呼吸。然而，因为可能在几分钟之上采集PET数据，所以患者典型地必须在PET采集期间呼吸数次，其导致图像伪影。

发明内容
在一个实施例中，提供了一种用于减少图像中的运动相关的成像伪影的方法。该方法包含获取感兴趣的区域的图像数据集，将多变量数据分析技术应用到图像数据集来生成运动信息，基于运动信息将图像数据集分类为多个箱(bin)，以及使用多个箱中的至少一个来生成感兴趣的区域的图像。在另一实施例中，提供了一种包含扫描仪和操作地耦合到该扫描仪的计算机的医学成像系统。编程计算机以使用扫描仪采集对象的成像数据集，使用图像数据集生成多个中间图像，将多变量数据分析技术应用到多个中间图像来生成运动信息，基于运动信息将中间图像分类为多个箱，以及使用多个箱的至少一个生成感兴趣的区域的图像。在另外的实施例中，提供了一种非暂时性计算机可读介质。提供了采用程序编码的计算机可读介质，该程序编程以指导计算机。为使用扫描仪来采集对象的成像数据集，使用图像数据集生成多个中间图像，将多变量数据分析技术应用到多个中间图像来生成运动信息，基于运动信息将中间图像分类为多个箱，以及使用多个箱的至少一个生成感兴趣的区域的图像。


图I是根据各种实施例形成的示范性成像系统的简化概略图示。图2是图示根据各种实施例的用于减少运动相关的成像伪影的方法的流程图。图3是图示根据各种实施例形成的示范性正弦图的简化框图。图4是表示根据各种实施例可以生成的平均信息的示范性图像。图5是表示根据各种实施例可以生成的运动信息的示范性图像。·图6是根据各种实施例可以生成的示范性运动信息。图7是表示根据各种实施例可以生成的运动信息的另一示范性图像。图8是根据各种实施例可以生成的示范性运动信号。图9是根据各种实施例可以生成的示范性图像。图10是表示根据各种实施例可以生成的运动信息的另一示范性图像。图11是根据各种实施例可以生成的示范性运动信息。图12是表示根据各种实施例可以生成的运动信息的另一示范性图像。图13是根据各种实施例可以生成的示范性运动信息。图14是根据各种实施例形成的示范性成像系统的图解图示。
具体实施例方式当结合附图来阅读时，将更好地理解上文概要以及下文各种实施例的详细描述。就图示各种实施例的功能块的图表的程度而言，功能块不必要指示硬件电路之间的划分。从而，例如，可用单件硬件(例如，通用信号处理器或随机存取存储器块、硬盘等)或多件硬件来实现一个或多个功能块(例如，处理器或存储器)。类似地，程序可以是独立程序，可以作为子程序并入操作系统中，可以是安装的软件包中的功能等。应理解，各种实施例不限于附图中所示的排列和工具。如本文所使用的，以单数叙述并冠以词语“一”的元件或步骤应该理解为不排除多个所述元件或步骤，除非明确声明这样的排除。此外，参照本发明的“一个实施例”并不旨在解释为排除也并入所叙述的特征的额外的实施例的存在。此外，除非明确地相反声明，否贝U “包括”或“具有”具有特定的性质的元件或多个元件的实施例可以包含不具有那种性质的额外的元件。还如本文所使用的，短语“重建图像”不旨在排除在其中生成表示图像而不是可视图像的数据的本发明的实施例。因此，如本文所使用的术语“图像”泛指可视图像和表示可视图像的数据两者。然而，许多实施例生成或配置为生成至少一个可视图像。图I是示范性正电子发射断层照相(PET)成像系统10的简化框图。PET成像系统10可以包含PET扫描仪12，可操作PET扫描仪12来生成患者16的四维(four-dimensional, 4D)发射数据集(emission dataset) 14。4D 发射数据集 14 表示三个空间维度和时间。4D发射数据集14可以存储在存储器装置20中。存储器装置20可以包含存储器电路的任何组合，包含电、磁性、和/或光学系统。存储器装置20可以包含例如只读存储器(ROM) 22、随机存取存储器(RAM) 24、和/或大容量存储器26。存储器装置20可在其中存储由处理器30执行的程序指令。程序指令可以以任何相配合适的计算机语言来书写，例如Matlab。处理器30可以是相配合适的处理系统的任何一个或其组合，其中例如微处理器、数字信号处理器、以及现场可编程逻辑阵列、等等。处理系统可以实施为任何相配合适的计算装置，例如计算机、个人数字助理(PDA)、膝上型计算机、笔记本计算机、基于硬驱动的装置，或者可以接收、送出和存储数据的任何装置。可以提供输入/输出(I/O)装置32用于接收用户输入和/或用于提供输出给用户。I/O装置32还可以是通信链路，其使数据能传送到其它装置和/或从其它装置传送，例如通过网络。I/O装置32可以包含用于显示PET图像的显示装置(未示出)。I/O装置32可以包含用于接收来自用户的输入的装置。例如，该装置可以包含鼠标、跟踪球、或触敏接口。成像系统10还包含运动特性化模块50，适配其以识别发射数据集14内的运动，并且利用所识别的运动来将发射数据集14分割为多个箱，其中箱的至少一个然后用来重建患者16的至少一个图像，如以下更详细地讨论的。在示范性实施例中，发射数据集14是4D发射数据集，并且运·动特性化模块安装在处理器30中。运动特性化模块可以实现为硬件装置、软件或其组合。图2是由在图I所示的PET成像系统10进行的示范性方法的框图。在示范性实施例中，可以使用运动特性化模块50来实现方法100。更具体地，方法100可以提供为具有在其上记录的指令的机器可读介质或多个介质，用于指引处理器30进行本文描述的方法的实施例。该介质或多个介质可以是任何类型的CD-R0M、DVD、软盘、硬盘、光盘、闪速RAM驱动，或其它类型的计算机可读介质的或其组合。方法100在发射数据集14上提供运动特性化、或运动减小来解决对象16的运动。在操作中，方法100识别对象16的运动并且重新组织发射数据集14来使对象16的图像能重建。应该认识到，虽然相对于从PET成像系统10获取的发射数据集14来描述方法100，但是方法100还可以应用到从CT成像系统获取的传送数据集。此外，方法100可以应用到使用本文讨论的任何成像模态获取的任何图像数据集或正弦图数据，并且发射数据集14仅是示范性的。在102，获取对象16的发射数据集14 (每个在图I中示出)。在示范性实施例中，使用PET成像系统10 (在图I中示出)获取发射数据集14。可以通过进行对象16的发射扫描而获取发射数据集14来产生发射数据集14。可选地，可以从在对象16的以前的扫描期间收集的数据来获取发射数据集14，其中发射数据集14已经存储在存储器中，例如存储器装置20 (在图I中示出)。发射数据集14可以以任何格式存储，例如列表模式数据集。可以在患者16的实时扫描期间获取发射数据集14。例如，可以在患者16的实时检查期间从PET扫描仪12接收发射数据时，对发射数据进行本文描述的方法。在104，利用发射数据集14来生成在时间上的多幅中间图像150。在示范性实施例中，中间图像150是正弦图150。可选地，中间图像可以是重建的图像。如图3中所示的，多个正弦图150可以包含在第一时间巾贞(例如300毫秒)之上米集的第一正弦图集150a。多个正弦图150可以包含在第一正弦图集150a后续的时间帧之上采集的第二正弦图集150b等。从而在扫描过程期间，对于每个时间帧(例如对于每个300ms的时间帧)生成正弦图集。生成的正弦图150的数量基于在成像系统10中安装的检测器的数量、期望的运动速度、在102进行的扫描的长度等。从而多个正弦图150可以包含在扫描过程中的连续时间帧之上采集的多个正弦图集(150a. . . 150η)，例如可以对于整个扫描过程之上的每个300毫秒采集单个正弦图集。可以通过以正弦图模式操作系统10而生成正弦图150。正弦图模式通常是指在其中可选地具有同一飞行时间(Time-of-Flight，TOF)的湮没事件以(距离轴的半径、角度)格式存储在正弦图中的采集模式。响应的阵列已知为正弦图。应该认识到，其它方法和/或装置可以用于数据存储，并且本文描述的正弦图表示存储数据的一个这样的示范性方法。在105，预处理在104生成的正弦图150来获取对应的集151。例如，可以通过这些来预处理正弦图150，即1)在正弦图的不同元素(即，数据尺寸和噪声减小)之上求平均；2)移除正弦图的某些元素(例如，在期望发生运动的地方或在患者数据贡献高空间信息的地方应用掩蔽，例如排除仅意外的或分散的重合的区域)；3)对于时间上的已知的改变(例如，由于放射性衰变)的校正；4)对数据应用逐元素(element-wise)操作来转换数据为近似高斯分布，优选用统一的方差。例如，对于泊松分布的数据，适当的变换包含Anscombe变换、Freeman-Tukey变换、或简单地取平方根。平均和掩蔽两者可以依靠外部信息,例如，以·前的扫描、和/或关于运动的主方向和位置的先验知识。例如，能由从以前CT扫描获取的图像确定患者轮廓或感兴趣的区域(例如膈膜)。然后可以使用PET扫描仪模型将此区域信息变换为正弦图。掩蔽还可以从平均数据自身(例如，基于阈值检测背景计数)导出。在106，利用预处理的正弦图集151a. . . 151η来生成时间平均信息152。例如，假设对于每个时间帧在104生成25个正弦图集150a... 150η，然后在105，生成25个预处理的正弦图150a. .. 150η，并且在106，对于150a中的每个正弦图中的每个元素计算25个对应元素的平均来生成平均信息152。图4图示可以从正弦图集150a. . . 150η生成的示范性平均信息152。对于正弦图集150a. . . 150η生成的平均信息的视觉表示指代为图4中的PC0。如以上所讨论的，通过对于每个预处理的正弦图集150 (在时间上)确定平均信息而采集平均信息152。例如，图4，例如PC0，图示通过扫描四个示范性统一对象160、162、164、以及166而采集平均信息，其中在此情况下，四个对象的每个的运动是在垂直方向上。如图4中所示的，四个对象160、162、164、以及166每个展现在每个对象的外边缘附近的平均像素强度的改变。此外，每个对象160、162、164、以及166分别的内部区域170、172、174、以及176，以及周围背景区180大体上是密实的，指示在时间上像素强度有很小的偏移或没有偏移。在108，对于每个正弦图集151a. .. 151η，从相应正弦图集150减去平均信息152以生成零均值(zero mean)信息，然后将其输入到PCA。因此，假设在示范性扫描期间有在时间上生成的25个正弦图集150，该方法包含对于25个正弦图集150的每个生成平均信息152，并且然后对于25个正弦图集150的每个生成零均值信息154。在110，正弦图集154然后用作对于PCA方法的输入(我们将称预处理和PCA的组合为逐掩蔽体积PCA(Masked Volume Wise PCA)或MVW-PCA)。PCA的描述是广泛可用的。总之，该方法找到数据的(估计)协方差矩阵的主特征向量。其输出是多个特征向量(每个对应于与集154a相同尺寸的正弦图集)和它们的对应特征值。有最大特征值的特征向量解释正弦图集154a... 154η中的最大变化(在时间上)。因为有最低特征值的特征向量对应于数据中的噪声，所以仅计算最初少数的特性向量是可能的并且可取的。这些特征向量通常称作“主成分”并且以降低特性值的顺序标记为PC1、PC2等。如图5中所示的，PCl图像表示在图4中使用的示范性数据中的第一主成分。在此例子中，PCl图像图示由于运动在时间上具有最大改变的区。例如，使用大体上密实的颜色(对应于非常接近于零的值)示出的对象160、162、164、以及166的中央区指示对象的此部分中的原始数据的值在时间上大体上不改变。此外，PCl使用表示运动的程度或量的阴影区域示出对象160、162、164、以及166(图5中所示)的边缘周围的存在一些运动。例如,PCl图示表示对象160的运动的阴影区域160a和160b。区域162a和162b表示对象162的运动。区域164a和164b表示对象164的运动，并且区域166a和166b表示对象166的运动。此外，阴影指示PCl中对应的元素中的值的符号。例如，160a中的阴影指示正值，而160b中的阴影指示负值。对象160、162、164、以及166本身示出为大体上密实的颜色，这是因为有在对象160、162、164、以及166的中心区域内的运动的相对较小的影响。此外，因为有在对象160、162、164、以及166的周围的运动的相对较小的影响，所以背景180也图示为大体上相同的颜色。再次参照图2，在110，利用对于每个正弦图集150生成的零均值信息154来生成·一个或多个权重因素156,在图6中不出。更具体地，对于每个分析的主成分，PCA输出单个权重因素、信号和/或轨迹。在此例子中，PCA输出对应第一主成分的权重因素156。权重因素156图示为图6中的点,每个时间巾贞为一个点。在示范性实施例中，通过将PCl的值与有零均值信息154a. . . 154η的每个正弦图集做逐元素相乘，并在正弦图集之上求和，从而对于正弦图集154a、154b. . . 154η的每个获取单个数字来生成权重因素156。更具体地，参照图5中示出的图像PCI，PCl指示在发射数据中有时间改变，发射数据由区域160a、160b、162a、162b、164a、164b、166a、以及166中的阴影图示。如果在特定时间帧中，对象160、162、164、以及166处于其最顶位置，则除了其将为例如正的区域160a、162a、164a、以及166a以外，零均值信息与PCl的相乘结果将在各处大体上为零。因此，此时间帧中的权重因素156的值将为正。相反地，如果对象处于其最底位置，则信号的值将为负(或实践中与处于最顶位置具有相反的符号)。因此，权重因素156的值可以是单个正数，其指示在第一方向上已观察到幅度改变。权重因素156的值可以是单个负数，其指示在相反的第二方向上已观察到幅度的负改变。权重因素156的值可以是零，其指示相对于对象的平均位置没有改变。此外,权重因素156的值可以在正值和负值的范围内。例如，低的正值可以指示相对于均值的对象的相对较小的移动。然而更高正值可以表示对象16的更大的移动。在示范性实施例中，MVW-PCA可以配置为更明确地识别零均值信息154。例如，图7图示图像PC2，其表示可以由PCA输出的第二成分。应该认识到，由MVW-PCA生成的PC图像的数量可以基于所希望的精确度来特征化零均值信息。例如，MVW-PCA可以仅利用PCl图像来生成运动信号，其将在下文中更详地讨论。对应于可以从PCA输出的不同的主成分的运动信号可以表示不同类型的运动，例如，呼吸运动、心脏运动、身体移动等，如将在下文中更详细地讨论的。再次参照图2，在112，利用在110导出的权重因素156来生成运动信号250。例如，图8图示使用MVW-PCA导出的示范性运动信号250，其中x轴表示时间并且y轴表示对于每个正弦图集150在110导出的权重因素156的幅度。在示范性实施例中，对于每个正弦图集150以时间顺序的方式导出权重因素156。例如，MVW-PCA对于在示范性实施例中作为在扫描期间生成的第一正弦图的正弦图集150a导出权重因素156a。然后MVW-PCA对于第二正弦图集150b导出权重因素156b，并且对于在扫描过程中采集的最后的正弦图集150η导出权重因素156η。如图8中所示的，以时间顺序方式安排多个权重因素156并且通过表示权重因素156的点绘制线。从权重因素156导出运动信号250。在一个实施例中，运动信号250与权重因素156大体相同。在另一实施例中，通过滤波权重因素156以例如减少噪声来获取运动信号250。在114，使用运动信号250 (在图8中示出)将发射数据集14 (例如正弦图150)分类为多个箱300。例如，图8图示示范性运动信号250，利用其将正弦图150映射为多个箱300。此外，图9图示编号为300. . . 310的多个箱，即η = 6箱。从而在示范性实施例中，将正弦图150分类为编号为300、302、304、306、308、以及310的6个箱300。然而，应该认识至IJ，在图8中图示的箱的数量是示范性的，并且在操作期间，可以利用少于6个箱或多于6个箱。因此，箱300、302、304、306、308、以及310的每个包含成像数据集14中的总信息的近似 1/6。
·
例如，假设在102进行的为采集发射数据的扫描的总长度为3分钟。此外，假设将成像数据集14分类为6个箱，其中每个相应箱包含近似30秒的信息。在示范性实施例中，基于上述运动信号250将正弦图150分类为相应箱。更具体地，基于如上述讨论的导出的运动信号250的值将正弦图150分类为具体的箱。因此，在示范性实施例中，每个箱包含表示在扫描过程期间处于相同的近似空间位置中的对象160、162、164、以及166的发射数据。例如，箱300可以包含表示在呼吸阶段的开始采集的对象160的信息，并且箱310可以包含表示在呼吸阶段的结束采集的成像数据。此外，每个介于其间的箱，例如302、304、306、以及308可以包含表示吸气和呼气之间的运动状态的成像数据。适配箱300、302、304、306、308、以及310的每个以接收在多个呼吸循环之上采集的成像数据。此外，适配箱300、302、304、306、308、以及310的每个以接收表示在患者的呼吸循环中的近似相同的点的成像数据。因此，箱300、302、304、306、308和310的每个包含表示患者16的某些运动状态的成像数据。在示范性实施例中，利用从运动信号250采集的信息将成像数据14分为6个大体上相等的部分，并将大体上相等的部分存储在相应箱300、302、304、306、308、以及310中。从而当在114处完成门控过程时，6个箱300、302、304、306、308、以及310的每个包含与相同的运动特性关联的发射数据。在示范性实施例中，在将发射数据分类为箱后，多个箱300、302、304、306、308、以及310可以重建为表示发射活动分布的图像，并且然后空间配准。例如，在一个实施例中，箱300可以选择为参照箱并且箱302、304、306、308、以及310可以配准到参照箱300。可以使用刚性或非刚性配准将箱302、304、306、308、以及310配准到参照箱300。可以由操作员手工地或由处理器30自动地进行刚性和非刚性配准。应该认识到，在很小或没有运动的某些区可以不要求配准。运动信息可以用来调整配准过程中的成本函数。一般而言，可以利用运动信息来识别受运动影响或未受运动影响的区来改善配准过程。在示范性实施例中，进行非刚性配准包含变换三维(3D)空间中的箱302、304、306,308,以及310内的信息来将箱302、304、306、308、以及310内的信息对准到参照箱300。例如，箱302中的图像可以相对于参照箱300中的图像轻微倾斜。因此，倾斜箱302内的图像来将该图像与参照箱300中的图像对准。还重新对准剩余的箱304、306、308以及310来大体上匹配参照箱300中的图像。在操作中，可以通过选择解剖的或其它的特征/点/界标，以及沿图像中检测的边缘或边界使用的这些特征或点对准的图像来实现刚性配准过程。备选地，不同的标记可以用来识别已知的解剖位置。刚性配准还可以基于弯曲的轮廓，例如图像内的骨。刚性配准还可以是基于体积或基于表面的。然而，应该意识到，可以进行任何的刚性配准过程，其包含最优化或计算某些可比较的准则或类似的措施。在另一实施例中，可以利用非刚性配准过程来进行对箱300、302、304、306、308、以及310内的信息的运动校正。在操作中，非刚性配准或弹性的配准包含非刚性变换。这些非刚性变换允许图像特性的局部扭曲并提供解决局部变形的配准。非刚性变换方法包含，例如多项式扭曲、基于平滑的函数(薄板样条或小波)的插值、以及物理连续体模型(粘性流体模型和大的变形微分同胚)。使用在箱300、302、304、306、308、以及310中分类的信息进行非刚性配准。非刚性配准可以包含例如点或界标的扭曲并且提供沿着有点或界标的插值和相关性的轮廓的最佳拟合。备选地，可以进行混合过程，其比较图像体素并且混合对应的区域。一般而言，局部非刚性配准包含任何类型的弹性·的变形模型，其允许不同图像集中的变化或移动。再次参照图2，在116，将在箱300、302、304、306、308、以及310中存储的2D图像，
即运动校正图像一起平均来生成单个运动校正图像。在另一实施例中，MVW-PCA用来构造多个运动信号，其对应不同类型的移动，例如呼吸、心脏的和/或整体身体移动。例如，图9图示通过扫描示范性统一对象504(类似于心肌层)而采集的两个时间帧500和502，其中在此情形下运动表示收缩/扩张和在垂直方向上以不同(较慢)速率的平移。图10图示在此例子中从在图9中所示的对象504获取的第一主成分PCI。MVW-PCA示出与图5中类似的特征。图11图示对应于此成分的权重因素510，如在以前的实施例中所构造的。在此例子中，其大部分由平移(或呼吸移动)确定。图12图示第二成分PC2，并且在图13中示出对应的权重因素520。在此例子中，PC2大部分由收缩性的(或心脏的)移动决定。更具体地，图12中示出的区域530大体上具有正值，而区域532大体上具有负值。因此，并且再次参照图13，当在特定时间帧中环面具有大的半径时，从PC2构造的权重因素520将如在540所示大体上为正。类似地，当环面具有小的半径时，信号将如在542所示大体上为负。类似于以前的实施例，然后权重因素用作双门控模块的输入来将发射数据分箱，其中每个门对应于特定的呼吸和心脏的状态。例如，如果识别6个呼吸状态和3个心脏状态，则将使用总共18个箱，其中信号510确定呼吸状态并且权重因素520确定心脏状态。在示范性实施例中，然后重建分箱的发射数据并且结果图像彼此配准。配准算法可以考虑到不同类型的运动，并且例如估计平移和收缩的组合运动。在另外的实施例中，分析主成分和/或其对应的信号以识别运动类型来允许取决于应用的进一步处理。作为例子，应用可以对于呼吸运动校正，而保持心脏的门控来确定射血分数(ejection fraction)。识别过程使用相应移动的已知的特性。例如，呼吸运动大部分发生在足至头的方向上，并且因此对应的主成分将在如与收缩性运动相反的结构的顶部和底部上大部分具有相反的符号的值。在一个实施例中，确定心脏的近似中心并且从主成分构造径向剖面。如果剖面在顶部和底部方向上大体上相同，则PC标记为对应于心脏，反之亦然。备选地或另外的，对于对应于主成分的信号进行频率分析。如果主频率在例如每秒1-2个循环的范围中，则PC标记为对应于心脏搏动，而例如每秒1/10到1/2个循环的范围标记为对应于呼吸移动。类似地，非常慢的移动或非常突然的移动可以标记为对应于全身的移动。在本实施例中，在识别后，根据其已知的特性处理PC或它们的信号来减少噪声和/或移除不同类型的运动之间的干扰是有利的。在给定的例子中，信号510，其一旦大部分标记为呼吸，可以用中通滤波器滤波来移除不对应于呼吸的频率，并且类似适当的滤波器用于其它信号。在一些应用中，使用已滤波的信号来确定门控是有利的。此外，可以利用备选的多变量数据分析技术来提取成分。这样的多变量数据分析技术可以包含，例如，独立成分分析(Independent Component Analysis, I CA)和正则化主成分分析(regularized Principal ComponentAnalysis,rPCA),尽管应该认识到可以利用其它的分析技术，但是其中的许多使用PCA作为初始步骤。这些技术中的一些允许关于期望的运动图案的先验信息的使用，其可以在嘈杂的情形中是有利的。作为简单例子，正弦图集150可以首先空间滤波和/或时间滤波来移除不想要的频率。例如，如果仅需要检测呼吸移动，则时间中通滤波器可以在使用PCA之前应用到数据。在示范性实施例中，利用逐掩蔽体积主成分分析(MVW-PCA)统计技术。在操作中，MVW-PCA分析发射数据集14来找到临·床有关的成分。在此情况下，利用MVW-PCA来识别运动相关的成分。更具体地，MVW-PCA是利用其来从诊断信息和噪声分离运动信息的数学算法。在示范性实施例中，MVW-PCA实现为可以存储在例如运动特性化模块50和/或存储器装置20上的指令。本文描述的各种实施例可以提供如图14中所示的多模态成像系统400。多模态成像系统400可以是任何类型的成像系统，例如，不同类型的医学成像系统，例如图I中示出的正电子发射断层照相(PET) 10、单光子发射计算机断层照相(SPECT)、计算机断层照相(CT)、超声系统、磁共振成像(MRI)或能够生成诊断图像的任何其它系统。各种实施例不限于多模态医学成像系统，而也可以用于单模态医学成像系统(例如独立PET成像系统、独立CT成像系统、磁共振成像(MRI)、以及光子发射计算机断层照相(SPECT))。此外，各种实施例不限于用于成像人类对象的医学成像系统，而也可以包含用于成像非人类对象的兽用或非医学系统等。参照图14，多模态成像系统400包含第一模态单元410和第二模态单元412。这两个模态单元使多模态成像系统400能使用第一模态单元410以第一模态扫描对象或患者16并且使用第二模态单元412以第二模态扫描患者16。多模态成像系统400允许对于以不同模态的多个扫描来在单个模态系统之上便于增加的诊断能力。在一个实施例中，多模态成像系统400是CT/PET成像系统400，例如第一模态410是CT成像系统并且第二模态412是PET成像系统10。成像系统400示出为包含与CT成像系统410关联的门架420和与PET成像系统10关联的门架422。在操作期间，使用例如机动台426将患者16定位在通过成像系统400限定的中央开口 424内。门架420包含X射线源430，其向门架420的相反侧上的检测器阵列432投射x射线光束。检测器阵列432由包含多个检测器元件的多个检测器行(未示出)形成，其一起感测穿过患者16的投射的X射线。每个检测器元件产生表示碰撞X射线光束的强度的电信号，并且因此如光束穿过患者16来允许光束的衰减的估计。在扫描来采集X射线衰减数据期间，门架420和装配在其上的部件关于转动中心转动。此外，PET成像系统包含配置为采集发射数据的检测器(未示出)。
在操作期间，机动台426将患者16移动入门架420和/或门架422的中央开口424以响应从操作员工作站440接收的一个或多个命令。然后工作站440操作第一和第二模态410和412来扫描患者16并采集患者16的衰减和/或发射数据两者。工作站440可以实施为个人计算机(PC)，其靠近成像系统400定位并且经由通信链路442硬接线到成像系统440。工作站440可以实施为便携式计算机，例如膝上计算机或手持式计算机，其传送信息到成像系统400并且从成像系统400接收信息。可选地，通信链路442可以是无线通信链路，其使信息能无线地传送到工作站440或从工作站440传送到成像系统400。在操作中，工作站440配置为实时控制成像系统400的操作。工作站440还编程为进行本文描述的医学图像诊断采集和重建过程。更具体地，工作站444可以包含运动特性化模块50，适配其以识别4D数据集14内的运动，并利用所识别的运动来将4D数据集14分割为多个箱，然后这些箱用来重建患者16的图像，如以上更详细地讨论的。操作员工作站440包含中央处理单元(CPU)或计算机444、显示器446、以及输入装置448。如本文所使用的，术语“计算机”可以包含任何基于处理器或基于微处理器的系统，包含使用微控制器、精简指令集计算机(RISC)、专用集成电路(ASIC)、现场可编程门阵列(FPGA)、逻辑电路、以及任何其它能够执行本文描述的功能的电路或处理器的系统。以上·例子仅是示范性的，因此并不旨在以任何方式限制术语“计算机”的定义和/或含义。在示范性实施例中，计算机444执行在一个或多个存储元件或存储器中存储的指令集，以便处理从第一和第二模态410和412接收的信息。存储元件还可以存储期望的或需要的数据或其它信息。存储元件可以以信息源或物理存储元件的形式位于计算机444内。指令集可能包含各种命令，其指导计算机444作为处理机器来进行例如本文描述的各种实施例的方法和过程的具体操作。指令集可以是以软件程序的形式。如本文所使用的，术语“软件”和“固件”是可互换的，并包含用于由计算机执行的在存储器中存储的任何计算机程序，包含RAM存储器、ROM存储器、EPROM存储器、EEPROM存储器、以及非易失性RAM(NVRAM)存储器。以上存储器类型只是示范性的，并且因此不作为对于计算机程序的存储可用的存储器类型的限制。软件可以是例如系统软件或应用软件的各种形式。另外，软件可以是单独程序的集合、在较大程序内的程序模块或部分程序模块的形式。软件还可以包含面向对象编程形式的模块化编程。由处理机器对输入数据的处理可以响应用户命令、或响应以前的处理结果、或响应由另一处理机器做出的请求。计算机444连接到通信链路442并从输入装置448接收输入，例如用户命令。输入装置448可以是例如键盘、鼠标、触摸屏面板、和/或声音识别系统等。通过输入装置448和关联的控制面板开关，操作员可以控制CT成像系统410和PET成像系统10的操作和用于扫描的患者16的定位。类似地，操作员可以控制显示器446上的结果图像的显示，并可以使用由计算机444执行的程序进行图像增强功能。一些本文描述的各种实施例的技术效果是为了改善图像质量并且对于呼吸运动、心脏运动、患者运动等补偿。更具体地，可以利用MVW-PCA来改善重要运动改变的检测和视觉化，并且基于所检测的改变将图像数据分类为箱以便改善图像质量并且捕获运动以便减少将出现的伪影，并且改善各种对象(例如肿瘤)的检测。可以使用例如PET成像系统、SPECT成像系统、MRI成像系统、超声成像系统、和/或MRI成像系统来采集图像数据。
本文描述的各种实施例提供有形的和非暂时性机器可读介质或多个介质，其具有用于处理器或计算机的记录在其上的指令以操作成像器件来进行本文描述的方法的实施例。介质或媒体可以是任何类型的CD-R0M、DVD、软盘、硬盘、光盘、闪速RAM驱动、或其它类型的计算机可读介质或其组合。各种实施例和/或部件(例如，监视器或显示器、或部件和其中的控制器)还可以实现为一个或多个计算机或处理器的部分。计算机或处理器可以包含计算装置、输入装置、显示单元和接口，例如，用于访问因特网。计算机或处理器可以包含微处理器。微处理器可以连接到通信总线。计算机或处理器还可以包含存储器。存储器可以包含随机存取存储器(RAM)和只读存储器(ROM)。计算机或处理器还可以包含存储装置，其可以是硬盘驱动或可移除存储驱动，例如软盘驱动、光盘驱动等。存储装置还可以是用于装载计算机程序或其它指令到计算机或处理器的其它的类似组件。应理解以上描述旨在说明性的而不是限制性的。例如，上述的实施例(和/或其各方面)可彼此组合使用。此外，可作出许多修改来使特定的情况或材料适配于各种实施·例的教导而没有背离它们的范围。虽然本文描述的材料的尺寸和类型旨在定义各种实施例的参数，但是它们绝对没有限制并且只是示范性的。回顾以上描述时，许多其它实施例对于本领域技术人员是明显的。因此，应该参照所附的权利要求来确定各种实施例的范围，连同这样的权利要求的等效的全部范围进行声明。在所附的权利要求中，术语“包含”和“在其中”用作相应术语“包括”和“其中”的简明英文等效。此外，在下文权利要求中，术语“第一”、“第二”和“第三”等只用作标签，并且并不旨在在它们的对象上施加数量要求。另外，下文权利要求的限制没有用组件加功能的格式书写并且并不旨在基于35U.S.C. § 112第六段来解释，除非并且直到权利要求限制清楚地使用后接没有进一步结构的功能的声明的短语“组件用于”。本书面描述使用包含最佳模式的例子来公开各种实施例，并且也使本领域的任何技术人员能实践各种实施例，包含制作并且使用任何装置或系统并且进行任何并入的方法。各种实施例的可专利范围由权利要求定义，并且可以包含本领域技术人员所想到的其它例子。如果这样的其它例子具有没有不同于权利要求的字面语言的结构元件，或这样的其它例子包含具有与权利要求的字面语言无实质差别的等效结构元件，则这样的其它例子旨在落入权利要求的范围内。
权利要求
1.一种医学成像系统(10)，包括扫描仪(I)和操作地耦合到所述扫描仪的计算机(30)，其中所述计算机编程为 使用所述扫描仪采集对象(16)的图像数据集(14)； 使用所述图像数据集生成多个中间图像(150)； 将多变量数据分析技术应用到所述多个中间图像来生成运动信息(250)； 基于所述运动信息将所述中间图像分类为多个箱(300);以及 使用所述多个箱中的至少一个来生成感兴趣的区域的图像(500)。
2.如权利要求I所述的医学成像系统(10)，其中所述计算机(30)还编程为获取四维(4D)正电子发射断层照相(PET)发射数据集和4D单光子发射计算机断层照相(SPECT)发射数据集中的至少一个、或4D计算机断层照相(CT)数据集。
3.如权利要求I所述的医学成像系统(10)，其中为进行所述多变量数据分析技术，所述计算机(30)还编程为进行所述多个中间图像的逐掩蔽体积主成分分析(MVW-PCA)来生成所述运动信息。
4.如权利要求I所述的医学成像系统(10)，其中所述计算机(30)还编程为 使用所述运动信息的频率分析来识别不同类型的运动。
5.如权利要求I所述的医学成像系统(10)，其中所述计算机(30)还编程为 应用所述多个中间图像(150)的所述多变量数据分析技术来生成多个权重因素(156);以及 基于所述权重因素将所述中间图像分类为多个箱(300)。
6.如权利要求I所述的医学成像系统(10)，其中所述计算机(30)还编程为 应用所述多个中间图像(150)的所述多变量数据分析技术来生成多个权重因素(156)； 使用所述多个权重因素生成运动信号(250);以及 基于所述运动信号将所述中间图像分类为多个箱(300)。
7.一种非暂时性计算机可读介质(50)，采用程序编码，所述程序编程来指导计算机(30)来 使用扫描仪⑴采集对象(16)的图像数据集(14)； 使用所述图像数据集生成多个中间图像(150)； 将多变量数据分析技术应用到所述多个中间图像来生成运动信息(250)； 基于所述运动信息将所述中间图像分类为多个箱(300);以及 使用所述多个箱的至少一个来生成感兴趣的区域的图像(500)。
8.如权利要求7所述的计算机可读介质(50)，其中为应用所述多变量数据分析技术，所述计算机可读介质还编程以指导计算机(30)进行所述多个中间图像(150)的逐掩蔽体积主成分分析(MVW-PCA)来生成所述运动信息(250)。
9.如权利要求7所述的计算机可读介质(50)，所述计算机可读介质还编程来指导计算机(30)以 应用所述多个中间图像(150)的所述多变量数据分析技术来生成多个权重因素(156);以及 基于所述权重因素将所述中间图像分类为所述多个箱。
全文摘要
本发明的名称是“用于运动校正的医学图像的方法和器件”。提供一种用于减少图像中的运动相关的成像伪影的方法，包含获取感兴趣的区域的图像数据集(14)，使用图像数据集生成多个中间图像(150)，将多变量数据分析技术应用到多个中间图像来生成运动信息(250)，基于运动信息将中间图像(150)分类为多个箱(300)，以及使用多个箱中的至少一个来生成感兴趣的区域的图像(500)。
文档编号A61B5/055GK102783966SQ20121026403
公开日2012年11月21日 申请日期2012年5月3日 优先权日2011年5月3日
发明者K·F·J·J·蒂勒曼斯, P·拉奇菲尔, S·P·S·拉索尔 申请人:通用电气公司