CN106062818A - 用于时间选通pet研究的基于选通特定的mr的衰减校正的方法和装置 - Google Patents

用于时间选通pet研究的基于选通特定的mr的衰减校正的方法和装置 Download PDF

Info

Publication number
CN106062818A
CN106062818A CN201480070827.8A CN201480070827A CN106062818A CN 106062818 A CN106062818 A CN 106062818A CN 201480070827 A CN201480070827 A CN 201480070827A CN 106062818 A CN106062818 A CN 106062818A
Authority
CN
China
Prior art keywords
image
baseline
images
pet
time
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.)
Pending
Application number
CN201480070827.8A
Other languages
English (en)
Inventor
C·W·斯蒂尔恩斯
S·D·沃伦韦伯
D·M·戈德哈伯
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.)
General Electric Co
Original Assignee
General Electric Co
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 General Electric Co filed Critical General Electric Co
Publication of CN106062818A publication Critical patent/CN106062818A/zh
Pending legal-status Critical Current

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • 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
    • A61B6/5282—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to scatter
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/0522—Magnetic induction tomography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
    • A61B5/085—Measuring impedance of respiratory organs or lung elasticity
    • 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]
    • 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/037—Emission tomography
    • 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/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
    • 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
    • A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
    • 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
    • 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]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/32—Determination of transform parameters for the alignment of images, i.e. image registration using correlation-based methods
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
    • G01T1/1611—Applications in the field of nuclear medicine, e.g. in vivo counting using both transmission and emission sources sequentially
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
    • G01T1/2914—Measurement of spatial distribution of radiation
    • G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10084—Hybrid tomography; Concurrent acquisition with multiple different tomographic modalities
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10088—Magnetic resonance imaging [MRI]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10104—Positron emission tomography [PET]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20048—Transform domain processing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20172—Image enhancement details
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing

Landscapes

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

Abstract

一种用于移动结构的多个时间选通的PET图像的衰减校正的方法包括：基于结构的基线MR图像建立ACF矩阵；展开多个图像变换，每一个图像变换将基线MR图像与移动结构的相应电影MR图像配准；将每一个图像变换应用于ACF矩阵以生成对应的电影ACF矩阵；将每一个电影MR图像及其对应的电影ACF矩阵时间关联到多个时间选通PET数据集之一；及从每一个时间选通PET数据集产生衰减校正的时间选通PET图像。

Description

用于时间选通PET研究的基于选通特定的MR的衰减校正的方 法和装置
技术领域
本发明的实施例一般涉及增强从PET-MR扫描器获得图像的质量。具体实施例涉及PET图像的衰减校正。
背景技术
正电子发射断层扫描(“PET”)机使用闪烁器或其他检测器的一个或多个环来从伽马射线(光子对)生成电信号，伽马射线从由示踪化合物中封装的放射性核素的衰变发出的正电子与在目标材料内的电子的重组产生。通常，重组事件在距离放射性核素衰变事件的约1mm内发生，并且重组光子在通常相反方向上发射以到达不同的检测器。将在检测窗口(通常小于几纳秒间隔)内发生的成对光子到达视为指示重组事件，并在此基础上将计算机断层扫描算法应用于闪烁器位置和检测数据，以便定位各种重组事件，从而产生目标材料内示踪剂处置的三维图像。
通常，目标材料是身体组织，示踪化合物是类似于生物流体的液体，并且利用生物流体将放射性核素主要设置在身体组织中。例如，常见形式的PET利用脱氧葡萄糖(18F)，其类似于葡萄糖，以18F放射性核素代替了通常组成葡萄糖的一个羟基。脑物质、肾脏和增殖细胞(例如转移性癌细胞)优先吸收葡萄糖和脱氧葡萄糖。因此，PET在肿瘤学研究中非常有用，用于定位特定器官并用于研究代谢过程。
在获得所希望的PET图像质量中的一个难题是在由正电子-电子相互作用产生的能量谱中的伽马射线容易被典型的身体组织衰减并且被不同身体组织不同衰减。变化的衰减将改变从患者检测到重组事件的可能性，从而扰乱从PET数据得到图像的过程。因此，非常希望提供用于衰减校正(“AC”)的部件。
例如，PET常常与计算机断层扫描(“CT”)结合，其使用移动的X射线源和检测器以获得内部结构的图像。X射线是光子，并且被衰减，非常类似于从正电子-电子重组事件产生的更高能量的光子。然而，与PET光子不同，CT光子的源强度是已知的。因而，CT图像数据是在源和检测器之间的光子衰减的直接测量。因此，CT图像数据为并行的PET成像提供了有用的基础。
组合的PET-CT扫描已经证明在呼吸选通的研究中特别有用，例如在正常呼吸期间为了诊断或评估肺癌而进行的研究。呼吸选通的图像可以减少运动模糊，它将会以另外方式在整个呼吸周期由PET图像获取产生。然而，呼吸选通的PET研究准确重构需要特定门的衰减校正信息。CT和PET数据的同时获取支持特定门的AC。
例如，考虑PET检测器的两个元件，其中，在检测器之间的响应线穿过患者的下胸部。在峰值吸气，该线可以通过下叶肺，其会相对少量地衰减511keV的光子。随着患者呼气，肝脏上移进入胸腔，使得在呼气结束，相同的响应线现在可以通过软组织，其对511keV光子的衰减比肺多的多。如果对PET数据进行选通，使得在峰值吸气和呼气结束获取的不同图像，每一个都将会从对应于在呼吸周期的该阶段期间衰减组织的实际分布的衰减校正获益。
逐渐地且出于包括寿命辐射剂量减少目标的各种原因，PET/CT扫描正被组合PET扫描与磁共振成像(“MR”)所代替。这个新组合(“PET-MR”)提出了新的技术问题。例如，尽管CT基于从源通过目标发射的X射线的检测来构成图像，但MR基于旋转的“弛豫”电磁场的检测来构成图像，该电磁场是响应于施加磁场波动而由具有奇数原子数量(即中子和质子的总数不能由二整除)的原子核在目标内产生的。因而，MR测量与由CT和PET检测的光子根本上不同的现象。
MR的一个优点是：磁场不会在人体组织中衰减，使得，可以仅根据所施加的磁场和响应场之间的频率偏移来确定(使用傅立叶分析)原子核的位置。另一个优点是：通过仔细选择脉冲序列，不同的组织或材料可以被高亮显示。因此，MR经常被用于区分患者内的组织类型，并且也用于识别细微的细节结构。通常，不同的脉冲序列被用于组织区分。例如，T1脉冲序列可以用于获得具有显现为较暗的水和较亮的脂肪的图像。另一方面，T2脉冲序列可用于获得具有较暗的脂肪和较亮的水的图像。
因此，结合PET和MR的单一装置(“PET-MR扫描器”)可以提供细微细节，组织区分和代谢数据。然而，因为MR信号与PET或CT信号衰减不同方式衰减，并且因为MR信号返回高度依赖于使用的脉冲序列的类型(其中每一个脉冲序列突出不同的材料)，而介入在重组事件和一对检测器之间的每一层的材料衰减PET信号，单扫描MR图像数据不一定为并行PET成像的AC提供可靠的基础。
MR成像具有以称为回顾性选通电影-MRI的方式获得呼吸周期上的大量图像的能力。然而，可以以电影方式获取的MR图像的类型被认为不适合于产生PET衰减校正值。常规视为对于产生PET衰减校正值有效的MR协议，例如2-点狄克逊水-脂肪扫描，通常需要超过10秒来获得，因此，不可用于在呼吸周期期间捕获多个图像。实际上，为了使得运动模糊和其他的呼吸伪影最小，这些MR图像数据通常是在长屏气期间获得的。然而，屏气的MR图像不示出中间肺位置，其是对于PET扫描的准确AC所需要的。
发明内容
本发明的实施例提供了一种用于移动结构的多个时间选通的PET图像的衰减校正的方法。该方法包括：基于结构的基线MR图像建立ACF矩阵；展开多个图像变换，每一个图像变换将基线MR图像与移动结构的相应电影MR图像配准；将每一个图像变换应用于ACF矩阵以生成对应的电影ACF矩阵；将每一个电影MR图像及其对应的电影ACF矩阵时间关联到多个时间选通PET数据集之一；及从每一个时间选通PET数据集产生衰减校正的时间选通PET图像。
本发明的方面提供了一种用于生成多个衰减校正的时间选通PET图像的装置。装置包括PET-MR扫描器；与PET-MR扫描器通信连接的控制器，及与控制器通信连接的图像处理器。控制器配置成获取结构的基线MR图像，获取结构的电影MR图像序列，及获取结构的多个时间选通PET数据集。图像处理器配置成基于基线MR图像生成基线ACF矩阵，将基线MR图像与电影MR图像的至少一些配准，基于基线MR图像与电影MR图像的配准生成电影ACF矩阵，及基于电影ACF矩阵衰减校正时间选通PET集。
本发明的其他方面提供了一种用于PET成像的衰减校正的方法。该方法包括：在第一MR获取中，获得适合于导出PET衰减校正因子的基线MR图像；及在第二MR获取中，与呼吸选通PET研究并行地获得呼吸周期期间的MR图像的电影序列。该方法还包括将基线MR图像与每一个电影MR图像非刚性配准，以便产生一组配准的MR图像，其对应于来自呼吸选通PET研究的PET数据集，且可以用于导出衰减校正因子。配准的MR图像随后用于为它们相应的时间选通PET数据集形成衰减校正因子。
附图说明
通过参考附图阅读下面非限制性实施例的描述，更好地理解本发明，在以下附图中：
图1示出用于基于并行CT图像的PET图像的衰减校正的常规方法。
图2示出根据本发明第一实施例的用于基于并行和基线MR图像的PET图像的衰减校正的发明方法。
图3示出根据本发明实施例的用于将基线MR图像与电影MR图像配准的示范性过程。
图4示出用于基线MR图像与电影MR图像的互信息的等式。
图5示出根据本发明第二实施例的用于基于并行和基线MR图像的PET图像的衰减校正的另一个方法。
图6示出根据本发明实施例的用于实现基于并行和基线MR图像的PET图像的衰减校正的的发明方法的装置。
具体实施方式
以下将详细参考本发明的示范性实施例，在附图中图示其示例。在可能的情况下，在全部附图中使用的相同参考字符指代相同或相似的部分，不再重复描述。相对于组合的PET-MR扫描器来描述本发明的示范性实施例，尽管实施例可以适于与其他成像系统一起使用。
本发明的方面涉及通过提供基于MR图像数据的PET图像的衰减校正(AC)来改进PET和MR图像的同时收集。
在PET/CT研究中，呼吸选通衰减的传统方法100(如图1所示)包括在与获取120时间选通PET数据集122的序列相同的解剖区上和在相同的时间段期间，获得110电影CT图像112的序列。可以在每秒一个图像的数量级的速率获取电影CT图像112的序列，允许在典型的呼吸周期期间6到10个图像的获取，使得CT图像的集合随后可以与PET原始数据文件的选通集合120相关联。CT图像110的序列用于导出130衰减校正因子矩阵(“ACF矩阵”)132的序列，并且随后每一个ACF矩阵与时间选通PET数据集122的一个时间关联140，以便重构150衰减校正(“AC”)的PET图像152的集合。
图2示出根据本发明的第一实施例的方法200，用于呼吸选通PET数据集222的序列的基于MR的衰减校正210，以产生衰减校正的PET图像212。除了获取220PET数据集222以外，方法200还包括获取230区分组织的或“基线”MR图像232，其适合于稍后生成240光子衰减校正因子的三维矩阵242(“基线ACF矩阵”)。在实施例中，使用患者屏气协议和适用于组织区分研究的扫描协议-例如2点狄克逊扫描，其清晰地高亮显示水和脂肪并需要超过十秒完成，作为基本上静止的图像获得230基线MR图像232。虽然，在某些实施例中，可以使用在正常呼吸期间的组织区分研究作为平均图像获得，或者在呼吸期间或在屏气期间作为复合图像从几个不同单一组织扫描编译基线MR图像232，但可以预计到平均图像或复合图像对于随后生成静止ACF矩阵242不是最佳的。
基于基线MR图像232中识别的组织，针对典型正电子重组能量进行标准光子衰减计算以便生成240基线ACF矩阵242。
除了生成240基线ACF矩阵242以外，方法200还包括获取250电影MR图像的252的序列，以及将每一个电影MR图像252配准260回静止MR图像232。在配准260的过程中，产生多个图像变换262，每一个图像变换都将相应的一个电影MR图像252与静止MR图像232配准。在一个实施例中，图像变换是非刚性的。例如，Vemuri，et al.，Med.Image Anal.7：1-20(2003)中描述了用于非刚性图像配准“水平集”算法。也可以使用其他基于互信息的算法(例如，如Woods(1992)或Viola和Wells(1994)所教导的)。在大多数实施例中，使用不同的(更快)用于获得基线MR图像232的扫描协议-例如，零-TE或超短-TE扫描协议获得电影MR图像252。而且，在某些实施例中，不同扫描协议可以用于不同的电影MR图像252。因此，配准260还可以包括像素强度缩放或移动以在不同扫描协议之间匹配组织类型。
例如，如图3所示，可以进行基于互信息的局部扰动(非刚性)刚性体配准260，其中，一个电影图像252表示为可微的三次样条函数fT(x)，将像素强度映射到线性像素阵列，其可以在组成图像的采样值的集合内的任何像素间位置内插。类似地，基线图像232表示为另一个可微的三次样条函数fR(x)。平滑联合直方图
μ＝arg min S(fR，fT og(·|μ))
(其将电影图像252有关回基线图像232)定义为两个样条函数的向量积。通过将变换参数的集合μ＝{γ，θ，φ，tx，ty，tz；δj}应用于基线图像232来变换平滑联合直方图，选择变换参数的集合以根据分级多分辨率优化方案优化互信息S(μ)。在变换参数μ的集合中，{γ，θ，φ}是体横摇-俯仰-偏航欧拉角，[tx，ty，tz]是体平移矢量，以及δj是变形系数的集合，每一个变形系数对应于多个控制点之一，控制点由在基线图像232上覆盖的稀疏网格定义。通常，稀疏网格是规则的网格。然而，在其他实施例中，可以根据像素强度或其他图像度量加权稀疏网格。对于非刚性变换参数的任何给定集合，基于分别对应于平滑联合直方图、用于电影图像252的边缘平滑直方图和用于基线图像232的边缘平滑直方图的概率分布P，pT和pR，计算互信息S(μ)(如图4所示)，以及计算互信息梯度概率分布被假定为与变换参数无关。图4中所示的参数τ和κ是分别用于电影图像252和基线图像232的像素强度直方图二进制索引(bin indexes)。改变直方图单元大小可以调整图像配准的质量。配准被描述为局部扰动刚性体，因为它将局部变形系数与刚性体变换{γ，θ，φ，[tx，ty，tz]}相结合。为了满意的S(μ)值，评估配合质量，配准过程260在修正的μ上迭代，或者返回μ(将基线图像232配准到特定电影图像252的图像变换262)。
随后，方法200包括将每一个图像变换262应用于静止ACF矩阵242以便生成270多个电影ACF矩阵272。每一个电影ACF矩阵272时间关联280到在呼吸选通PET序列220内的多个PET数据集222之一，并应用210到对应的PET数据集222以便产生重构和衰减校正的PET图像212。在某些实施例中，每一个PET数据集222时间关联到多于一个电影ACF矩阵，并且在将多个电影ACF矩阵272应用210于PET数据集222中使用简单平均、时间加权平均、速度加权平均或时间-速度加权平均中的一个。
尽管将静止MR图像获取230示为在电影MR序列获取250之前，但以其他顺序同样可以完成这些步骤。在某些实施例中，可以复制电影MR序列(例如，在屏气静止图像之前和之后)或可以复制屏气静止图像(例如，在电影序列之前和之后)。类似地，配准260可以与获取250并行完成，或在获取了所有图像之后。
尽管将多个配准的ACF矩阵272描述为通过将图像变换262应用于基线ACF矩阵242而获得，但配准ACF矩阵272也可以从基线MR图像232计算，如同按照多个图像变换262所变换的。因而，根据本发明其他实施例的如图5中所示的备选的方法500包括第一MR获取230，以获得静止(屏气)MR图像232，其适合于导出240基线PET衰减校正因子242；及第二MR获取250，以与呼吸选通PET研究220并行地获得呼吸周期期间获取的电影MR图像252的序列。将静止MR图像232与每一个电影MR图像252非刚性配准，以便产生560一组变换的MR图像562，其对应于来自呼吸选通PET研究的PET数据集222，且可以用于导出570衰减校正因子572。配准的MR图像随后用于形成570衰减校正因子572，它们是时间关联280的，用于从它们相应的时间选通PET数据集222产生210衰减校正的PET图像212。
图6示意性地示出根据本发明实施例的装置600，用于完成如上所述的方法200。装置600包括组合的PET-MR扫描器610，其与控制器620通信连接，控制器又与图像处理器630通信连接。控制器和图像处理器二者都可以实现为单一计算机中单独的软件或者单独在与PET-MR扫描器610共同定位的计算机中；或者控制器620和图像处理器630的任意一个或二者可以与PET-MR扫描器和/或彼此非现场(远离)定位。控制器620配置成完成至少获取220PET数据集222，获取230基线MR图像232，及获取250电影MR图像252的序列。图像处理器630配置成完成至少生成240基线ACF矩阵242，将基线MR图像232与电影MR图像252的序列的每一个配准260，生成270多个电影ACF矩阵272，及基于多个电影ACF矩阵272衰减校正210时间选通PET数据集222。控制器620或图像处理器630可以配置成完成将电影MR图像252与时间选通PET数据集222时间关联280。备选地，控制器620和图像处理器630可以集成到单一软件/单一处理器中；或者可以在云(SaaS)范例中实现。
因而，本发明的实施例提供了基于不是在单一选通期间内获取的衰减校正数据的用于时间选通PET数据的选通特定的衰减校正因子。
例如，本发明的一些实施例提供了一种用于移动结构的多个时间选通PET图像的衰减校正的方法。方法包括基于结构的基线MR图像建立ACF矩阵；展开多个图像变换，每一个图像变换将基线MR图像与移动结构的相应电影MR图像配准；将每一个图像变换应用于ACF矩阵以生成对应的电影ACF矩阵；将每一个电影MR图像及其对应的电影ACF矩阵时间关联到多个时间选通PET数据集之一；及从每一个时间选通PET数据集产生衰减校正的时间选通PET图像。在某些实施例中，基线MR图像是复合图像。在一些其他实施例中，基线MR图像是平均图像。在一些实施例中，在结构基本上是静止的同时获得基线MR图像。在一些实施例中，使用组织区分扫描协议获得基线MR图像。在一些实施例中，每一个电影MR图像使用与每一个其他电影MR图像相同的扫描协议。在某些实施例中，展开多个图像变换包括像素强度缩放或像素强度移动。在某些实施例中，展开多个图像变换包括迭代计算从基线MR图像与从电影MR图像的中间变换构成的联合像素强度直方图的互信息。在某些实施例中，展开多个图像变换包括为每一个电影MR图像建立图像变换μ，其包括刚性变换参数和变形系数{δj}，用于将基线MR图像配准到那个电影MR图像。在某些实施例中，产生衰减校正的时间选通PET图像包括将与那个PET数据集时间关联的至少一个电影ACF矩阵应用于多个时间选通PET数据集之一。在某些实施例中产生衰减校正的时间选通PET图像包括将与那个PET数据集时间关联的多个电影ACF矩阵的简单平均、时间加权平均、速度加权平均或时间-速度加权平均中的一个应用于多个时间选通PET数据集之一。
本发明的其他方面提供了一种用于生成多个衰减校正的时间选通PET图像的装置。该装置包括PET-MR扫描器；与PET-MR扫描器通信连接的控制器，及与控制器通信连接的图像处理器。控制器配置成获取结构的基线MR图像，获取结构的电影MR图像序列，及获取结构的多个时间选通PET数据集。图像处理器配置成基于基线MR图像生成基线ACF矩阵，将基线MR图像与电影MR图像的至少一些配准，基于基线MR图像与电影MR图像的配准生成电影ACF矩阵，及基于电影ACF矩阵衰减校正时间选通PET数据集。在某些方案中，控制器配置成将每一个电影ACF矩阵时间关联到一个时间选通PET数据集。在某些方面中，图像处理器配置成将每一个电影ACF矩阵时间关联到一个时间选通PET数据集。在某些方案中，图像处理器配置成使用局部扰动刚性配准将基线MR图像配准到每一个电影MR图像。在某些方面中，图像处理器配置成使用互信息算法将基线MR图像配准到每一个电影MR图像。在某些方面中，控制器配置成在结构基本上是静止的同时获得基线MR图像。在某些方面中，控制器配置成在结构移动的同时并行获得电影MR图像和PET数据集。在某些方面中，图像处理器配置成使用与那个PET数据集时间关联的多个电影ACF矩阵的简单平均、时间加权平均、速度加权平均或时间-速度加权平均中的一个衰减校正多个时间选通PET数据集至少之一。
本发明的其他方面提供了一种用于PET成像的衰减校正的方法。该方法包括：在第一MR获取中，获得适合于导出PET衰减校正因子的基线MR图像；及在第二MR获取中，与呼吸选通PET研究并行地获得呼吸周期期间的MR图像的电影序列。该方法还包括将基线MR图像与每一个电影MR图像非刚性配准，以便产生一组配准的MR图像，其对应于来自呼吸选通PET研究的PET数据集，且可以用于导出衰减校正因子。配准的MR图像随后用于为它们相应的时间选通PET数据集形成衰减校正因子。
要理解，上面描述意图是说明性的而不是限制性的。例如，上述实施例(和/或其方面)可彼此结合使用。另外，可进行多种修改以使特定情形或材料适合本发明主题的教导，而没有背离其范围。虽然本文所述材料的尺寸和类型意图定义本发明主题的参数，但是它们决不是限制性的，而是示范性的实施例。在审查上面描述时，许多其他实施例对于本领域的技术人员将是显而易见的。因此，应参考所附权利要求书连同这类权利要求书所被赋予的等同物的全部范围来确定本发明主题的范围。在所附权利要求书中，术语“包含”和“其中”用作相应术语“包括”和“其中”的易懂英语等同物。此外，在下面权利要求书中，术语“第一”、“第二”和“第三”等只用作标记，而不是意图对其对象强加数字要求。此外，下面的权利要求书的限制没有以方法加功能形式来书写并且不意图基于35U.S.C.§112(f)第六段，除非并且直到这类权利要求限制确切地使用后面是缺乏进一步结构的功能陈述的短语“用于…的部件”。
本书面描述使用示例来公开本发明主题的若干实施例，并且还使本领域的任何技术人员能够实施本发明主题的实施例，包含制作和使用任何装置或系统以及执行任何包含的方法。本发明主题的可取得专利的范围由权利要求书限定，并且可包含本领域的技术人员想到的其他示例。如果这类其他示例具有没有不同于权利要求书的文字语言的结构元件，或者如果它们包含具有与权利要求书的文字语言的无实质差异的等效结构元件，则它们意图处于权利要求书的范围之内。
如本文所使用的，以单数阐述并且以单词“一”或“一个”进行的元件或步骤应当被理解为并不排除多个所述元件或步骤，除非明确规定这种排除。此外，对本发明主题的“一个实施例”的提及并不意图被解释为排除也结合了所阐述特征的附加实施例的存在。此外，除非相反的明确规定，否则，“包括”、“包含”或“具有”具有特定性质的元件或者多个元件的实施例可包含没有那种性质的附加的这类元件。
由于在不脱离本发明涉及的精神和范围的情况下上述实施例中可以做出某些变化，意图上面描述或附图中所示的全部主题都应仅解释为图示本发明的概念的示例，而不应理解为限制本发明。

Claims (20)

1.一种用于移动结构的多个时间选通的PET图像的衰减校正的方法，包括：
基于所述结构的基线MR图像建立ACF矩阵；
展开多个图像变换，所述图像变换的每一个将所述基线MR图像与所述移动结构的相应电影MR图像配准；
将所述图像变换的每一个应用于所述ACF矩阵以生成对应的电影ACF矩阵；
将所述电影MR图像的每一个及其对应的电影ACF矩阵时间关联到所述多个时间选通PET数据集之一；以及
从所述时间选通PET数据集的每一个产生衰减校正的时间选通PET图像。
2.根据权利要求1所述的方法，其中，所述基线MR图像是复合图像。
3.根据权利要求1所述的方法，其中，所述基线MR图像是平均图像。
4.根据权利要求1所述的方法，进一步包括在所述结构基本上是静止的同时获得所述基线MR图像的步骤。
5.根据权利要求1所述的方法，进一步包括使用组织区分扫描协议获得所述基线MR图像的步骤。
6.根据权利要求1所述的方法，其中，所述电影MR图像的每一个使用与所述电影MR图像的每一个其他相同的扫描协议。
7.根据权利要求1所述的方法，其中，展开所述多个图像变换包括像素强度缩放或像素强度移动。
8.根据权利要求1所述的方法，其中，展开所述多个图像变换包括，对于所述电影MR图像的每一个，迭代计算从所述基线MR图像与从所述电影MR图像的中间变换构成的联合像素强度直方图的互信息。
9.根据权利要求1所述的方法，其中，展开所述多个图像变换包括为每一个电影MR图像建立图像变换μ，其包括刚性变换参数{γ，θ，φ，[tx，ty，tz]}和变形系数{δj}，用于将所述基线MR图像配准到那个电影MR图像。
10.根据权利要求1所述的方法，其中，产生衰减校正的时间选通PET图像包括将与那个PET数据集时间关联的至少一个电影ACF矩阵应用于所述多个时间选通PET数据集的一个。
11.根据权利要求10所述的方法，其中，产生衰减校正的时间选通PET图像包括将与那个PET数据集时间关联的多个电影ACF矩阵的简单平均、时间加权平均、速度-加权平均或时间-速度加权平均中的一个应用于所述多个时间选通PET数据集的一个。
12.一种用于生成多个衰减校正的时间选通PET图像的装置，所述装置包括：
PET-MR扫描器；
控制器，所述控制器与所述PET-MR扫描器通信连接，其中所述控制器配置成获取结构的基线MR图像，获取所述结构的电影MR图像序列，以及获取所述结构的多个时间选通PET数据集；以及
图像处理器，所述图像处理器与所述控制器通信连接，其中，所述图像处理器配置成基于所述基线MR图像生成基线ACF矩阵，将所述基线MR图像与所述电影MR图像的至少一些配准，基于所述基线MR图像与所述电影MR图像的配准生成电影ACF矩阵，以及基于所述电影ACF矩阵衰减校正所述时间选通PET数据集。
13.根据权利要求12所述的装置，其中，所述控制器配置成将所述电影ACF矩阵的每一个时间关联到所述时间选通PET数据集的一个。
14.根据权利要求12所述的装置，其中，所述图像处理器配置成将所述电影ACF矩阵的每一个时间关联到所述时间选通PET数据集的一个。
15.根据权利要求12所述的装置，其中，所述图像处理器配置成使用局部扰动刚性配准将所述基线MR图像配准到所述电影MR图像的每一个。
16.根据权利要求12所述的装置，其中，所述图像处理器配置成使用互信息算法将所述基线MR图像配准到所述电影MR图像的每一个。
17.根据权利要求12所述的装置，其中，所述控制器配置成在所述结构基本上是静止的同时获得所述基线MR图像。
18.根据权利要求12所述的装置，其中，所述控制器配置成在所述结构移动的同时并行获得所述电影MR图像和所述PET数据集。
19.根据权利要求12所述的装置，其中，所述图像处理器配置成使用与那个PET数据集时间关联的多个电影ACF矩阵的简单平均、时间加权平均、速度-加权平均或时间-速度加权平均中的一个衰减校正所述时间选通PET数据集的至少一个。
20.一种用于呼吸选通PET研究中的PET成像的衰减校正的方法，包括：
在第一MR获取中，获得适合于导出PET衰减校正因子的基线MR图像；
在第二MR获取中，与所述呼吸选通PET研究并行地获得呼吸周期期间的MR图像的电影序列；以及
将所述基线MR图像与所述电影MR图像的每一个非刚性配准，以便产生一组配准的MR图像，其对应于来自所述呼吸选通PET研究的PET数据集，并且可以用于导出衰减校正因子，
其中，配准的MR图像随后用于为它们相应的时间选通PET数据集形成衰减校正因子。
CN201480070827.8A 2013-12-23 2014-08-06 用于时间选通pet研究的基于选通特定的mr的衰减校正的方法和装置 Pending CN106062818A (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US14/138,311 US9168015B2 (en) 2013-12-23 2013-12-23 Method and apparatus for gate specific MR-based attenuation correction of time-gated PET studies
US14/138311 2013-12-23
PCT/US2014/049885 WO2015099837A1 (en) 2013-12-23 2014-08-06 Method and apparatus for gate specific mr-based attenuation correction of timed-gated pet studies

Publications (1)

Publication Number Publication Date
CN106062818A true CN106062818A (zh) 2016-10-26

Family

ID=51392410

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201480070827.8A Pending CN106062818A (zh) 2013-12-23 2014-08-06 用于时间选通pet研究的基于选通特定的mr的衰减校正的方法和装置

Country Status (4)

Country Link
US (2) US9168015B2 (zh)
CN (1) CN106062818A (zh)
DE (1) DE112014005888T5 (zh)
WO (1) WO2015099837A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106691487A (zh) * 2017-01-05 2017-05-24 沈阳东软医疗系统有限公司 成像方法和成像系统
CN109961419A (zh) * 2019-03-26 2019-07-02 江苏赛诺格兰医疗科技有限公司 对pet活度分布图像进行衰减校正的校正信息获取方法

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9168015B2 (en) * 2013-12-23 2015-10-27 General Electric Corporation Method and apparatus for gate specific MR-based attenuation correction of time-gated PET studies
DE102019216648B3 (de) * 2019-10-29 2021-04-22 Siemens Healthcare Gmbh Verfahren zur Aufnahme eines PET-Bilddatensatzes mit einer Magnetresonanz-PET-Einrichtung, Magnetresonanz-PET-Einrichtung, Computerprogramm und elektronisch lesbarer Datenträger

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130320973A1 (en) * 2012-05-30 2013-12-05 Siemens Medical Solutions Usa, Inc. Method of attenuation correction of positron emission tomography data and combined positron emission tomography and magnetic resonance tomography system

Family Cites Families (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5251128A (en) * 1990-11-19 1993-10-05 General Electric Company Motion artifact reduction in projection imaging
US6377833B1 (en) * 1999-01-25 2002-04-23 Douglas Albert System and method for computer input of dynamic mental information
WO2004042546A1 (en) * 2002-11-04 2004-05-21 V-Target Technologies Ltd. Apparatus and methods for imaging and attenuation correction
US7180074B1 (en) * 2001-06-27 2007-02-20 Crosetto Dario B Method and apparatus for whole-body, three-dimensional, dynamic PET/CT examination
US6990368B2 (en) * 2002-04-04 2006-01-24 Surgical Navigation Technologies, Inc. Method and apparatus for virtual digital subtraction angiography
US7117026B2 (en) * 2002-06-12 2006-10-03 Koninklijke Philips Electronics N.V. Physiological model based non-rigid image registration
CN101080746B (zh) * 2004-12-15 2012-07-04 皇家飞利浦电子股份有限公司 多模态图像的配准
US7787675B2 (en) * 2005-03-17 2010-08-31 General Electric Company System, program product, and methods for attenuation correction of emission data on PET/CT and SPEC/CT
WO2008151202A2 (en) * 2007-06-03 2008-12-11 The Regents Of The University Of California Elastic deformable heart and torso phantom for nuclear imaging of heart with realistic modes of cardiac and respiratory motion
US8098916B2 (en) * 2007-10-30 2012-01-17 General Electric Company System and method for image-based attenuation correction of PET/SPECT images
JP2011503570A (ja) * 2007-11-09 2011-01-27 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 減衰マップを形成するための装置及び方法
US8208709B2 (en) * 2008-04-17 2012-06-26 The Ohio State University Research Foundation System and method for improved real-time cine imaging
US8218848B2 (en) * 2008-07-23 2012-07-10 Siemens Aktiengesellschaft System and method for the generation of attenuation correction maps from MR images
US8724875B2 (en) * 2008-08-15 2014-05-13 Koninklijke Philips N.V. Attenuation correction for PET or SPECT nuclear imaging systems using magnetic resonance spectroscopic image data
US9453922B2 (en) * 2008-09-09 2016-09-27 Multi-Magnetics Incorporated System and method for correcting attenuation in hybrid medical imaging
DE102009030714A1 (de) * 2009-06-26 2010-12-30 Siemens Aktiengesellschaft Verfahren zur Absorptionskorrektur von PET-Daten und MR-PET-Anlage
US8600132B2 (en) * 2011-05-03 2013-12-03 General Electric Company Method and apparatus for motion correcting medical images
WO2013177661A1 (en) * 2012-05-29 2013-12-05 University Of Manitoba Systems and methods for improving the quality of images in a pet scan
US8942445B2 (en) * 2012-09-14 2015-01-27 General Electric Company Method and system for correction of lung density variation in positron emission tomography using magnetic resonance imaging
US9168015B2 (en) * 2013-12-23 2015-10-27 General Electric Corporation Method and apparatus for gate specific MR-based attenuation correction of time-gated PET studies

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130320973A1 (en) * 2012-05-30 2013-12-05 Siemens Medical Solutions Usa, Inc. Method of attenuation correction of positron emission tomography data and combined positron emission tomography and magnetic resonance tomography system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CHRISTIAN BUERGER等: "Investigation of MR-Based Attenuation Correction and Motion Compensation for Hybrid PET/MR", 《IEEE TRANSACTIONS ON NUCLEAR SCIENCE》 *
DAVID MATTES 等: "PET-CT Image Registration in the Chest Using Free-form Deformations", 《IEEE TRANSACTIONS ON MEDICAL IMAGING》 *
FREDERIK MAES 等: "Medical Image Registration Using Mutual Information", 《PROCEEDINGS OF THE IEEE》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106691487A (zh) * 2017-01-05 2017-05-24 沈阳东软医疗系统有限公司 成像方法和成像系统
CN109961419A (zh) * 2019-03-26 2019-07-02 江苏赛诺格兰医疗科技有限公司 对pet活度分布图像进行衰减校正的校正信息获取方法

Also Published As

Publication number Publication date
DE112014005888T5 (de) 2016-09-15
US20150178906A1 (en) 2015-06-25
US20160029990A1 (en) 2016-02-04
US9168015B2 (en) 2015-10-27
US9427206B2 (en) 2016-08-30
WO2015099837A1 (en) 2015-07-02

Similar Documents

Publication Publication Date Title
US10810768B2 (en) System and method for segmenting medical image
CN103620437B (zh) 使用具有fid和多梯度回波采集的超短回波时间脉冲序列以及水‑脂肪分离处理的骨mri
US11020077B2 (en) Simultaneous CT-MRI image reconstruction
Li et al. Model‐based image reconstruction for four‐dimensional PET
CN105263407B (zh) 使用磁共振成像的对骨组织的探测
CN109741284A (zh) 用于校正pet成像中由呼吸运动引起的失配的系统和方法
CN109035355A (zh) 用于pet图像重建的系统和方法
Paganelli et al. Liver 4DMRI: a retrospective image‐based sorting method
CN109009200A (zh) 用于正电子发射断层成像图像重建的系统和方法
CN107133549B (zh) Ect运动门控信号获取方法及ect图像重建方法
Lamare et al. Evaluation of respiratory and cardiac motion correction schemes in dual gated PET/CT cardiac imaging
Rausch et al. Technical and instrumentational foundations of PET/MRI
Liu et al. Four dimensional magnetic resonance imaging with retrospective k‐space reordering: a feasibility study
WO2018006419A1 (en) System and method for generating attenuation map
CN109961834A (zh) 影像诊断报告的生成方法及设备
CN109009199A (zh) 用于正电子发射断层扫描中的图像数据处理的系统和方法
Petibon et al. Impact of motion and partial volume effects correction on PET myocardial perfusion imaging using simultaneous PET-MR
BR112018010953B1 (pt) Sistema de imageamento por ressonância magnética, produto de programa de computador e método para operar um sistema de imageamento por ressonância magnética
Fayad et al. The use of a generalized reconstruction by inversion of coupled systems (GRICS) approach for generic respiratory motion correction in PET/MR imaging
CN105579861B (zh) 在具有dixon脉冲序列的pet/mr成像中的基于mr的衰减校正
CN106062818A (zh) 用于时间选通pet研究的基于选通特定的mr的衰减校正的方法和装置
Feng et al. Dual respiratory and cardiac motion estimation in PET imaging: methods design and quantitative evaluation
EP2907107B1 (en) Multi-structure atlas and/or use thereof
Zhang et al. A multisource adaptive magnetic resonance image fusion technique for versatile contrast magnetic resonance imaging
Liu et al. Myocardial perfusion SPECT using a rotating multi‐segment slant‐hole collimator

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20161026