用于发射断层摄影的噪声模型选择制造技术

对于发射断层摄影的许多应用，例如，具有可信度的运动建模或ＳＵＶ量化，准确的误差估计是有益的。由于影响ＰＥＴ图像的噪声特性的参数的多样性，对于所有成像情况和数据处理配置而言，使用单一误差模型导致不准确的误差估计。本发明专利技术通过提供包含用于不同成像情况的多个预定噪声模型的数据库而避免了这个问题。于是依据给定的成像情况，能够手动地或自动地选择最适当的噪声模型。因此，只需为每个模型执行一次提取正确噪声模型的费时的过程，例如，通过使用自举方法或通过分析重复测量，并且该过程可以由采集系统的厂商来完成，从而使得临床医师能够立即访问来自数据库的最佳误差模型。

【技术实现步骤摘要】
【国外来华专利技术】用于发射断层摄影的噪声模型选择本专利技术涉及一种系统，设备和方法，用以为发射断层摄影成像情况和数据处理配置(set-up)提供噪声模型。在发射断层摄影中，所获得的数据是基于放射性核衰变的，它是 一种随机过程。已经充分了解了该过程的本质和统计学特性。然而， 对于从所获得的数据中重构得到的图像的统计学特性而言还不能这 么说。从探测器自身通过获得电子信号开始，到数据内插、校正和重 构方法，图像处理流水线的每一部分都影响所得到的图像的噪声特 性。因此，不可能定义一个简单的、适合于所有成像情况和所有数据 处理配置的噪声模型。临床医师和研究人员必须为每种配置个别地和独立地开发噪声 模型。在许多情况下，有可利用的可靠的预定义噪声模型来帮助确定 准确的误差估计是有利的。例如，在动态数据的运动建模中，误差估 计被用作调整过程的加权。较好的误差估计会导致更好的配合。此外， 在静态成像中，误差估计被用来向临床医师提供生物参数的可信度， 例如对于肿瘤学应用中的SUV量化。由于影响PET图像噪声特性的参数的多样性，对于所有成像情 况和数据处理配置仅使用一个误差模型是行不通的。在文献中可以找 到的对误差估计的分析方法是针对特定问题的，因此仅适于所有成像 配置的一小类子集，例如见H. H. Barrett等人，A^/m /VopeW/as 0/7& 五M爿/gonY/wr /.TTzeoo^ Phys. Med. Bio.， 39, pp.833-46， 1994，其全部 内容合并于此作为参考。自举(bootstrap)方法是用于噪声特性估计 的另一禾中方法，例如见，I,Buavt，爿7Vow-P(3rame/77'c 5oofWra/ ^p; raac/z ^wa/戸'"g iStof/Wca/尸ra/ eWas q/";SPEC7"a"(i尸五7Vmages， Phys. Med. Bio., 47， pp.1761-75, 2002,其全部内容合并于此作为参考。与分 析方法相比，自举方法具有这样的优势能够产生数据副本，并且可 以为几乎全部成像配置针对该数据副本进行统计学分析。整个图像采集和处理链被视为黑箱，无需与其要素有关的任何详细信息作为该自 举分析的输入。然而，这个方法存在极其费时的缺陷，这阻碍了其在 大多数临床情况下被单个临床医师和研究人员使用。现在参考图l，示出了用于误差分析的典型过程，其中，在步骤101中获得数据，在步骤102中将数据校正并重构为图像。随后在步 骤103中为进一步的图像分析使用了简化的噪声模型(如果其确实存 在的话)。总的来说，目前没有足够灵活、能够覆盖临床和研究关注的所有 情况的可利用的噪声模型。在当前标准图像分析软件中使用的模型在 大多数情况下不能够反映正确的噪声特性。本专利技术的系统、设备和方法提供了一种有效且高效率的方法，为 通常使用的可预定义的配置提供了可重复使用的预定的噪声模型，从 而消除了临床医师和研究人员为这些配置独立地、个体地执行噪声模 型确定的需要。本专利技术的系统、设备和方法为临床医师/研究人员提供了一种工 具，其从此类模型的数据库中自动选择噪声模型，且该噪声模型适用 于特定成像情况。在优选实施例中，数据库被填充了噪声模型(根据 成像情况及其相应的配置而编制索引)，这些噪声模型是预先使用用 于多个不同但典型的成像情况和配置的PET图像的自举分析来提取 的。从而，能够克服上述方法的这两个主要缺陷1. 分析型噪声模型速度快但其灵活性不足以覆盖全部相关成像 情况，例如不同重构方法、散射校正(scattercorrection)方法等的使 用。2. 统计学自举分析(或类似的方法，例如重复测量)用来为每种 成像情况确定噪声特性，因此非常灵活，但这个方法极其费时。对于可靠的误差/噪声估计的问题， 一种快速且灵活的解决方案 是，为临床医师/研究人员提供用于标准成像情况的预定噪声模型的 数据库以及用于易于基于特定应用的特性而为该应用选择(或自动选 择)适当模型的工具。图2A示出了本专利技术的优选实施例。在步骤101 中的数据的采集、校正和图像重构之后，基于该配置以及采集、数据校正和重构的参数，从数据库201中选择适当的噪声模型202。特定 应用的噪声模型202允许改进的噪声估计，其得到了更好的图像，并 为临床医师/研究人员提供了有价值的信息。附图说明图1示出了采集和成像过程的典型过程； 图2A示出了根据本专利技术而改进的图1的过程；图2B示出了实现图2A过程的图像分析设备；图3示出了影像(phantom)数据集重构和统计学分析；图4示出了相对于计数速度的平方根而绘制的标准偏差；以及 图5示出了包含图2B的图像分析设备的一种成像系统。 本领域普通技术人员会理解提供以下说明的目的是为了解释而 不是为了限制。技术人员会理解在本专利技术的精神和所附的权利要求的 范围之内存在许多变化。从当前说明中省略了已知的功能和操作的不 必要的细节，以免使得本专利技术含糊不清。 本专利技术的一个优选实施例如下1. 为一组预定的临床上相关的成像情况/配置产生方差图像。例 如，这些方差图像可以用自举方法或通过使用重复测量来产生。可以 使用的另一种技术是蒙特卡罗模拟(Monte-Carlo simulation)。2. 每个配置的噪声特性按以下所述被参数化，这些参数与定义 特定成像情况的参数一起被存储在数据库中，例如"脑部成像、基于 CT的衰减校正、非一致性散射校正(non-uniform scatter correction)、 采用2次迭代和8个子集的OSEM重构，以及2 X 2 X 4mm3三维像素 尺寸"。3. 由采集系统的厂商提供该数据库，但它也能够被更新以包括 由用户产生的新参数集。4. 在步骤101中用户象通常一样获取并处理数据。依据成像流 水线的设定，从数据库201中选择适当的噪声模型(手动或自动的)(见以上步骤2中的实例)。然后，在步骤203中，这个模型被用于 对重构的图像102进行进一步分析，例如，运动建模、具有可信等级 的SUY量化等。以下的实例示出了步骤2的优选实施例。现在参照图3，在此所示的是分别以滤波反投影(filtered back projection) (a)和迭代行活 动最大似然算法2d (RAMLA2D) (iterative row action maximum likelihood algorithm) (c)(例如见J. A. Browne禾卩A. R. De Pierro， ^/"五mz'wz'ow 7bmogra/^_y， IEEE Transaction on Medical Imaging,Vo1.15， pp. 687-699， 1996)重构的影像数据集。以自举方法产生相应的方差 图像(b)禾n (d)。自举方法用于产生图3所示的方差图像，其基于 相同的数据，但以不同的重构方法来重构。自举方法是基于计算机的统计学方法，用于确定依据实验数据而 估计的统计值9 (例如中间数)的准确度(例如见Efron和Tibshirani, 爿w /"的c/wc"'ow to f/ze丑oo加ra; , New York: Chapman and Hall, 1993 )。它需要实验样本X气X,本文档来自技高网...

【技术保护点】
一种用于为发射断层摄影成像流水线选择噪声模型的方法，包括以下步骤：为发射断层摄影提供具有至少一个预定噪声模型的数据库，通过用相应成像情况／配置的至少一个噪声特性和至少一个参数作为搜索项，对所述预定噪声模型编制索引；借助于至少 一个噪声特性和至少一个参数表征所述成像流水线；将所提供的数据库的搜索项与所述成像流水线的至少一个噪声特性和至少一个参数相比较，以进行匹配；并且当匹配时，选择相应成像情况／配置的预定噪声模型作为所述成像流水线的噪声模型。

【技术特征摘要】
【国外来华专利技术】...

【专利技术属性】
技术研发人员：M布施，R布林克斯，
申请(专利权)人：皇家飞利浦电子股份有限公司，
类型：发明
国别省市：NL[荷兰]