CN101223553B - 带有粗疏角和层重排的三维飞行时间pet - Google Patents

带有粗疏角和层重排的三维飞行时间pet Download PDF

Info

Publication number
CN101223553B
CN101223553B CN2006800120670A CN200680012067A CN101223553B CN 101223553 B CN101223553 B CN 101223553B CN 2006800120670 A CN2006800120670 A CN 2006800120670A CN 200680012067 A CN200680012067 A CN 200680012067A CN 101223553 B CN101223553 B CN 101223553B
Authority
CN
China
Prior art keywords
flight time
response line
time response
line data
angle
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.)
Expired - Fee Related
Application number
CN2006800120670A
Other languages
English (en)
Other versions
CN101223553A (zh
Inventor
S·J·C·范登伯格
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
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 CN101223553A publication Critical patent/CN101223553A/zh
Application granted granted Critical
Publication of CN101223553B publication Critical patent/CN101223553B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus 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/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/508—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for non-human patients

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Heart & Thoracic Surgery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • General Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine (AREA)
  • Processing And Handling Of Plastics And Other Materials For Molding In General (AREA)
  • Separation, Recovery Or Treatment Of Waste Materials Containing Plastics (AREA)
  • Refuse Collection And Transfer (AREA)

Abstract

在飞行时间正电子发射断层扫描(TOF-PET)成像方法中，采集三维飞行时间反应线(TOF-LOR)数据。每个TOF-LOR对应于带有飞行时间空间定位的反应线。在飞行时间空间定位的基础上将TOF-LOR数据按层排在多个二维TOF-LOR数据集中。层排的TOF-LOR数据中的至少一些对应于对于二维数据集是斜的反应线。将TOF-LOR数据粗疏角重排到多个粗疏角仓中，每个仓具有至少约10°的角跨距。对粗疏角排列的TOF-LOR数据进行重建，以产生图像层。

Description

带有粗疏角和层重排的三维飞行时间PET
技术领域
本发明涉及辐射检测领域。本发明尤其涉及飞行时间正电子发射断层扫描(time-of-flight positron emission tomography，TOF-PET)成像，并特别参考这种技术对本发明进行描述。不过，本发明大体上更多地涉及利用同时产生辐射粒子发射事件的放射性示踪剂和成像剂的成像、光谱术和类似的技术。
背景技术
在PET中，将放射性药剂给予患者或其它成像对象。放射性药剂产生辐射衰变事件，这些事件发射正电子，这些正电子在电子-正电子湮没事件中与周围的成像对象的电子快速相互作用以产生两条相对指向的伽马射线之前传播非常短的距离。这些伽马射线由成像对象周围的辐射检测器检测为两个基本上同时发生的辐射检测事件，这些事件在它们之间限定反应线(line-of-response，LOR)。
虽然对这些伽马射线“基本上同时”进行检测，但若所涉及的两个辐射检测器中的一个比另一个检测器更接近于电子-正电子湮没事件，那么就会在这两个辐射检测事件之间产生小的时间差。由于伽马射线以光速传播，所以这些检测之间的时间差通常在约几纳秒或更少。在TOF-PET中，辐射检测器以足够高的速度运行，以能够产生这种小的飞行时间差的测量值，然后，这种测量值用于将电子-正电子湮没事件定位于LOR的区段。所采集的LOR数据限定投射数据，这些投射数据可由经过过滤的背投、迭代背投(iterativebackprojection)或其它算法进行重建，以产生重建图像。
重建速度和精度中的一个因素是角排列(angular binning)分辨率。典型的PET扫描仪具有一个或多个检测器环，每个检测器环包括600至700个检测器像素。因此，在检测器环的横向平面中采集的LOR数据的角分辨率(所采集的角排列尺寸)约在0.5°至1°。通过将LOR数据重排(rebinning)到较少的角仓中来提高重建速度，每个角仓跨越较大的角仓尺寸或间隔。不过，通常认为角重排可将图像的分辨率降低到一定的程度，这种程度相当于角重排的LOR数据的降低的角分辨率。
在三维PET中，辐射检测器的多个环围绕成像对象以向多个层提供成像数据。相对指向的伽马射线通常可由在相同或不同环上的辐射检测器进行检测。由相同环的辐射检测器检测的伽马射线对从位于包含该环的平面中的正电子-电子湮没事件发出；不过，由不同环的辐射检测器检测的伽马射线对从位于由检测这些伽马射线的环所界定的体积中的正电子-电子湮没事件发出。在一种重建方式中，仅对与由相同的环所检测的重合伽马射线对应的LOR数据进行处理。这种方式允许对每个层单独进行重建，但导致大量LOR数据的废弃。为了利用包括使用两个不同的检测器环采集的LOR投射的所有LOR数据，以三维方式进行背投或其它重建，这种方式属于计算密集型并要求大量的数据储存。
本发明构思出克服前面所描述的限制和其它缺点的改进设备和方法。
发明内容
根据一个方面，提供一种成像方法。采集三维正电子发射断层扫描飞行时间反应线(TOF-LOR)数据。每个TOF-LOR对应于带有飞行时间空间定位的反应线。在飞行时间空间定位的基础上将TOF-LOR数据按层排在多个二维TOF-LOR数据集中。这些层排的TOF-LOR数据中的至少一些对应于对于二维数据集是斜的反应线。将二维TOF-LOR数据集重建成限定三维重建图像的对应图像层。其中，每个二维飞行时间反应线数据集的所述重建步骤包括：
将所述二维飞行时间反应线数据集的飞行时间反应线数据角排列成多个角仓；以及
将所述角排列的飞行时间反应线数据重建以产生所述图像层，
并且其中所述角排列步骤包括：选择用于所述角仓的至少约ΔΦ＝atan(2g/ξ)的均匀角跨距ΔΦ，式中，g等于所涉及的空间区域的半径除以所述飞行时间反应线的统计宽度，且ξ等于跨越所涉及的空间区域的像素数量。
根据另一个方面，公开一种正电子发射断层(PET)扫描仪，这种扫描仪用于实现发明内容的第一段中所提出的方法。
根据再一个方面，公开一种飞行时间正电子发射断层扫描仪，包括：
用于采集三维正电子发射断层扫描飞行时间反应线(TOF-LOR)数据的采集装置，每个TOF-LOR对应于带有飞行时间空间定位的反应线；
用于在飞行时间空间定位的基础上将TOF-LOR数据按层排在多个二维TOF-LOR数据集中的层排装置，这些层排的TOF-LOR数据中的至少一些对应于反应线，对于二维数据集来讲，这些反应线是斜的；以及
用于将二维TOF-LOR数据集重建成限定三维重建图像的对应图像层的重建装置，其中所述重建装置包括用于将所述二维飞行时间反应线数据集粗疏角地排列成多个粗疏角仓的粗疏角排列装置，
其中，对每个二维飞行时间反应线数据集而言，所述重建装置适于：
通过所述粗疏角排列装置将所述二维飞行时间反应线数据集的飞行时间反应线数据角排列成多个角仓；以及
对所述角排列的飞行时间反应线数据进行重建以产生所述图像层，
并且其中所述粗疏角排列装置适于：
选择用于所述角仓的至少约ΔΦ＝atan(2g/ξ)的均匀角跨距ΔΦ，式中，g等于所涉及的空间区域的半径除以所述飞行时间反应线的统计宽度，且ξ等于跨越所涉及的空间区域的像素数量。
本发明的一个优点在于提供快速TOF-PET图像重建。
另一个优点在于提供高分辨率TOF-PET图像重建。
再一个优点在于使TOF-PET图像重建在分层的基础上进行，而并不废弃由不同的检测器环所采集的LOR数据。
再一个优点在于利用粗疏角重排提供快速TOF-PET图像重建，而并无伴随的图像分辨率的损失。
在阅读了下面的详细描述之后，本领域中熟练的技术人员就会明白多种另外的优点和好处。
附图说明
本发明通过各种器件和器件的布置以及各种程序操作和程序操作的布置成形。附图仅用于示出优选实施例且并不解释为对本发明进行限制。
图1图示出采用粗疏角重排且极可能是层重排(slice rebinning)的TOF-PET系统。
图2图示出在用于飞行时间定位的基础上将三维TOF-LOR数据集排到多个二维TOF-LOR数据集中的层排。
图3图示出用于二维TOF-LOR数据的粗疏角重排的最小仓尺寸的偏差。
图4图示出一种粗疏角重排方式。
图5图示出适于最可能由图1中TOF-PET系统的层重排处理器和二维重建处理器进行的集成三维TOF-PET重建过程的框图。
具体实施方式
参看图1，飞行时间正电子发射断层(TOF-PET)扫描仪8包括多个辐射检测器环10，这些环10布置成观察成像区域12，成像对象安排在成像区域12中。在图1中，成像对象是患者13的头部；不过，也可形成其它解剖区域的图像，而且成像对象可以不是人类。为了便于说明，在图1中示出了五个辐射检测器环10，这五个辐射检测器环分别用r1、r2、r3、r4和r5表示；不过，通常包括20至60个或更多的检测器环，以用于高分辨率三维成像。而且，将会理解，用图形示出了这些辐射检测器环10；这些辐射检测器环通常容纳在扫描仪8的外壳14中，因此从外部看不到。每个辐射检测器环10通常包括数百个辐射检测器，例如，限定每个环600至700个的辐射检测器像素。
在TOF-PET成像的起动之前将适当的放射性药剂给予患者或其它成像对象。放射性药剂包括放射性物质，这种放射性物质经受放射衰变事件，这些事件发射正电子。这些正电子迅速与附近的成像对象的电子湮灭。图1图示出一个示例性正电子-电子湮灭事件16。每个正电子-电子湮灭事件产生具有511keV的能量的两条相对指向的伽马射线。该示例性正电子-电子湮灭事件产生伽马射线17、18。还构思出其它的放射性药剂，这些放射性药剂导致引起多个粒子产生的放射衰变，将沿着反向的或其它的可预见的相对轨道发射这些粒子。伽马射线以光速传播，即～3×108米/秒。由于成像区域12通常具有约为2米或更小的直径或其它性能尺寸，所以从正电子-电子湮灭事件的位置至多个辐射检测器10中的一个检测器的飞行时间约几纳秒或更少。因此，两条相对指向的伽马射线基本上同时撞击这些辐射检测器中的两个。
这两条伽马射线可由在相同的辐射检测器环上的辐射检测器进行检测。在这种情形中，对应的LOR位于辐射检测器环的平面中，且可容易地明白，正电子-电子湮灭事件位于由辐射检测器环所限定的空间层中。
不过，这两条伽马射线更普遍地是由两个不同的辐射检测器环检测。例如，所示出的伽马射线17、18由辐射检测器环10中的两个不同的检测器进行检测。特别地，伽马射线17由辐射检测器环r1检测，而伽马射线18由辐射检测器环r4检测。在这种情形中，对应的LOR并不与辐射检测器环10中的任何一个共面，且正电子-电子湮灭事件16可位于由检测器环r1、r2、r3、r4所限定的四个空间层中的任何一个中。
这些辐射检测器例如利用Anger逻辑或其它处理将每个粒子检测定位，以确定粒子的能量和检测的精确位置。可将这种逻辑集成在辐射检测器电子装置中，或可容纳在PET外壳14中或可位于离开扫描仪的位置。将合成的辐射检测事件储存在事件清单缓冲器22中，且每个检测事件由检测时间加上标签或附上说明，以便于TOF处理。伽马射线对检测电路24对这些辐射检测事件进行处理，以识别属于对应的电子-正电子湮灭事件的基本上同时发生的伽马射线检测对。这种处理可包括如能量加窗(即废弃以约511keV设置的已选能量过滤窗以外的辐射检测事件)和重合检测电路(即废弃以大于已选时间过滤间隔在时间上相互分开的辐射检测事件对)。
在识别伽马射线对时，反应线(LOR)处理器28处理与两个伽马射线检测事件有关的空间信息，以识别连接这两个伽马射线检测的空间反应线(LOR)。由于沿着完全相对的路径发射这两条伽马射线，所以分辨出在LOR上的某个位置已发生了电子-正电子湮灭事件。
在TOF-PET中，检测器环10的辐射检测器具有足够高的瞬时分辨率，以检测两个“基本上同时发生的”伽马射线检测之间的飞行时间差。飞行时间处理器30对这两个伽马射线检测事件的时间之间的差异进行分析，以沿着该LOR将正电子-电子湮灭事件定位。在飞行时间信息的基础上确定飞行时间反应线(TOF-LOR)，如组织投射。在图1中示出了对应于示例性正电子-电子湮灭事件16的组织投射31。由于与检测伽马射线17的辐射检测器相比，示例性正电子-电子湮灭事件16更接近于检测伽马射线18的辐射检测器，所以伽马射线18的飞行时间短于伽马射线17的飞行时间，因此，沿着该LOR朝向检测伽马射线18的辐射检测器相应地移动TOF-LOR或组织投射31。
TOF-LOR31具有沿着该LOR的有限长度，该长度非零，因为这些辐射检测器具有有限的时间分辨率。在一些实施例中，每个TOF-LOR体现或代表表示沿着LOR的正电子-电子湮灭事件的位置的高斯或其它可能性分布。在其它实施例中，该TOF-LOR体现或代表整个统计宽度的不变值，如沿着LOR的可能性分布的半极大处全宽度(FWHM)。TOF-LOR可在TOF信息的基础上采用沿着LOR的正电子-电子湮灭事件的位置的其它适当表示。将累积用于大量的正电子-电子湮灭事件的合成TOF-LOR数据作为三维TOF-LOR数据32集储存。
继续参看图1并进一步参看图2，最可能的层排处理器40在这些TOF-LOR的飞行时间定位的基础上将三维TOF-LOR数据集32排列在多个二维TOF-LOR数据集42中。层排将TOF-LOR数据32适当地排在由辐射检测器环所限定的层中。对于所示出的带有五个检测器环10的实施例来讲，最可能的层排处理器40将三维TOF-LOR数据32排列在分别对应于检测器环r1、r2、r3、r4和r5的五个二维TOF-LOR数据集42中。
层排将每个TOF-LOR适当地排到含有TOF-LOR的至少中心部分的空间层。若每个TOF-LOR代表表示基础正电子-电子湮灭事件的高斯或其它可能性分布，那么就直接确定表示正电子-电子湮灭事件的最有可能的空间位置的可能性分布的高斯峰值或其它中间值或平均值。然后，层排处理器40将每个TOF-LOR适当地排到空间层，该空间层含有对应于该TOF-LOR的正电子-电子湮灭事件的最有可
z ML = ( z 1 + z 2 ) 2 - c Δt 2 ( z 2 - z 1 ) d - - - ( 1 )
能的空间位置。
在图2中，可以看出将示例性TOF-LOR31适当地排到由辐射检测器环r3所限定的层，该层含有大多数或所有的TOF-LOR31。在一些实施例中，层装包括将每个TOF-LOR投射到层中，将该TOF-LOR层装到该层中。这在图2中示出，在图2中，将示例性TOF-LOR31投射到由检测器环r3所限定的层中，以限定位于由检测器环r3所限定的层中的TOF-LOR31′。在所构思出的其它实施例中，将TOF-LOR适当地旋转到该层中，这样就在二维数据集中产生通常较大的TOF-LOR。或者，空间TOF-LOR排到两个或更多的层中，且TOF-LOR对按照在每个层中发生的事件的可能性权重的每个层由作用。
下面是一种用于层装的方式的一个数量示例。在此示例中，每个LOR具有由在时间t1发生的辐射检测事件所确定的第一端点(x1、y1、z1)个由在时间t2发生的辐射检测事件所确定的第一端点(x2、y2、z2)。飞行时间信息由时间差Δt＝t2-ti量化。将最有可能的层(ZML)适当地确定为：
式中，“c”表示光速，d表示LOR的长度，这种长度由下面的等式给出：
d = ( x 2 - x 1 ) 2 + ( y 2 - y 1 ) 2 + ( z 2 - z 1 ) 2 - - - ( 2 )
ZML的第一项是LOR的第一和第二端点之间的平均层位置，而ZML的第二项是说明TOF信息的相关或调节项。
继续参看图2，用θ表示LOR所排到的图像层的各自的LOR的角度。在将高斯核表示用于飞行时间信息时，就可通过用光速(c＝3.0×1011毫米/秒)乘时间上的飞行时间偏差计算飞行时间(FWHMTOF)的空间上的偏差。将由最有可能的层排处理器40所进行的层的误置的可能性的估计计算为：
FWHMslice＝sin(θ)×FWHMTOF    (3)
对于在带有15°的轴向开度角(θ(max)＝7.5°)和FWHMTOF＝300皮秒(90毫米)的系统中具有最斜LOR的TOF-PET来讲：
FWHMslice＝sin(7.5°)×(90毫米)＝12毫米    (4)。
这是最大的层误置。对于不太斜的LOR(即θ＜7.5°)来讲，这种误差会更小。有利的是，层误置误差仅略微取决于PET扫描仪的视野内的TOF-LOR的位置，并且基本上由飞行时间分辨率确定。
二维TOF-LOR数据集42中的每一个由二维重建处理器50重建，以产生对应的图像层。这些重建图像层共同限定存储在图像存储器52中的三维重建图像。可将这种三维重建图像显示在用户界面54上，并且可打印、储存或通过内联网或国际互联网传送或以其它方式使用。在所示出的实施例中，用户界面54还能够使放射线医生或其它用户控制TOF-PET扫描仪8；在其它的实施例中，可设有分立的控制器或控制计算机。
继续参看图1，二维重建处理器50包括粗疏角重排处理器60，粗疏角重排处理器60将二维TOF-LOR数据集重排到多个粗疏角仓中。在一些实施例中，每个角仓具有至少约10°的角跨距。在一些实施例中，角仓的数量少于约二十个仓，这些仓共同跨越180°的角跨距。背投射仪62将粗疏角重排的TOF-LOR数据进行背投，以产生重建图像层。在一些实施例中，背投射仪62是在单一迭代中产生重建图像的过滤背投射仪(filtered backprojector)。在其它实施例中，背投射仪62采用二维TOF-LOR数据的迭代背投来产生图像层。其它的重建算法也可用于重建粗疏角重排的TOF-LOR数据，以产生重建图像层。
参看图3，粗疏重排，例如利用带有至少约10°的角跨距的仓的粗疏重排，并不在实质上降低重建图像的分辨率。对于高斯飞行时间核(Gaussian time-of-flight kernels)来讲，将空间飞行时间分辨率适当地表示为FWHMTOF。考虑半径为R的所涉及的区域，ξ个像素跨越这个区域(即Δx·R＝ξ，其中，Δx是如图3所示的像素尺寸)。在限定g＝R/FWHMTOF且R是所涉及的区域的半径时，若角重排所导致的误差小于像素尺寸Δx，那么分辨率就不会降低。对于用于角仓的ΔΦ的均匀角跨距来讲，这通过选择用于角仓的至少约ΔΦ＝atan(2g/ξ)的均匀角跨距来实现。这推导用于代表对称的高斯可能性分布的TOF-LOR，对称的高斯可能性分布与LOR的角度无关。
图4示出了用于粗疏角重排的一种方式。TOF-LOR 31′在图4中示出，图4还示出了用实线画出的作为高斯轮廓的TOF信息。TOF-LOR31′在层平面中处于角Φ。将经过角重排的TOF-LOR32″重排到用标记n识别的角仓Φn中。用于经过角重排的TOF-LOR31′′的TOF信息示出为用虚线画出的高斯轮廓。示于图4中的坐标S是TOF-LOR31′的反应线与扫描仪的中心之间的径向距离。坐标I确定TOF-LOR31′沿着反应线的位置。坐标S′和I′是经过角重排的TOF-LOR31″的对应径向和TOF坐标。
参看图5，可改变由排列处理器40、60进行的TOF-LOR排列操作的排序。图5示出了由最有可能的层重排处理器40适当进行的集成三维TOF-PET重建过程的框图，而且还示出了二维重建处理器50。将带有接近的取向的TOF-LOR数据加到用于每个角度的不同的原像中。原像的生成用于每个角度，且将利用时间差计算的用于每个“关闭的”TOF-LOR的对应于每个TOF-LOR的正电子-电子湮灭的最有可能的点加到该原像。
参看图1，熟练的技术人员会理解，可将包括粗疏角重排的二维重建处理器50直接应用于二维TOF-LOR数据集。例如，具有单检测器环的TOF-PET扫描仪直接采集二维TOF-LOR数据集。可将这种二维TOF-LOR数据集直接输入到二维重建处理器50以提供利用粗疏角重排的快速二维重建，而并无层排，因为具有单检测器环的PET扫描仪仅从单层采集数据。
已通过参看优选实施例对本发明进行了描述。显然，通过阅读和理解前面的详细描述，其他人可对本发明进行修改和变化。本发明旨在解释为包括所有的这些修改和变化，只要这些修改和变化在所附的权利要求书或等同的描述的范围之内。

Claims (15)

1.一种成像方法，所述成像方法包括以下步骤：
采集三维正电子发射断层扫描飞行时间反应线数据，每个飞行时间反应线对应于带有飞行时间空间定位的反应线；
在所述飞行时间空间定位的基础上将所述飞行时间反应线数据层排在多个二维飞行时间反应线数据集中，所述层排的飞行时间反应线数据中的至少一些对应于这样的反应线，这些反应线对于所述二维数据集而言是斜的；以及
将所述二维飞行时间反应线数据集重建成限定三维重建图像的对应图像层；
其中，每个二维飞行时间反应线数据集的所述重建步骤包括：
将所述二维飞行时间反应线数据集的飞行时间反应线数据角排列成多个角仓；以及
将所述角排列的飞行时间反应线数据重建以产生所述图像层，
并且其中所述角排列步骤包括：
选择用于所述角仓的至少ΔΦ＝atan(2g/ξ)的均匀角跨距(ΔΦ)，式中，g等于所涉及的空间区域的半径除以所述飞行时间反应线的统计宽度，且ξ等于跨越所涉及的空间区域的像素数量。
2.如权利要求1所述的成像方法，其特征在于：所述采集飞行时间反应线数据的步骤采用多个辐射检测器环，且所述层排步骤包括：
将所述飞行时间反应线数据排到对应于由所述辐射检测器环所限定的层的多个二维飞行时间反应线数据集中。
3.如权利要求2所述的成像方法，其特征在于：所述层排步骤进一步包括：
将每个斜的飞行时间反应线投射到其中所述飞行时间反应线被层排的层中。
4.如权利要求1所述的成像方法，其特征在于：所述采集飞行时间反应线数据的步骤采用多个辐射检测器环，每个辐射检测器环限定空间层，且所述层排步骤包括：
将每个斜的飞行时间反应线排列到含有所述飞行时间反应线的至少中心部分的空间层。
5.如权利要求1所述的成像方法，其特征在于：所述采集飞行时间反应线数据的步骤采用多个辐射检测器环，每个辐射检测器环限定空间层，且所述层排步骤包括：
对于每个斜的飞行时间反应线来讲，确定对应于所述飞行时间反应线的正电子-电子湮灭事件的最有可能的空间位置；以及
将每个斜的飞行时间反应线排到含有对应于所述斜的飞行时间反应线的正电子-电子湮灭事件的最有可能的空间位置的空间层。
6.如权利要求1所述的成像方法，其特征在于：每个二维飞行时间反应线数据集的所述重建包括：
将所述飞行时间反应线数据过滤背投以产生所述图像层。
7.如权利要求1所述的成像方法，其特征在于：每个二维飞行时间反应线数据集的所述重建步骤包括：
采用所述飞行时间反应线数据的迭代背投以产生所述图像层。
8.如权利要求1所述的成像方法，其特征在于：每个二维飞行时间反应线数据集的所述重建步骤包括：
将近似图像限定为由所述二维飞行时间反应线数据集的飞行时间反应线限定的空间点的集合；以及
对所述近似的图像进行逆向飞行时间过滤以产生所述图像层。
9.如权利要求8所述的成像方法，其特征在于：将所述近似图像限定为由所述二维飞行时间反应线数据集的飞行时间反应线限定的空间点的集合步骤包括：
将每个空间点限定为对应于所述飞行时间反应线的正电子-电子湮灭事件的最有可能的空间位置。
10.如权利要求1所述的成像方法，其特征在于：所述角仓的数量少于二十个仓，这些仓共同跨越180°的角跨距。
11.如权利要求1所述的成像方法，其特征在于：每个角仓的角跨距至少为10°。
12.如权利要求1所述的成像方法，其特征在于：
将所述二维飞行时间反应线数据集的飞行时间反应线数据角排列成多个角仓的步骤进一步包括：将所述飞行时间反应线数据粗疏角地排列成多个粗疏角仓，每个仓具有至少10°的角跨距；以及
将所述角排列的飞行时间反应线数据重建以产生所述图像层的步骤进一步包括：对所述粗疏角排列的飞行时间反应线数据进行重建，以产生所述图像层。
13.一种用于实现权利要求1所提出的成像方法的正电子发射断层(PET)扫描仪(8)。
14.一种飞行时间正电子发射断层扫描仪(8)，包括：
用于采集三维正电子发射断层扫描飞行时间反应线数据的采集装置(10、22、24、28、30)，每个飞行时间反应线对应于带有飞行时间空间定位的反应线；
用于在所述飞行时间空间定位的基础上将所述飞行时间反应线数据层排在多个二维飞行时间反应线数据集中的层排装置(40)，所述层排的飞行时间反应线数据中的至少一些对应于反应线，对于所述二维数据集来讲，这些反应线是斜的；以及
用于将所述二维飞行时间反应线数据集重建成限定三维重建图像的对应图像层的重建装置(50)，其中所述重建装置(50)包括用于将所述二维飞行时间反应线数据集粗疏角地排列成多个粗疏角仓的粗疏角排列装置(60)，
其中，对每个二维飞行时间反应线数据集而言，所述重建装置(50)适于：
通过所述粗疏角排列装置(60)将所述二维飞行时间反应线数据集的飞行时间反应线数据角排列成多个角仓；以及
对所述角排列的飞行时间反应线数据进行重建以产生所述图像层，
并且其中所述粗疏角排列装置(60)适于：
选择用于所述角仓的至少ΔΦ＝atan(2g/ξ)的均匀角跨距(ΔΦ)，式中，g等于所涉及的空间区域的半径除以所述飞行时间反应线的统计宽度，且ξ等于跨越所涉及的空间区域的像素数量。
15.如权利要求14所述的扫描仪(8)，其特征在于：所述飞行时间反应线数据采集装置包括：
多个检测器环(10)，所述检测器环限定对应于所述二维飞行时间反应线数据集的层，每条斜的反应线由两个不同的环中的检测器来检测。
CN2006800120670A 2005-04-14 2006-03-28 带有粗疏角和层重排的三维飞行时间pet Expired - Fee Related CN101223553B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US67115805P 2005-04-14 2005-04-14
US60/671,158 2005-04-14
PCT/IB2006/050942 WO2006109203A1 (en) 2005-04-14 2006-03-28 Three-dimensional time-of-flight pet with course angular and slice rebinning

Publications (2)

Publication Number Publication Date
CN101223553A CN101223553A (zh) 2008-07-16
CN101223553B true CN101223553B (zh) 2012-11-28

Family

ID=36652473

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2006800120670A Expired - Fee Related CN101223553B (zh) 2005-04-14 2006-03-28 带有粗疏角和层重排的三维飞行时间pet

Country Status (8)

Country Link
US (1) US7897926B2 (zh)
EP (1) EP1875437B1 (zh)
JP (1) JP5378787B2 (zh)
CN (1) CN101223553B (zh)
AT (1) ATE460715T1 (zh)
DE (1) DE602006012820D1 (zh)
RU (1) RU2007142016A (zh)
WO (1) WO2006109203A1 (zh)

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008146186A2 (en) 2007-05-30 2008-12-04 Koninklijke Philips Electronics, N.V. Pet local tomography
CN101809463B (zh) * 2007-09-24 2013-05-22 皇家飞利浦电子股份有限公司 临床前飞行时间pet成像
DE102007052035A1 (de) * 2007-10-30 2009-05-07 Forschungszentrum Jülich GmbH Verfahren zur Positronen-Emissions-Tomographie sowie PET-Scanner
WO2010048325A2 (en) * 2008-10-21 2010-04-29 University Of Southern California Exact and approximate rebinning of time-of-flight pet positron emission tomography data
US8698087B2 (en) 2008-11-03 2014-04-15 The Trustees Of The University Of Pennsylvania Limited angle tomography with time-of-flight PET
WO2010143082A1 (en) * 2009-06-08 2010-12-16 Koninklijke Philips Electronics N.V. Time-of-flight positron emission tomography reconstruction using image content generated event-by-event based on time-of-flight information
US9111381B2 (en) * 2010-01-27 2015-08-18 Koninklijke Philips N.V. Shift-varying line projection using graphics hardware
US20140107476A1 (en) * 2011-06-16 2014-04-17 Koninklijke Philips N.V. Spatial sampling improvement for list-mode pet acquisition using planned table/gantry movement
US8903152B2 (en) * 2012-06-29 2014-12-02 General Electric Company Methods and systems for enhanced tomographic imaging
KR101912715B1 (ko) * 2012-11-20 2018-10-29 삼성전자주식회사 방사선이 방출된 위치의 분포를 추정하는 방법 및 장치
US11356549B2 (en) * 2014-01-07 2022-06-07 Brian Way System and method for discouraging inappropriate use of a mobile device
US9466132B2 (en) * 2014-12-08 2016-10-11 General Electric Company Systems and methods for motion mitigation determinations
US10339676B2 (en) * 2015-09-07 2019-07-02 Uniwersytet Jagiellonski Method for reconstructing multi-tracer metabolic and morphometric images and tomography system for multi-tracer metabolic and morphometric imaging
US10977841B2 (en) * 2016-08-03 2021-04-13 Koninklijke Philips N.V. Time-of-flight (TOF) PET image reconstruction using locally modified TOF kernels
CN109009197B (zh) * 2018-08-06 2020-07-07 南京航空航天大学 一种用于pet检测的双晶体条穿越符合响应线检测系统和方法
CN109350099A (zh) * 2018-09-13 2019-02-19 中山市明峰医疗器械有限公司 一种应用于临床pet系统的随机事件去除处理方法
CN110400360B (zh) * 2019-07-25 2021-03-19 北京航空航天大学 一种基于全卷积神经网络的声波渡越时间检测方法
KR102444195B1 (ko) * 2020-08-24 2022-09-19 주식회사 메디트 3차원 스캐너 디바이스, 동작 방법, 및 그 동작 방법을 실행하기 위한 프로그램이 저장된 컴퓨터 판독 가능 저장 매체

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1203668A (zh) * 1995-10-12 1998-12-30 阿达克实验室 用于双头伽马照相机的分辨率增强
CN1416781A (zh) * 2002-12-30 2003-05-14 北京质子科技开发有限公司 正电子断层扫描中热源高分辨快速图像迭代重建方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4559597A (en) * 1982-07-07 1985-12-17 Clayton Foundation For Research Three-dimensional time-of-flight positron emission camera system
JPS6173083A (ja) * 1984-09-17 1986-04-15 クレイトン　フアウンデイシヨン　フオ−　リサ−チ 時経過陽電子放射カメラによる放射線検出方法及びその実施に使用するカメラ
JPH06347555A (ja) * 1993-06-10 1994-12-22 Hamamatsu Photonics Kk ポジトロンイメージング装置
JP2535762B2 (ja) * 1993-10-14 1996-09-18 東北大学長 陽電子断層撮影装置におけるγ線吸収体による散乱同時計数測定法及び陽電子断層撮影装置
US5744802A (en) 1995-10-25 1998-04-28 Adac Laboratories Image generation from limited projections in positron emission tomography using multi-slice rebinning
DE69613711T2 (de) 1995-11-22 2002-05-08 Picker Internat Inc Pet-echtzeitabbildungsverfahren für einzelphotonen-gammakamera
DE19854438B4 (de) * 1998-11-25 2011-04-14 Siemens Ag Verfahren zur Bildrekonstruktion für einen Computertomographen
US6462342B1 (en) 2000-06-22 2002-10-08 Ge Medical Systems Global Technology Co. Llc Method and system for pet image reconstruction
US6915004B2 (en) * 2002-02-28 2005-07-05 Cti Pet Systems, Inc. Continuous tomography bed motion data processing apparatus and method
US6804325B1 (en) 2002-10-25 2004-10-12 Southeastern Universities Research Assn. Method for position emission mammography image reconstruction
FR2859299B1 (fr) * 2003-08-28 2006-02-17 Ge Med Sys Global Tech Co Llc Procede de reconstruction tomographique par rectification
US7057178B1 (en) * 2004-11-15 2006-06-06 General Electric Company Method and system for imaging using a filter for Time-of-Flight PET
US7417231B2 (en) * 2004-12-22 2008-08-26 Cti Pet Systems, Inc. Fourier re-binning of time-of-flight positron emission tomography data

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1203668A (zh) * 1995-10-12 1998-12-30 阿达克实验室 用于双头伽马照相机的分辨率增强
CN1416781A (zh) * 2002-12-30 2003-05-14 北京质子科技开发有限公司 正电子断层扫描中热源高分辨快速图像迭代重建方法

Also Published As

Publication number Publication date
CN101223553A (zh) 2008-07-16
JP2008537113A (ja) 2008-09-11
RU2007142016A (ru) 2009-05-27
EP1875437A1 (en) 2008-01-09
US7897926B2 (en) 2011-03-01
EP1875437B1 (en) 2010-03-10
US20090124900A1 (en) 2009-05-14
ATE460715T1 (de) 2010-03-15
JP5378787B2 (ja) 2013-12-25
DE602006012820D1 (de) 2010-04-22
WO2006109203A1 (en) 2006-10-19

Similar Documents

Publication Publication Date Title
CN101223553B (zh) 带有粗疏角和层重排的三维飞行时间pet
EP3224801B1 (en) Multi-modality imaging system and method
US7408149B2 (en) Detector head position correction for hybrid SPECT/CT imaging apparatus
CN105473072B (zh) 具有晶体或探测器单元间距的pet系统
US8116426B2 (en) Computed tomography device and method using circular-pixel position-adaptive interpolation
US20070106154A1 (en) Method for reducing an electronic time coincidence window in positron emission tomography
US20120141001A1 (en) System and method for triangular interpolation in image reconstruction for pet
WO2007020543A1 (en) Positron emission tomography time-of-flight list mode reconstruction with detector response function
US6590213B2 (en) Method and system for estimating scatter in a pet scanner
CN106255994A (zh) 针对正电子发射断层摄影(pet)列表模式迭代重建的重建中滤波
JP2005315887A (ja) 陽電子放出断層システムの規格化のための方法及びシステム
EP0976109B1 (en) Direct tomographic reconstruction
WO2014074666A1 (en) Inter-detector scatter enhanced emission tomography
CN103596504A (zh) 使用规划的患者台/扫描架运动对列表模式pet采集进行空间采样改进
US10222490B2 (en) PET scanner with emission and transmission structures in a checkerboard configuration
Casey et al. Siemens biograph vision 600
Groiselle et al. 3D PET list-mode iterative reconstruction using time-of-flight information
Koeppe et al. Instrumentation for positron emission tomography: tomographs and data processing and display systems
US6373059B1 (en) PET scanner septa
Sajedi et al. Development and preliminary results of Xtrim-PET, a modular cost-effective preclinical scanner
US6423971B1 (en) Emission computed tomography through the detection of paired gamma rays
Virador et al. 3-D reconstruction in PET cameras with irregular sampling and depth of interaction
Pellot-Barakat et al. Detection of motion in hybrid PET/SPECT imaging based on the correlation of partial sinograms
Andreyev et al. Spatial resolution of the multiple coincidences Compton camera
Vidal et al. Reconstruction in positron emission tomography

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

Granted publication date: 20121128

Termination date: 20150328

EXPY Termination of patent right or utility model