CN113039582A - 用于步进式拍摄和运动补偿正电子发射断层摄影(pet)研究的边缘保留惩罚重建 - Google Patents
用于步进式拍摄和运动补偿正电子发射断层摄影(pet)研究的边缘保留惩罚重建 Download PDFInfo
- Publication number
- CN113039582A CN113039582A CN201980075271.4A CN201980075271A CN113039582A CN 113039582 A CN113039582 A CN 113039582A CN 201980075271 A CN201980075271 A CN 201980075271A CN 113039582 A CN113039582 A CN 113039582A
- Authority
- CN
- China
- Prior art keywords
- volume
- penalty function
- voxel
- local
- image
- 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
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/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/037—Emission tomography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/424—Iterative
Abstract
一种存储能够由工作站(18)读取和执行的指令的非瞬态计算机可读介质,所述工作站包括至少一个电子处理器(20)以执行成像方法(100)。该方法包括:针对沿着轴向方向的帧在逐帧的基础上接收成像数据,其中,相邻帧沿着所述轴向方向交叠,其中,所述帧至少包括体积(k)和沿着所述轴向方向至少部分地与所述体积(k)交叠的后继体积(k+1);并且使用迭代图像重建过程来生成所述体积(k)的图像,其中,所述迭代图像重建过程的迭代包括:计算用于抑制所述体积(k)上的噪声的局部惩罚函数,包括减小交叠区域中的所述局部惩罚函数的值;使用来自所述体积(k)的成像数据并且还使用所述局部惩罚函数来生成所述体积(k)的更新图像。
Description
技术领域
以下总体上涉及医学成像技术、正电子发射断层摄影(PET)成像和图像重建技术、正则化迭代图像重建技术以及相关技术。
背景技术
正电子发射断层摄影术(PET)中的图像质量强烈地受噪声主导。使用边缘保留或对比度保留先验的迭代最大后验惩罚重建(PR)是一种抑制噪声同时保留可能具有临床意义的小特征(例如肿瘤)的方法。一些边缘保留或对比度保留惩罚包括:(1)相对差异惩罚(RDP)(例如参见J.Nuyts等人,“A concave prior penalizing relative differencesfor maximum-a-posteriori reconstruction in emission tomography”IEEETrans.Nucl.Sci.,2002年);(2)各向异性扩散惩罚(ADP)(例如,参见H.Zhu,“Imagereconstruction for positron emission tomography using fuzzy nonlinearanisotropic diffusion penalty”,Med Bio Eng Comput,2006年),等等。在大多数保留边缘惩罚中,调谐参数建立对重建图像中边缘过渡的灵敏度,并且低于灵敏度阈值的成像数据由强大的低通滤波器进行平滑。
由于PET探测器环或阵列的有限轴向范围,商业PET相机具有有限的轴向视场(FOV)。因此,通常在多个床位置处以系列“步进式拍摄”模式采集“全身”患者扫描,在床位置之间有特定交叠。出于实际原因,通常不等待下一个床位置的PET数据被完全采集而分开地重建每个床位置,并且将重建的个体帧图像编织为单幅全身图像。每个床位置交叠,以确保在轴向方向上平滑的信号转换。通常,对于全身研究,超过50%的扫描轴向范围将在交叠的切片中覆盖。因此,在每个重建图像体积中,将存在由从上一个床位置或下一个床位置采集的计数所覆盖的一个或两个子体积(注意,如本文所讨论的,假设在所有情况下考虑为床位置的数量至少大于1)。
代替于分开地重建每个床位置,并且然后在图像空间中将图像编织在一起,在备选方法中,在重建之前将计数预先组合(以正弦图格式或列表模式),而不会使成像工作流程过度复杂(参见,例如,Z.Sun等人,“Reconstruction and combination of PET multi-bed image”,美国公开US 2017/0103551)。然而,先前的工作示出,至少有序子集期望最大化(OSEM)迭代图像重建算法说明拟线性行为,并且在重建之前将计数预组合有很少实际好处,除了超低计数研究(参见,例如,S.Ross等人,“A method of overlap correction forfully 3D OSEM reconstruction of PET data”,MIC,2004)。在分别重建每个床位置的情况下,由于边缘处的数据的有效“丢失”,每幅床位置图像在边缘切片中能够具有较嘈杂的外观。但是,通常可以在图像编织流程期间(即,当将个体床位置图像组合成单幅全身图像时)补偿噪声。
以下公开了解决这些问题的新的和改进的系统和方法。
发明内容
在一个公开的方面中,一种非瞬态计算机可读介质存储指令,所述指令能够由工作站读取和执行,所述工作站包括至少一个电子处理器以执行成像方法。所述方法包括:针对沿着轴向方向的帧在逐帧的基础上接收成像数据,其中,相邻帧沿着所述轴向方向交叠,其中,所述帧至少包括体积(k)和沿着所述轴向方向至少部分地与所述体积(k)交叠的后继体积(k+1);并且使用迭代图像重建过程来生成所述体积(k)的图像,其中,所述迭代图像重建过程的迭代包括:计算用于抑制所述体积(k)上的噪声的局部惩罚函数,包括减小交叠区域中的所述局部惩罚函数的值;使用来自所述体积(k)的成像数据并且还使用所述局部惩罚函数来生成所述体积(k)的更新图像。
在另一个公开的方面中,一种成像系统,包括:正电子发射断层摄影(PET)成像设备。至少一个电子处理器被编程为:针对沿着轴向方向的帧在逐帧的基础上接收成像数据,其中,相邻帧沿着所述轴向方向交叠,其中,所述帧至少包括体积(k)和沿着所述轴向方向至少部分地与所述体积(k)交叠的后继体积(k+1);并且使用迭代图像重建过程来生成所述体积(k)的图像,其中,所述迭代图像重建过程的迭代包括:计算用于抑制所述体积(k)上的噪声的局部惩罚函数,包括减小交叠区域中的所述局部惩罚函数的值,所述局部惩罚函数取决于穿过被惩罚的所述体素的总符合计数,并且所述局部惩罚函数的值在所述交叠区域中减小补偿穿过在所述后继体积k+1中惩罚的所述体素的所述总符合计数的量;使用来自所述体积(k)的成像数据并且还使用所述局部惩罚函数来生成所述体积(k)的更新图像。所述重建包括:将来自所述成像数据中的所述体积(k)的计数的数量组合到总体积中;并且将所述总体积重建为重建图像。
在另一个公开的方面中,一种针对沿着轴向方向的帧在逐帧的基础上对成像数据执行的成像方法,其中,相邻帧沿着所述轴向方向交叠,其中,所述帧至少包括体积(k)和沿着所述轴向方向至少部分地与所述体积(k)交叠的后继体积(k+1),所述方法包括:使用迭代图像重建过程来生成所述体积(k)的图像,其中,所述迭代图像重建过程的迭代包括:计算用于抑制所述体积(k)上的噪声的局部惩罚函数,包括减小交叠区域中的所述局部惩罚函数的值,所述局部惩罚函数取决于所述PET成像设备的局部几何灵敏度,并且所述局部惩罚函数的值在所述交叠区域中减小补偿由于所述体积(k)的轴向截断而造成的几何灵敏度损失的量;使用来自所述体积(k)的成像数据并且还使用所述局部惩罚函数来生成所述体积(k)的更新图像。所述重建包括:与所述成像数据的其他体积(k)分开地将所述成像数据的所述体积(k)重建为重建图像。
在另一个公开方面中,一种针对沿着轴向方向的帧在逐帧的基础上对成像数据执行的成像方法,其中,相邻帧沿着所述轴向方向交叠,其中,所述帧至少包括体积(k)和沿着所述轴向方向至少部分地与所述体积(k)交叠的后继体积(k+1),所述方法包括:使用迭代图像重建过程来生成所述体积(k)的图像,其中,所述迭代图像重建过程的迭代包括:使用两个或更多个门相位和在所述门相位中检测到的运动利用呼吸门控来调节所述成像数据;计算用于抑制所述体积(k)上的噪声的局部惩罚函数,包括减小交叠区域中的所述局部惩罚函数的值;使用来自所述体积(k)的成像数据并且还使用所述局部惩罚函数来生成所述体积(k)的更新图像。
一个优点在于利用来自PET成像设备的所有可用计数来适当地调节重建惩罚函数的参数。
另一个优点在于提供一种保留边缘的惩罚重建过程,其避免了导致潜在丢失图像特征的最终图像的过度平滑。
另一个优点在于提供了一种保留边缘的惩罚重建过程,其避免了导致图像中的潜在噪声的最终图像的平滑不足。
另一个优点在于,当在图像重建之前将成像数据的交叠区域中的数据进行组合时,与通过组合交叠区域中的相邻图像而引入的有效灵敏度增加成比例地减小惩罚函数,以保持重建图像的空间分辨率。
另一个优点在于,当在图像重建之前不组合成像数据的交叠区域中的数据时,减小惩罚函数以保留重建图像的空间分辨率。
另一个优点在于利用两个成像帧之间的采集持续时间比率来缩放惩罚函数。
另一个优点在于利用计数端点调节的运动校正呼吸或心脏门控或动态研究动态调节惩罚函数的惩罚参数。
给定的实施例可以不提供任何一个前述优点,可以提供一个、两个、更多或所有前述优点,和/或可以提供其他优点,如对于本领域的普通技术人员在阅读和理解本公开后将变得显而易见的。
附图说明
本公开可以采取各种部件和各部件布置的形式,并且可以采取各种步骤和各步骤安排的形式。附图仅出于图示优选实施例的目的,并且不应被解释为对本公开的限制。
图1示意性地示出了根据一个方面的成像系统;
图2和图3示出了图1的系统的示例性流程图操作;并且
图4a-c示出了标准重建图像与利用图2的流程图操作重建的图像的比较。
具体实施方式
迭代惩罚PET图像重建采用了一种更新函数,其包括一个惩罚项,被指定为惩罚(即抑制)噪声,同时保留边缘。例如,相对差异惩罚(RDP)采用惩罚P∝1/γ,其中,项γ是边缘保留参数。较大的γ值产生较小的惩罚P,其增强了对特征(即“边缘”)(例如小病变)的保留,但也导致保留了更多噪声;而较小的γ值增加惩罚,从而改进噪声抑制,但也会导致去除临床上重要的结构(例如“边缘”)(例如小病变)。在变型方法中,例如通过使用对惩罚进行局部调谐,其中,λi是由i索引的体素的值,而si是与穿过体素i的响应线(LOR)的数量成比例的体素i处的几何灵敏度。此公式的前提是以下预期:如果更多的LOR穿过体素i(因此更高si),则附近应该有更多的计数,并且因此体素i处的噪声应该更低,使得可以使用较小的惩罚P增强边缘,同时仍然抑制(相对低的)噪声。在备选调谐中,将几何灵敏度si替换为量化体素i周围的局部邻域中的总计数的项。再次,局部邻域中较高的总计数倾向于与较低的噪声相对应,使得可以使用较低的惩罚P。
在本文中认识到,在步进式拍摄PET成像的交叠区域的情况下,这种惩罚框架可能不太有效。在步进式拍摄方法(也称为多床成像或类似命名法)中,使患者步进通过扫描器以对患者的连续轴向部分(“床位置”)成像。因此,在每个床位置处采集的数据在对应于该床位置处的PET扫描器最外侧探测器环的位置的边缘切片处被轴向截断。通常,在每个床位置处采集的PET计数是分别重建的(但可能使用来自相邻床位置的信息来估计FOV外散射),并且通过对床位置图像的交叠区域中的体素值进行加权平均在交叠处将得到的床位置图像编织在一起。然而,在这种方法中,使用针对该床位置的几何灵敏度si计算惩罚的局部调谐。对于床位置的外部切片,灵敏度si变小,因为终止于床位置轴向FOV之外的LOR不贡献于灵敏度si。对于极端边缘切片,在切片的平面中的仅LOR贡献于si。这意味着惩罚在外部切片处变大,导致交叠区域中的过度噪声抑制,并且增强抑制临床相关的图像特征(例如小肿瘤)的潜力。随后在交叠处将相邻的床位置图像编织在一起不能恢复这些抑制的特征,因为两个交叠区域将在惩罚方面具有该相同的“边缘效应”,从而导致特征的丢失。si替换为局部邻域中的总计数的备选公式化也会类似地受到影响,因为针对在床位置的截断轴向FOV之外终止的LOR,计数没有采集。
在一些公开的实施例中,减小交叠区域中的局部惩罚以便考虑上述效果。
在本文公开的一些实施例中,减小与如果将截断的体积k继续到相邻的交叠的截断的体积k+1中则将实现的有效灵敏度增加(或有效总邻域计数增加)成比例。在一种采用几何灵敏度以调节惩罚的方法中,这相当于在交叠区域中用代替其中,si,k+1项考虑不是体积k的部分(并且因此没有贡献于几何灵敏度si,k)但是是交叠体积k+1的部分(并且因此贡献于几何灵敏度si,k+1)的那些LOR。如果通采用来自交叠体积k和k+1两者的体素i的邻域中的总计数的加和,使用局部邻域中的总计数对γ进行局部调节,则可以使用类似的公式化。
在上面的实施例中,每个体积k仍然分别被重建。在备选联合重建实施例中,来自所有体积的PET计数首先被组合,并且然后一起重建。在这种情况下,再次使用此处,合适的公式化可以是其中,λi是通过组合数据集的(单个)联合重建产生的体素i的值。
在变型实施例中,如果各个床位置的采集时间不同,则灵敏度si,k和si,k+1与采集时间成比例地缩放。这是由于较短的采集时间导致较少的计数被采集和相应较低的灵敏度。因此,如果对于相应交叠床位置k和k+1,采集时间是Tk和Tk+1,则缩放为和
在本文公开的其他实施例中,类似方法可以应用于使用(例如呼吸)门控的运动补偿的PET成像。在这种情况下,考虑两个相继的门相位φk和φk+1,其中,在两个相位之间发生由运动向量Δr表示的一些运动。然后,经调节的边缘保留参数的可能公式化为
参考图1,示出了说明性医学成像系统10。如图1所示,系统10包括成像或图像采集设备12。在一个示例中,图像采集设备12可以包括PET成像设备,包括PET机架13和辐射探测器14的阵列(在图1中图解地指示;通常,PET机架的辐射探测器被布置为一系列PET探测器环,其被布置为跨越轴向FOV)。说明性图像采集设备12是PET/CT扫描器,其还包括计算机断层摄影(CT)机架15。患者台(或床)16被布置为将患者加载到成像设备12的检查区域17中,例如到CT机架15的膛中或PET机架13的膛中。通常,CT机架15被用于采集对象的CT图像,该CT图像被转换为对象的衰减图。PET机架13从CT图像采集使用衰减图重建的PET成像数据,以在PET图像重建中考虑患者内的辐射吸收。
系统10还包括具有典型部件的计算机或工作站或其他电子数据处理设备18,例如至少一个电子处理器20、至少一个用户输入设备(例如,鼠标、键盘、轨迹球、触控板和和/或类似物)22和显示设备24(例如,LCD显示器、OLED显示器、等离子显示器等)。在一些实施例中,显示设备24可以是与计算机18分开的部件,和/或可以包括两个或更多个显示器。工作站18还可以包括一个或多个数据库或非瞬态存储介质26(例如磁盘、RAID或其他磁存储介质;固态驱动器、闪存驱动器、电可擦除只读存储器(EEROM)或其他电子存储器;光盘或其他光存储设备;它们的各种组合;等等)。显示设备24被配置为显示由成像系统10采集的图像,并且通常还显示包括各种用户对话框的图形用户接口(GUI)28,例如,每个具有一个或多个字段、径向选择按钮等,以接收来自用户输入设备22的用户输入。
至少一个电子处理器20与一个或多个数据库26操作地连接,该数据库存储由至少一个电子处理器20可读和可执行以执行所公开的操作的指令,所公开的操作包括执行成像方法或过程100。在一些示例中,成像方法或过程100可以至少部分地通过云处理来执行。
参考图2,图解地示出了成像方法100的说明性实施例作为流程图。在102处,对至少一个电子处理器20进行编程以接收,以其他方式操作或控制CT/PET成像设备12(并且更具体地,PET机架13),以针对沿轴向方向的帧在逐帧的基础上采集PET成像数据,其中,相邻帧沿着轴向方向交叠,其中,帧至少包括体积(k)和沿着轴向方向与该体积(k)至少部分交叠的后继体积(k+1)。在104处,至少一个电子处理器20被编程为使用包括操作106和108的迭代图像重建过程来生成体积(k)的图像。在106处,计算局部惩罚函数以抑制体积(k)上的噪声,包括减小交叠区域中的局部惩罚函数的值。在108处,使用来自体积(k)的成像数据并且还使用局部惩罚函数来生成体积(k)的更新图像。尽管未图示,但是将理解,迭代图像重建104可以使用从使用CT机架15采集的对象的CT图像生成的衰减图来考虑对象的辐射吸收。
在本文公开的一个实施例中,局部惩罚函数取决于PET成像设备12的局部几何灵敏度,并且局部补偿函数的值在交叠区域中减小补偿由于体积(k)的轴向截断而造成的几何灵敏度损失的量。接下来描述说明性示例,其中,边缘保留惩罚函数是相对差异先验(RDP),并且使用有序子集期望最大化(OSEM)重建。对于迭代n的索引为n+1的更新图像适当地写为:
其中,λi是由i索引的体素处的估计活动,并且是局部惩罚加权因子,ajj′是系统矩阵值,si,k是体素i处的几何灵敏度(如果不同的床位置具有不同的采集时间,则优选地按采集时间缩放),yj是数据投影分箱,并且参数γi>0是边缘保留参数并且操纵RDP先验。先验是在体素i周围的局部图像邻域Ni上估计的。通常,较大的γi值产生较大的边缘保留(即,减少惩罚)。值γi=0将完全消除边缘保留,并且RDP变为二次先验。
如本文所公开的,在体积k和k+1的床位置交叠的交叠区域中(以及如果当前床位置k在两端交叠,则体积k-1和k的床位置交叠的交叠区域中),局部惩罚函数的值减小。在一些合适的实施例中,局部惩罚函数取决于PET成像设备12对于在体积k的床位置的图像体积中由i索引的体素的局部几何灵敏度si,k。体素i处的几何灵敏度可以视为穿过该体素i的响应线(LOR)的数量。在体积k中与相邻体积k-1或k+1都不交叠的部分中,可以使用标准灵敏度,并且对于第n次迭代,可以使用值其中,γ0是常数。然而,在交叠区域中,惩罚减小。这是通过将交叠区域中的边缘保留参数γ的值增加一量来补偿的,该量补偿由于体积k的轴向截断而造成的几何灵敏度的损失。因此,使用“联合”边缘保留参数其可以适当地写为其中,si,k+1是下一个相邻体积k+1中(相同)体素i的几何灵敏度。再次强调,该调节增加边缘保留参数(即),其由于γ位于等式(1)的局部RDP惩罚的分母中而导致减小交叠区域中的局部RDP惩罚。RDP惩罚的这种减小考虑在体积k+1中的额外数据,这意味着当在交叠区域中将两个重建体积k和k+1编织在一起时噪声较低,使得当重建体积k以实现期望的降噪时可以在该交叠区域中采用较低的惩罚,并且因此由于交叠区域中的减小的惩罚,更好地保留特征(边缘)。
在另一个实施例中,局部惩罚函数取决于穿过被惩罚的体素的总符合计数,并且局部补偿函数的值在交叠区域中减小补偿穿过在后继体积k+1中被惩罚的体素的总符合计数的量。惩罚可以根据实际的局部符合计数统计来制定,而不是根据几何灵敏度。这需要用实际计数替换以上表达式中的几何灵敏度s。在这种情况下,(无交叠调节),其中,Ci,k是数据集中针对其LOR穿过体素i的床位置k的符合事件的计数。针对飞行时间(TOF)PET,可以将这公式化为其中,索引v遍历针对其LOR穿过体素i的床位置k的数据集的Ci,k符合计数,并且wv是源自于体素i的由v索引的计数的飞行时间概率。对于非TOF情况,联合边缘保留参数(即,提供交叠调节)为其中,再次类似地,项Ci,k+1是数据集中针对其LOR通过(相同)体素i的床位置k+1的符合事件计数。在这种公式化中,甚至更明确的是,减小的惩罚(其再次由分母项的增加值提供)考虑了体积k+1中的额外数据Ci,k+1,这意味着,当在交叠区域中将两个重建体积k和k+1编织在一起时噪音较低。TOF的等效公式化为其中,加和位于体积k+1的数据集中。
在以上示例中,重建过程106、108包括:将成像数据的每个体积(k)与成像数据的其他体积(k)分开地重建为重建床位置图像,并且然后在图像空间中将床位置图像编织在一起以生成最终图像。在其他实施例中,重建过程106、108可以包括将来自成像数据中的所有体积(k)的计数的数量组合到总体积中;并且将总体积重建为重建图像。在后一个实施例中,局部惩罚函数取决于
在任一实施例中,如果连续床位置的采集时间不是全部相同,则还基于体积(k)和后继体积k+1的相应采集时间的比率适当地缩放交叠区域中的局部惩罚函数的值的减小。
参考图3,图解地示出了呼吸门控成像方法200的说明性实施例作为流程图。这是不同的应用,但是所公开的在(现在是时间的,而不是空间/轴向的)交叠区域中调节惩罚的方法是相似的。在202处,至少一个电子处理器20被编程为操作或控制PET成像设备12以采集呼吸患者的成像数据。与PET成像数据采集同时,由呼吸监测器201采集呼吸门控信号,呼吸监测器201提供呼吸相位与时间的信息。例如,呼吸门控信号可以由呼吸监测器201提供,呼吸监测器201包括呼吸监测带,该呼吸监测带测量呼吸期间胸部的扩张和收缩;或者呼吸门控信号可以由呼吸监测器201提供,该呼吸监测器201包括流量传感器,该流量传感器测量流入和流出口腔和/或鼻子的空气流量;等等。在203处,使用呼吸门控信号通过呼吸相位将PET数据进行分箱。通过类似于图2的多床位置成像示例,每个呼吸相位分箱由索引k指代,并且两个相邻的呼吸相位分箱k和k+1在呼吸序列中相邻,即在每个呼吸间隔中呼吸相位k+1在时间上紧随呼吸相位k之后。应当理解,按呼吸相位的分箱将导致在连续呼吸中在相同呼吸相位k上采集的PET数据被组合以形成用于呼吸相位分箱k的PET数据集;同样,在连续呼吸中在相同呼吸相位k+1上采集的PET数据将被组合以形成呼吸相位分箱k+1的PET数据集;等等。此外,分箱在时间上交叠。在204处,至少一个电子处理器20被编程为使用迭代图像重建过程来生成呼吸相位(k)的图像,该迭代图像重建过程对在多次呼吸中累积在呼吸相位分箱(k)中的所采集的PET成像数据进行操作。迭代重建204包括操作206、208和210。在206处,针对在门相位中检测到的呼吸运动调节每个分箱的成像数据。这需要根据导出的运动向量来移位LOR端点的计数。对于每个体素i,都有一个相关联的运动向量Δri,指示体素i从其参考位置(例如,可以任意指定其在分箱k=0处的位置,使得对于所有i,Δri,k=0=0)开始的运动。在208处,计算局部惩罚函数以抑制呼吸相位分箱(k)中的噪声,包括减小与呼吸相位分箱k+1的时间交叠区域中的局部惩罚函数的值。在210处,使用来自呼吸相位分箱(k)的成像数据并且还使用局部惩罚函数来生成呼吸相位分箱(k)的更新图像。
在一些实施例中,根据以下公式执行在操作206处的调节:
其中,γ是边缘保留参数,i是体素,λ是体素的灰度值,s是穿过体素的几何灵敏度,并且Δr是两个门相位φk和φk+1之间的运动向量。
图4示出了使用均匀体模的重建图像的示例,其具有位于交叠区域中心的六边形图案的六个6mm直径的球体。标记为(a)的图是编织后通过标准惩罚性重建重建的图像。标记为(b)的图是通过成像方法100的联合惩罚重建而重建的图像。模型设置的侧视图示于标记为(c)的图中。示出了圆柱体模型50,其在第一床位置54和第二床位置56上具有多个球体52。
已经参考优选实施例描述了本公开。在阅读和理解了前面的详细描述之后,其他人可能会想到修改和变更。本公开旨在解释为包括所有这样的修改和改变,只要它们落入所附权利要求或其等同物的范围之内。
Claims (20)
1.一种存储能够由工作站(18)读取和执行的指令的非瞬态计算机可读介质,所述工作站包括至少一个电子处理器(20)以执行成像方法(100),所述方法包括:
针对沿着轴向方向的帧在逐帧的基础上接收成像数据,其中,相邻帧沿着所述轴向方向交叠,其中,所述帧至少包括体积(k)和沿着所述轴向方向至少部分地与所述体积(k)交叠的后继体积(k+1);并且
使用迭代图像重建过程来生成所述体积(k)的图像,其中,所述迭代图像重建过程的迭代包括:
计算用于抑制所述体积(k)上的噪声的局部惩罚函数,包括减小交叠区域中的所述局部惩罚函数的值;
使用来自所述体积(k)的成像数据并且还使用所述局部惩罚函数来生成所述体积(k)的更新图像。
2.根据权利要求1所述的非瞬态计算机可读介质,其中,所述局部惩罚函数取决于所述PET成像设备(12)的局部几何灵敏度,并且所述局部惩罚函数的所述值在所述交叠区域中减小补偿由于所述体积(k)的轴向截断而造成的几何灵敏度损失的量。
4.根据权利要求1所述的非瞬态计算机可读介质,其中,所述局部惩罚函数取决于穿过被惩罚的所述体素的总符合计数,并且所述局部惩罚函数的所述值在所述交叠区域中减小补偿穿过在所述后继体积k+1中惩罚的所述体素的所述总符合计数的量。
7.根据权利要求1-6中的任一项所述的非瞬态计算机可读介质,其中,所述重建包括:
与所述成像数据的其他体积(k)分开地将所述成像数据的所述体积(k)重建为重建图像。
8.根据权利要求1-6中的任一项所述的非瞬态计算机可读介质,其中,所述重建包括:
将来自所述成像数据中的所述体积(k)的计数的数量组合到总体积中;并且
将所述总体积重建为重建图像。
10.根据权利要求1-8中的任一项所述的非瞬态计算机可读介质,其中,基于针对所述体积(k)与所述后继体积k+1的相应采集时间的比率来缩放所述交叠区域中的所述局部惩罚函数的所述值的所述减小。
11.一种成像系统(10),包括:
正电子发射断层摄影(PET)成像设备(12);以及
至少一个电子处理器(20),其被编程为:
针对沿着轴向方向的帧在逐帧的基础上接收成像数据,其中,相邻帧沿着所述轴向方向交叠,其中,所述帧至少包括体积(k)和沿着所述轴向方向至少部分地与所述体积(k)交叠的后继体积(k+1);并且
使用迭代图像重建过程来生成所述体积(k)的图像,其中,所述迭代图像重建过程的迭代包括:
计算用于抑制所述体积(k)上的噪声的局部惩罚函数,包括减小交叠区域中的所述局部惩罚函数的值,所述局部惩罚函数取决于穿过被惩罚的所述体素的总符合计数,并且所述局部惩罚函数的值在所述交叠区域中减小补偿穿过在所述后继体积k+1中惩罚的所述体素的所述总符合计数的量;
使用来自所述体积(k)的成像数据并且还使用所述局部惩罚函数来生成所述体积(k)的更新图像;
其中,所述重建包括:
将来自所述成像数据中的所述体积(k)的计数的数量组合到总体积中;并且
将所述总体积重建为重建图像。
15.根据权利要求11-14中的任一项所述的系统(10),其中,基于针对所述体积(k)与所述后继体积k+1的相应采集时间的比率来缩放所述交叠区域中的所述局部惩罚函数的值的所述减小。
16.一种针对沿着轴向方向的帧在逐帧的基础上对成像数据执行的成像方法(100),其中,相邻帧沿着所述轴向方向交叠,其中,所述帧至少包括体积(k)和沿着所述轴向方向至少部分地与所述体积(k)交叠的后继体积(k+1),所述方法包括:
使用迭代图像重建过程来生成所述体积(k)的图像,其中,所述迭代图像重建过程的迭代包括:
计算用于抑制所述体积(k)上的噪声的局部惩罚函数,包括减小交叠区域中的所述局部惩罚函数的值,所述局部惩罚函数取决于所述PET成像设备的局部几何灵敏度,并且所述局部惩罚函数的值在所述交叠区域中减小补偿由于所述体积(k)的轴向截断而造成的几何灵敏度损失的量;
使用来自所述体积(k)的成像数据并且还使用所述局部惩罚函数来生成所述体积(k)的更新图像;
其中,所述重建包括:
与所述成像数据的其他体积(k)分开地将所述成像数据的所述体积(k)重建为重建图像。
18.根据权利要求16和17中的任一项所述的成像方法(100),其中,基于针对所述体积(k)与所述后继体积k+1的相应采集时间的比率来缩放所述交叠区域中的所述局部惩罚函数的所述值的所述减小。
19.一种针对沿着轴向方向的帧在逐帧的基础上对成像数据执行的成像方法(200),其中,相邻帧沿着所述轴向方向交叠,其中,所述帧至少包括体积(k)和沿着所述轴向方向至少部分地与所述体积(k)交叠的后继体积(k+1),所述方法包括:
使用迭代图像重建过程来生成所述体积(k)的图像,其中,所述迭代图像重建过程的迭代包括:
使用两个或更多个门相位和在所述门相位中检测到的运动利用呼吸门控来调节所述成像数据;
计算用于抑制所述体积(k)上的噪声的局部惩罚函数,包括减小交叠区域中的所述局部惩罚函数的所述值;
使用来自所述体积(k)的成像数据并且还使用所述局部惩罚函数来生成所述体积(k)的更新图像。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201862768131P | 2018-11-16 | 2018-11-16 | |
US62/768,131 | 2018-11-16 | ||
PCT/EP2019/081515 WO2020099649A1 (en) | 2018-11-16 | 2019-11-15 | Edge preserving penalized reconstruction for step-and-shoot and motion compensated positron emission tomography (pet) studies |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113039582A true CN113039582A (zh) | 2021-06-25 |
Family
ID=68583433
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201980075271.4A Pending CN113039582A (zh) | 2018-11-16 | 2019-11-15 | 用于步进式拍摄和运动补偿正电子发射断层摄影(pet)研究的边缘保留惩罚重建 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20220012928A1 (zh) |
EP (1) | EP3881290A1 (zh) |
JP (1) | JP7404362B2 (zh) |
CN (1) | CN113039582A (zh) |
WO (1) | WO2020099649A1 (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102324089A (zh) * | 2011-07-13 | 2012-01-18 | 南方医科大学 | 基于广义熵与mr先验的pet图像最大后验重建方法 |
CN102846326A (zh) * | 2011-06-03 | 2013-01-02 | 通用电气公司 | 用于处理选通图像数据的方法和系统 |
CN103299345A (zh) * | 2011-01-10 | 2013-09-11 | 皇家飞利浦电子股份有限公司 | 双能量断层摄影成像系统 |
CN103559728A (zh) * | 2013-10-31 | 2014-02-05 | 南方医科大学 | 基于解剖功能联合先验模型的pet图像最大后验重建方法 |
US20160063741A1 (en) * | 2014-08-28 | 2016-03-03 | Kabushiki Kaisha Toshiba | Method and Apparatus for Estimating Scatter in a Positron Emission Tomography Scan at Multiple Bed Positions |
CN106466188A (zh) * | 2015-08-20 | 2017-03-01 | 通用电气公司 | 用于发射断层显像定量的系统和方法 |
US20170061629A1 (en) * | 2015-08-25 | 2017-03-02 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image calibration |
US20170103551A1 (en) * | 2015-10-13 | 2017-04-13 | Shenyang Neusoft Medical Systems Co., Ltd. | Reconstruction and combination of pet multi-bed image |
CN107004259A (zh) * | 2014-12-10 | 2017-08-01 | 皇家飞利浦有限公司 | 多对比度成像中的统计加权正则化 |
US20170249758A1 (en) * | 2016-02-26 | 2017-08-31 | Wisconsin Alumni Research Foundation | System And Method For Resolving Artifacts In Four-Dimensional Angiographic Data |
CN108010093A (zh) * | 2016-10-31 | 2018-05-08 | 上海东软医疗科技有限公司 | 一种pet图像重建方法和装置 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7173248B2 (en) * | 2004-10-20 | 2007-02-06 | General Electric Company | Methods and systems for positron emission tomography data correction |
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 |
JP5804681B2 (ja) | 2010-08-30 | 2015-11-04 | 株式会社東芝 | 放射線イメージング装置及び画像処理装置 |
US9478049B2 (en) * | 2013-09-30 | 2016-10-25 | Koninklijke Philips N.V. | Method for local adjustment of regularization parameters for image quality optimization in fully 3D iterative CT reconstruction |
US9031300B1 (en) * | 2013-10-25 | 2015-05-12 | General Electric Company | System and method reconstructing a nuclear medicine image using deformed attenuation image |
US10022097B2 (en) * | 2014-06-16 | 2018-07-17 | University Of Southern California | Direct patlak estimation from list-mode PET data |
US9474495B2 (en) * | 2014-12-22 | 2016-10-25 | General Electric Company | System and method for joint estimation of attenuation and activity information |
-
2019
- 2019-11-15 EP EP19805275.5A patent/EP3881290A1/en active Pending
- 2019-11-15 WO PCT/EP2019/081515 patent/WO2020099649A1/en unknown
- 2019-11-15 JP JP2021525553A patent/JP7404362B2/ja active Active
- 2019-11-15 US US17/293,958 patent/US20220012928A1/en active Pending
- 2019-11-15 CN CN201980075271.4A patent/CN113039582A/zh active Pending
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103299345A (zh) * | 2011-01-10 | 2013-09-11 | 皇家飞利浦电子股份有限公司 | 双能量断层摄影成像系统 |
CN102846326A (zh) * | 2011-06-03 | 2013-01-02 | 通用电气公司 | 用于处理选通图像数据的方法和系统 |
CN102324089A (zh) * | 2011-07-13 | 2012-01-18 | 南方医科大学 | 基于广义熵与mr先验的pet图像最大后验重建方法 |
CN103559728A (zh) * | 2013-10-31 | 2014-02-05 | 南方医科大学 | 基于解剖功能联合先验模型的pet图像最大后验重建方法 |
US20160063741A1 (en) * | 2014-08-28 | 2016-03-03 | Kabushiki Kaisha Toshiba | Method and Apparatus for Estimating Scatter in a Positron Emission Tomography Scan at Multiple Bed Positions |
CN107004259A (zh) * | 2014-12-10 | 2017-08-01 | 皇家飞利浦有限公司 | 多对比度成像中的统计加权正则化 |
CN106466188A (zh) * | 2015-08-20 | 2017-03-01 | 通用电气公司 | 用于发射断层显像定量的系统和方法 |
US20170061629A1 (en) * | 2015-08-25 | 2017-03-02 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image calibration |
US20170103551A1 (en) * | 2015-10-13 | 2017-04-13 | Shenyang Neusoft Medical Systems Co., Ltd. | Reconstruction and combination of pet multi-bed image |
US20170249758A1 (en) * | 2016-02-26 | 2017-08-31 | Wisconsin Alumni Research Foundation | System And Method For Resolving Artifacts In Four-Dimensional Angiographic Data |
CN108010093A (zh) * | 2016-10-31 | 2018-05-08 | 上海东软医疗科技有限公司 | 一种pet图像重建方法和装置 |
Non-Patent Citations (1)
Title |
---|
马小波;桂志国;张銮;: "正电子发射断层图像重建中基于有序子集惩罚最小二乘算法的应用", 中国组织工程研究与临床康复, no. 43 * |
Also Published As
Publication number | Publication date |
---|---|
EP3881290A1 (en) | 2021-09-22 |
JP7404362B2 (ja) | 2023-12-25 |
WO2020099649A1 (en) | 2020-05-22 |
JP2022507135A (ja) | 2022-01-18 |
US20220012928A1 (en) | 2022-01-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8351671B2 (en) | Motion correction in nuclear imaging | |
US8331639B2 (en) | Radiological imaging incorporating local motion monitoring, correction, and assessment | |
EP2707853B1 (en) | List mode dynamic image reconstruction | |
US8600132B2 (en) | Method and apparatus for motion correcting medical images | |
US8675936B2 (en) | Multimodal image reconstruction | |
US11410349B2 (en) | Methods for data driven respiratory motion estimation | |
US11210820B2 (en) | Iterative image reconstruction with dynamic suppression of formation of noise-induced artifacts | |
US11064965B2 (en) | Noise-robust real-time extraction of the respiratory motion signal from PET list-data | |
CN113039582A (zh) | 用于步进式拍摄和运动补偿正电子发射断层摄影(pet)研究的边缘保留惩罚重建 | |
US11009615B2 (en) | Time-of-flight resolution-adaptive image regularization and filtering in positron emission tomography | |
CA3117317A1 (en) | Activity image reconstruction using anatomy data | |
US20220262048A1 (en) | Determination of motion frames based on image data | |
Jin et al. | Four-dimensional gated cardiac SPECT reconstruction and evaluation study | |
Bhagalia et al. | Non-rigid Registration-based Coronary Artery Motion Correction |
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 |