CN103415779B - 利用混合式核成像/mr中的核发射数据进行mr分割的系统及方法 - Google Patents

利用混合式核成像/mr中的核发射数据进行mr分割的系统及方法 Download PDF

Info

Publication number
CN103415779B
CN103415779B CN201280011934.4A CN201280011934A CN103415779B CN 103415779 B CN103415779 B CN 103415779B CN 201280011934 A CN201280011934 A CN 201280011934A CN 103415779 B CN103415779 B CN 103415779B
Authority
CN
China
Prior art keywords
image
bone
air
fuzzy
cutting
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
CN201280011934.4A
Other languages
English (en)
Other versions
CN103415779A (zh
Inventor
Z·胡
N·欧嘉
C-H·通
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
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 Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN103415779A publication Critical patent/CN103415779A/zh
Application granted granted Critical
Publication of CN103415779B publication Critical patent/CN103415779B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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
    • 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/1603—Measuring radiation intensity with a combination of at least two different types of detector
    • 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)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Molecular Biology (AREA)
  • Theoretical Computer Science (AREA)
  • Pulmonology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Nuclear Medicine (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

在生成磁共振（MR）衰减图（39）时，分割MR图像以识别患者的身体轮廓、软组织结构和包括骨和/或空气的模糊结构。为了在模糊结构中区分骨和空气，分割同一患者或感兴趣区域的核发射图像（例如PET）。将经分割的功能图像数据与经分割的MR图像数据相关，以在模糊结构中区分骨和空气。分别向经分割的MR图像中的空气体素和骨体素分配适当的辐射衰减值，从增强的经分割的MR图像生成MR衰减图，其中已经对空气和骨之间的模糊进行了分辨。使用MR衰减图以生成衰减校正的核图像，向用户显示该衰减校正的核图像。

Description

利用混合式核成像/MR中的核发射数据进行MR分割的系统及方法
技术领域
本申请具体应用于组合式PET/MR医疗成像系统。不过，要认识到，所描述的技术也可以应用于其他诊断系统、其他成像情形或其他诊断技术中。
背景技术
在核成像中，穿过受试者的某些伽玛辐射与患者的身体相互作用并被患者的身体吸收。骨组织和致密组织吸收的伽玛辐射比软组织多，而软组织吸收的伽玛辐射比空气多(例如，在肺中)。如果未经校正，这些不同量的辐射衰减可能会引起所得图像中观察到的相对放射性浓度错误地增加或减少。校正这种衰减通常基于衰减图，即，受试者的图或图像，其中每个体素值表示患者体内相应组织的相对或绝对辐射衰减属性。
有时将当前市场上诸如正电子发射断层摄影(PET)扫描机或单光子发射计算机断层摄影(SPECT)扫描机的核扫描机与计算断层摄影(CT)扫描机相组合。在CT图像中，体素值表示患者体内相应组织的辐射衰减属性。因此，CT扫描机有时用于为核扫描机生成衰减图。
近来，已提出了组合式核成像/磁共振(MR)成像。MR扫描机对共振偶极子，尤其是水和脂类的1H偶极子进行成像。MR图像的体素值和患者体内相应组织的辐射衰减属性之间没有直接的关联。通过识别与MR体素相关联的组织以及所识别组织的辐射衰减属性，MR图像可以用来提供衰减图。但是，在标定中会有模糊，这可能会导致次优的衰减图。例如，常规MR图像将骨密质显示为暗体素，但是也将空气显示为暗体素，这就使得仅基于图像强度来区分骨密质和空气出现困难。骨密质具有高辐射衰减属性而空气具有很低的辐射衰减属性。不正确地分辨这种空气/骨模糊可能会导致衰减图和最终的PET或SPECT图像中出现潜在的显著误差。
本申请涉及新的和改进的系统和方法，其便于区分MR衰减图中呈现的骨密质与其他组织或(一个或多个)对象，这克服了上述问题和其他问题。
发明内容
根据一个方面，一种便于对磁共振(MR)图像或衰减图中模糊进行分辨的系统包括：MR重建处理器，其从在MR扫描机扫描受试者期间采集的原始MR数据生成MR图像；以及功能图像重建处理器，其从在功能图像扫描机扫描受试者期间采集的功能图像数据生成功能图像。该系统还包括处理器，所述处理器被编程为：分割所述MR图像以生成经分割的MR图像，所述经分割的MR图像具有空气/骨模糊区域，其中，空气体素和骨体素之间的模糊未分辨；分割所述功能图像以生成具有骨区域和其他区域的经分割的功能图像，以及比较所述经分割的功能图像中的骨区域与所述经分割的MR图像中的骨/空气模糊区域，以对所述经分割的MR图像中对应于骨和空气的骨/空气区域中体素之间的模糊进行分辨。所述处理器还被编程为：在所述经分割的MR图像中向识别出的骨体素分配与骨相符的辐射衰减值，向识别出的空气体素分配与空气相符的辐射衰减值；以及利用分配的辐射衰减值生成MR衰减图。
根据另一方面，一种对MR图像或衰减图中的模糊进行分辨的方法包括：从在MR扫描机扫描受试者期间采集的原始MR数据生成MR图像；从在功能图像扫描机扫描受试者期间采集的功能图像数据生成功能图像；以及分割所述MR图像以生成经分割的MR图像，所述经分割的MR图像具有空气/骨模糊区域，其中，空气体素和骨体素之间的模糊未分辨。该方法还包括：分割所述功能图像以生成具有骨区域的经分割的功能图像；将所述经分割的MR图像中的骨/空气模糊区域与所述分割的功能图像中的骨区域中的至少一个进行比较，以对所述经分割的MR图像中与骨和空气相对应的骨/空气区域中体素之间的模糊进行分辨；以及在所述经分割的MR图像中向识别出的骨体素分配与骨相符的辐射衰减值，向识别出的空气体素分配与空气相符的辐射衰减值。此外，该方法包括：利用分配的辐射衰减值生成MR衰减图，利用所述MR衰减图将所述功能图像数据重建成图像以校正所述功能图像数据中的衰减，以及在显示器上向用户显示所述图像。
根据另一方面，一种利用核发射数据增强MR衰减图的方法包括：分割MR图像以识别MR图像中的软组织结构和模糊结构，所述模糊结构包括骨和空气中的一种或多种；分割核图像以识别所述核图像中的骨结构；以及在所述经分割的MR图像上叠加分割的核图像。该方法还包括：通过将所述模糊结构与所述重叠分割核图像中的骨结构进行比较，将所述经分割的MR图像中的模糊结构识别为骨或空气；在所述经分割的MR图像中向识别出的骨体素分配第一辐射衰减值，向经分割的MR图像中识别出的空气体素分配第二辐射衰减值；以及从所述经分割的MR图像生成MR衰减图。
一个优点是衰减图的可靠性得以改进。
另一个优点在于分辨了受试者中的骨组织和气穴之间的模糊。
本领域的普通技术人员在阅读和理解以下详细描述之后将理解本创新的更多优点。
附图说明
本创新可以采取各种部件和部件布置以及各种步骤和步骤布置的形式。附图仅仅为了例示各方面，不应被解释为限制本发明。
图1示出了利用单示踪剂或多示踪剂功能(例如，PET，包括它的变型，例如飞行时间PET(TOF-PET)等)成像以便于生成改进的衰减图的系统，其中，从解剖学成像模态，尤其是MR，导出衰减图。
图2示出了利用单示踪剂或多示踪剂功能成像来生成改进的衰减图的方法，其中，从诸如MR的解剖学成像模态中导出衰减图。
图3示出了18F氟化钠PET全身解剖图像的对比说明。
具体实施方式
本创新通过以下方式克服上述问题：使用用于分割骨结构的核发射数据并将分割的核发射数据与MR图像分割相组合以导出精确的衰减图。首先生成解剖结构的基于MR的图像，从中可以识别出躯体、肺和各种软组织的轮廓。尽管可以使用生理学知识来区分肺和骨骼，但是区分气穴和骨结构，尤其是骨结构中的气穴，仍具有模糊。在本申请中，重建和分割核发射图像来识别骨结构，将其与经分割的MR图像中模糊的气-骨结构相比较以消除与区分空气和骨相关的模糊。可以结合单示踪剂或多示踪剂核成像协议来实施本创新，其中，从诸如MR的解剖学成像模态中导出衰减图。
图1示出了利用单示踪剂或多示踪剂核(例如，PET，包括它的变型，例如飞行时间PET(TOF-PET)等)成像以便于生成改进的衰减图的系统10，其中从解剖学成像模态，尤其是MR中导出衰减图。系统10包括MR扫描机11和功能图像扫描机或核扫描机12。在一个实施例中，MR扫描机和核扫描机是独立的扫描机。在另一个实施例中，将MR扫描机11和核扫描机12容纳到单个单元中，例如多模态MR/核扫描机13中。使用MR扫描机扫描受试者以生成存储在MR扫描数据存储器或缓存14中的原始MR扫描数据。在一个实施例中，在核扫描之前但是在给患者注射放射性示踪团剂(例如18F放射性示踪剂)之后采集MR扫描数据。可以在核扫描之前或之后执行MR扫描。也使用核扫描机12扫描受试者，以采集存储在核原始数据存储器或缓存15中的原始核扫描数据。由MR重建处理器16来重建原始的MR数据，由功能图像重建处理器17来重建原始的核数据。任选地，功能图像重建处理器17在所采集的核扫描数据上执行散射校正以改进从中重建的功能图像。该系统还包括图像数据库18，图像数据库18存储了由MR重建处理器生成的MR图像20以及由核重建处理器生成的功能图像21。
系统10包括工作站22，工作站22包括输入设备24(例如键盘、鼠标、指示笔、触摸屏、方向垫、麦克风或任何其他适当的输入设备)，经由输入设备，用户向系统中输入信息。工作站还包括存储器26和处理器28，存储器26存储并且处理器28执行，计算机可执行指令(例如例程、程序、算法、软件代码等)，用于执行本文中描述的各种功能、方法、流程等。该系统还包括显示器30，在其上向用户呈现信息。此外，本文中所使用的“模块”表示一组计算机可执行指令、软件代码、程序、例程或用于执行所描述的功能的其他计算机可执行单元等，本领域技术人员将会理解。额外地或可选地，可以手动地执行本文中关于模块所述的一种或多种功能。
存储器可以是存储控制程序计算机可读介质，例如，磁盘、硬盘驱动器等。常见形式的非暂态计算机可读介质包括，例如，软盘、软碟、硬盘、磁带、或任何其他磁存储介质，CD-ROM、DVD、或任何其他光学介质，RAM、ROM、PROM、EPROM、FLASH-EPROM、它们的变型，其他存储器芯片或存储器盒，或任何其他处理器能够从中读取和执行的有形介质。在该背景下，本文中描述的系统可以实现在或实现为：一个或多个通用计算机；(一个或多个)专用计算机；编程的微处理器或微控制器以及外围集成电路元件；ASIC或其他集成电路；数字信号处理器；硬布线电子或逻辑电路，如分立元件电路；可编程逻辑器件，如PLD、PLA、FPGA、图形卡CPU(GPU)或PAL；等等。
存储器存储由处理器执行的MR分割模块32以分割患者的MR图像20，以便于识别不同的组织类型，例如软组织、骨密质、气穴等，并生成一个或多个经分割的MR图像33。分割模块32还可以包括一个或多个附加处理器以提高处理速度。在一个实施例中，这种分割是自动的。在另一个实施例中，使用输入设备24和显示器30手动地或半自动地执行该分割。不过，当试图区分MR图像中的骨密质和空气时，经分割的MR图像可能会留下未分辨的模糊。因此，存储器存储有助于分割相同受试者的功能核图像21的功能图像分割模块34。可以使用输入设备24和显示器30全自动地、半自动地或手动地执行核图像分割程序。
功能图像分割模块有助于生成一个或多个分割的核发射图像35，其中，可以将骨，尤其是骨密质与包括空气在内的组织的其余部分，区分开来。模糊分辨模块36将核图像中分割的骨区域与经分割的MR图像中的空气/骨区域相比较，以区分MR图像中的骨和空气，并将适当的辐射衰减值μ分配给对应的体素。可以用各种方式对分割的核图像中的骨区域和经分割的MR图像中的骨/空气区域进行比较。为了简化比较，可以使用肺的基本衰减属性来识别属于空气类别的肺体积。比较骨与骨/空气模糊区域的一种方法是将骨区域和骨/空气模糊区域叠加，并查看两个区域中相对应的每对体素。如果在核骨区域和MR空气/骨模糊区域中都有体素，它就被指定为骨。在另一种方法中，将核骨区域和MR空气/骨模糊区域相减。例如，可以通过将空气/骨模糊区域用作掩模来分割核图像，以消除在MR图像中不模糊的核图像的任何体素，即，分割对应于MR图像空气/骨模糊区域的核图像的第一区域。在空气区域和骨区域之间很容易分割第一区域。从MR空气/骨模糊区域减核图像骨区域后就只剩下了对应于MR图像中的空气的体素。例如，可以为骨密质分配1.8的μ值(假定水的μ值为1)，可以为空气分配0.0的μ值。通过这种方式，就提升了在经分割的MR图像中识别出的组织μ值的精确度。衰减图生成模块38使用精确度提高的μ值信息来生成一个或多个核发射数据增强的MR衰减图39。要认识到，可以自动地、半自动地或手动地执行MR图像和/或PET图像的分割。
一旦生成了MR衰减图，就重建所采集的原始功能图像数据以生成校正的衰减功能图像，其可以在显示器30上显示出来。在另一实施例中，生成校正的组合式MR和衰减功能图像，并将其在显示器上呈现给用户。
在分割MR和/或功能核图像时，可以采用多种技术中的一种或多种。例如，存储器26存储由处理器28执行的阈值化模块40，以识别出具有比预定阈值大(或小)的值或处于预定窗口内的值的体素。额外地或可选地，存储器存储由处理器执行的区域增长模块42以便于图像分割。另外，存储器存储基于图集的分割模块44，其有助于使用解剖图像分割的图集，通过使模型拟合当前图像来分割MR和/或核图像。在另一实施例中，存储器存储由处理器执行的基于模型的分割模块46，以便于分割MR和/或核图像。要认识到，可以自动地、半自动地或者手动地执行所述分割技术中的任何一种技术。
在一个实施例中，功能扫描机是PET扫描机。现有技术中已知，当电子和正电子相遇时，它们就会湮灭，根据动量守恒原理，发射出两条相反方向的511keV的伽玛射线。在PET数据采集中，假定两个基本同时的511keV伽玛射线探测事件源于相同的正负电子湮灭事件，因此该事件发生于将两个基本同时的511keV的伽玛射线探测事件连接起来的“响应线”(LOR)的某处。这种响应线有时也被称为投影，所收集的PET数据就被称为投影数据。在常规PET中，将基本同时的511keV伽玛射线探测事件定义为两个在选定的短时间窗口内，例如彼此在一纳秒之内发生的511keV的伽玛射线探测事件。由于探测器元件湮没位置相对于探测器元件的变化，基本同时的伽玛光子探测事件之间就会发生小的(例如，亚纳秒)时间差。
在另一实施例中，PET扫描机是飞行时间(TOF)PET扫描机。TOFPET成像利用这一小的时间差，沿响应线对正负电子湮灭事件进行进一步定位。通常，湮灭事件发生于沿投影线更接近于首先发生伽马射线探测事件的点。如果两个伽玛射线探测事件在探测器的时间分辨率范围内同时发生，那么湮灭事件就发生在投影的中点处。
在另一实施例中，功能扫描机是单光子发射计算机断层摄影(SPECT)扫描机。在SPECT成像中，给被成像的受试者施予一种或多种放射性药物或放射性同位素，使得如上所述地从其发出发射辐射。使安装在旋转扫描架上的探测头围绕受试者旋转以从多个方向探测辐射。探测头可以以360度旋转围绕正在成像的受试者转动，并且随着旋转在多个离散位置进行扫描。另外，探测头可以沿较小的弧旋转或围绕受试者进行多圈旋转。重建由探测头接收到的发射投影数据或所测正弦图以生成SPECT图像。所采用的重建发射数据的重建技术可以包括但不限于，迭代重建、基于傅里叶变换的重建、滤波后向投影或其他一些合适的重建技术。
在另一实施例中，以列表模式采集核扫描数据。列表模式投影数据通常包括一列探测到的事件，列表中每个条目都包括诸如探测到事件的时间以及对应LOR的位置和方向这样的信息。具有TOF能力的扫描机而言，也提供了沿LOR对的湮没位置的估计。另外，可以对所采集的数据分类或将其分箱(bin)到正弦图分箱或投影分箱中。
MR扫描机能够生成软组织的详细图像。在MR成像中，利用组织内发现的各种化合物的特定属性来生成图像，例如最常将水用于该目的。在受到强的外部磁场的作用时，质子与外磁场对齐，导致净磁矩。由射频RF脉冲激发后，这种磁化生成可被探测到的射频(RF)信号。RF信号的特征在于其与磁场强度相关的频率。因此，利用磁场梯度来编码用于从探测到的信号中重建图像的空间信息。
图2图示了利用单示踪剂或多示踪剂核成像来生成增强衰减图的方法，其中从解剖学成像模态，尤其是MR中导出衰减图。在80，将MR扫描数据采集和重建为MR图像，用于：(a)分割，以及(b)与PET图像组合以生成同时表示患者解剖和功能信息的组合图像。在一个实施例中，在给患者注射放射性示踪剂团剂(例如18F放射性示踪剂或放射性药物)之后采集MR扫描数据。在82，分割MR图像以识别身体轮廓、肺和软组织等，空气/骨模糊区域可以是骨或空气，例如骨密质结构中的气穴。在84，利用注射的单示踪剂或多示踪剂采集相同受试者的发射数据，并将其重建以生成一阶发射图像，即，没有衰减校正的发射图像。在86，分割一阶发射图像以分割骨，尤其是基于18F示踪剂的摄取的骨密质结构。在88，通过比较在发射图像分割中识别出的骨区域的体素和MR图像的空气/骨模糊区域中相对应的体素，来对MR图像中空气和骨体素之间的模糊进行分辨。在90，为经分割的MR图像中识别出的解剖结构分配和/或校正辐射衰减值μ，以导出精确的衰减图。例如，诸如脂肪、肌肉、内脏等的软组织，可以具有分配的大约1.0的μ值，并且该值从经分割的MR图像确定。骨密质具有大约1.8的μ值，而空气具有大约0.0的μ值。不过，由于使用经分割的MR图像不能准确地区分空气和骨密质，因而采用所分割的发射图像来对它们之间的模糊进行分辨，然后在MR图像中给发射图像中识别出的骨密质组织分配1.8的μ值，而给MR图像中的空气体素分配0.0的μ值。在92，正如从分割的发射图像所确定地，使用校正的μ值生成改进的MR衰减图。在94，使用衰减图重建原始的PET数据，并将其显示给用户。在96，任选地生成PET/MR叠加图像并将其显示给用户，叠加图像包括使用衰减图重建的PET图像和已经分辨了骨/空气模糊的MR图像。
在一个实施例中，在82的MR图像分割和/或在86的PET图像分割可以使用多种技术中的任何一种，包括但并不限于，阈值化、区域增长、基于图集的方法、模型适配或者以上一种或多种方法与经验推理的组合等。
在另一实施例中，为了在86导出更好的骨密质结构识别，使用从发射数据重建的发射图像，其中，发射数据是用注射的增强骨结构的诸如18F氟化钠的放射性示踪剂采集的。在另一实施例中，采用多示踪剂注射。在另一实施例中，使用吸入示踪剂(例如，氙气等)。
在另一实施例中，迭代地执行该方法。例如，可以采集和采用多个核图像和/或多个MR图像来改进MR衰减图，其中，每一个后继的衰减图都相对于上一个得到改进，并采用其在PET图像的一个或多个图像中校正衰减。
图3示出了18F氟化钠PET解剖图像的对比说明。左边示出了原始的PET图像100。根据本文中所述的一个或多个方面，右边是显示了使用阈值分割法已识别出的分割骨结构的改进图像110。可以采用分割PET图像110来准确区分经分割的MR图像中的骨和空气。在一个实施例中，将分割PET图像110叠加到经分割的MR图像上，使分割PET图像中识别的骨体素与经分割的MR图像中模糊体素对准。通过这种方式，将MR图像中对应的模糊体素分辨为骨体素并给其分配适当的μ值。将剩余的模糊MR体素判定为空气体素。
要认识到，尽管结合使用核发射数据分辨MR图像分割中骨密质和空气之间模糊以增强MR衰减图精确度描述了以上系统和方法，但它们不限于此。而是，可以采用描述的系统和方法，利用设计成由模糊组织类型吸收的适当放射示踪剂来分辨任何组织类型之间的模糊。
已经参考几个实施例描述了该创新。他人在阅读和理解以上详细描述之后可能想到修改和变更。应当将本创新解释为包括所有这样的修改和变更，只要它们在权利要求或其等价方案的范围之内。

Claims (15)

1.一种便于对磁共振(MR)图像或衰减图中的模糊进行分辨的系统(10)，包括：
MR重建处理器(16)，其从在MR扫描机(11)扫描受试者期间采集的原始MR数据(14)生成MR图像(20)；
功能图像重建处理器(17)，其从在功能图像扫描机(12)扫描受试者期间采集的功能图像数据(15)生成功能图像(21)；以及
处理器(28)，其被编程为：
分割所述MR图像(20)以生成经分割的MR图像(33)，所述经分割的MR图像具有空气/骨模糊区域，在所述空气/骨模糊区域中空气体素和骨体素之间的模糊未分辨；
分割所述功能图像(21)以生成具有骨区域的经分割的功能图像(35)；
将所述经分割的功能图像(35)中的骨区域与所述经分割的MR图像(33)中的所述空气/骨模糊区域进行比较，以对所述经分割的MR图像(33)中与骨和空气相对应的所述空气/骨模糊区域中体素之间的模糊进行分辨；
在所述经分割的MR图像(33)中，向识别出的骨体素分配与骨相符的辐射衰减值，并且向识别出的空气体素分配与空气相符的辐射衰减值；并且
利用所分配的辐射衰减值生成MR衰减图(39)；
其中，所述功能图像数据包括来自PET扫描机的发射数据并且是利用单放射示踪剂采集的。
2.根据权利要求1所述的系统，其中，所述放射示踪剂包括18F同位素氟化物。
3.根据权利要求1-2中任一项所述的系统，其中，所述处理器被编程为：
将所述经分割的功能图像(35)的至少一个骨区域与所述经分割的MR图像(33)中的至少一个解剖学上相对应的空气/骨模糊区域进行叠加；
识别所述骨区域和所述空气/骨模糊区域中的对应体素对；并且
将所述经分割的MR图像中在所述经分割的功能图像的所述骨区域中具有对应体素的所述空气/骨模糊区域中的体素识别为所述经分割的MR图像中的骨体素。
4.根据权利要求1-2中任一项所述的系统，其中，所述处理器被编程为：
将所述经分割的功能图像(35)的至少一个空气区域与所述经分割的MR图像(33)中的至少一个解剖学上相对应的空气/骨模糊区域进行叠加；
识别所述空气区域和所述空气/骨模糊区域中的对应体素对；
从所述空气/骨模糊区域中的所述对应体素减去所述空气区域中的所述对应体素；并且
将所述空气/骨模糊区域中的剩余体素识别为骨体素。
5.根据权利要求1-2中任一项所述的系统，其中，所述功能图像扫描机(12)是以下中的至少一种：
正电子发射断层摄影(PET)扫描机，其采集所述功能图像数据；以及
单光子发射计算机断层摄影(SPECT)扫描机，其采集所述功能图像数据。
6.根据权利要求5所述的系统，其中，所述PET扫描机和所述MR扫描机(11)包括在单个多模PET/MR扫描设备(13)中。
7.根据权利要求1-2中任一项所述的系统，其中，所述处理器还被编程为：
执行阈值化模块(40)以分割所述MR图像(20)和所述功能图像(21)中的至少一个。
8.根据权利要求1-2中任一项所述的系统，其中，所述处理器还被编程为：
执行区域增长模块(42)以分割所述MR图像(20)和所述功能图像(21)中的至少一个。
9.根据权利要求1-2中任一项所述的系统，其中，所述处理器还被编程为：
执行基于图集的分割模块(44)以分割所述MR图像(20)和所述功能图像(21)中的至少一个。
10.根据权利要求1-2中任一项所述的系统，其中，所述处理器还被编程为：
执行基于模型的适配模块(46)以分割所述MR图像(20)和所述功能图像(21)中的至少一个。
11.一种对磁共振(MR)图像或衰减图中的模糊进行分辨的方法，包括：
从在MR扫描机(11)扫描受试者期间采集的原始MR数据(14)生成MR图像(20)；
从在功能图像扫描机(12)扫描受试者期间采集的功能图像数据(15)生成功能图像(21)；
分割所述MR图像(20)以生成经分割的MR图像(33)，所述经分割的MR图像具有空气/骨模糊区域，在所述空气/骨模糊区域中，空气体素和骨体素之间的模糊未分辨；
分割所述功能图像(21)以生成具有骨区域和其他组织区域的经分割的功能图像(35)；
将所述经分割的功能图像(35)中的所述骨区域与所述经分割的MR图像(33)中的所述空气/骨模糊区域进行比较，以对所述经分割的MR图像(33)中与骨和空气相对应的所述空气/骨模糊区域中体素之间的模糊进行分辨；
在所述经分割的MR图像(33)中，向识别出的骨体素分配与骨相符的辐射衰减值，向识别出的空气体素分配与空气相符的辐射衰减值；
利用所分配的辐射衰减值生成MR衰减图(39)；
利用所述MR衰减图将所述功能图像数据重建为图像，以校正所述功能图像数据中的衰减；并且
在显示器上向用户显示所述图像；
其中，所述功能图像数据包括来自PET扫描机的发射数据并且是利用单放射示踪剂采集的。
12.根据权利要求11所述的方法，其中，所述放射示踪剂包括18F同位素氟化物。
13.根据权利要求11-12中任一项所述的方法，还包括：
将所述经分割的功能图像(35)的至少一个骨区域与所述经分割的MR图像(33)中的至少一个解剖学上相对应的空气/骨模糊区域进行叠加；
识别所述骨区域和所述空气/骨模糊区域中的对应体素对；并且
将所述经分割的MR图像中在所述经分割的功能图像的所述骨区域中具有对应体素的所述空气/骨模糊区域中的体素识别为所述经分割的MR图像中的骨体素。
14.根据权利要求11-12中任一项所述的方法，还包括：
将所述经分割的功能图像(35)的至少一个空气区域与所述经分割的MR图像(33)中的至少一个解剖学上相对应的空气/骨模糊区域进行叠加；
识别所述空气区域和所述空气/骨模糊区域中的对应体素对；
从所述空气/骨模糊区域中的所述对应体素减去所述空气区域中的所述对应体素；并且
将所述空气/骨模糊区域中的剩余体素识别为骨体素。
15.根据权利要求11-12中任一项所述的方法，其中，所述功能图像扫描机(12)是以下中的至少一种：
正电子发射断层摄影(PET)扫描机，其采集所述功能图像数据(15)；以及
单光子发射计算机断层摄影(SPECT)扫描机，其采集所述功能图像数据。
CN201280011934.4A 2011-03-07 2012-03-02 利用混合式核成像/mr中的核发射数据进行mr分割的系统及方法 Active CN103415779B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201161449810P 2011-03-07 2011-03-07
US61/449,810 2011-03-07
PCT/IB2012/050998 WO2012120422A1 (en) 2011-03-07 2012-03-02 Mr segmentation using nuclear emission data in hybrid nuclear imaging/mr

Publications (2)

Publication Number Publication Date
CN103415779A CN103415779A (zh) 2013-11-27
CN103415779B true CN103415779B (zh) 2016-05-11

Family

ID=45937462

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201280011934.4A Active CN103415779B (zh) 2011-03-07 2012-03-02 利用混合式核成像/mr中的核发射数据进行mr分割的系统及方法

Country Status (5)

Country Link
US (1) US9271652B2 (zh)
EP (1) EP2684066B1 (zh)
CN (1) CN103415779B (zh)
RU (1) RU2595808C2 (zh)
WO (1) WO2012120422A1 (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9271652B2 (en) * 2011-03-07 2016-03-01 Koninklijke Philips N.V. MR segmentation using nuclear emission data in hybrid nuclear imaging/MR
EP2685899B1 (en) * 2011-03-17 2015-05-20 Koninklijke Philips N.V. Multiple modality cardiac imaging
EP2747658B1 (en) * 2011-11-16 2020-01-15 Koninklijke Philips N.V. Method to compute and present brain amyloid in gray matter
US9536303B2 (en) * 2012-03-29 2017-01-03 Koninklijke Philips N.V. MRI method for assigning individual pixels or voxels tissue-specific PET attenuation values
US20130303899A1 (en) * 2012-05-14 2013-11-14 University Of Maryland, Baltimore Quantitative evaluation of fractional regional ventilation using four-dimensional computed tomography
US8923592B2 (en) * 2012-05-29 2014-12-30 General Electric Company Methods and systems for performing attenuation correction
GB2510632B (en) * 2013-02-12 2015-03-11 Siemens Medical Solutions Residual activity correction at reconstruction
JP5966112B1 (ja) * 2013-04-02 2016-08-10 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 磁気共鳴画像法を使用した骨組織の検出
WO2015092622A2 (en) * 2013-12-20 2015-06-25 Koninklijke Philips N.V. Density guided attenuation map generation in pet/mr systems
US9943703B2 (en) 2014-07-28 2018-04-17 The University Of Maryland, Baltimore System and method for irradiation therapy using voxel based functional measurements of organs at risk
US9530206B2 (en) * 2015-03-31 2016-12-27 Sony Corporation Automatic 3D segmentation and cortical surfaces reconstruction from T1 MRI
CN108351396B (zh) * 2015-11-02 2021-09-14 皇家飞利浦有限公司 用于组织分类的方法、计算机程序产品和磁共振成像系统
US10210634B2 (en) 2016-07-20 2019-02-19 Shanghai United Imaging Healthcare Co., Ltd. System and method for segmenting medical image
CN107123095B (zh) * 2017-04-01 2020-03-31 上海联影医疗科技有限公司 一种pet图像重建方法、成像系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101077303A (zh) * 2006-04-04 2007-11-28 西门子公司 测定身体部位的正电子发射测量信息的方法及所属装置
CN101231745A (zh) * 2007-01-24 2008-07-30 中国科学院自动化研究所 一种对图像初始分割边界进行优化的自动分割方法
CN101390754A (zh) * 2007-09-20 2009-03-25 西门子公司 用于确定患者的正电子发射断层造影数据的衰减值的方法
CN101470180A (zh) * 2007-12-29 2009-07-01 西门子（中国）有限公司 磁共振成像中失真校准的方法和装置
CN101953693A (zh) * 2009-07-16 2011-01-26 西门子公司 完整或部分正电子发射断层造影衰减图的基于模型的估计
CN101976435A (zh) * 2010-10-07 2011-02-16 西安电子科技大学 基于对偶约束的联合学习超分辨方法

Family Cites Families (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6823205B1 (en) * 2001-02-08 2004-11-23 Boston University Radiology Associates Synthetic images for a magnetic resonance imaging scanner using linear combination of source images to generate contrast and spatial navigation
US6927406B2 (en) * 2002-10-22 2005-08-09 Iso-Science Laboratories, Inc. Multimodal imaging sources
US7831088B2 (en) * 2003-06-13 2010-11-09 Georgia Tech Research Corporation Data reconstruction using directional interpolation techniques
US7567696B2 (en) * 2004-03-02 2009-07-28 Siemens Medical Solutions Usa, Inc. System and method for detecting the aortic valve using a model-based segmentation technique
US20080021502A1 (en) * 2004-06-21 2008-01-24 The Trustees Of Columbia University In The City Of New York Systems and methods for automatic symmetry identification and for quantification of asymmetry for analytic, diagnostic and therapeutic purposes
PT103200B (pt) * 2004-09-30 2006-08-24 Taguspark-Soc. Prom.Desenv.Parq.Ci.Tec.Area Lisboa Sistema de tomografia por emissão de positrões (pet)
WO2007050437A2 (en) * 2005-10-21 2007-05-03 The General Hospital Corporation Methods and apparatus for segmentation and reconstruction for endovascular and endoluminal anatomical structures
DE602006020618D1 (de) * 2005-12-22 2011-04-21 Visen Medical Inc Kombiniertes röntgen- und optisches tomographie-bildgebungssystem
US7936922B2 (en) * 2006-11-22 2011-05-03 Adobe Systems Incorporated Method and apparatus for segmenting images
US7907766B2 (en) * 2007-01-02 2011-03-15 General Electric Company Automatic coronary artery calcium detection and labeling system
DE102008029897B4 (de) * 2008-06-24 2017-12-28 Siemens Healthcare Gmbh Verfahren zur Erfassung von MR-Daten eines Messobjektes bei einer MR-Untersuchung in einer Magnetresonanzanlage und entsprechend ausgestaltete Magnetresonanzanlage
US8218848B2 (en) 2008-07-23 2012-07-10 Siemens Aktiengesellschaft System and method for the generation of attenuation correction maps from MR images
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
US9453922B2 (en) * 2008-09-09 2016-09-27 Multi-Magnetics Incorporated System and method for correcting attenuation in hybrid medical imaging
US8406495B2 (en) * 2008-09-17 2013-03-26 Koninklijke Philips Electronics N.V. MR segmentation using transmission data in hybrid nuclear/MR imaging
EP2221627A1 (en) * 2009-02-20 2010-08-25 IBBT vzw Method and assembly for correcting a relaxation map for medical imaging applications
US9213074B2 (en) 2009-04-14 2015-12-15 The General Hospital Corporation Stem and method for acquiring MRI data from bone and soft tissues
US9271652B2 (en) * 2011-03-07 2016-03-01 Koninklijke Philips N.V. MR segmentation using nuclear emission data in hybrid nuclear imaging/MR

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101077303A (zh) * 2006-04-04 2007-11-28 西门子公司 测定身体部位的正电子发射测量信息的方法及所属装置
CN101231745A (zh) * 2007-01-24 2008-07-30 中国科学院自动化研究所 一种对图像初始分割边界进行优化的自动分割方法
CN101390754A (zh) * 2007-09-20 2009-03-25 西门子公司 用于确定患者的正电子发射断层造影数据的衰减值的方法
CN101470180A (zh) * 2007-12-29 2009-07-01 西门子（中国）有限公司 磁共振成像中失真校准的方法和装置
CN101953693A (zh) * 2009-07-16 2011-01-26 西门子公司 完整或部分正电子发射断层造影衰减图的基于模型的估计
CN101976435A (zh) * 2010-10-07 2011-02-16 西安电子科技大学 基于对偶约束的联合学习超分辨方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Advances in Attenuation Correction Techniques in PET;Habib zaidi et .al;《PET CLINICS》;20070401;第2卷(第2期);第199页左栏第3段-右栏第3段 *

Also Published As

Publication number Publication date
EP2684066A1 (en) 2014-01-15
US20130336564A1 (en) 2013-12-19
WO2012120422A1 (en) 2012-09-13
RU2013144764A (ru) 2015-04-20
CN103415779A (zh) 2013-11-27
EP2684066B1 (en) 2020-08-19
RU2595808C2 (ru) 2016-08-27
US9271652B2 (en) 2016-03-01

Similar Documents

Publication Publication Date Title
CN103415779B (zh) 利用混合式核成像/mr中的核发射数据进行mr分割的系统及方法
US9474495B2 (en) System and method for joint estimation of attenuation and activity information
EP2441050B1 (en) Time-of-flight positron emission tomography reconstruction using image content generated event-by-event based on time-of-flight information
CN103315760B (zh) 用于核医疗成像中基于发射数据进行衰减补偿的系统和方法
CN102763138B (zh) 辐射治疗中的运动校正
US8218848B2 (en) System and method for the generation of attenuation correction maps from MR images
JP5385394B2 (ja) ハイブリッド型の核／ｍｒ画像形成における透過データを使用したｍｒセグメンテーション
CN102812492B (zh) 用于心脏成像的方法及系统
CN103339652A (zh) 靠近伪影源的诊断图像特征
US20150065854A1 (en) Joint estimation of attenuation and activity information using emission data
US20040066909A1 (en) Multi modality imaging methods and apparatus
US9684973B2 (en) Systems and methods for selecting imaging data for principle components analysis
EP2245592B1 (en) Image registration alignment metric
US20140010428A1 (en) Method for extraction of a dataset from a medical image dataset and also medical imaging device
JP2009528139A (ja) リストモードデータに基づく局所動き補償
US8415630B2 (en) Apparatus and methods for determining a boundary of an object for positron emission tomography scatter correction
CN107111867A (zh) 多模态成像系统及方法
CN105849775A (zh) Pet/mr系统中的密度引导的衰减图生成
US20170084060A1 (en) Image Construction With Multiple Clustering Realizations
CN110536640B (zh) 从pet列表数据中的呼吸运动信号的噪声鲁棒的实时提取
CN104736062A (zh) Mr接收线圈定位和基于mr的衰减校正
CN106548473B (zh) 一种构建相位图像的方法及装置
US20230022425A1 (en) Apparatus, system, method and computer probram for providing a nuclear image of a region of interest of a patient
CN106548464B (zh) 一种构建相位图像的方法及装置
Tang et al. Effect of MR truncation compensation on quantitative PET image reconstruction for whole-body PET/MR

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