CN102934143A - 用于在pet-mr中产生衰减图的方法 - Google Patents
用于在pet-mr中产生衰减图的方法 Download PDFInfo
- Publication number
- CN102934143A CN102934143A CN2009801364695A CN200980136469A CN102934143A CN 102934143 A CN102934143 A CN 102934143A CN 2009801364695 A CN2009801364695 A CN 2009801364695A CN 200980136469 A CN200980136469 A CN 200980136469A CN 102934143 A CN102934143 A CN 102934143A
- Authority
- CN
- China
- Prior art keywords
- under inspection
- person under
- image
- geometric areas
- pad value
- 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.)
- Granted
Links
Images
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]
-
- 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
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/424—Iterative
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pulmonology (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Nuclear Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
磁共振(MR)图像分割处理器(32),其配置成使用受检者的MR图像识别所述受检者的一个或多个几何区域。发射数据重建处理器(40),其配置成通过向所述受检者的几何区域分配初始衰减值(52)来产生受检者的衰减图(54),并且配置成:(i)处理(56)从所述受检者采集的发射数据以产生所述受检者的发射图像(58),所述处理采用所述受检者的衰减图;(ii)基于使用所述受检者的发射图像计算的校正来更新(60)所述衰减图;并且(iii)迭代操作(i)和(ii)以迭代地产生所述受检者的重建的发射图像。
Description
下文涉及成像技术、发射断层摄影技术、正电子发射断层摄影(PET)技术、单光子发射计算断层摄影(SPECT)技术、磁共振技术和相关技术。具体参考PET/MR系统和协同式PET/MR成像应用对本发明进行了描述,但本发明也将适用于PET成像、SPECT成像和其他基于辐射发射的成像模态和用于这些的重建技术中。
通过使用适当的衰减图计入被成像受检者体内的吸收,来改进通过发射断层摄影,例如PET或SPECT进行的成像。结合透射计算断层摄影(CT)执行的发射断层摄影有利地得益于由CT模态提供的辐射衰减数据的可用性。重建的CT图像实质上是被成像受检者对产生CT图像数据期间使用的X射线辐射的衰减图。尽管CT中使用的X射线辐射一般与在PET中测量的511keV辐射或在SPECT或其他发射断层摄影技术中测量的发射是不同的,但已知能够通过适当缩放CT灰度级水平以计入辐射类型的差异,来从重建的透射CT图像产生用于发射断层摄影的衰减图。
人们对与磁共振(MR)成像协同组合来执行发射断层摄影也是感兴趣的。在这里,没有要从其产生衰减图的透射CT图像。使用MR图像作为衰减图来分析PET或SPECT数据的尝试已经失败了。这是因为MR中的对比度机制与PET或SPECT(或就此而论,CT)是根本不同的。因此,不能说“暗的”MR像素一定对应于高的或低的衰减。例如,对于一些MR成像模式而言骨骼组织和空气具有相似的灰度级亮度,但骨骼对发射辐射的衰减比空气的衰减要高得多。
在这种环境中,已经想过使用各种技术产生衰减图。在一种预期的方式中,采用典型受检者,例如典型人类受检者的衰减“图集”。衰减图集标识出典型受检者各部分或区域的衰减。不过,实际受检者,例如实际人类受检者,变动很大,将衰减图集适配到特定受检者并不是那么简单直接的。
预期的另一种方法是采用典型衰减图案的机器学习来构建一种将MR图像变换成适用于PET中的衰减图的变换算法。这种方法难以实施,并且机器学习方法的经验性质可能导致难以预测或估计的错误。
另一种方法是将衰减图的产生与发射数据重建结合起来。例如,参见Hawkins等人的美国专利No.6,310,968和Nuyts等人的以下文献:“Simultaneous Maximum A Posteriori Reconstruction of Attenuation andActivity Distributions from Emission Sinograms”,IEEE Trans.on MedicalImaging,vol.18,No.5,pp.393-403(1999)。有利的是,这些方法不需要使用外部提供的衰减图。不过,它们也不能利用可从对应MR图像获得的信息。有时会在使用这些单模态技术重建的图像中发现人为噪声。
下文提供了能够克服上述和其他问题的新的改进的设备和方法。
根据公开的一个方面,公开了一种结合受检者的磁共振(MR)图像和从受检者采集的发射数据进行工作的方法,所述方法包括:使用所述MR图像识别所述受检者的一个或多个几何区域;以及利用由所述几何区域界定的衰减图迭代地重建所述发射数据以产生所述受检者的发射图像,其中,所述衰减图具有基于所述迭代重建进行近似得到的对于所述几何区域的衰减值。
根据公开的另一方面,公开了一种发射数据重建处理器,其被配置成执行上一段所述的方法。
根据公开的另一方面,公开了一种设备,包括重建处理器,所述重建处理器被配置成使用受检者的MR图像识别所述受检者的一个或多个几何区域,通过向所述受检者的所述几何区域分配初始衰减值产生所述受检者的衰减图,并且:(i)处理从所述受检者采集的发射数据以产生所述受检者的发射图像,所述处理采用所述受检者的衰减图;(ii)基于使用所述受检者的发射图像计算的校正来更新所述衰减图;并且(iii)迭代操作(i)和(ii)以迭代地产生所述受检者的重建的发射图像。
根据公开的另一方面,公开了一种系统,包括:磁共振(MR)扫描器和MR重建处理器,二者配置成协作地产生受检者的MR图像;发射辐射检测器,其配置成采集所述受检者的发射数据;MR图像分割处理器,其配置成使用所述MR图像识别所述受检者的一个或多个几何区域;以及发射数据重建处理器,其配置成执行所述发射数据的迭代重建以产生所述受检者的发射图像,所述迭代重建利用随着所述迭代重建的迭代而迭代地更新的衰减图,至少一次迭代衰减图更新受到由所述MR图像分割处理器识别的所述受检者的几何区域约束。
根据公开的另一方面,公开了一种校正发射图像的衰减的方法,所述方法包括:分割磁共振(MR)图像以形成几何区域;向所述几何区域分配衰减值;基于所述几何区域的衰减值重建发射数据;以及基于所述重建发射数据校正所述几何区域的衰减值。
一个优点在于PET、SPECT或其他发射数据的更精确图像重建,获得很多新颖诊断应用所需的定量图像信息。
另一个优点在于将MR数据与PET、SPECT或其他发射数据协同耦合以增进后者的图像重建。
在阅读并理解以下详细描述之后,对于本领域普通技术人员而言,更多优点将是显而易见的。
图1图解示出了混合式MR/PET成像系统,包括PET重建处理器,PET重建处理器利用具有从对应MR图像导出的几何结构和从PET发射数据迭代导出的衰减值的衰减图。
图2图解示出了由PET重建处理器执行的处理。
参考图1,混合式扫描器包括磁共振(MR)扫描器10和正电子发射断层摄影(PET)检测器12的集成环,在图示实施例中磁共振扫描器为水平膛磁共振扫描器。MR扫描器10和PET检测器12的集成环配置成从设置于MR扫描器10内膛中的感兴趣区域14采集图像数据。(在概略的图1中,示出了具有集成PET检测器12的水平膛MR扫描器10,其被切掉了MR膛的一半和PET检测器12的环的对应一半,以显露出膛的内部和PET检测器12集成环的剩余一半)。MR扫描器10包括的部件例如有,用于产生静态(B0)磁场的主磁体,用于叠加磁场梯度的磁场梯度线圈,以及一个或多个用于激励和检测磁共振的射频线圈,为简单起见在图1中未示出这样的部件。PET检测器12的环包括电子器件的主干16,用于执行选定的“机上”操作,例如任选地执行辐射检测事件导致的电信号的模拟到数字转换,任选地执行辐射检测事件的数字时间戳操作等。备选地,可以由远程电子器件(未示出)执行这些操作中的一些。
MR扫描器10使用选定的空间编码采集磁共振(MR)数据,例如k空间样本,并且在MR成像数据缓冲器20中存储采集的MR数据。MR重建处理器22使用与所选空间编码相称的重建技术处理MR数据。例如,如果空间编码是在激励期间采用切片选择的梯度且在磁共振信号衰落和读出期间分别采用相位和读出编码梯度的常规笛卡尔编码,那么MR重建处理器22适当地采用基于傅里叶变换的重建技术。MR重建处理器22的输出是MR图像,该MR图像存储于在MR图像存储器24中并任选地显示于计算机26或其他显示装置上或以其他方式对其进行利用。
PET检测器12和关联的机上电子器件16和/或远程电子器件(未示出)检测辐射检测事件并执行能量和时间开窗操作,以识别表示电子-正电子湮没事件的基本同时的511keV检测事件。每一对基本同时的511keV检测事件界定连接两个511keV检测事件的投射或响应线。这一组数据被称为PET数据,或更一般地称为发射数据,其被存储在PET成像数据缓冲器30中。
在一些实施例中,PET检测器12具有充分高的时间分辨率以分辨源于相同电子-正电子湮没事件的两个511keV伽马粒子的基本同时检测之间的飞行时间时间差(或没有时间差),机上和/或远程电子器件还被配置成沿着响应线定位电子-正电子湮没事件,从而产生飞行时间PET数据。
图示的混合式扫描器10、12、16是范例。更一般地,这里公开的用于执行发射数据重建的技术适当地与任何PET扫描器相结合来实践,或者更一般地,与产生发射数据的任何扫描器相结合来实践,例如图示的PET检测器12、独立的PET检测器、产生单光子发射计算断层摄影(SPECT)数据的伽马射线摄影机等。进一步地,这里公开的执行发射数据重建的技术适当地与从其采集发射数据的相同受检者的任何MR图像相结合来实践,这种MR图像是由图示的混合式扫描器10、12、16或由独立的MR扫描器采集的。
继续参考图1,PET重建处理器32采用迭代重建技术重建由PET系统12、16(或更一般地,由诸如PET、SPECT等辐射发射成像系统)采集的发射数据,以形成重建的发射图像,重建的发射图像被存储于PET图像存储器34中并任选地显示于计算机26或另一种显示装置上或以其他方式被利用。
为了执行发射图像重建,PET重建处理器32采用计入受检者体内对发射辐射再吸收导致的发射损耗的衰减图。使用由MR图像提供的几何信息与从发射数据集自身导出的衰减值信息的组合产生适当的衰减图。
为此,由MR图像分割处理器40分割MR图像以识别受检者的一个或多个几何区域,并在MR分割图像存储器42中适当存储这种几何信息。分割处理器40能够采用任何适当的分割过程刻画受检者解剖学上的不同区域(假设人类或动物受检者具有解剖结构)。一些适当的分割过程包括,例如:基于阈值的分割方法;基于索贝尔(Sobel)算子的分割方法;区域生长分割方法;分水岭分割方法;基于模型的分割等。任选地,在计算机26上显示由分割处理器40识别的受检者的一个或多个几何区域供操作人员查看,操作人员可以选择使用一个或多个计算机26的接口装置来调节或以其他方式修改自动识别的受检者一个或多个几何区域。由分割处理器40从MR图像识别所得的受检者的一个或多个几何区域充当对PET重建处理器32的输入,并且PET重建处理器32在构建发射数据重建中使用的衰减图时使用由所识别的受检者的一个或多个几何区域表示的关于受检者的几何信息。
继续参考图1,并进一步参考图2,描述由PET重建处理器32采用的适当重建过程。由PET重建处理器32执行初始衰减图产生操作50以产生初始衰减图,用于在迭代发射数据重建过程的初始迭代中使用。通过向由MR图像分割识别的受检者的一个或多个几何区域分配初始衰减值52来构建初始衰减图。能够通过各种方式获得初始衰减值52。在一种方式中,受检者的一个或多个几何区域界定受检者的外部轮廓,一开始,为受检者轮廓之内的每个体素或像素分配默认衰减值,例如,与水的衰减值对应的衰减值,并向设置于受检者轮廓之外的体素或像素分配零衰减值(假设周围环境是空气,对于感兴趣的大部分发射辐射而言空气具有低吸收)。在这样的实施例中,受检者的一个或多个几何区域可以是界定受检者轮廓的单一区域,或者可以是一起共同界定受检者轮廓的多个几何区域。
在其他实施例中,初始衰减值52包括基于几何区域的几何结构选择的对于每个几何区域的默认衰减值。例如,基于一些或所有几何区域的形状或位置可以识别那些几何区域的组织标识。通过这种方式,例如,可以将不同区域标记为骨骼区域、脂肪区域、肌肉区域等,并可以基于其组织类型向每个这样被标记的区域分配适当的默认衰减值。
也可以使用这些方式的组合来选择初始衰减值52,例如,基于组织类型来标记一些几何区域并基于组织类型标记来分配衰减值,然后向不能基于不可识别区域的形状或位置判断其组织类型的任何区域分配默认衰减值,例如水的衰减值。
通过衰减值分配操作50产生衰减图54的初始迭代。由PET重建处理器32执行的发射数据重建迭代56使用衰减图54来产生重建图像58。迭代地进行由重建处理器32执行的发射数据重建过程。在图2中,用迭代索引n来标注迭代发射数据重建过程的当前迭代,并将第n次发射数据重建迭代56中使用的衰减图54标注为第n衰减图54,将由第n次发射数据重建迭代56输出的重建图像58标注为第n重建图像58。
基于第n重建图像58,在操作60中计算衰减值校正,并使用这些校正来产生更新的衰减图,亦即,第(n+1)衰减图62,其用于执行迭代重建过程的第(n+1)次迭代,等等。继续进行迭代处理,直到符合停止标准,例如相继重建图像之间的相似性度量高于阈值。将重建过程的最后迭代产生的重建图像输出到PET图像存储器34作为最终的重建PET图像。尽管未示出,但还预期存储、显示或以其他方式利用最后产生的衰减图。
基本可以采用适于重建发射数据的任何迭代重建算法用于迭代56,例如最大似然期望最大化(MLEM)算法。衰减值校正60也能够采用任何适当的校正算法,例如在如下文献中公开的那些算法之一:Nuyts等人的“Simultaneous Maximum A Posteriori Reconstruction of Attenuation andActivity Distributions from Emission Sinograms”,IEEE Trans.on MedicalImaging,vol.18,no.5,pp.393-403(1999),在此通过引用将其全文并入;或Hawkins等人的美国专利No.6,310,968,在此也通过引用将其全文并入。用于第n次发射数据重建迭代56、采用MLEM的适当方法如下:
其中:n表示已经描述的重建迭代序号;j是感兴趣区域中的体素的索引;i是检测基本同时的511keV辐射的PET检测器对12的索引(对于SPECT而言,该索引i适当地标注单一检测器);λj表示索引为j的体素中的估计活动;yi表示在索引为i的检测器对中测量的并发光子;cij表示系统矩阵项目,该系统矩阵项目表示索引为i的检测器对对于索引为j的体素的灵敏度;并
且
其中:cik表示系统矩阵项目,该系统矩阵项目表示索引为i的检测器对对于索引为k的体素的灵敏度;lik表示索引为i的检测器对的响应线和索引为k的体素之间的交互长度;λk表示索引为k的体素中的估计活动;而μk表示索引为k的体素的衰减系数,对于第n次迭代而言该衰减系数是由第n衰减图54提供的。对于飞行时间PET数据,方程(1)如下适当地并入电子-正电子湮没事件的飞行时间定位:
用于执行衰减值校正60的适当方法是正向投射更新的活动图像以为检测器对i产生投射数据或响应线yi,并根据下式计算校正场Kj:
其中:N表示重建体积的直径;而αp表示松弛因子。根据下式使用校正场Kj更新衰减图:
其中:Rr表示索引为r的几何区域,其中,索引r遍布存储器42中存储的由MR图像分割识别的区域;表示对于第n衰减图54的区域特异性衰减系数;表示对于第(n+1)衰减图62的区域特异性衰减系数;而方程(6)中的分数项表示对于区域Rr的平均衰减值校正。
方程(6)的衰减图更新逐个对每个几何区域进行操作,每个几何区域Rr被分配以对于该区域所有像素或体素的单一更新衰减值方程(6)的衰减图更新向每个几何区域分配大致不同的衰减值。“大致不同”的意思是衰减图更新向每个区域分配衰减值,该衰减值一般与分配给其他区域的衰减值不同,但在有些情况下可以与分配给另一区域或多个区域的衰减值相同。例如,如果两个不同区域对应于公共组织类型(例如,均表示骨骼区域),衰减图更新可能向对应于相同组织类型的这两个不同区域分配相同或类似的衰减值。
预期的另一种变体是针对前几次迭代根据方程(6)执行每个区域的衰减图更新,然后针对后续迭代根据切换到每个像素或每个体素的衰减图更新。这种方式允许初始收敛到相对于区域优化的衰减图,其需要优化相对少的参数并应当给出与实际受检者衰减的很好近似;然后继之以从相对于区域优化的(接近正确的)衰减图开始的计算强度更大的逐个像素或逐个体素的优化。从每个区域衰减图更新切换到每个像素或每个体素衰减图更新的时间选择可以基于各种标准。在一简单方式中,使用每个区域的衰减图更新进行固定数目的第一迭代,例如用于前三次迭代,然后使用每个像素或每个体素的更新进行剩余迭代,例如用于第四和后面的迭代。在一种更复杂的方法中,能够监测切换标准,例如计算从一次迭代到下一次之间任何区域衰减值的最大变化。在这个最大变化下降到阈值以下时,则表示每个区域的衰减图优化接近收敛,切换到每个像素或每个体素的更新。
在预期的另一种变体中,能够针对被认为具有基本均匀吸收的区域中的像素或体素执行每个区域的衰减值更新,而能够针对预计具有更多不均匀吸收的区域中的像素或体素执行每个像素或每个体素的衰减值更新。将区域指派为均匀的还是不均匀的可以基于对该区域可能的组织类型或器官特征的先验知识,或者可以基于该区域上的定量度量,例如校正场Kj的变化。在后一种方式中,例如,逐个像素或逐个体素地更新校正场Kj变化超过阈值量的任何区域,假设整个区域上校正场Kj的大变化表示该区域具有不均匀吸收。
再者,能够使用除方程(5)和(6)的更新或更新之外的更新。例如,在一些实施例中,如这里所述那样来计算值或但基于所计算的或值从一组标准值中选择分配给衰减图的区域或体素或像素的衰减值。例如,标准值可以包括分别针对空气、骨骼、脂肪、组织和血液的那么,对于要更新的每个区域或体素或像素,分配标准值中最接近计算值或的衰减值。
在上文中,假设MR图像和发射数据采用相同的受检者坐标系。对于图示的混合式扫描器10、12、16可能是这种情况。在由相应的分立或独立成像系统采集MR图像和发射数据时,能够使用设置于受检者身上且在MR和发射系统中都可见的基准标记来实现MR和发射坐标系的对应。或者,一开始在不考虑吸收的情况下对PET数据进行任选的“粗略”重建,并将所得的“粗略”图像在空间上与MR图像对准或配准,以确定用于MR图像和发射数据的公共受检者坐标系。然后,能够将参考图2公开的迭代发射数据重建算法应用于空间上配准的MR图像和发射数据。
在一些预期的实施例中,对衰减图进行迭代更新还包括调节几何区域的尺寸和/或形状。在一种适当方式中,针对几何区域中的至少一个,例如使用基于多边形或样条的表面模型定义形状模型。能够将形状模型拟合到经分割MR图像的相关图段,以定义初始形状模型参数。之后,当在每次迭代中更新衰减图时,使用衰减值校正60来更新形状模型。例如,可以预计在边界处校正场Kj是大的,在边界处,当前的表面模型与发射数据中表示的区域形状不一致。然后能够根据从迭代产生的校正场Kj的信息迭代地调节形状模型。例如,如果要适应MR和PET数据采集之间患者的变化,这样的几何校正是有利的。已知在人类受检者的腹部区域和某些其他区域中这样的变化相对频繁。例如,能够使用校正场Kj来调整表示胃、腹部外轮廓等的区域形状来适应数据采集之间可能发生的典型变化。
可以通过各种方式实现这里公开的各种计算部件22、32、40,例如,通过计算机或包括数字处理器并被编程控制或包括固件以执行所公开的处理的其他装置,或通过配置成执行所公开处理或其部分等的混合式或模拟电路。在一些实施例中,可以由图示的具有适当固件或程序设计的计算机26实现计算部件22、32、40。这里公开的MR和发射数据和图像处理方法可以通过这样的处理器或其他硬件来实现,和/或可以具体化为存储指令的存储介质,在由这样的处理器或其他硬件执行该指令时,执行所公开的方法。可以通过一种或多种存储媒质实现这样的存储介质,例如如下一种或多种:磁盘;光盘;闪速存储器或其他静电存储器;随机存取存储器(RAM);只读存储器(ROM);等等。
已经参考优选实施例描述了本发明。在阅读并理解了前述详细说明的同时,本领域技术人员可以想到修改和变化。这意味着,应当将本发明推断为包括所有此类落在权利要求及其等同要件的范围内的修改和变化。
Claims (25)
1.一种结合受检者的磁共振(MR)图像和从所述受检者采集的发射数据工作的方法,所述方法包括:
使用所述MR图像识别所述受检者的一个或多个几何区域;以及
利用由所述几何区域界定的衰减图迭代地重建所述发射数据以产生所述受检者的发射图像,所述衰减图具有基于所述迭代重建进行近似得到的对于所述几何区域的衰减值。
2.根据权利要求1所述的方法,还包括:
显示所述发射图像。
3.根据权利要求1-2中的任一项所述的方法,还包括:
使用正电子发射断层摄影(PET)和单光子发射计算断层摄影(SPECT)之一从所述受检者采集所述发射数据。
4.根据权利要求1-3中的任一项所述的方法,还包括:
基于通过所述迭代重建产生的所述发射图像的迭代更新迭代地更新对于所述几何区域的所述衰减值。
5.根据权利要求4所述的方法,还包括:
通过产生对于所述几何区域的初始衰减值启动所述迭代重建,所述初始衰减值是如下之一:
(i)对于所有几何区域的相同衰减值,以及
(ii)基于所述几何区域的形状、位置和尺寸中的至少一种选择的对于每个几何区域的默认衰减值。
6.根据权利要求4-5中的任一项所述的方法,其中,对于所述几何区域的所述衰减值的所述迭代更新包括:
对于所述迭代重建中的至少一些迭代逐个像素或逐个体素地更新衰减值。
7.根据权利要求4-6中的任一项所述的方法,其中,对于所述几何区域的所述衰减值的所述迭代更新包括:
对于所述迭代重建中的至少一些迭代逐个区域地更新衰减值,其中,逐个区域更新衰减值包括向给定区域的所有像素或体素分配所述相同衰减值。
8.根据权利要求4-7中的任一项所述的方法,其中,对于所述几何区域的所述衰减值的所述迭代更新包括:
在所述迭代重建的每次迭代之后,基于从由所述迭代重建的最近迭代产生的经更新的发射图像导出的衰减信息,更新对于所述区域的衰减值。
9.根据权利要求8所述的方法,其中,所述衰减值的所述更新包括:
正向投射由所述迭代重建的所述最近迭代产生的所述经更新的发射图像以产生正向投射数据;以及
基于所述正向投射数据和从所述受检者采集的所述发射数据的比较计算衰减值校正。
10.根据权利要求1-9中的任一项所述的方法,其中,所述受检者为动物或人类受检者,且使用所述MR图像识别所述受检者的几何区域包括:
分割所述MR图像以识别多个解剖学上的不同区域。
11.根据权利要求1-9中的任一项所述的方法,其中,所述受检者为动物或人类受检者,且使用所述MR图像识别所述受检者的一个或多个几何区域包括:
识别所述MR图像中包括所述受检者的外轮廓的单一区域。
12.根据权利要求1-10中的任一项所述的方法,其中,在所述迭代重建中,由基于所述迭代重建更新的形状模型表示所述几何区域中的至少一个的几何结构。
13.一种配置成执行根据权利要求1和4-12中的任一项所述的方法的发射数据重建处理器(32,40)。
14.一种设备,包括:
重建处理器(32,40),其配置成使用受检者的MR图像识别所述受检者的一个或多个几何区域,通过向所述受检者的所述几何区域分配初始衰减值(52)产生所述受检者的衰减图(54),并且所述重建处理器配置成:
(i)处理(56)从所述受检者采集的发射数据以产生所述受检者的发射图像(58),所述处理采用所述受检者的所述衰减图,
(ii)基于使用所述受检者的所述发射图像计算的校正来更新(60)所述衰减图,并且
(iii)迭代操作(i)和(ii)以迭代地产生所述受检者的重建的发射图像。
15.根据权利要求14所述的设备,还包括:
显示器(26),其配置成显示所述受检者的所述重建的发射图像。
16.根据权利要求14-15中的任一项所述的设备,其中,所述更新操作(ii)(60)包括:
(ii)(a)正向投射所述受检者的所述发射图像(58)以产生正向投射数据,以及
(ii)(b)基于所述正向投射数据和从所述受检者采集的所述发射数据的比较来计算所述衰减图(54)的校正。
17.根据权利要求14-16中的任一项所述的设备,还包括:
MR扫描器(10),其配置成产生所述受检者的所述MR图像;以及
单光子发射计算断层摄影(SPECT)或正电子发射断层摄影(PET)检测器(12),其配置成采集从所述受检者采集的所述发射数据。
18.根据权利要求14-17中的任一项所述的设备,其中,所述重建处理器(32,40)包括MR图像分割处理器(40),所述MR图像分割处理器配置成对所述受检者的所述MR图像进行分割以识别所述受检者的所述一个或多个几何区域。
19.根据权利要求14-18中的任一项所述的设备,其中,所述重建处理器(32,40)配置成使用所述受检者的所述MR图像识别所述受检者的多个几何区域,并且所述衰减图更新操作(ii)(60)向所述几何区域中的每个分配大致不同的衰减值。
20.根据权利要求14-18中的任一项所述的设备,其中,所述衰减图更新操作(ii)(60)向所述衰减图(54)的每个像素或体素分配大致不同的衰减值。
21.根据权利要求14-18中的任一项所述的设备,其中,所述重建处理器(32,40)配置成使用所述受检者的所述MR图像识别所述受检者的多个几何区域,并且所述衰减图更新操作(ii)(60):
为最初一次或多次所述迭代(iii)逐个几何区域地分配衰减值,以及
为后续一次或多次所述迭代(iii)逐个像素或逐个体素地分配衰减值。
22.一种系统,包括:
磁共振(MR)扫描器(10)和MR重建处理器(22),二者配置成协作地产生受检者的MR图像;
发射辐射检测器(12),其配置成采集所述受检者的发射数据;
MR图像分割处理器(40),其配置成使用所述MR图像识别所述受检者的一个或多个几何区域;以及
发射数据重建处理器(32),其配置成执行所述发射数据的迭代重建以产生所述受检者的发射图像,所述迭代重建利用随着所述迭代重建的迭代而迭代地更新的衰减图(54),至少一次迭代衰减图更新受到由所述MR图像分割处理器识别的所述受检者的所述几何区域约束。
23.根据权利要求22所述的系统,其中,约束初始迭代衰减图更新以向由所述MR图像分割处理器识别的选定几何区域的所有像素或体素分配奇异衰减值。
24.一种校正发射图像的衰减的方法,所述方法包括:
分割磁共振(MR)图像以形成几何区域;
向所述几何区域分配衰减值;
基于所述几何区域的所述衰减值重建发射数据;以及
基于所重建的发射数据校正所述几何区域的所述衰减值。
25.根据权利要求24所述的方法,其中,所述校正包括:
正向投射所重建的发射数据以产生正向投射数据;以及
基于所述正向投射数据和所述发射数据的比较来校正所述几何区域的所述衰减值。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US9825408P | 2008-09-19 | 2008-09-19 | |
US61/098,254 | 2008-09-19 | ||
PCT/IB2009/053946 WO2010032168A2 (en) | 2008-09-19 | 2009-09-09 | Method for generation of attenuation map in pet-mr |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102934143A true CN102934143A (zh) | 2013-02-13 |
CN102934143B CN102934143B (zh) | 2016-04-27 |
Family
ID=42039964
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200980136469.5A Expired - Fee Related CN102934143B (zh) | 2008-09-19 | 2009-09-09 | 用于在pet-mr中产生衰减图的方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US8600136B2 (zh) |
EP (1) | EP2340524B1 (zh) |
JP (1) | JP5530446B2 (zh) |
CN (1) | CN102934143B (zh) |
WO (1) | WO2010032168A2 (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105308470A (zh) * | 2013-06-07 | 2016-02-03 | 皇家飞利浦有限公司 | 用于磁共振协议的视觉预扫描患者信息 |
CN105339804A (zh) * | 2013-06-20 | 2016-02-17 | 皇家飞利浦有限公司 | 根据MR Dixon数据的皮质骨分割 |
CN106846430A (zh) * | 2014-11-21 | 2017-06-13 | 上海联影医疗科技有限公司 | 一种图像重建方法 |
CN107644421A (zh) * | 2016-07-20 | 2018-01-30 | 上海联影医疗科技有限公司 | 医学图像分割方法及系统 |
CN107871331A (zh) * | 2016-09-28 | 2018-04-03 | 通用电气公司 | 用于重构发射活动图像的系统和方法 |
US11543482B2 (en) * | 2018-10-16 | 2023-01-03 | Koninklijke Philips N.V. | Magnetic resonance imaging using motion-compensated image reconstruction |
Families Citing this family (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2161594B1 (en) * | 2008-09-09 | 2017-11-08 | Multi Magnetics Incorporated | System and method for correcting attenuation in hybrid medical imaging |
DE102008047840B4 (de) * | 2008-09-18 | 2018-08-30 | Siemens Healthcare Gmbh | Verfahren zur wenigstens teilweisen Ermittlung und/oder Anpassung einer zur Schwächungskorrektur von Positronenemissionstomographie-Bilddatensätzen genutzten Schwächungskarte in einem kombinierten Magnetresonanz-Positronenemissionstomographie-Gerät |
US8242777B2 (en) * | 2009-09-28 | 2012-08-14 | Siemens Aktiengesellschaft | Calibration of an emission tomography subsystem |
WO2012025855A1 (en) | 2010-08-25 | 2012-03-01 | Koninklijke Philips Electronics N.V. | Dual modality imaging including quality metrics |
CN103299345B (zh) * | 2011-01-10 | 2016-10-19 | 皇家飞利浦电子股份有限公司 | 双能量断层摄影成像系统 |
US9336613B2 (en) * | 2011-05-24 | 2016-05-10 | Koninklijke Philips N.V. | Apparatus for generating assignments between image regions of an image and element classes |
EP2715664A1 (en) * | 2011-05-24 | 2014-04-09 | Koninklijke Philips N.V. | Apparatus and method for generating an attenuation correction map |
US8913810B2 (en) * | 2011-07-26 | 2014-12-16 | Siemens Medical Solutions Usa, Inc. | Simultaneous reconstruction of emission activity and attenuation coefficient distribution from TOF data, acquired with external shell source |
EP2833163B1 (en) * | 2012-03-28 | 2018-07-25 | National Institutes for Quantum and Radiological Science and Technology | Method for generating pet absorption-corrected image from mr image and computer program |
WO2013144799A1 (en) | 2012-03-29 | 2013-10-03 | Koninklijke Philips N.V. | Mri method for assigning individual pixels or voxels tissue - specific pet attenuation values |
KR101331504B1 (ko) * | 2012-05-22 | 2013-11-26 | 서울대학교산학협력단 | 핵의학 방출 영상 시스템 및 단일 투과 영상을 이용한 감쇠보정 방법 |
US8958660B2 (en) * | 2012-06-22 | 2015-02-17 | General Electric Company | Method and apparatus for iterative reconstruction |
US20150065854A1 (en) * | 2012-10-31 | 2015-03-05 | General Electric Company | Joint estimation of attenuation and activity information using emission data |
US9928617B2 (en) * | 2012-11-08 | 2018-03-27 | The General Hospital Corporation | System and method for multi-modality time-of-flight attenuation correction |
CN103908278B (zh) * | 2013-01-06 | 2017-07-28 | 上海联影医疗科技有限公司 | 图像重建方法及装置、医疗成像系统 |
US11922604B2 (en) | 2014-10-20 | 2024-03-05 | Koninklijke Philips N.V. | Classified truncation compensation |
US9474495B2 (en) | 2014-12-22 | 2016-10-25 | General Electric Company | System and method for joint estimation of attenuation and activity information |
US20160327622A1 (en) * | 2015-05-05 | 2016-11-10 | General Electric Company | Joint reconstruction of activity and attenuation in emission tomography using magnetic-resonance-based priors |
US10078889B2 (en) | 2015-08-25 | 2018-09-18 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image calibration |
WO2018006419A1 (en) | 2016-07-08 | 2018-01-11 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for generating attenuation map |
US10593071B2 (en) * | 2017-04-14 | 2020-03-17 | Siemens Medical Solutions Usa, Inc. | Network training and architecture for medical imaging |
WO2018236748A1 (en) | 2017-06-19 | 2018-12-27 | Washington University | DIGITAL ASSISTED IMAGE RECONSTRUCTION FOR TOMOGRAPHIC IMAGING |
KR101958099B1 (ko) * | 2017-07-31 | 2019-07-04 | 인제대학교 산학협력단 | 양전자 방출 단층 촬영기의 비행시간 정보를 이용한 감쇠보정 방법 |
US20210012546A1 (en) * | 2018-03-26 | 2021-01-14 | Koninklijke Philips N.V. | Automatic fault detection in hybrid imaging |
KR102542505B1 (ko) * | 2020-11-17 | 2023-06-13 | (주)자비스 | 학습데이터 생성 방법, 장치, 프로그램 및 이를 이용한 이물질 검출 방법 |
US11854126B2 (en) * | 2021-07-07 | 2023-12-26 | Siemens Medical Solutions Usa, Inc. | Methods and apparatus for deep learning based image attenuation correction |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6310968B1 (en) * | 1998-11-24 | 2001-10-30 | Picker International, Inc. | Source-assisted attenuation correction for emission computed tomography |
US20030146388A1 (en) * | 1999-06-06 | 2003-08-07 | Naor Wainer | Pet and spect systems with attenuation correction |
US20040161140A1 (en) * | 2003-02-19 | 2004-08-19 | Isabella Yuzefovich | Nuclear imaging |
CN101156780A (zh) * | 2007-09-04 | 2008-04-09 | 陈武凡 | 在pet成像中最大后验优化图像重建方法 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5376795A (en) * | 1990-07-09 | 1994-12-27 | Regents Of The University Of California | Emission-transmission imaging system using single energy and dual energy transmission and radionuclide emission data |
US5672877A (en) * | 1996-03-27 | 1997-09-30 | Adac Laboratories | Coregistration of multi-modality data in a medical imaging system |
US6384416B1 (en) | 1998-04-27 | 2002-05-07 | Duke University | Transmission scanning technique for gamma-camera coincidence imaging |
DE10231061A1 (de) * | 2002-07-10 | 2004-01-22 | Philips Intellectual Property & Standards Gmbh | Verfahren und System zur Verbesserung des Informationsgehaltes in einem Bild |
AU2003276658A1 (en) * | 2002-11-04 | 2004-06-07 | V-Target Technologies Ltd. | Apparatus and methods for imaging and attenuation correction |
US20060004274A1 (en) * | 2004-06-30 | 2006-01-05 | Hawman Eric G | Fusing nuclear medical images with a second imaging modality |
US8108024B2 (en) * | 2004-12-15 | 2012-01-31 | Koninklijke Philips Electronics N.V. | Registration of multi-modality images |
DE602005020844D1 (de) * | 2004-12-17 | 2010-06-02 | Koninkl Philips Electronics Nv | Kompensationsalgorithmus bei abschneidung für iterative rekonstruktion |
US7402807B2 (en) * | 2005-11-02 | 2008-07-22 | Siemens Medical Solutions Usa, Inc. | Method for reducing an electronic time coincidence window in positron emission tomography |
US7348564B2 (en) | 2005-12-12 | 2008-03-25 | General Electric Company | Multi modality imaging methods and apparatus |
WO2008064319A2 (en) | 2006-11-22 | 2008-05-29 | The General Hospital Corporation | Attenuation correction of pet image using image data acquired with an mri system |
DE102007034953B4 (de) * | 2007-07-26 | 2016-09-22 | Siemens Healthcare Gmbh | Verfahren zur Bewegungsvorgänge berücksichtigenden Aufnahme von Messdaten eines Patienten und zugehörige medizinische Einrichtung |
US8098916B2 (en) * | 2007-10-30 | 2012-01-17 | General Electric Company | System and method for image-based attenuation correction of PET/SPECT images |
-
2009
- 2009-09-09 WO PCT/IB2009/053946 patent/WO2010032168A2/en active Application Filing
- 2009-09-09 US US13/062,224 patent/US8600136B2/en not_active Expired - Fee Related
- 2009-09-09 JP JP2011527435A patent/JP5530446B2/ja not_active Expired - Fee Related
- 2009-09-09 CN CN200980136469.5A patent/CN102934143B/zh not_active Expired - Fee Related
- 2009-09-09 EP EP09787150.3A patent/EP2340524B1/en not_active Not-in-force
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6310968B1 (en) * | 1998-11-24 | 2001-10-30 | Picker International, Inc. | Source-assisted attenuation correction for emission computed tomography |
US20030146388A1 (en) * | 1999-06-06 | 2003-08-07 | Naor Wainer | Pet and spect systems with attenuation correction |
US20040161140A1 (en) * | 2003-02-19 | 2004-08-19 | Isabella Yuzefovich | Nuclear imaging |
CN101156780A (zh) * | 2007-09-04 | 2008-04-09 | 陈武凡 | 在pet成像中最大后验优化图像重建方法 |
Non-Patent Citations (2)
Title |
---|
ALVARO DE PIERRO ET AL.: "Simultaneous Activity Attenuation Reconstruction in Positron Emission Tomography via Maximum Likelihood and Iterative Methods", 《IEEE INTERNATIONAL CONFERENCE ON IMAGE PROCESSING,2005》 * |
HABIB ZAIDI ET AL.: "Strategies for attenuation compensation in neurological PET studies", 《NEUROIMAGE》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105308470A (zh) * | 2013-06-07 | 2016-02-03 | 皇家飞利浦有限公司 | 用于磁共振协议的视觉预扫描患者信息 |
CN105308470B (zh) * | 2013-06-07 | 2019-11-05 | 皇家飞利浦有限公司 | 用于磁共振协议的视觉预扫描患者信息 |
CN105339804A (zh) * | 2013-06-20 | 2016-02-17 | 皇家飞利浦有限公司 | 根据MR Dixon数据的皮质骨分割 |
CN105339804B (zh) * | 2013-06-20 | 2019-03-01 | 皇家飞利浦有限公司 | 根据MR Dixon数据的皮质骨分割 |
CN106846430A (zh) * | 2014-11-21 | 2017-06-13 | 上海联影医疗科技有限公司 | 一种图像重建方法 |
CN106846430B (zh) * | 2014-11-21 | 2020-06-26 | 上海联影医疗科技有限公司 | 一种图像重建方法 |
CN107644421A (zh) * | 2016-07-20 | 2018-01-30 | 上海联影医疗科技有限公司 | 医学图像分割方法及系统 |
CN107644421B (zh) * | 2016-07-20 | 2021-03-09 | 上海联影医疗科技股份有限公司 | 医学图像分割方法及系统 |
CN107871331A (zh) * | 2016-09-28 | 2018-04-03 | 通用电气公司 | 用于重构发射活动图像的系统和方法 |
CN107871331B (zh) * | 2016-09-28 | 2023-07-14 | 通用电气公司 | 用于重构发射活动图像的系统和方法 |
US11543482B2 (en) * | 2018-10-16 | 2023-01-03 | Koninklijke Philips N.V. | Magnetic resonance imaging using motion-compensated image reconstruction |
Also Published As
Publication number | Publication date |
---|---|
WO2010032168A3 (en) | 2014-04-24 |
EP2340524B1 (en) | 2019-06-05 |
WO2010032168A2 (en) | 2010-03-25 |
US8600136B2 (en) | 2013-12-03 |
CN102934143B (zh) | 2016-04-27 |
EP2340524A2 (en) | 2011-07-06 |
US20110158497A1 (en) | 2011-06-30 |
JP2012506530A (ja) | 2012-03-15 |
JP5530446B2 (ja) | 2014-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102934143B (zh) | 用于在pet-mr中产生衰减图的方法 | |
CN102067176B (zh) | 结合局部运动监测、校正和评估的辐射成像 | |
CN107644421B (zh) | 医学图像分割方法及系统 | |
CN102483852B (zh) | 利用基于飞行时间信息逐个事件生成的图像内容的飞行时间正电子发射断层摄影重建 | |
US10943349B2 (en) | Positron emission tomography system and image reconstruction method using the same | |
EP3226766B1 (en) | System and method for image calibration | |
CN101959454B (zh) | 通过磁共振图像的自动分割进行ct替代 | |
CN1809841B (zh) | 运动补偿的重建方法、设备与系统 | |
CN101849247B (zh) | 用于生成衰减图的设备和方法 | |
CN107133996A (zh) | 产生用于pet数据重建的衰减图的方法及pet/ct系统 | |
CN104700438B (zh) | 图像重建方法及装置 | |
CN106920246A (zh) | 在存在金属伪影的情况下用于分割的不确定性图 | |
CN110678906B (zh) | 用于定量分子成像的准确混合数据集的生成 | |
CN106725570A (zh) | 成像方法及系统 | |
CN101305297A (zh) | 使用解剖列表模式掩模的pet成像 | |
CN105147312A (zh) | Pet图像获取方法及系统 | |
US20190090826A1 (en) | Systems and methods for improved diagnostics for nuclear medicine imaging | |
CN109961419B (zh) | 对pet活度分布图像进行衰减校正的校正信息获取方法 | |
US20210398329A1 (en) | Artificial intelligence (ai)-based standardized uptake vaule (suv) correction and variation assessment for positron emission tomography (pet) | |
CN111373445A (zh) | 对具有重叠的全身正电子发射断层摄影(pet)扫描重建图像,并且改变个体床位的曝光时间 | |
CN110070588A (zh) | Pet图像重建方法、系统、可读存储介质和设备 | |
US20230056685A1 (en) | Methods and apparatus for deep learning based image attenuation correction | |
WO2024126321A1 (en) | Bilateral cross filtering of semantic class probabilities |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160427 Termination date: 20200909 |