CN101842806A - 脏同位素pet重建 - Google Patents
脏同位素pet重建 Download PDFInfo
- Publication number
- CN101842806A CN101842806A CN200780010626A CN200780010626A CN101842806A CN 101842806 A CN101842806 A CN 101842806A CN 200780010626 A CN200780010626 A CN 200780010626A CN 200780010626 A CN200780010626 A CN 200780010626A CN 101842806 A CN101842806 A CN 101842806A
- Authority
- CN
- China
- Prior art keywords
- data
- correction
- false
- indication
- generate
- 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
- 238000012937 correction Methods 0.000 claims abstract description 94
- 238000000034 method Methods 0.000 claims abstract description 43
- 238000003384 imaging method Methods 0.000 claims abstract description 20
- 238000009826 distribution Methods 0.000 claims description 32
- 230000005855 radiation Effects 0.000 claims description 22
- 238000002591 computed tomography Methods 0.000 claims description 17
- 230000000694 effects Effects 0.000 claims description 11
- 238000007689 inspection Methods 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 238000002600 positron emission tomography Methods 0.000 description 31
- 238000005516 engineering process Methods 0.000 description 14
- 230000006870 function Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000000155 isotopic effect Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000002285 radioactive effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 238000007519 figuring Methods 0.000 description 2
- 239000000700 radioactive tracer Substances 0.000 description 2
- 239000012217 radiopharmaceutical Substances 0.000 description 2
- 229940121896 radiopharmaceutical Drugs 0.000 description 2
- 230000002799 radiopharmaceutical effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000010523 cascade reaction Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005264 electron capture Effects 0.000 description 1
- 230000005263 electron capture decay Effects 0.000 description 1
- 238000000295 emission spectrum Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- 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/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic 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)
-
- 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
-
- 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)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Nuclear Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
一种用于脏同位素正电子成像中的方法使用关于对象(118)的测量特性的信息生成假符合校正。将假成像校正应用于来自对象的正电子成像检查的数据。使用校正数据产生指示对象的人可读图像。
Description
本发明涉及正电子发射断层摄影的领域,具体而言涉及对从脏同位素(dirty isotope)的使用中产生的假符合的校正。其特别应用于医学成像以及希望降低假符合影响的其他应用。
正电子发射断层摄影(PET)是核医学的一个分支,其中正电子发射的放射性药物被引入患者或其他受检查对象的身体中。常规的PET同位素仅仅通过正电子发射衰变。正电子在被称为正电子湮灭事件中与电子反应,由此产生沿着响应线(LOR)朝相反方向传播的511keV伽玛射线的符合对。在符合时间内探测到的伽玛射线对通常由PET扫描器作为湮灭事件进行记录。
然而,实际上,康普顿散射可以影响511keV光子之一或两者。这样的散射可以使得符合事件错误发生在受检查的对象内部,甚至在散射介质外部,从而导致图像质量的降低。因此,单散射模拟(SSS)技术已用于校正从常规放射性同位素的衰变产生的511keV符合对的康普顿散射。参见Accorsi等人的“Implementation of a Single Scatter Simulation Algorithm for3D PET:Application to Emission and Transmission Scanning”,Proceedings ofIEEE Nuclear Science Symposium and Medical Imaging Conference(2002)。
根据以下关系产生校正投影数据:
方程1
P校正=P测量-a×S康普顿-b×CF
其中,P校正是校正投影数据,P测量是测量投影数据,S康普顿是使用SSS技术获得的康普顿散射估计值,并且CF是用于说明诸如随机的误差的常数或均匀校正因子。定标因子a和b已经使用最佳拟合算法获得。校正数据P校正用于重建指示受检查对象中放射性核素分布的图像。
近来的趋势是PET技术的应用范围,尤其是在医学应用中快速增长。在许多情况下,新的PET示踪剂的发展受到生物化学要求和使用从单光子发射计算机断层摄影(SPECT)或其他方式已知的、具有相对良好特征的化学性质的示踪剂的愿望推动。结果,所谓的脏同位素,例如76Br、86Y、89Zr、94mTc、124I、和66Ga开始引起注意。与常规PET同位素相比,脏同位素交替地通过正电子发射和电子俘获衰变,从而导致从正电子发射事件产生511keV湮灭光子、从电子俘获事件产生各种能量的光子、以及从正电子发射事件产生与湮灭光子级联的各种能量的所谓瞬发光子。例如,碘-124产生511keV湮灭光子,而且产生603keV单光子(两者来自独立的核衰变和来自单核衰变的级联反应)。
不幸的是,脏同位素比常规同位素更难以成像。它们通常更低的正电子丰度通常导致计数的数量减少。而且,单光子会导致假符合,其中单光子不仅彼此,而且与从正电子湮灭产生的511keV光子在时间上符合。这样的假符合产生错误的LOR,由此降低结果图像的质量。
用于减小瞬发光子的影响、因此减少假符合的一种技术是基于被探测光子的能量进行鉴别。然而根据特定同位素的发射光谱和特定扫描器的能量分辨率,可能难以在瞬发光子和511keV光子之间进行区分。
假符合同样已经用级数展开的二阶进行说明,其中零阶和一阶分别提供了对均匀和线性背景贡献的校正。参见Kull等人的“Quantitative Imagingof Yttrium-86PET With the ECAT EXACT HR+in2D Mode”,CancerBiotherapy and Radiopharmaceuticals,第19卷,Number4,2004年。然而,Kull技术假设脏同位素背景分布是二阶或n2函数。同样它并未说明可以影响假符合背景的患者或对象特异性变化。
因此,人们希望提供一种用于减少正电子成像数据中假符合的影响的改进技术。
本发明的各个方面解决了这些和其他问题。
根据本发明的第一方面,一种正电子成像方法包括使用指示对象的衰减分布的信息来生成假符合校正,将假符合校正应用于来自对象的正电子成像检查的数据,以及生成指示校正数据的人可读图像。
根据本发明的另一方面,一种正电子成像装置包括用于获得指示受检查对象中放射性核素衰变的数据的器件。所述衰变包括导致发射时间符合光子对的正电子衰变和导致生成单光子的衰变。所述装置还包括用于测量对象的物理特性的器件,使用测得的物理特性校正被探测到符合的光子对中的单光子和单光子的器件,和用于生成指示校正数据的人可读图像的器件。
根据本发明的另一方面,一种包含指令的计算机可读存储介质,所述指令在由计算机执行时使得计算机执行一种方法,所述方法包括获得指示受检查对象中正电子湮灭的投影数据,获得指示对象的物理特性的对象特异性数据,生成假符合校正,使用假符合校正来校正投影数据,和生成指示校正投影数据的体积图像数据。所述校正是对象特异性数据的函数。
根据另一方面,一种成像方法包括获得指示受检查对象中正电子湮灭的数据,将假符合校正应用于所述数据以减小所述数据中假符合的影响,和生成指示校正数据的人可读图像。所述假符合校正是对象特异性物理特性的函数。
根据另一方面,一种方法包括获得指示受检查对象中正电子湮灭的测量投影数据,使用对象特异性校正函数来校正测量投影数据中的康普顿散射和假符合,和生成指示校正投影数据的人可读图像。
根据另一方面,一种装置包括围绕检查区域布置的多个辐射敏感探测器,用于生成指示由所述辐射敏感探测器探测到的时间符合光子的数据的符合探测器件,用于根据对象特异性校正函数校正所述数据中的散射和假符合以生成校正数据的器件,和用于生成指示校正数据的人可读图像的器件。
根据本发明的另一方面,一种利用正电子湮灭的方法包括获得指示受检查对象中多个正电子湮灭的数据,和将假符合和衰减校正应用于所述数据以生成校正数据。所述假符合校正与下列成比例:
其中,x是离辐射敏感探测器的环的中心的距离。所述方法还包括生成指示校正数据的人可读图像。
根据另一方面,一种正电子成像装置包括围绕检查区域在半径为R的弧中布置的多个辐射敏感探测器,和用于生成指示由所述辐射敏感探测器探测到的时间符合光子的数据的符合探测器件。所述数据包括从脏同位素的衰变产生的假符合。所述装置还包括用于将假符合和衰减校正应用于所述数据以生成校正数据的器件,和用于生成指示校正数据的人可读图像的器件(129,128)。所述假符合校正与下列成比例:
本领域的技术人员在阅读并理解附图和说明书的基础上将领会到本发明的另外其他方面。
本发明通过举例方式进行说明,但是并不局限于各幅附图,在附图中相同的附图标记指示相似的元件,在附图中:
图1描绘了组合式PET/计算机断层摄影(CT)系统;
图2描绘了使用嵌入在圆柱形散射介质中的线源获得的PET正弦图的信号分布;
图3描绘了脏同位素PET过程的二维模型;
图4描绘了分析导出的假符合分布;
图5描绘了脏同位素PET过程的二维模型;
图6描绘了相对于横向平面的不对称假符合源位置;
图7描绘了稍微不对称源的假符合分布;
图8描绘了高度不对称源的假符合分布;
图9描绘了一种用于生成假符合校正的技术;
图10描绘了一种用于校正正电子发射数据的技术。
参照图1,组合式PET/CT系统100包括PET机架部分102和CT机架部分104。PET机架部分102包括由辐射敏感探测器106组成的一个或多个轴向环,所述轴向环围绕检查区域108。探测器106探测在PET检查区域108内发生的正电子湮灭事件的伽玛辐射特性。
CT部分104包括辐射源110,例如围绕CT检查区域112旋转的x射线管。辐射敏感探测器114探测由x射线源发射的辐射,所述辐射横穿检查区域112。
PET机架部分102和CT机架部分104优选地邻近它们各自的检查区域108、112进行定位,所述检查区域沿着公共纵轴或z轴进行布置。对象支架116支撑待成像的对象118,例如人类患者。对象支架116优选地可与PET/CT系统100的操作配合在纵向上移动,使得对象118可以由PET和CT机架部分102、104两者在多个纵向位置处进行扫描。
在静态调强或基于帧的模式中,对象支架116移动到第一纵向位置。在采集完期望数据之后,支架116移动到第二纵向位置,并且获得期望数据。重复该过程直到已经将对象118的期望纵向范围扫描完为止。在连续模式中,对象支架116在扫描期间基本连续地移动。
CT数据采集系统122处理来自CT探测器114的信号以生成指示沿着通过检查区域112的多条线或射线的辐射衰减的数据。CT重建器126使用合适的重建算法重建数据以生成指示对象118的辐射衰减的体积图像数据。
PET数据采集系统120生成投影数据,所述投影数据包括指示时间符合光子的信息以及因此指示由探测器106探测到的正电子湮灭事件的信息。具体而言,投影数据提供每个事件的有关LOR的信息,例如LOR的横向和纵向位置、它的横向和方位角、以及在扫描器具有TOF能力情况下的飞行时间(TOF)信息。所述数据还可以被重新分箱(rebinned)成一个或多个正弦图或投影箱。当用于成像脏同位素时,所述数据可以包括从其衰变产生的假符合。
优选地包括一个或多个计算机处理器的PET重建器129典型地使用迭代重建算法生成指示对象118中放射性核素的分布的体积图像数据。另外,PET重建器129使用来自CT系统的信息将衰减和其他期望校正应用于PET数据。而且,正如下面将要详细描述的,PET重建器129还校正正电子发射数据中的假符合。
使得重建器129执行重建和相关校正的计算机可读指令优选地搭载在一个或多个计算机可读介质上,例如计算机盘片、易失性或非易失性存储器、或(各)处理器可访问的类似存储装置。所述指令还可以通过合适的通信网络,例如因特网被传送到(各)处理器可访问的存储介质。
工作站计算机用作操作者控制台128。控制台128包括诸如监视器或显示器的人可读输出设备,和诸如键盘和鼠标的输入设备。驻留在控制台128上的软件允许操作者观察或以其他方式操纵由PET和CT重建器129、126生成的体积图像数据。驻留在控制台128上的软件还允许操作者通过建立期望扫描方案、启动和终止扫描、和以其他方式与扫描器100相互作用来控制系统100的操作。
同样可以预见到有关系统100的变更。这样,例如,扫描器的CT部分可以被省略、位于远离PET机架部分102的位置、或者用另一成像设备例如磁共振(MR)扫描器代替。也可以使用与PET机架部分102相关联的透射源生成或者从另一源获得衰减或解剖信息。
如上所述,脏同位素的特征在于511keV湮灭光子,以及由瞬发光子生成的附加发射线。图2描绘了使用嵌入在圆柱形散射介质中并且沿着在PET扫描器探测器环106的中心的z轴布置的线源获得的PET正弦图的信号分布。横坐标表示与探测器环106的中心的横向或径向偏移,纵坐标表示计数的数量。
分布202示出了使用含氟-18的体模获得的分布,而分布204示出了使用含碘-124的体模获得的分布。两个分布被归一化为实际计数峰值中的相同整数。第三分布206通过从第二、碘-124分布204减去第一、氟-18分布202获得。第三分布206指示假符合并且通常为碗形,假符合在探测器106环的中心处最低并且通常作为离那里的距离的非线性函数增加。
假符合分布也可以使用脏同位素PET过程的简化二维模型分析导出。参照图3,PET扫描器探测器106被布置在半径为R的圆302中。脏同位素304被建模为位于探测器环的中心的理想点源。在碘-124的情况下,511keV伽玛306和602keV伽玛308被建模为彼此成随机角α发射。它们各自与探测器环302的交点限定附加的瞬间符合或假符合背景LOR310。应当注意模型已经省略了物理效应,例如不完美探测效率、散射和衰减。
光子对306、308可以被建模为在x轴之上和之下成角β=α/2(0<β<π)对称地被发射。得到的LOR310垂直于x轴,因而表示瞬间符合背景在x轴上的投影。由于几何形状的旋转对称,因此成其他角的投影与图3中所示的情况相同。
LOR310的x坐标可以被表示为:
方程2
x(β)=R·cos(β)
因此:
方程3
表示为β和x的函数的LOR310的密度分布g(β)和h(x)的关系如下:
方程4
由于作为β的函数的密度分布g是常数(即,对于每个无穷小的角步长dβ,加入平行于y的另一LOR301),从方程3得出LOR310密度h与导数dβ/dx成比例,x在范围(-R<x<+R)内:
方程5
应当注意:为了便于本分析,比例常数被设置为1以示出瞬间符合背景的形状。
这也是瞬间符合背景LOR的分布所需的最终结果(由于仅对湮灭伽玛和瞬发伽玛中的一个进行探测),原因是任何投影角的线密度由于如上面解释的对称原因与h相等。
图4示出了在R=1的情况下在探测器环302的整个直径上的假符合分布h(x)的绘图。在实际PET扫描器中,跨轴视场仅仅在环直径的一半上延伸,大致对应于图4中的范围(-0.5<x+0.5)。所以在x=-R和x=+R出现的奇点实际上是不重要的。
如应当领会到的,方程5的LOR分布从位于探测器环302的中心的点源304导出。可以将该导数扩展以确定不对称定位源(即,在横向平面中不居中的源)的影响。
现在参见图5的模型,湮灭光子306和附加的瞬发伽玛308可以被建模为彼此成随机角从不对称源位置304a被发射。在数值模拟中,第一光子的发射角α可以在外部程序循环中以步长Δα从0°递增到360°。对于每个α值,第二光子β的方向然后可以在内部循环中以步长Δβ从0°递增到360°。得到的与PET探测器圆302的交点限定出瞬间符合或假LOR310。LOR310根据其离原点的距离r和LOR310与x轴之间的角θ被累积成正弦图。正弦图然后可以在所有角上求和以作为离横向平面的中心的距离的函数显示线密度分布曲线。
图6描绘了在具有半径R=1的探测器环302的扫描器的横向平面中位于坐标(0.1,0)的轻微不对称源304s和位于坐标(0.5,0.5)的高度不对称点源304h。图7和8分别描绘了轻微304s和高度304h不对称点源的假符合分布。
如应当领会到的,图7和8的假符合分布非常类似于方程5的解析表达式。而且,分布的形状对点源304、602、604的位置相对不敏感。结果,在没有衰减和散射的情况下,方程5产生的假符合分布近似于扩展的脏同位素活度分布。
由于活度分布和诸如衰减和散射的物理效应典型地对于不同患者有所不同,甚至在特定患者的连续扫描之间也有所不同,因此人们期望说明一个或多个这些因素的技术产生假符合的更加实际的分布,当用在如随后所述的校正方案中时,所述技术导致假符合对图像质量的影响进一步减小。在图9中示出了用于产生对象特异性假符合校正的技术,该技术解决了这样的影响。
在902,典型地从受检查患者或其他对象的扫描中获得未校正的发射和透射数据及图像。
在904,使用未校正的发射和透射图像估计同位素的活度分布和受检查对象的衰减分布。
在906,将多个假符合对(即,一个511keV湮灭光子和一个瞬发伽玛光子)的源点随机分布在衰减体积中。增加源点的数量通常增加校正的精度,而减少源点的数量通常减少处理时间。
在908,选择LOR。
在910,选择源点。
在912,使用活度分布估计、立体角、衰减分布和同位素的衰变图确定由源点产生的沿着LOR的假事件的数量。
在914,将选定源点的贡献加入选定的LOR。
在916,过程返回步骤910直到选择完所有源点为止。
在918,过程返回步骤908直到已选择完所有LOR为止。
在920,解释LOR空间以获得假背景正弦图。
应当注意可以回溯地确定校正。所述校正方案可以一次地或以迭代方式应用于未校正数据。迭代用法将利用与图9中所示相同的步骤,但是从前一迭代步骤的校正图像数据集开始以估计活度和衰减分布。迭代的数量可以由用户预先给出或根据停止准则进行确定。
假符合校正可以用于减小测量投影数据中假符合的影响。具体而言,可以根据以下关系产生校正投影数据P校正:
方程6
P校正=P测量-a×S康普顿-b×C-c×S假
其中,P校正是校正投影数据,P测量是测量投影数据,S康普顿是优选使用已知的SSS技术执行的康普顿散射校正,CF是用于说明诸如随机的误差的常数或均匀校正因子,而S假是用于假符合的校正。假符合校正S假可以有利地基于方程6的函数或上面结合图9所述的技术获得。定标因子a、b和c使用最佳拟合法获得。
校正数据P校正的产生和指示放射性核素图像数据的体积图像数据的重建将进一步参照图10进行描述。
在1002,将合适的脏同位素引入受检查对象中并且使用PET扫描器获得测量投影数据P测量。
在1004,定标康普顿散射校正S康普顿并且将其应用于测量投影数据P测量。
在1006,定标和应用均匀校正系数CF。
在1008,定标和应用假散射校正S假以产生校正投影数据P校正。
在1010,重建校正投影数据P校正以产生指示脏同位素分布的体积数据。
在1012,例如在控制台128或其他合适的显示器上、屏幕上、或类似物上,以人可读的形式显示体积数据。
在这点上,应当注意到,可以按任何期望的顺序应用校正1004、1006、1008。可以省略一个或多个校正,也可以加入其他的。可以使用患者特异性的或是另外方式的其他校正技术。应当注意到,定标因子a、b和c也可以在单拟合过程中获得。一些或所有校正也可以外部地应用于PET重建器129,例如在PET数据采集系统120中。也可以预见到,也可以应用合适的衰减校正以校正受检查对象中的衰减效应。
当然,他人在阅读和理解说明书的基础上将容易想到修改和变更。本发明旨在解释为当所有这样的修改和变更落在权利要求书或其等价物的范围内时包括所有这样的修改和变更。
Claims (40)
1.一种正电子成像方法,包括:
使用指示对象(118)的衰减分布的信息生成假符合校正;
将所述假符合校正应用于来自所述对象的正电子成像检查的数据;
生成指示所述校正数据的人可读图像。
2.根据权利要求1所述的方法,包括使用指示所述对象中活度分布的信息生成所述假符合校正。
3.根据权利要求2所述的方法,包括使用来自所述正电子成像检查的数据生成指示所述活度分布的所述信息。
4.根据权利要求1所述的方法,包括使用来自所述对象的CT扫描的数据生成指示所述衰减分布的所述信息。
5.根据权利要求1所述的方法,其中,使用指示所述对象的衰减分布的信息生成所述假符合校正的所述步骤包括:
选择LOR;
选择假符合源位置(304);
确定由所述源位置产生的沿着所述LOR的假事件的数量;
针对多条LOR的每条重复确定所述事件数量的所述步骤。
6.根据权利要求5所述的方法,其中,使用指示所述对象的衰减分布的信息生成所述假符合校正的所述步骤包括针对多个源位置的每一个重复所述选择步骤和确定所述事件数量的所述步骤。
7.根据权利要求5所述的方法,包括使用指示所述衰减分布的信息确定所述事件数量。
8.根据权利要求7所述的方法,包括使用指示活度分布的信息确定所述事件数量。
9.根据权利要求1所述的方法,包括将康普顿散射校正和随机校正中的至少一个应用于来自所述正电子成像检查的数据。
10.根据权利要求9所述的方法,包括根据以下函数校正来自所述正电子成像检查的数据:
P校正=P测量-a×S康普顿-b×CF-c×S假
其中,P校正是所述校正数据,P测量是来自所述正电子成像检查的数据,S康普顿是所述康普顿散射校正,CF是均匀校正,S假是所述假符合校正,而a、b和c是定标因子。
11.一种正电子成像装置(100),包括:
用于获得指示受检查对象(118)中放射性核素衰变的数据的器件(102),其中,所述衰变包括导致发射时间符合光子对的正电子衰变和导致产生单光子的衰变;
用于测量所述对象的物理特性的器件(104);
使用所述测得的物理特性校正被探测到时间符合的光子对中的单光子和单光子的器件(129);
用于产生指示所述校正数据的人可读图像的器件(129)。
12.根据权利要求11所述的装置,其中,所述测得的物理特性包括辐射衰减分布。
13.根据权利要求11所述的装置,其中,所述测得的物理特性包括活度分布。
14.一种包含指令的计算机可读存储介质,所述指令在由计算机执行时使得所述计算机执行一种方法,所述方法包括:
获得指示受检查对象(118)中正电子湮灭的投影数据,所述投影数据包括从脏同位素的衰变产生的假符合;
获得指示所述对象的物理特性的对象特异性数据;
生成假符合校正,所述校正是所述对象特异性数据的函数;
使用所述假符合校正来校正所述投影数据;
生成指示所述校正投影数据的体积图像数据。
15.根据权利要求14所述的计算机可读存储介质,其中,所述假符合校正是在多个假符合源位置(304)处生成的假符合的函数。
16.根据权利要求14所述的计算机可读存储介质,其中,所述假符合校正是对象活度分布的函数。
17.根据权利要求14所述的计算机可读存储介质,其中,生成所述假符合校正的所述步骤包括:
确定由假符合源位置产生的沿着LOR的事件的数量;
针对多个源位置重复所述确定步骤。
18.根据权利要求17所述的计算机可读存储介质,其中,所述方法包括确定沿着多条LOR的每条的假符合分布。
19.根据权利要求14所述的计算机可读存储介质,其中,所述物理特性包括衰减分布。
20.根据权利要求14所述的计算机可读存储介质,其中,使用所述假符合校正的所述步骤包括将定标因子应用于所述假符合校正。
21.一种成像方法,包括:
获得指示受检查对象(118)中正电子湮灭的数据,所述数据包括假符合;
将假符合校正应用于所述数据以减小所述假符合的影响,其中,所述假符合校正是对象特异性物理特性的函数;
生成指示所述校正数据的人可读图像。
22.根据权利要求21所述的方法,其中,所述对象特异性物理特性包括空间变化特性。
23.根据权利要求22所述的方法,其中,所述特性包括辐射衰减。
24.根据权利要求22所述的方法,包括:
选择LOR;
作为所述空间变化特性的函数确定由所述源位置产生的沿着LOR的假事件的数量;
针对多条LOR的每条重复所述选择和确定的步骤。
25.根据权利要求21所述的方法,包括使用所述数据估计所述对象特异性物理特性。
26.一种方法,包括:
获得指示受检查对象(118)中正电子湮灭的测量投影数据;
使用对象特异性校正函数来进行所述测量投影数据中的康普顿散射和假符合校正;
生成指示所述校正投影数据的人可读图像。
27.根据权利要求26所述的方法,包括进行所述测量投影数据中的随机校正。
28.根据权利要求26所述的方法,其中,所述对象特异性校正函数采用以下形式:
P校正=P测量-a×S康普顿-b×CF-c×S假
其中,P校正是所述校正投影数据,P测量是所述测量投影数据,S康普顿是康普顿散射校正,CF是空间均匀校正,S假是对象特异性假符合校正,而a、b和c是定标因子。
29.根据权利要求26所述的方法,包括使用指示所述对象的内部特性的数据校正所述假符合。
30.根据权利要求29所述的方法,其中,所述指示特性的数据包括x射线计算机断层摄影数据。
31.根据权利要求26所述的方法,包括重复以下步骤:获得测量投影数据、使用对象特异性校正函数、以及针对多个对象的每一个生成人可读图像。
32.根据权利要求26所述的方法,其中,获得测量投影数据的所述步骤包括获得使用PET扫描器采集的测量投影数据,所述扫描器包括半径为R的辐射敏感探测器的环,和应用与下列成比例的假符合校正:
其中,x是离辐射敏感探测器的所述环的中心的距离。
33.一种装置,包括:
围绕检查区域(108)布置的多个辐射敏感探测器(106);
用于生成指示由所述辐射敏感探测器探测到的时间符合光子的数据的符合探测器件(120);
用于根据对象特异性校正函数校正所述数据中的散射和假符合以生成校正数据的器件(129);
用于生成指示所述校正数据的人可读图像的器件(129,128)。
34.根据权利要求33所述的装置,其中,辐射敏感探测器被布置在半径为R的弧中,并且其中,用于校正的所述器件应用与下列成比例的假符合校正:
其中,x是离辐射敏感探测器的所述弧的中心的距离。
35.根据权利要求33所述的装置,其中,用于校正的所述器件应用假符合校正,所述假符合校正是所述对象的辐射衰减的函数。
36.根据权利要求35所述的装置,其中,所述假符合校正是所述对象中的活度的函数。
37.一种利用使用PET扫描器采集的正电子湮灭数据的方法,所述扫描器包括半径为R的辐射敏感探测器(106)的环,所述方法包括:
获得指示受检查对象(118)中多个正电子湮灭的数据,所述数据包括从脏同位素的衰变产生的假符合;
将假符合和衰减校正应用于所述数据以生成校正数据,其中,所述假符合校正与下列成比例:
其中,x是离辐射敏感探测器的所述环的中心的距离;
生成指示所述校正数据的人可读图像。
38.根据权利要求37所述的方法,其中,根据以下函数产生所述校正数据:
P校正=P测量-a×S康普顿-b×CF-c×S假
其中,P校正是所述校正数据,P测量是所述数据,S康普顿是康普顿散射校正,CF是均匀校正,S假是所述假符合校正,而a、b和c是定标因子。
39.一种正电子成像装置,包括:
围绕检查区域(108)在半径为R的弧中布置的多个辐射敏感探测器(106);
用于生成指示由所述辐射敏感探测器探测到的时间符合光子的数据的符合探测器件(120),所述数据包括从脏同位素的衰变产生的假符合;
用于将假符合和衰减校正应用于所述数据以生成校正数据的器件(129),其中,所述假符合校正与下列成比例:
其中,x是离所述辐射敏感探测器环的中心的距离;
用于生成指示所述校正数据的人可读图像的器件(129,128)。
40.根据权利要求39所述的装置,还包括CT扫描器(104),所述扫描器采集指示所述对象的辐射衰减的信息。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US76741306P | 2006-03-27 | 2006-03-27 | |
US60/767,413 | 2006-03-27 | ||
PCT/US2007/063561 WO2007117801A2 (en) | 2006-03-27 | 2007-03-08 | Dirty isotope pet reconstruction |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101842806A true CN101842806A (zh) | 2010-09-22 |
CN101842806B CN101842806B (zh) | 2014-03-05 |
Family
ID=38537106
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200780010626.9A Expired - Fee Related CN101842806B (zh) | 2006-03-27 | 2007-03-08 | 脏同位素pet重建 |
Country Status (5)
Country | Link |
---|---|
US (1) | US7777189B2 (zh) |
EP (1) | EP2080166A2 (zh) |
JP (1) | JP5290145B2 (zh) |
CN (1) | CN101842806B (zh) |
WO (1) | WO2007117801A2 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104173070A (zh) * | 2013-05-21 | 2014-12-03 | 上海联影医疗科技有限公司 | 确定pet飞行时间的方法和装置 |
CN105894525A (zh) * | 2016-04-29 | 2016-08-24 | 上海联影医疗科技有限公司 | 正电子发射断层成像系统探测器校正方法 |
WO2017032297A1 (en) * | 2015-08-25 | 2017-03-02 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image calibration |
CN110189387A (zh) * | 2019-05-27 | 2019-08-30 | 上海联影医疗科技有限公司 | Pet图像校正方法方法、装置、计算机设备和存储介质 |
WO2023051719A1 (en) * | 2021-09-30 | 2023-04-06 | Shanghai United Imaging Healthcare Co., Ltd. | Methods and systems for attenuation map generation |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008084438A2 (en) | 2007-01-11 | 2008-07-17 | Koninklijke Philips Electronics N.V. | Pet/mr scanners for simultaneous pet and mr imaging |
US7894652B2 (en) * | 2007-04-27 | 2011-02-22 | Siemens Medical Solutions Usa, Inc. | Prompt gamma correction for non-standard isotopes in a PET scanner |
US9069089B2 (en) | 2011-03-30 | 2015-06-30 | The Board Of Trustees Of The Leland Stanford Junior University | Methods and systems for increasing the sensitivity of simultaneous multi-isotope positron emission tomography |
WO2012149552A2 (en) | 2011-04-29 | 2012-11-01 | The Board Of Trustees Of The Leland Stanford Junior University | Dual-isotope positron emitting tomography for disease evaluation |
WO2015006123A1 (en) * | 2013-07-08 | 2015-01-15 | Massachusetts Institute Of Technology | System and method to improve image quality of emission tomography when using advanced radionuclides |
PL227658B1 (pl) * | 2013-08-30 | 2018-01-31 | Uniwersytet Jagiellonski | Tomograf TOF-PET i sposób obrazowania za pomocą tomografu TOF-PET w oparciu o prawdopodobieństwo produkcji i czas życia pozytonium |
PL228483B1 (pl) | 2013-08-30 | 2018-04-30 | Univ Jagiellonski | Tomograf hybrydowy TOF-PET/ MRI |
KR102259856B1 (ko) * | 2019-07-01 | 2021-06-04 | 가천대학교 산학협력단 | 정량적 pet 영상 정밀도 향상을 위한 경동맥 내 방사선 검출 시스템 및 이를 이용한 혈관 내 방사선량 보정 방법 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6008493A (en) | 1997-05-30 | 1999-12-28 | Adac Laboratories | Method and apparatus for performing correction of emission contamination and deadtime loss in a medical imaging system |
US6198104B1 (en) | 1998-10-23 | 2001-03-06 | Adac Laboratories | Randoms correction using artificial trigger pulses in a gamma camera system |
US6490476B1 (en) * | 1999-10-14 | 2002-12-03 | Cti Pet Systems, Inc. | Combined PET and X-ray CT tomograph and method for using same |
AU2001266579A1 (en) * | 2000-05-16 | 2001-11-26 | Dario B. Crosetto | Method and apparatus for anatomical and functional medical imaging |
US7312455B2 (en) * | 2005-01-12 | 2007-12-25 | The General Electric Company | Method and system for scatter correction in a positron emission tomography system |
CN101300601B (zh) * | 2005-11-01 | 2012-03-21 | 皇家飞利浦电子股份有限公司 | 使用部分事件数据进行pet图像重建的方法和系统 |
-
2007
- 2007-03-08 EP EP07758142A patent/EP2080166A2/en not_active Withdrawn
- 2007-03-08 JP JP2009503123A patent/JP5290145B2/ja not_active Expired - Fee Related
- 2007-03-08 CN CN200780010626.9A patent/CN101842806B/zh not_active Expired - Fee Related
- 2007-03-08 US US12/293,089 patent/US7777189B2/en not_active Expired - Fee Related
- 2007-03-08 WO PCT/US2007/063561 patent/WO2007117801A2/en active Application Filing
Non-Patent Citations (2)
Title |
---|
BRADLEY J. BEATTIE. ETC: "Quantitative imaging of bromine-76 and yttrium-86 with PET: A method for the removal of spurious activity introduced by cascade gamma rays", 《MED. PHYS.》 * |
ROBERTO ACCORSI, LARS-ERIC ADAM, MATTHEW E. WERNER: "Implementation of a Single Scatter Simulation Algorithm for 3D PET:Application to Emission and Transmission Scanning", 《IEEE》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104173070A (zh) * | 2013-05-21 | 2014-12-03 | 上海联影医疗科技有限公司 | 确定pet飞行时间的方法和装置 |
CN104173070B (zh) * | 2013-05-21 | 2017-11-28 | 上海联影医疗科技有限公司 | 确定pet飞行时间的方法和装置 |
WO2017032297A1 (en) * | 2015-08-25 | 2017-03-02 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image calibration |
CN105894525A (zh) * | 2016-04-29 | 2016-08-24 | 上海联影医疗科技有限公司 | 正电子发射断层成像系统探测器校正方法 |
CN105894525B (zh) * | 2016-04-29 | 2019-03-05 | 上海联影医疗科技有限公司 | 正电子发射断层成像系统探测器校正方法 |
CN110189387A (zh) * | 2019-05-27 | 2019-08-30 | 上海联影医疗科技有限公司 | Pet图像校正方法方法、装置、计算机设备和存储介质 |
WO2023051719A1 (en) * | 2021-09-30 | 2023-04-06 | Shanghai United Imaging Healthcare Co., Ltd. | Methods and systems for attenuation map generation |
Also Published As
Publication number | Publication date |
---|---|
US7777189B2 (en) | 2010-08-17 |
CN101842806B (zh) | 2014-03-05 |
WO2007117801A3 (en) | 2010-05-27 |
JP2009534631A (ja) | 2009-09-24 |
EP2080166A2 (en) | 2009-07-22 |
JP5290145B2 (ja) | 2013-09-18 |
US20090057561A1 (en) | 2009-03-05 |
WO2007117801A2 (en) | 2007-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101842806B (zh) | 脏同位素pet重建 | |
US9619905B2 (en) | Apparatus and method for generation of attenuation map | |
EP3224801B1 (en) | Multi-modality imaging system and method | |
US9747701B2 (en) | Systems and methods for emission tomography quantitation | |
Zimmerman et al. | Multi-centre evaluation of accuracy and reproducibility of planar and SPECT image quantification: an IAEA phantom study | |
CN103329168A (zh) | 针对spect/ct系统的迭代锥形射束ct重建的截断补偿 | |
CN101305297A (zh) | 使用解剖列表模式掩模的pet成像 | |
US8712124B2 (en) | Artifact removal in nuclear images | |
CN106659452B (zh) | 在定量单光子发射计算机断层扫描中利用多个光电峰的重构 | |
US11642093B2 (en) | Calibration of radiation dose estimation in medical nuclear imaging | |
Buvat et al. | Unified description and validation of Monte Carlo simulators in PET | |
US7129497B2 (en) | Method and system for normalization of a positron emission tomography system | |
D′ Ambrosio et al. | Attenuation correction for small animal PET images: a comparison of two methods | |
Ljungberg | Quantitative SPECT imaging | |
US7564035B2 (en) | Generating detector efficiency estimates for a pet scanner | |
Saha et al. | Performance characteristics of PET scanners | |
Hutton et al. | SPECT and SPECT/CT | |
Turkington | PET imaging basics | |
CN116324515A (zh) | 作为剂量计的微型单光子发射计算机断层摄影 | |
Dahlbom | PET Calibration, Acceptance Testing, and Quality Control | |
Strother et al. | Testing quantitation in PET | |
Crystals et al. | Physics of Positrons, 364 Production of Positron Emission Tomography Radiotracers, 365 Coincidence Detection in Positron Emission Tomography, 366 | |
Buvat | Methodologies for quantitative SPECT | |
Azhari et al. | Positron Emission Tomography (PET) | |
Saha et al. | Data Acquisition and Corrections |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
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: 20140305 Termination date: 20160308 |