CN114930384A - Pet图像数据的实况显示 - Google Patents
Pet图像数据的实况显示 Download PDFInfo
- Publication number
- CN114930384A CN114930384A CN202080093115.3A CN202080093115A CN114930384A CN 114930384 A CN114930384 A CN 114930384A CN 202080093115 A CN202080093115 A CN 202080093115A CN 114930384 A CN114930384 A CN 114930384A
- Authority
- CN
- China
- Prior art keywords
- image volume
- cartesian
- data
- aerial image
- frame
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000002600 positron emission tomography Methods 0.000 claims abstract description 91
- 238000000034 method Methods 0.000 claims abstract description 35
- 238000003384 imaging method Methods 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims description 7
- 230000015654 memory Effects 0.000 claims description 5
- 238000012986 modification Methods 0.000 claims description 4
- 230000004048 modification Effects 0.000 claims description 4
- 230000004044 response Effects 0.000 claims description 3
- 230000033001 locomotion Effects 0.000 description 10
- 239000000700 radioactive tracer Substances 0.000 description 9
- 230000000875 corresponding effect Effects 0.000 description 7
- 238000002347 injection Methods 0.000 description 7
- 239000007924 injection Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000013170 computed tomography imaging Methods 0.000 description 5
- 238000012879 PET imaging Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 4
- 230000004807 localization Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- AOYNUTHNTBLRMT-SLPGGIOYSA-N 2-deoxy-2-fluoro-aldehydo-D-glucose Chemical compound OC[C@@H](O)[C@@H](O)[C@H](O)[C@@H](F)C=O AOYNUTHNTBLRMT-SLPGGIOYSA-N 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000013152 interventional procedure Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000031018 biological processes and functions Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- PNDPGZBMCMUPRI-UHFFFAOYSA-N iodine Chemical compound II PNDPGZBMCMUPRI-UHFFFAOYSA-N 0.000 description 1
- 230000002107 myocardial effect Effects 0.000 description 1
- 230000005258 radioactive decay Effects 0.000 description 1
- 229910052701 rubidium Inorganic materials 0.000 description 1
- IGLNJRXAVVLDKE-UHFFFAOYSA-N rubidium atom Chemical compound [Rb] IGLNJRXAVVLDKE-UHFFFAOYSA-N 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 210000003813 thumb Anatomy 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- 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/48—Diagnostic techniques
- A61B6/486—Diagnostic techniques involving generating temporal series of image data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5223—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data generating planar views from image data, e.g. extracting a coronal view from a 3D image
-
- 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
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/10—Geometric effects
- G06T15/20—Perspective computation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/04—Indexing scheme for image data processing or generation, in general involving 3D image data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10104—Positron emission tomography [PET]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Optics & Photonics (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Geometry (AREA)
- Computer Graphics (AREA)
- Computing Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Nuclear Medicine (AREA)
Abstract
一种系统和方法,包括:由成像设备获取的正电子发射断层摄影数据的第一帧到笛卡尔数据的第一帧的定位;基于笛卡尔数据的第一帧的对第一笛卡尔图像体积的生成;第一笛卡尔图像体积的显示;由成像设备获取的正电子发射断层摄影数据的第二帧到笛卡尔数据的第二帧的定位;基于笛卡尔数据的第二帧的对第二笛卡尔图像体积的生成;以及组合的笛卡尔图像体积的显示。
Description
背景技术
根据常规的正电子发射断层摄影(PET)成像,放射性同位素示踪剂被初始地注射到患者身体中。在身体内部,放射性同位素示踪剂发射正电子,该正电子与电子湮灭以产生伽马射线。位于身体外部的检测器系统检测所发射的伽马射线,并且将相关联的湮灭事件记录在正弦图中或者将其记录为列表模式数据。然后,可以基于该正弦图/列表模式数据来重建三维笛卡尔图像。
重建是耗时和耗资源的。所重建的图像仅可用于在PET扫描完成之后良好地查看。这种延迟限制了结合介入过程来使用PET成像的能力。此外,当前的技术无法提供时间示踪剂分布的合适指示。系统被期望以改进基于PET数据生成的图像的及时性和有用性。
附图说明
图1图示了根据一些实施例的从PET数据生成三维图像体积的过程;
图2图示了根据一些实施例的事件从检测器空间到笛卡尔坐标的定位;
图3是根据一些实施例的PET/CT成像系统的框图;
图4包括根据一些实施例的从PET数据生成三维图像体积的过程的流程图;
图5包括根据一些实施例的从在连续床运动期间获取的PET数据来生成三维图像体积的过程的流程图;
图6图示了根据一些实施例的PET扫描的执行;
图7A至7H包括根据一些实施例的从PET数据生成的三维图像体积的二维描绘;
图8包括根据一些实施例的生成从PET数据生成的过去和当前图像体积的加权组合的过程的流程图;
图9A至9D包括根据一些实施例的从由PET数据生成的过去和当前图像体积所生成的三维图像体积的二维描绘;以及
图10A至10D包括根据一些实施例的从由PET数据生成的过去和当前图像体积所生成的三维图像体积的二维描绘。
具体实施方式
以下描述被提供以使得本领域任何人员能够制作和使用所描述的实施例,并且阐述了预期用于实施所描述的实施例的最佳模式。然而,各种修改对本领域技术人员来说将仍然是明显的。
通常,一些实施例提供了PET数据到笛卡尔空间位置的直接定位,以便生成笛卡尔图像数据以供显示。附加PET数据被获取,并且也被投影到笛卡尔空间位置以便生成附加笛卡尔图像数据。该附加笛卡尔图像数据与已经显示的笛卡尔图像数据组合,并且图像基于该组合的图像数据来显示。
因此,一些实施例可以在PET扫描期间提供接近实时的图像。这种图像可以允许基于同时获取的PET数据的实况(live)介入过程。
根据一些实施例,该附加PET数据是从与从其中获取了先前PET数据的区域不同的身体区域获取的。该附加笛卡尔图像数据与已经显示的笛卡尔图像数据的组合因此可以允许显示整个移动扫描的视野,其中随着扫描的进行,该视野的区域被逐渐地“填充”。
在一些实施例中,笛卡尔图像数据集合的组合可以是时间相关的。例如,在当前显示的图像数据与新生成的笛卡尔图像数据的上述组合中,当前显示的图像数据可以与新生成的笛卡尔图像数据相比被更少地加权。随着新PET数据帧被获取并且从其中生成的笛卡尔图像数据被与当前显示的图像组合,该加权用于降低每个随后显示的组合帧中的较早获取的PET数据的显著性(prominence)。
笛卡尔图像数据集合的组合也可以或者可以替代地是空间相关的。例如,该组合可以对某些区域(例如,远在感兴趣区域外部的区域)的图像数据进行比其他区域的图像数据更少的加权。
图1图示了根据一些实施例的三维图像体积从PET数据的生成。沿着时间轴t前进,PET数据102至132表示PET数据的接续获取的帧。
一般而言,并且如本领域所已知,被注射到身体化合物中的PET示踪剂的放射性衰变生成正电子,该正电子最终遇到电子并且由此被湮灭。湮灭产生了在大致相反的方向上行进的两个伽马光子。当被设置在身体的相对侧的两个检测器在特定重合(coincidence)时间窗口内检测到两个相反地行进的伽马光子的到达时,湮灭事件被标识。
因为两个伽马光子在大致相反的方向上行进,所以两个检测器的位置确定了沿着其发生该湮灭事件的响应线(LOR)。飞行时间(TOF)PET附加地测量从每个湮灭事件产生的两个伽马光子的到达时间之间的差异。现代PET扫描器能够以足够的准确度来测量每个光子的到达时间中的差异,以提供沿着LOR的其中发生了正电子湮灭的空间位置的指示。
PET数据102至132的每一帧可以包括描述了在与该帧相关联的特定时间段内检测到的每个重合事件或LOR的列表模式数据、连同每个事件的TOF信息。常规地,当事件被检测到时,每个事件的LOR和TOF被置于列表模式流中,并且每当所定义的帧持续时间过去时,新的数据“帧”被建立。
在一些实施例中,PET数据102至132被存储在称为正弦图的数据阵列中。正弦图指示每个LOR的角度对比位移。每个正弦图存储了每个重合事件的LOR的位置,以使得穿过该体积中的单个点的所有LOR在正弦图中描绘出正弦曲线。每个正弦图包括包含针对特定方位角的LOR的一个行。沿着同一LOR检测到的多个事件被直方图化(histogrammed),以保留接续的事件。
如图1中所描绘,帧102至132中的每一个经受笛卡尔定位以生成笛卡尔数据106至136的对应帧。从列表模式TOF PET数据到笛卡尔数据的定位可以根据任何已知或变得已知的过程来进行。图2图示了根据一些实施例的事件从检测器空间到笛卡尔坐标的定位。
图2图示了由伽马光子检测器组成的检测器环210。将假定的是,湮灭事件沿着LOR220发生,从而导致生成两个相反地行进的伽马光子,这两个伽马光子被位于检测器坐标x a ,y a 和x b ,y b 处的检测器检测到。由于伽马光子在检测器处的相应到达时间,可以确定的是,该湮灭事件在TOF窗口230内沿着LOR 220发生。根据一些实施例,假定该湮灭事件在沿着LOR 220定位的体素处并且在TOF窗口230的中心处发生。
TOF列表模式帧到笛卡尔帧的定位可以包括将每个湮灭事件定位到笛卡尔坐标。例如,湮灭事件的LOR可以被认为是在空间上以和定位的一对晶体、以及描述了该湮灭事件与LOR中心的空间偏移量的TOF值。然后,所估计的湮灭事件的空间位置可以被计算为:
其中L 3D 是两个晶体a与b之间的距离,
这可以针对帧的每个LOR被连续地重复,以在LOR被处理时累积活动分布的体积表示。可以通过由针对给定LOR的CT导出的校正因子将每个事件对笛卡尔数据的贡献进行加权来校正每个事件以用于衰减。
上述过程针对PET数据的每个随后获取的帧而发生。然而,如图1中所示,所显示的笛卡尔数据118是所显示的笛卡尔数据108和基于PET数据112生成的笛卡尔数据116的组合。类似地,所显示的笛卡尔数据128是所显示的笛卡尔数据128和基于PET数据122生成的笛卡尔数据126的组合。如上所描述,用于组合笛卡尔坐标系(Cartesian)的技术可以在实施例当中变化,从而产生各种结果。
根据一些实施例,笛卡尔数据106至136的每一帧可以被认为是在帧开始时间处被初始化为0的直方图。当新的检测器对事件在PET数据帧中被收集时,每个事件被定位到在TOF窗口的中心处沿着LOR的笛卡尔坐标。该笛卡尔坐标在笛卡尔帧(即,直方图)中被递增。在帧结束时间处,该笛卡尔帧与当前显示的笛卡尔体积组合,并且被显示。
组合的笛卡尔数据可以被显示为可导航的定量笛卡尔体积(即,可手动选择的横轴(transaxial)、矢状和冠状切片)、被显示为矢状和冠状最大强度投影(MIP)视图、或者被显示为单个旋转MIP。
图3是执行本文中描述的一个或多个过程的PET/CT成像系统300的框图。实施例不限于系统300。
系统300包括定义了孔洞312的台架310。如本领域已知的,台架310容纳用于获取PET图像数据的PET成像组件和用于获取CT图像数据的CT成像组件。CT成像组件可以包括一个或多个x射线管和一个或多个对应的x射线检测器,如本领域已知的那样。
PET成像组件可以包括采用任何配置的任何数量或类型的检测器(例如,硅光电倍增器(SiPM)),如本领域已知的那样。检测器与切片厚度(空间分辨率)相关联,使得组件能够独立地对以大于或等于切片厚度的距离而分隔的两个切片进行成像。切片厚度(例如,2.0mm)对应于检测器的分辨率。
注射系统318可以操作以在PET扫描之前和/或期间向患者递送铷、氟脱氧葡萄糖(FDG)、碘或其他放射性同位素的经校准的注射。在一些实施例中,注射系统318被并入台架310中。注射系统318可以支持与控制系统320的有线或无线通信链路,以用于接收指定了剂量、注射方案和扫描延迟的信息。
床315和基座316可操作以在成像之前、期间和之后将躺在床315上的患者移动进入孔洞312和离开孔洞312。在一些实施例中,床315被配置成在基座316上平移,并且在其他实施例中,基座316可与床315一起移动、或者替代地从床315移动。
患者进入和离开孔洞312的移动可以允许使用台架310的CT成像元件和PET成像元件来扫描患者。这种扫描可以基于扫描参数(诸如,扫描范围和对应的扫描速度)来进行。根据一些实施例,床315和基座316可以在这种扫描期间提供连续床运动和/或步进式(step-and-shoot)运动。
控制系统320可以包括任何通用或专用计算系统。因此,控制系统320包括被配置成执行处理器可执行程序代码以使系统320如本文中描述的那样操作的一个或多个处理单元322、以及用于存储该程序代码的存储设备330。存储设备330可以包括安装在对应接口(例如,USB端口)中的一个或多个固定盘、固态随机存取存储器和/或可移除介质(例如,拇指驱动器)。
存储设备330存储控制程序331的程序代码。一个或多个处理单元322可以执行控制程序331以结合PET系统接口323、床接口325和注射接口327来控制硬件元件,以将放射性同位素注射到患者中,将患者移动到孔洞312中从而经过台架310的PET检测器,并且检测发生在患者内的湮灭事件。检测到的事件可以作为PET数据333被存储在存储器330中,该PET数据333可以包括列表模式数据和/或正弦图。
一个或多个处理单元322还可以执行控制程序331以结合CT系统接口324使台架310内的辐射源从不同的投影角度朝向孔洞312内的身体发射辐射,并且控制对应的检测器以获取二维CT数据。CT数据可以与如上所描述的PET数据基本同时地获取,并且可以被存储为CT数据334。这种CT数据334可以用于同时获取的PET数据333的衰减校正,如本领域已知的那样。在这点上,控制程序331还可以被执行以使用任何已知或变得已知的重建算法将PET扫描的PET数据333重建成三维切片。
一个或多个处理单元322可以执行笛卡尔定位程序332,以帮助基于PET数据333的相应帧来生成笛卡尔图像体积335。笛卡尔图像体积335可以包括如本文中描述的组合的图像体积。
PET图像、CT图像和/或图像体积335可以经由终端接口326被传输到终端340。终端340可以包括耦合到系统320的显示设备和输入设备。终端340可以显示PET图像、CT图像和/或图像体积335。终端340可以接收用于控制数据的显示、系统300的操作和/或本文中描述的处理的用户输入。在一些实施例中,终端340是分离的计算设备,诸如但不限于台式计算机、膝上型计算机、平板计算机和智能电话。
系统300的每个组件可以包括对于其操作所必要的其他元件、以及用于提供除了本文中描述的那些功能之外的功能的附加元件。本文中描述的每个功能组件可以在计算机硬件、程序代码和/或执行这种程序代码的一个或多个计算系统中实现,如本领域已知的那样。这种计算系统可以包括执行被存储在存储器系统中的处理器可执行程序代码的一个或多个处理单元。
图4包括根据一些实施例的从PET数据生成三维图像体积的过程400的流程图。流程图400和本文中描述的其他过程可以使用硬件和软件的任何合适的组合来执行。体现了这些过程的软件程序代码可以由任何非暂时性有形介质来存储,该介质包括固定盘、易失性或非易失性随机存取存储器、DVD、闪存驱动器和磁带。实施例不限于下面描述的示例。
在S405处,初始地获取PET数据的第一帧。在放射性同位素示踪剂到主体体积(例如,患者)中的注射之后,可以通过常规的静态PET扫描来获取PET数据,如本领域已知的那样。根据一些实施例,由PET扫描器获取的数据是如上所描述的列表模式数据。
接下来,在S410处,将PET数据定位到笛卡尔坐标空间中以生成笛卡尔数据的帧。如上所描述,特定事件到三维笛卡尔空间的定位基于吸收了该事件的重合伽马光子的检测器的位置、以及该事件的TOF数据。笛卡尔数据的帧可以包括直方图,使得与特定笛卡尔坐标相关联的值可以与定位到该坐标的PET事件的数量直接相关。
在S415处,基于笛卡尔数据的帧来显示笛卡尔图像数据。如上所提到的,该图像数据可以使用用于显示三维图像数据的任何已知或变得已知的技术来显示。在S420处,获取PET数据的下一帧。如下面将更详细描述的,PET数据的下一帧可以表示与先前PET帧的成像区域不同的成像区域(即,身体的不同区域)。在S425处,将PET数据的下一帧定位到笛卡尔坐标空间中以生成笛卡尔数据的下一帧,如上面关于S410所描述的那样。
接下来,在S430处,基于在S415处显示的笛卡尔图像数据和在S425处生成的笛卡尔数据来生成组合的笛卡尔图像数据。S430处的组合可以简单地包括基于在S425处生成的笛卡尔数据来生成图像数据,并且将所生成的图像数据添加到在S415处显示的笛卡尔图像数据。在一些实施例中,可以对图像数据的组合进行加权,其中与特定图像数据值相关联的权重取决于获取时间、笛卡尔坐标和/或示踪剂衰变速率。
在S435处,显示组合的笛卡尔图像数据。然后,在S440处确定扫描是否完成。如果否,则流程返回到S420并且如上所描述的那样继续。特别地,在S420处获取PET数据的下一帧,并且在S425处将其定位到笛卡尔坐标空间中以生成笛卡尔数据的下一帧。
在S430处,基于最近显示的笛卡尔图像数据和笛卡尔数据的最近生成的帧来生成组合的笛卡尔图像数据。在S435处,显示笛卡尔图像数据的新组合。流程继续到S440,并且通过S420、S425、S430、S435和S440循环,直到然后确定扫描完成。
根据过程400的一些实施例,每个随后显示的图像包括每个先前显示的图像的图像数据。如下所描述,在一些实施例中,随着越来越多的图像数据被生成和显示,先前显示的图像的图像数据可能会逐渐消失(fade)或者最终被完全移除。
图5包括根据一些实施例的从在连续床运动期间获取的PET数据来生成三维图像体积的过程的流程图。
在S505处,初始地从视野的第一部分获取PET数据的第一帧。PET数据可以在连续床运动(例如,1mm/s)期间在S505处被获取,如本领域已知的那样,在其中该视野是在扫描期间经过检测器的身体的整个区域。在S510处,将PET数据定位到笛卡尔坐标空间中以生成笛卡尔数据的帧,并且在S515处,基于笛卡尔数据的帧来显示笛卡尔图像数据。
在S420处,从该视野的第二部分获取PET数据的下一帧。由于连续床运动,因此第二部分可能不同于从其中获取了第一PET数据的第一部分。在一些实施例中,第一部分和第二部分可以在空间上重叠。
图6图示了根据一些实施例的在连续床运动期间的PET数据的获取。箭头表示在PET数据的获取期间的时间流逝(passage of time),该获取是关于在该时间流逝期间从其中获取PET数据的身体600的部分的。在所图示的示例中,当身体600的下部部分被移动从而经过PET检测器(凭借使床315移动穿过孔洞312)时,PET数据的获取开始,并且该获取继续直到头部已经被移动从而经过检测器为止。过程500中使用的PET数据的每一帧可以包括在该运动期间的指定时间段内获取的数据。
在S525处,将表示该视野的第二部分的PET数据的下一帧定位到笛卡尔坐标空间中以生成笛卡尔数据的下一帧。接下来,在S530处,基于在S515处显示的笛卡尔图像数据和在S525处生成的笛卡尔数据的帧来生成组合的笛卡尔图像数据。由于图像数据表示整个轴向视野的不同部分,因此S530处的组合可以简单地包括基于在S525处生成的笛卡尔数据来生成图像数据,并且将所生成的图像数据添加到在S515处显示的笛卡尔图像数据。
在S535处,显示组合的笛卡尔图像数据,并且流程返回到S520以如上所描述的那样继续,直到扫描完成。在这种循环期间,表示该视野的接续部分的笛卡尔图像数据继续被添加到所显示的图像数据,从而“完成”该视野的图像。根据一些实施例,在将图像数据与当时当前显示的图像数据进行组合之前,同位素相关的衰变校正函数被在运行中(on-the-fly)应用于新生成的图像数据。同位素相关的衰变校正函数基于像素值的相关联的获取时间来增加像素值的强度,以补偿同位素衰变。
图7A至7H包括根据一些实施例的基于过程500生成的三维图像的二维描绘。尽管仅示出了八个图像,但是应当理解,全身PET扫描可以生成数百个PET数据帧,这些PET数据帧可以被定位到笛卡尔数据的对应帧。
图7A可以描绘基于PET数据的第一帧所生成的第一笛卡尔图像数据。图7B进而可以描绘基于PET数据的另一个帧所生成并且被添加到图7A的第一笛卡尔图像数据的笛卡尔图像。图7B的附加图像数据可以占据(inhabit)与图7A图像数据相同的笛卡尔区域中的一些。在这点上,在过程500期间,PET数据到笛卡尔数据的定位计及了沿着z轴的床移动。例如,在其处获取PET数据的全身视野中的z位置被用于确定从其中定位的笛卡尔数据的z坐标。
如所示出,随着附加PET数据帧被获取并且被定位到笛卡尔空间中,所显示的图像在扫描的z方向上“增长”。如上所描述,对应于稍后获取的帧的图像像素的强度可以基于预期示踪剂衰变而增加,以便归一化所显示的图像内的强度。还如上所描述,连续改变的所显示的图像可以包括例如可手动选择的横轴、矢状和冠状切片、矢状和冠状MIP视图、或单个旋转MIP。
图8包括根据一些实施例的生成从PET数据生成的过去和当前图像体积的加权组合的过程800的流程图。步骤S805至S825可以与过程400的步骤S405至S425类似地进行。因此,在S825处,基于所获取的PET数据帧而生成的笛卡尔图像数据当前被显示,并且笛卡尔数据的下一帧已经基于下一个所获取的PET数据帧而被定位。
接下来,在S830处,基于当前显示的图像数据和笛卡尔数据的新生成的(即,“下一个”)帧来生成组合的笛卡尔图像数据。如下面详细描述的,S830处的生成将不同的权重应用于所显示的图像数据和从笛卡尔数据的下一帧新生成的笛卡尔图像数据。
在根据一些实施例的S830的以下描述中,组合图像被表示为g t ,当前显示的图像被表示为g t-1,并且新生成的笛卡尔图像数据被表示为f t 。在S830处生成的组合图像g t 是新生成的笛卡尔图像数据f t 与当前显示的图像g t-1的一小部分(fraction)的总和,使得g t =αf t +βg t-1。α(x, y, z, t)是控制图像数据f t 对全局体积的贡献程度的加权因子。β(x, y, z,t)是确定事件从接续的组合图像g t 中被多么迅速地移除的持续性参数。
例如,如果β=1,则所有先前的图像数据被累积在组合图像中。如果β=0,则组合图像仅包括最新生成的笛卡尔图像数据f t 。α和β可以保持恒定,或者可以是在空间上和/或在时间上变化的,以计及已知的示踪剂分布、示踪剂衰变、特定临床方案(例如,心肌血流研究)和良好表征的生物过程。最优参数可以通过分析来计算或者经由深度学习方法来计算。
在S835处,显示组合图像,并且然后流程如上所描述的那样从S840循环到S820,直到扫描完成。
图9A至9D包括根据一些实施例的在S830处生成并且在S835处显示的三维图像体积的二维描绘。图9A至9D图示了其中β=1的实现方式。因此,每个接续图像包括先前组合图像的所有图像数据、以及新生成的笛卡尔图像数据。在这点上,图9A至9D的实现方式类似于图7A至7H中描绘的实现方式。
图10A至10D包括根据一些实施例的在S830处生成并且在S835处显示的三维图像体积的二维描绘。特别地,图10A至10D图示了其中β=0.9的实现方式。因此,每个接续图像包括先前组合图像的所有90%的图像数据、以及新生成的笛卡尔图像数据。由于β到接续图像的应用,特定图像数据在组合图像中出现的程度取决于自生成了该特定图像数据以来已经生成的组合图像的数量。例如,对于α=.1和β=0.9,第三组合图像可以被生成为0.9((0.9 xImage1)+Image2) + Image3 = 0.81Image1 + 0.9Image2 +Image3。
前述示图表示根据一些实施例的逻辑架构和过程,并且实际实现方式可以包括以其他方式布置的更多或不同的组件和/或步骤。此外,本文中描述的每个组件或设备可以通过经由任何数量的其他公共和/或专用网络通信的任何数量的设备来实现。这种计算设备中的两个或更多个可以被定位成彼此远离,并且可以经由(一个或多个)网络和/或专用连接的任何已知方式来彼此通信。每个组件或设备可以包括适合于提供本文中描述的功能以及任何其他功能的任何数量的硬件和/或软件元件。例如,在根据一些实施例的系统的实现方式中使用的任何计算设备可以包括处理器以执行程序代码,使得计算设备如本文中描述的那样进行操作。
本文中讨论的所有系统和过程可以体现在被存储在一个或多个非暂时性计算机可读介质上的程序代码中。这种介质可以包括例如硬盘、DVD-ROM、闪存驱动器、磁带和固态随机存取存储器(RAM)或只读存储器(ROM)存储单元。实施例不限于硬件和软件的任何特定组合。
本领域技术人员将领会,在不脱离权利要求的情况下,可以配置对上面描述的实施例的各种适配和修改。因此,要理解的是,权利要求可以以除了本文中具体描述的方式之外的方式来实践。
Claims (20)
1.一种系统,包括:
成像设备,用于:
获取正电子发射断层摄影数据的帧;
处理系统,用于:
将由所述成像设备获取的正电子发射断层摄影数据的第一帧定位到空间数据的第一帧;
基于空间数据的第一帧来生成第一空间图像体积;
显示第一空间图像体积;
将由所述成像设备获取的正电子发射断层摄影数据的第二帧定位到空间数据的第二帧;
基于空间数据的第二帧来生成第二空间图像体积;
基于第一空间图像体积和第二空间图像体积来生成组合的空间图像体积;以及
显示所述组合的空间图像体积。
2.根据权利要求1所述的系统,其中第一空间图像体积包括视野的第一部分,并且其中第二空间图像体积包括所述视野的第二部分。
3.根据权利要求2所述的系统,其中所述组合的空间图像体积的生成包括:基于放射性同位素的衰变概况来修改第二空间图像体积、以及经修改的第二空间图像体积到第一空间图像体积的添加。
4.根据权利要求1所述的系统,其中所述组合的空间图像体积的生成包括:权重到第一空间图像体积的应用、以及经加权的第一空间图像体积到第二空间图像体积的添加。
5.根据权利要求4所述的系统,其中所述权重是在时间上变化的。
6.根据权利要求5所述的系统,其中所述权重是在空间上变化的。
7.根据权利要求6所述的系统,其中所述组合的空间图像体积的生成包括:基于放射性同位素的衰变概况的对第二空间图像体积的修改、以及经加权的第一空间图像体积到经修改的第二空间图像体积的添加。
8.一种方法,包括:
将由成像设备获取的正电子发射断层摄影数据的第一帧定位到笛卡尔数据的第一帧;
基于笛卡尔数据的第一帧来生成第一笛卡尔图像体积;
显示第一笛卡尔图像体积;
将由所述成像设备获取的正电子发射断层摄影数据的第二帧投影到笛卡尔数据的第二帧;
基于笛卡尔数据的第二帧来生成第二笛卡尔图像体积;
基于第一笛卡尔图像体积和第二笛卡尔图像体积来生成组合的笛卡尔图像体积;以及
显示所述组合的笛卡尔图像体积。
9.根据权利要求8所述的方法,其中第一笛卡尔图像体积包括视野的第一部分,并且其中第二笛卡尔图像体积包括所述视野的第二部分。
10.根据权利要求9所述的方法,其中生成所述组合的笛卡尔图像体积包括:基于放射性同位素的衰变概况来修改第二笛卡尔图像体积、以及将经修改的第二笛卡尔图像体积添加到第一笛卡尔图像体积。
11.根据权利要求8所述的方法,其中生成所述组合的笛卡尔图像体积包括:将权重应用于第一笛卡尔图像体积、以及将经加权的第一笛卡尔图像体积添加到第二笛卡尔图像体积。
12.根据权利要求11所述的方法,其中所述权重是在时间上变化的。
13.根据权利要求12所述的方法,其中所述权重是在空间上变化的。
14.根据权利要求13所述的方法,其中生成所述组合的笛卡尔图像体积包括:基于放射性同位素的衰变概况来修改第二笛卡尔图像体积、以及将经加权的第一笛卡尔图像体积添加到经修改的第二笛卡尔图像体积。
15.一种计算系统,包括:
存储器,其存储处理器可执行的过程步骤;
处理单元,用于执行所述处理器可执行的过程步骤,以便:
确定与正电子发射断层摄影数据的第一帧的每个响应线相关联的第一空间图像体积的坐标;
显示第一空间图像体积;
确定与正电子发射断层摄影数据的第二帧的每个响应线相关联的第二空间图像体积的坐标;
基于第一空间图像体积和第二空间图像体积来生成组合的空间图像体积;以及
显示所述组合的空间图像体积。
16.根据权利要求15所述的系统,其中第一空间图像体积包括视野的第一部分,并且其中第二空间图像体积包括所述视野的第二部分。
17.根据权利要求16所述的系统,其中所述组合的空间图像体积的生成包括:基于放射性同位素的衰变概况来修改第二空间图像体积、以及经修改的第二空间图像体积到第一空间图像体积的添加。
18.根据权利要求16所述的系统,其中所述组合的空间图像体积的生成包括:权重到第一空间图像体积的应用、以及经加权的第一空间图像体积到第二空间图像体积的添加。
19.根据权利要求18所述的系统,其中所述权重是在时间上变化的和在空间上变化的。
20.根据权利要求19所述的系统,其中所述组合的空间图像体积的生成包括:基于放射性同位素的衰变概况的对第二空间图像体积的修改、以及经加权的第一空间图像体积到经修改的第二空间图像体积的添加。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/US2020/013423 WO2021145856A1 (en) | 2020-01-14 | 2020-01-14 | Live display of pet image data |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114930384A true CN114930384A (zh) | 2022-08-19 |
Family
ID=69529026
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202080093115.3A Pending CN114930384A (zh) | 2020-01-14 | 2020-01-14 | Pet图像数据的实况显示 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20220265237A1 (zh) |
EP (1) | EP4073743A1 (zh) |
CN (1) | CN114930384A (zh) |
WO (1) | WO2021145856A1 (zh) |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006111869A2 (en) * | 2005-04-22 | 2006-10-26 | Koninklijke Philips Electronics N.V. | Pet/mr scanner with time-of-flight capability |
US7983465B2 (en) * | 2007-05-09 | 2011-07-19 | Société De Commercialisation Des Produits De La Recherche Appliquée - Socpra Sciences Santé Et Humaines, S.E.C. | Image reconstruction methods based on block circulant system matrices |
-
2020
- 2020-01-14 US US17/634,577 patent/US20220265237A1/en active Pending
- 2020-01-14 CN CN202080093115.3A patent/CN114930384A/zh active Pending
- 2020-01-14 WO PCT/US2020/013423 patent/WO2021145856A1/en unknown
- 2020-01-14 EP EP20704706.9A patent/EP4073743A1/en active Pending
Also Published As
Publication number | Publication date |
---|---|
WO2021145856A1 (en) | 2021-07-22 |
EP4073743A1 (en) | 2022-10-19 |
US20220265237A1 (en) | 2022-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9990741B2 (en) | Motion correction in a projection domain in time of flight positron emission tomography | |
EP2441050B1 (en) | Time-of-flight positron emission tomography reconstruction using image content generated event-by-event based on time-of-flight information | |
JP5172347B2 (ja) | 反復的制約デコンボリューションによる核医学2次元平面画像の復元 | |
US9053569B2 (en) | Generating attenuation correction maps for combined modality imaging studies and improving generated attenuation correction maps using MLAA and DCC algorithms | |
US10925554B2 (en) | Outside-FOV activity estimation using surview and prior patient data in positron emission tomography | |
US20120278055A1 (en) | Motion correction in radiation therapy | |
US20220015721A1 (en) | Determination of metabolic rate from static pet scan | |
JP2014517277A (ja) | リストモードダイナミックイメージ再構成 | |
US7750304B2 (en) | Concurrent reconstruction using multiple bed frames or continuous bed motion | |
EP2984631A1 (en) | Method for modeling and accounting for cascade gammas in images | |
JP7242708B2 (ja) | リストモード再構成における定量化精度向上のための補正方法 | |
JP2020511659A (ja) | Petリストデータからの呼吸運動信号のノイズロバストリアルタイム抽出 | |
US11524178B2 (en) | Internal dose tomography | |
CN114930384A (zh) | Pet图像数据的实况显示 | |
EP4026103A1 (en) | Continuous bed motion acquisition with axially short phantom for pet imaging system setup and quality control | |
Bal et al. | Studying effects of missing data for clinical TOF PET | |
US11622737B2 (en) | Radiolabeled cell tracking and imaging | |
US11663758B2 (en) | Systems and methods for motion estimation in PET imaging using AI image reconstructions | |
Verra | Feasibility and Quality Assessment of Model-based Respiratory Motion Compensation in Positron Emission Tomography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |