CN101953693B - 完整或部分正电子发射断层造影衰减图的基于模型的估计 - Google Patents

完整或部分正电子发射断层造影衰减图的基于模型的估计 Download PDF

Info

Publication number
CN101953693B
CN101953693B CN201010232204.XA CN201010232204A CN101953693B CN 101953693 B CN101953693 B CN 101953693B CN 201010232204 A CN201010232204 A CN 201010232204A CN 101953693 B CN101953693 B CN 101953693B
Authority
CN
China
Prior art keywords
pet
model
decay
parameter
decay pattern
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.)
Active
Application number
CN201010232204.XA
Other languages
English (en)
Other versions
CN101953693A (zh
Inventor
马塞厄斯·芬谢尔
拉尔夫·雷德贝克
克里斯琴·J·米歇尔
查尔斯·C·沃森
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
Siemens Medical Solutions USA Inc
Original Assignee
Siemens AG
Siemens Medical Solutions USA Inc
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, Siemens Medical Solutions USA Inc filed Critical Siemens AG
Publication of CN101953693A publication Critical patent/CN101953693A/zh
Application granted granted Critical
Publication of CN101953693B publication Critical patent/CN101953693B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/4808—Multimodal MR, e.g. MR combined with positron emission tomography [PET], MR combined with ultrasound or MR combined with computed tomography [CT]
    • G01R33/481—MR combined with positron emission tomography [PET] or single photon emission computed tomography [SPECT]
    • 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
    • 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/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • A61B6/5247—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from an ionising-radiation diagnostic technique and a non-ionising radiation diagnostic technique, e.g. X-ray and ultrasound
    • 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/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • 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/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
    • A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels

Landscapes

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

Abstract

本发明涉及完整或部分正电子发射断层造影衰减图的基于模型的估计。示例性实施例涉及一种在磁共振(MR)扫描器和正电子发射断层造影(PET)单元中校正衰减的方法。所述方法包括采集在PET单元的视野内部的对象的PET正弦图数据。所述方法还包括：基于参数化的模型实例和PET正弦图数据的最大似然性期望最大化(MLEM)，产生衰减图。

Description

完整或部分正电子发射断层造影衰减图的基于模型的估计
技术领域
示例性实施例涉及一种用于估计在正电子发射断层造影和磁共振断层造影(MR-PET)中的衰减图的方法。
背景技术
正电子发射断层造影(PET)与磁共振断层造影(MR)一起被应用于医疗诊断。MR是用于表示身体内部的结构和层的成像方法，而PET允许显示和量化在活体中的代谢活动。
PET利用了正电子发射体和正电子湮灭的特别的特征，以量化地确定器官和/或细胞区域的功能。利用该技术，在检查之前对患者给予具有放射性核的合适的放射性药物标记物。随着它们衰减，放射性核放射出正电子，后者在一个短的距离之后与电子相互作用，导致湮灭发生。这产生两个伽马量子，它们在相对的方向(偏移180)上飞开。在一个特定的时间窗之内由两个相对的PET探测器模块探测到伽马量子(相合测量)，其结果是，在所述两个探测器模块的直线连接线上的一个位置上定位湮灭位置。
在PET情况下，为探测的目的，探测器模块一般地覆盖机架弧形长度的大部分。探测器模块被划分为具有数毫米的边长的探测器元件。当探测到伽马量子时，每个探测器元件产生指定了时间和探测位置的事件记录。该信息被传递到快速逻辑单元并且被比较。如果两个事件在一个最大的时间间隔之内同时发生，则假定，在两个相关的探测器元件之间的连接线上正在发生伽马衰减过程。使用断层造影算法、例如反投影，来重建PET图像。
在PET系统中，诸如MR-PET系统中，伽马量子通过位于各个伽马量子的起源位置(site of origin)和PET探测器之间的任何东西而被衰减。为了防止图像伪影，在PET图像的重建中必须考虑该衰减。对象，诸如在患者身体内的组织、空气和MR/PET系统本身的一部分(例如患者定位台)，位于在患者身体中的伽马量子的起源位置和起作用的PET探测器之间。在伽马量子的起源位置和起作用的PET探测器之间的对象的衰减值被考虑并且被编辑到衰减图(μ图)中。
衰减图包含被检查的体积的每个体积元素(体素)的衰减值。因此，例如，可以产生对于患者定位台的衰减图。例如，同样适用为了MR检查而附加到患者的局部线圈。为了产生衰减图，确定并结合衰减值。可以通过例如各个组件的CT记录或PET传输测量来确定这些衰减值。可以仅一次地测量这种衰减图，因为衰减值贯穿各个组件的寿命都不改变。
公知如下的方法，在该方法中，从解剖结构的MR图像可以确定患者身体的衰减值，并且可以被添加到衰减图。在这种情况下，使用特定的MR序列，通过其可以识别例如不同的衰减的组织类别(例如肺组织)。借助于MR图像就可以根据衰减的组织类别的位置，来将合适的衰减值指定给衰减图。
然而，横向的(transaxial)MR视野通常小于PET视野。因此，待检查的对象的仅一部分在PET视野中。因此，位于MR视野之外的衰减值的获得变得困难。
通过将MR图像分割为不同的组织类型并且给不同的组织类型指定相应的衰减值来进行PET衰减图的基于MR的估计。然而，该方法没有涉及位于MR视野外部的扫描区域。
近来，最大似然性期望最大化(maximum-likelihood expectation maximization，MLEM算法)已经被用来从PET正弦图(sinogram)数据同时重建发射和衰减图。PET正弦图数据可以涉及PET原始数据、PET计数或PET计数数据。术语“图像”是从PET正弦图数据重建的图像。从基于MR的分割或其它公知方法所得到的衰减图可以被用来初始化MLEM算法。
用于基于MR的衰减校正的其它方法包括使用具有已知的衰减的图集、模型或参考图像，诸如互相配准(coregistered)的相应的CT、PET传输图像或从光学3D扫描得到的身体轮廓。然后，将实际的MR图像与具有已知的衰减的图集或参考配准，并且从配准信息和附加的后处理方法得到实际的衰减图。
发明内容
示例性实施例针对使用MLEM对完整或部分PET衰减图的基于模型的估计。
至少一个示例性实施例公开了一种在MR扫描器和PET单元中校正衰减的方法。所述方法包括采集在PET单元的视野内部的对象的PET正弦图数据，和基于参数化的模型实例(model instance)的和PET正弦图数据的最大似然性期望最大化(MLEM)，产生衰减图。
至少另一个示例性实施例提供了一种在MR扫描器和PET单元中校正衰减的方法。该方法包括采集在PET单元的视野内部的对象的PET正弦图数据，和采集在MR扫描器的视野内部的对象的MR数据。基于参数化的模型实例的和PET正弦图和MR数据的最大似然性期望最大化(MLEM)，产生衰减图。MLEM被参数化的模型实例的模型参数限制。
另一个示例性实施例提供一种包括正电子发射断层造影(PET)单元的设备，该PET单元具有多个探测单元并且配置为采集在PET单元的视野内部的对象的PET正弦图数据。磁共振(MR)扫描器被配置为采集在MR扫描器的视野内部的对象的MR数据。计算机被配置为基于参数化的模型实例的和采所集的PET正弦图和MR数据的最大似然性期望最大化(MLEM)产生衰减图，其中，MLEM被参数化的模型实例的模型参数限制。
附图说明
从以下借助附图的详细描述中，示例性实施例将变得更加易理解。图1-4表示如在此描述的非限制性的示例性实施例。
图1示出了可以用于示例性实施例中的用于重叠的MR和PET图像显示的设备；
图2示出了按照一种示例性实施例使用用于重建PET图像的MLEM来估计完整的PET衰减图的方法；
图3示出了按照一种示例性实施例使用基于MR的衰减图和用于重建PET图像的MLEM来估计完整的PET衰减图的方法；以及
图4示出了按照一种示例性实施例使用MLEM来细化初始衰减图的方法。
具体实施方式
现在参考附图详细描述不同的示例性实施例，在附图中示出了一些示例性实施例。在附图中，为了清楚起见可能夸张层和区域的厚度。
相应地，示例性实施例可以有不同的修改和替换方式，通过在附图中的例子示出其实施例并且在此详细描述。然而应当理解，不是旨在限制示例性实施例公开的具体方式，而是相反，示例性实施例要覆盖落入示例性实施例的范围之内的所有的修改、等价和替换情况。贯穿附图的描述相似的标号指代相似的元件。
可以理解的是，尽管术语第一、第二等在此可能用来描述不同的元素，这些元素不应该被这些术语限制。这些术语仅用来区分一个元素与另一个。例如，第一元件可以被称为第二元件，并且，类似的，第二元件可以被称为第一元件，而不脱离示例性实施例的范围。如此处使用的那样，术语“和/或”包括一个或多个相关列出的项目的任何和所有组合。
可以理解的是，当提到一个元件被“连接”或“耦接”到另一个元件，可以是直接连接或耦接到另一个元件或也可以存在中间元件。相反，当提到一个元件是“直接连接”或“直接耦接”到另一个元件时，没有中间元件。用来描述在元件之间的关系的其它词语应该被解释为以类似方式(例如“在...之间”相对“直接地在...之间”，“相邻”相对“直接相邻”等。)
在此使用的术语为了描述具体实施例的目的仅并且不是旨在限制示例性实施例。如在此使用的，单数形式“一个”和“这个”是旨在也包括复数形式，除非上下文中清楚地指出。还将理解的是，当在此使用术语“包括”“包含”时，具体指出所述特征、整数、步骤、操作、元件和/或组件的出现，但是不排除出现或附加一个或多个其它特征、整数、步骤、操作、元件、组件和/或他们的组
有关空间的术语，例如“在...之下”、“在...之下”、“低于”、“在...之上”“高于”等，在此可以被用来便于描述如在附图中示出的元件或在元件或特征和其它元件或特征之间的关系。可以理解的是，有关空间的术语旨在包含除了在附图中示出的取向之外使用或操作中的设备的不同的取向。例如，如果要翻转在附图中的设备，被描述为“低于”其它元件或特征或“在其下面”的元件就将被取向为“在其它元件或特征的上面”。因此，例如，术语“在...之下”可以包括在其之上以及在其之下两种取向。设备可以另外取向(旋转90度或在其它方向看或参考)并且有关空间的描述在此使用的应该相应地被解释。
还应该注意的是，在一些替换的实施方式中，功能/动作提到的可以不是按照附图中标注的顺序发生。例如，连续示出的两个附图实际上可以同时并列进行或有时候可以按照相反顺序进行，这取决于涉及的功能/动作。
除非另有定义，这里使用的所有术语(包括技术和科学术语)具有与示例性实施例所属的领域中的技术人员通常理解的相同含义。可以进一步理解的是，诸如通常使用的字典中定义的那些术语，应当被解释为具有与它们在相关技术和本公开的上下文中的含义一致的含义，而不用理想化或以过于正式的含义来解释它们，除非在此这样明确地定义。
以软件、或算法和对在计算机存储器之内的数据比特的操作的符号表示呈现示例性实施例的部分和相应的详细描述。这些描述和表示是本领域普通技术人员向本领域其它普通技术人员有效地转达其工作的实质的那些。算法，如此处使用的术语，并且如其通常使用的，被理解为导致一个期望的结果的自相一致的步骤顺序。步骤是需要物理量的物理操作的那些。通常，尽管不是必须的，这些量采取能够被存储、传输、组合、比较和其它操作的光、电或磁信号的形式。有时被证明是方便的是，原则上为了通常使用的原因，将这些信号作为比特、值、元素、符号、特征、项、号等。
在以下的描述中，将参考动作和操作的符号表示(例如以流程图形式)描述示例性实施例，所述操作是可以被作为包括用于执行特定任务或执行特定抽象数据类型的例程、程序、对象、组件、数据结构等被执行并且可以使用现有硬件程序模块或功能过程被执行的。这样的现有硬件可以包括一个或多个中央处理单元(CPU)、数字信号处理器(DSP)、特殊用途集成电路、现场可编程门阵列(FPGA)计算机等。
然而，应该记住的是，所有这些和类似术语与合适的物理量相关并且仅仅是应用于这些量的方便的标注。除非另有说明，或者从讨论中明显得出，诸如“处理”或“计算”或“确定”或“显示”等的术语指处理和转换在计算机系统的寄存器和存储器内部的代表物理量、电子量的数据的计算机系统或类似的电子计算设备的动作或处理。
还要注意，程序存储介质可以是磁的(例如软盘或硬驱)或光的(例如压缩盘只读存储器，或“CD ROM”)，并且可以是只读或随机存取的。类似地，传输介质可以是绞合线对、同轴电缆、光纤，或者一些其它对于现有技术所公知的合适的传输介质。示例性实施例不限于任何所给出的实施的这些方面。
术语模型可以指对象的衰减图集、解剖结构衰减模型、衰减参考图像的任意种类，或用来估计PET衰减图的任何其它参考。此外，可以实现可形变模型，其可以捕获衰减对象的所有合理的形状并且捕获在各个空间位置的所有可能的衰减值。可形变模型可以在形状和外观上变化。
术语模型的衰减外观模型指的是，在对象内部的一个空间位置上所有可能的衰减系数的集合，从而每个空间位置可以具有可能发生的一组潜在的衰减系数。衰减外观模型的一个实例是衰减系数的特别的设置，对于对象的每个空间位置一个。
图1示出了在示例性实施例中可以使用的用于重叠的MR和PET图像显示的设备1。设备1包括公知的MR扫描器2。MR扫描器2定义与图1的图平面正交地延伸的纵轴z。
如在图1中所示，在MR扫描器2内部同轴地设置具有多个围绕纵轴z成对地相对设置的PET探测单元3的PET单元。PET探测单元3优选地由具有在前连接的、由LSO晶体4组成的阵列的APD光电二极管阵列5以及电气放大电路(AMP)6组成。但是，本发明不限于具有APD光电二极管阵列5和在前连接的硅酸镥(LSO)晶体4的阵列的PET探测单元3，而是可以使用不同形成的光电二极管、晶体和装置用于探测目的。
由计算机7执行用于重叠的MR和PET图像显示的图像处理。
MR扫描器2沿着其纵轴z定义圆柱形的第一视野。多个PET探测单元3沿着纵轴z定义圆柱形的第二视野。根据本发明，PET探测单元3的第二视野基本上与MR扫描器2的第一视野一致。通过相应地匹配PET探测单元3沿着纵轴z的设置密度来实现这一点。
图2示出了使用用于重建PET图像的MLEM来估计完整PET衰减图的方法。MLEM可以是任何公知的MLEM。图2的方法可以在任何PET设备或诸如图1所示的设备1的具有PET模态的混合设备中执行。
如图2所示，在S200产生统计学模型。虽然为解释的目的使用统计学模型，然而应该理解的是，可被参数化的任何其它模型都可以被用在其它示例性实施例中。
通过执行由数据组的互相配准(coregistrations)产生的形变区域(deformation field)和衰减图的主成分分析可以构建统计学模型。可以从对多个个体的扫描以及或者简单的衰减图、或者MR图像数据和从其可以得到衰减图的图像的相应的对(例如来自每个个体的MR和CT图像对)中获得数据组。通过互相配准数据组，可以捕获形状(例如手臂)的统计学变化和衰减值。主成分分析允许待展开的参数空间的更紧凑的表示。
从诸如形变参数和衰减参数的输入数据的协方差矩阵的主轴变换(principal axis transformations)中可以获得主成分。主成分是输入数据的协方差矩阵的主本征矢量(principal Eigen vector)。将输入数据变换到主轴产生输入数据的紧凑的线性表示，从这些输入数据中通过线性组合可以产生模型的新实例。改变主成分的线性组合系数产生统计学模型的其它实例。
应该理解的是，不同于主成分分析的统计学分析方法，诸如聚类分析(clustering analysis)，可以被用来减少维数。
统计学模型可以是整个身体或任意身体部分。例如，统计学模型可以是与人的手臂的衰减图结合的运动学手臂模型或例如整个身体的统计学图集和统计学衰减图。统计学模型可以与在Rueckert et al.″Automatic Construction of 3DStatistical Deformation Models Using Non-rigid Registration.″Lecture Notes inComputer Science，vol.2208(2001)，77-84或者Fenchel et al.″Automatic Labelingof Anatomical Structures in MR FastView Images Using a Statistical Atlas.″LectureNotes in Computer Science，vol.5241(2008)，576-84中描述的模型类似，只是这些模型是基于灰度值图像而不是衰减图。
通过对于每个实例i的形变参数di和衰减外观参数(attenuation appearanceparameters)ai来参数化统计学模型。形变参数di参数化对象的形状。衰减外观参数ai参数化在空间位置上的衰减系数。衰减外观参数ai的例子是不同组织类型在其各个空间位置上在511keV的衰减值，例如肺组织衰减0.018/cm。
从互相配准的数据组获得形变参数di和衰减外观参数ai。然后计算在所有输入实例上的协方差矩阵。从参数的协方差矩阵提取主成分。然后可以通过线性模型描述统计学模型的一个实例：
μ = Σ i = 0 ( p i * w i ) - - - ( 1 )
其中，μ是统计学模型的实例，pi是主成分并且wi是对于线性方程中第i个主成分的系数。一般地，从主要值的三个求和(sigma)的区间选择系数。系数wi可以是对于l＜i＜m的形变参数di，并且wi可以是对于m+l＜i＜n的衰减外观参数ai。
因此，通过指定不同的系数可以创建任意实例。虽然上面描述了产生和参数化统计学模型的示例性实施例，然而应该理解的是，其它方法可以被用于其它模型。
包括诸如旋转和平移的空间变换参数的仿射参数(affine parameter)Ai，可以被用来在空间任意对准(align)和缩放统计学模型。而且，Ai可以被用来建立矩阵M和平移矢量t，通过它们每个空间位置在形变之后可以被转换为：
A(x)＝M*x+t    (2)
其中，x是空间位置的矢量。
应该理解的是，替代或除了形变参数、衰减外观参数和仿射参数之外，可以通过其它参数来参数化统计学模型。
基于统计学模型，在S210创建通过模型实例μ给出的PET衰减图(例如平均模型)。可以由计算机分析通过模型实例μ给出的PET衰减图。更具体地，用于模型矢量μ的衰减图可以是di、ai和Ai的函数并且定义为：
μ(di，ai，Ai)    (3)
L是发射图像(发射体分布)的对数似然性(log likelihood)L(λ，u)，其中λ是发射图像(正电子发射的空间分布)。发射图像λ基于从PET正弦图数据、例如通过反投影计算的初始发射体分布图像。如上所述，通过模型实例μ给出的PET衰减图是形变参数di和衰减外观参数ai和仿射参数Ai的函数。因此，
(λ，μ(di，ai，Ai))＝arg max(L(λ，μ(di，ai，Ai)))    (4)成为用于最大似然性的参数设置。此外，应该理解的是，发射图像λ还可以通过模型、例如统计学模型被参数化。此外，应当理解的是，其它测量可以被集成到扩展的似然性中。例如，如果模型参数的统计学分布是已知的或者可以被估计，则模型实例本身的似然性可以被集成到似然性测量中。
数据组的数量越大，则统计学模型越全面并且因此统计学模型越概括。应当理解的是，统计学模型是可形变模型(deformable model)的可能的实施例，并且统计学模型的每个实例i是对于该实例的模型参数的函数。例如，形变参数di和衰减外观参数ai和放射参数Ai是模型参数。
在S215采集在PET单元的视野内的对象的PET正弦图数据。可以通过在图1中示出的PET单元采集PET正弦图数据。基于PET正弦图数据，与模型参数同时地计算发射图像λ(PET图像)。替换地，可以交替地通过首先保持发射图像λ为固定并且更新统计学模型，然后保持统计学模型为固定并且更新发射图像，来计算发射图像λ和统计学模型。
在S220基于PET正弦图数据优化发射图像λ和对于该实例的模型参数。以迭代方式优化发射图像λ和对于该实例的模型参数。在优化过程中，计算发射图像λ，更新统计学模型并且重新计算发射图像λ，直到达到了优化。模型参数可以是形变参数di、衰减外观参数ai和仿射参数Ai。
在S220，基于MLEM函数同时地重建通过模型实例μ给出的PET衰减图和发射图像λ。图1所示的计算机可以重建发射图像λ以及通过模型实例μ给出的PET衰减图。可以基于如下定义的(未截取)所测量的PET发射的对数似然性重建和优化发射图像λ以及通过模型实例μ给出的PET衰减图：
L = ( λ , μ ( d i , a i , A i ) ) = Σ i = 0 ( y i * log ( y i ′ ) - y i ′ ) - - - ( 5 )
其中，yi是所测量的PET正弦图数据，yi′是所估计的yi值。通过前向投影发射图像λ以及校正通过模型给出的衰减图的衰减来估计yi′。然后可以作为参数di、ai、Ai和λ的函数来优化L。在此，将MLEM约束到模型的参数空间是用来估计整个衰减图。参数空间是指参数化的统计学衰减模型的所有可能的值。在图2所示的示例性实施例中，参数空间可以包括所有模型参数，例如，所有形变、衰减外观和仿射变换参数。
以迭代方式优化模型参数，直到达到了最大似然性。可以通过任意通常的优化算法求得最佳值。最佳值是这样的参数设置，对于该参数设置，最大似然性达到最大值。最佳值定义了用于衰减和发射体图像的参数化模型的最可能的实例。
当在S220优化算法收敛时，已经达到了最佳的最大似然性。然后，在S225直接从模型实例获得衰减图并且从发射图像获得PET图像。应该理解的是，当使用衰减图触发了另一个PET图像重建时，可以舍弃在S225获得的发射图像。
图3示出了按照一种示例性实施例通过使用基于MR的衰减图和用于重建PET图像的MLEM来估计完整的PET衰减图的方法。可以在任意PET设备或诸如图1所示的设备1的具有PET模态的混合设备中执行图3的方法。
如图3所示，在S300产生统计学模型。S300与S200相同。因此，为清楚和简洁起见省略S300的详细描述。
在S305a，在MR单元的视野内通过MR扫描对象以采集MR数据。图1所示的MR扫描器可以被用来采集MR数据。
一旦扫描了对象，则在S305b产生基于MR的衰减图。基于MR的衰减图可以通过产生基于MR的衰减图的任意已知方法来产生并且可以通过图1所示的计算机来产生。
根据基于MR的衰减图和统计学模型，在S310创建对于模型实例μ的PET衰减图。S310与S210相同，只是用于模型实例μ的PET衰减图通过基于MR的衰减图限制。因此，为清楚和简洁起见省略S310的详细描述。
在S315采集在PET单元的视野内的对象的PET正弦图数据。S315与S215相同。
在S320，优化用于实例的模型参数。S320与S220相同。优化模型参数，直到已经达到最大的似然性。
一旦已经达到了对于模型实例μ的PET衰减图和发射图像λ的最大似然性，在S325对于模型实例μ的衰减图就定义了对于模型实例μ的最佳衰减图。重建的PET图像也是在S325产生，但是当对于模型实例μ的最佳衰减图可以被用在用于获得PET图像的另一个重建过程时，可以舍弃重建的PET图像。组合基于MR的衰减图和PET衰减图的方法是现有技术中公知的。因此，为清楚和简洁起见不再讨论。
图4示出了按照一种示例性实施例通过使用MLEM来细化(refining)初始衰减图的方法。可以在任意PET设备或诸如图1所示的设备1的具有PET模态的混合设备中执行图4的方法。如图4所示，在S400以图2和3中相同的方式产生统计学模型。
在S402，将被参数化的统计学模型适应于初始衰减图。可以事先从例如低分辨率的MR图像或传输扫描(transmission scan)中产生初始衰减图。在另一个示例性实施例中，包括频繁地延伸到MR视野外部的解剖结构的参数化的模型(例如人手臂的动力学模型)可以被添加到初始衰减图或可以被用来完善初始衰减图。
在S405，通过MR单元扫描在MR单元的视野内部的对象以产生MR图像。在S410创建模型实例。关于最小平方途经基于初始衰减图的最佳适应(best fit)来创建模型实例。对于统计学模型，可以通过执行例如初始衰减图到统计学模型的线性空间的正交投影来计算所述最佳适应。因此，所创建的模型是缩放到初始衰减图的统计学模型的平均实例。
在S415，扫描PET正弦图数据然后在S420优化对于模型实例μ的衰减图的模型参数。S415和S420与S215和S220相同。因此，为清楚和简洁起见，不详细描述S415和S420。在S420，通过使用MLEM细化初始衰减图。
在S425，基于对于模型实例μ的衰减图的最大似然性产生细化的衰减图。也产生PET图像。
如上所述，这些方法可以被用来基于使用MLEM重建的模型，来估计对象的完整衰减图，和/或被用来完善在通过模型方式之前利用其它方法计算的衰减图的缺失部分。此外，示例性实施例可以被用来细化(包括模型的初始化和对数据的细化)衰减图，诸如，从基于MR的衰减图计算方法计算的衰减图。示例性实施例帮助避免局部最大并且在其参数空间、图集或参考图像内产生有效和有意义的模型实例。
上面描述了示例性实施例，很显然可以按照多种方式改变上述示例性实施例。这样的变化不被认为是脱离了所述示例性实施例的精神和范围，并且对于本领域技术人员是明显的所有这样的修改旨在被包括在所述示例性实施例的范围之内。

Claims (11)

1.一种在磁共振(MR)扫描器和正电子发射断层造影(PET)单元中校正衰减的方法，该方法包括：
采集在PET单元的视野内部的对象的PET正弦图数据；
采集在MR扫描器的视野内部的对象的MR数据；
基于参数化的模型实例的最大似然性期望最大化(MLEM)以及PET正弦图和MR数据，产生衰减图，所述MLEM通过参数化的模型实例的模型参数而限制。
2.根据权利要求1所述的方法，其中，产生衰减图包括产生PET衰减图。
3.根据权利要求2所述的方法，其中，产生PET衰减图包括
产生模型，
通过模型参数来参数化所述模型，和
基于参数化的模型创建参数化的模型实例。
4.根据权利要求3所述的方法，其中，所述模型是统计学模型。
5.根据权利要求4所述的方法，其中，所述模型参数是形变参数和衰减外观参数。
6.根据权利要求5所述的方法，其中，所述模型参数还包括仿射参数。
7.根据权利要求2所述的方法，还包括：
产生MR衰减图。
8.根据权利要求1所述的方法，其中，产生衰减图包括，
产生模型，和
将所述模型适应于初始衰减图，基于所适应的模型创建参数化的模型实例。
9.根据权利要求1所述的方法，其中，所述PET单元被配置为采集PET正弦图数据。
10.根据权利要求1所述的方法，其中，连接到所述MR扫描器和所述PET单元的计算机被配置为产生所述衰减图。
11.根据权利要求1所述的方法，其中，所述MR扫描器被配置为采集MR数据。
CN201010232204.XA 2009-07-16 2010-07-16 完整或部分正电子发射断层造影衰减图的基于模型的估计 Active CN101953693B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US12/458,575 US8299438B2 (en) 2009-07-16 2009-07-16 Model based estimation of a complete or partial positron emission tomography attenuation map using maximum likelihood expectation maximization
US12/458,575 2009-07-16

Publications (2)

Publication Number Publication Date
CN101953693A CN101953693A (zh) 2011-01-26
CN101953693B true CN101953693B (zh) 2015-02-18

Family

ID=43465884

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010232204.XA Active CN101953693B (zh) 2009-07-16 2010-07-16 完整或部分正电子发射断层造影衰减图的基于模型的估计

Country Status (2)

Country Link
US (1) US8299438B2 (zh)
CN (1) CN101953693B (zh)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102010004384B4 (de) * 2010-01-12 2012-03-08 Siemens Aktiengesellschaft Verfahren zur Ermittlung von der Berechnung eines Bestrahlungsplans zugrunde zu legenden Informationen und kombinierte Magnetresonanz-PET-Vorrichtung
DE102010005287A1 (de) * 2010-01-21 2011-07-28 Siemens Aktiengesellschaft, 80333 Verfahren zur Aufnahme und Auswertung von mit einer kombinierten Magnetresonanz-PET-Vorrichtung gleichzeitig mit Magnetresonanzdaten aufgenommen PET-Daten und kombinierte Magnetresonanz-PET-Vorrichtung
CN103415779B (zh) * 2011-03-07 2016-05-11 皇家飞利浦有限公司 利用混合式核成像/mr中的核发射数据进行mr分割的系统及方法
CN103186882B (zh) * 2011-12-30 2015-09-30 沈阳东软派斯通医疗系统有限公司 Pet系统中图像的衰减校正方法及装置
DE102012201412B4 (de) * 2012-02-01 2023-01-12 Siemens Healthcare Gmbh Verfahren zum Berechnen eines Wertes eines Absorptionsparameters der Positronen-Emissions-Tomographie, Verfahren zur Positronen-Emissions-Tomographie, Magnetresonanzanlage und Positronen-Emissions-Tomograph
CN104271045B (zh) * 2012-05-04 2018-10-02 皇家飞利浦有限公司 具有正电子发射断层摄影中的散射的符合的衰减图
WO2014074148A1 (en) * 2012-11-08 2014-05-15 The General Hospital Corporation System and method for multi-modality time-of-flight attenuation correction
US9082193B2 (en) * 2013-04-19 2015-07-14 Siemens Medical Solutions Usa, Inc. Shape-based image segmentation
EP2994776B1 (en) 2013-05-08 2019-11-13 Koninklijke Philips N.V. Apparatus and method for the evaluation of gamma radiation events
CN106963406B (zh) * 2013-10-30 2020-06-02 上海联影医疗科技有限公司 衰减图及正电子发射断层图像生成装置
EP3125765A4 (en) 2014-03-04 2018-04-18 The Trustees of Columbia University in the City of New York Regularization of images
US9684973B2 (en) 2014-12-08 2017-06-20 General Electric Company Systems and methods for selecting imaging data for principle components analysis
US9466132B2 (en) 2014-12-08 2016-10-11 General Electric Company Systems and methods for motion mitigation determinations
WO2016197127A1 (en) * 2015-06-04 2016-12-08 Rensselaer Polytechnic Institute Attenuation map reconstruction from tof pet data
US9858691B2 (en) 2015-07-31 2018-01-02 The Board Of Trustees Of The Leland Stanford Junior University Statistical data acquisition model for GPU based MLEM joint estimation of tissue activity distribution and photon attenuation map from PET data
WO2018220686A1 (ja) * 2017-05-29 2018-12-06 株式会社島津製作所 吸収係数画像推定方法、吸収係数画像推定プログラム並びにそれを搭載したポジトロンｃｔ装置
CN107610198B (zh) * 2017-09-20 2021-05-14 赛诺联合医疗科技（北京）有限公司 Pet图像衰减校正方法及装置

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5903008A (en) * 1997-07-02 1999-05-11 General Electric Company Scatter correction methods and systems in single photon emission computed tomography
CN1220009A (zh) * 1996-03-27 1999-06-16 阿达克实验室 医学成象系统中多形态数据的互补重合
US6310968B1 (en) * 1998-11-24 2001-10-30 Picker International, Inc. Source-assisted attenuation correction for emission computed tomography
CN1445560A (zh) * 2002-03-20 2003-10-01 株式会社日立制作所 放射成象装置、放射成象方法和放射成象支持方法
US6740883B1 (en) * 1998-08-14 2004-05-25 Robert Z. Stodilka Application of scatter and attenuation correction to emission tomography images using inferred anatomy from atlas
CN1903128A (zh) * 2005-07-25 2007-01-31 株式会社岛津制作所 正电子ct设备的3d图像重构方法和正电子ct设备
CN101152084A (zh) * 2006-09-29 2008-04-02 西门子公司 组合式磁共振-正电子发射断层成像设备的磁场发生单元
CN101352348A (zh) * 2007-07-26 2009-01-28 西门子公司 考虑运动过程的患者测量数据拍摄方法及有关医疗设备

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE60020214T2 (de) * 1999-01-11 2006-03-23 British Nuclear Fuels Plc, Warrington Vorrichtung und verfahren zur untersuchung radioaktiver strahlenquellen in einer probe
US7105824B2 (en) * 2002-05-09 2006-09-12 Neurologica, Corp. High resolution photon emission computed tomographic imaging tool
AU2003276658A1 (en) * 2002-11-04 2004-06-07 V-Target Technologies Ltd. Apparatus and methods for imaging and attenuation correction
US7599540B2 (en) * 2003-06-18 2009-10-06 Koninklijke Philips Electronics N.V. Motion compensated reconstruction technique
US8423125B2 (en) * 2004-11-09 2013-04-16 Spectrum Dynamics Llc Radioimaging
US7534418B2 (en) * 2004-12-10 2009-05-19 The Regents Of The University Of Michigan Imaging agents
JP2008523871A (ja) * 2004-12-15 2008-07-10 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ マルチモダリティ画像のレジストレーション
ATE465425T1 (de) * 2004-12-17 2010-05-15 Koninkl Philips Electronics Nv Kompensationsalgorithmus bei abschneidung für iterative rekonstruktion
EP1861824B1 (en) * 2005-03-16 2018-09-26 Koninklijke Philips N.V. Method and device for the iterative reconstruction of tomographic images
RU2427037C2 (ru) * 2005-10-05 2011-08-20 Конинклейке Филипс Электроникс Н.В. Способ и система для реконструкции рет-изображения с использованием суррогатного изображения
WO2007054843A1 (en) * 2005-11-10 2007-05-18 Koninklijke Philips Electronics, N.V. Pet imaging using anatomic list mode mask
US7348564B2 (en) * 2005-12-12 2008-03-25 General Electric Company Multi modality imaging methods and apparatus
US7756310B2 (en) * 2006-09-14 2010-07-13 General Electric Company System and method for segmentation
US7813783B2 (en) * 2006-11-02 2010-10-12 General Electric Company Methods and systems for attenuation correction in medical imaging
US20080135769A1 (en) 2006-11-22 2008-06-12 Rosen Bruce R Attenuation correction of pet image using image data acquired with an mri system
CN101681520B (zh) * 2007-05-30 2013-09-25 皇家飞利浦电子股份有限公司 Pet局部断层摄影
EP2316043A1 (en) * 2008-08-15 2011-05-04 Koninklijke Philips Electronics N.V. Attenuation correction for pet or spect nuclear imaging systems using magnetic resonance spectroscopic image data
US7888632B2 (en) * 2008-09-22 2011-02-15 Siemens Medical Solutions Usa, Inc. Co-registering attenuation data and emission data in combined magnetic resonance/positron emission tomography (MR/PET) imaging apparatus
DE102010004384B4 (de) * 2010-01-12 2012-03-08 Siemens Aktiengesellschaft Verfahren zur Ermittlung von der Berechnung eines Bestrahlungsplans zugrunde zu legenden Informationen und kombinierte Magnetresonanz-PET-Vorrichtung

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1220009A (zh) * 1996-03-27 1999-06-16 阿达克实验室 医学成象系统中多形态数据的互补重合
US5903008A (en) * 1997-07-02 1999-05-11 General Electric Company Scatter correction methods and systems in single photon emission computed tomography
US6740883B1 (en) * 1998-08-14 2004-05-25 Robert Z. Stodilka Application of scatter and attenuation correction to emission tomography images using inferred anatomy from atlas
US6310968B1 (en) * 1998-11-24 2001-10-30 Picker International, Inc. Source-assisted attenuation correction for emission computed tomography
CN1445560A (zh) * 2002-03-20 2003-10-01 株式会社日立制作所 放射成象装置、放射成象方法和放射成象支持方法
CN1903128A (zh) * 2005-07-25 2007-01-31 株式会社岛津制作所 正电子ct设备的3d图像重构方法和正电子ct设备
CN101152084A (zh) * 2006-09-29 2008-04-02 西门子公司 组合式磁共振-正电子发射断层成像设备的磁场发生单元
CN101352348A (zh) * 2007-07-26 2009-01-28 西门子公司 考虑运动过程的患者测量数据拍摄方法及有关医疗设备

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Corrections for Accidental Coincidences and Attenuation in Maximum-Likelihood Image Reconstruction for Positron-Emission Tomography;David G.Politte等;《IEEE TRANSACTIONS ON MEDICAL IMAGING》;19910331;第10卷(第1期);第82-89页 *
Maximum-Likelihood Expectation-Maximization Reconstruction of Sinograms with arbitrary Noise Distribution Using NEC-Transformations;J.Nuyts等;《IEEE TRANSACTIONS ON MEDICAL IMAGING》;20010531;第20卷(第5期);第365-375页 *
Penalized-Likelihood Sinogram Restoration for Computed Tomography;Patrick J.La Riviere等;《IEEE TRANSACTIONS ON MEDICAL IMAGING》;20060831;第25卷(第8期);第1022-1036页 *
基于非线性统计估计技术的PET散射校正;吴朝霞等;《生物医学工程学杂志》;20051231;第22卷(第1期);第74-77页 *

Also Published As

Publication number Publication date
US8299438B2 (en) 2012-10-30
US20110015904A1 (en) 2011-01-20
CN101953693A (zh) 2011-01-26

Similar Documents

Publication Publication Date Title
CN101953693B (zh) 完整或部分正电子发射断层造影衰减图的基于模型的估计
Arabi et al. The promise of artificial intelligence and deep learning in PET and SPECT imaging
CN101681520B (zh) Pet局部断层摄影
US8620053B2 (en) Completion of truncated attenuation maps using maximum likelihood estimation of attenuation and activity (MLAA)
EP3226766B1 (en) System and method for image calibration
CN109961834B (zh) 影像诊断报告的生成方法及设备
US7948503B2 (en) Techniques for 3-D elastic spatial registration of multiple modes of measuring a body
JP5486182B2 (ja) 情報処理装置および情報処理方法
RU2413245C2 (ru) Обработка изображений позитронно-эмиссионной томографии с использованием анатомической маски режима списка
US7349563B2 (en) System and method for polyp visualization
US20100027861A1 (en) Segmentation of regions in measurements of a body based on a deformable model
CN102460509B (zh) 基于图像信息建立结构的轮廓
US7304644B2 (en) System and method for performing a virtual endoscopy
CN101739708B (zh) 医用图像处理装置、超声波图像取得装置和医用图像处理方法
CN106466188A (zh) 用于发射断层显像定量的系统和方法
CN101204326B (zh) 用于运行混合医学成像单元的方法
EP1506742A1 (en) Fusion of structural images in PET-computer tomography
CN101017569A (zh) 用于减小断层造影图像数据组中的噪声的方法
CN106575441A (zh) 重构内运动校正
CN101352348A (zh) 考虑运动过程的患者测量数据拍摄方法及有关医疗设备
Apostolopoulos et al. Applications of generative adversarial networks (GANs) in positron emission tomography (PET) imaging: A review
CN108553121A (zh) 一种pet延迟扫描的方法和装置
US7907756B2 (en) System and method for validating an image segmentation algorithm
US20110081054A1 (en) Medical image analysis system for displaying anatomical images subject to deformation and related methods
CN102028484B (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
TR01 Transfer of patent right

Effective date of registration: 20220130

Address after: Erlangen

Patentee after: Siemens Healthineers AG

Patentee after: SIEMENS medical solutions company

Address before: Munich, Germany

Patentee before: SIEMENS AG

Patentee before: SIEMENS medical solutions company

TR01 Transfer of patent right