Zhou Yun, Huang Sung-Cheng, Bergsneider Marvin, Wong Dean F
The Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins Medical Institutions, Baltimore, MD 21287, USA.
Neuroimage. 2002 Mar;15(3):697-707. doi: 10.1006/nimg.2001.1021.
The value of parametric images that represent both spatial distribution and quantification of the physiological parameters of tracer kinetics has long been recognized. However, the inherent high noise level of pixel kinetics of dynamic PET makes it unsuitable to generate parametric images of the microparameters of tracer kinetic model by conventional weighted nonlinear least squares (WNLS) fitting. Based on the concept that both spatial and temporal information should be integrated to improve parametric image quality, a nonlinear ridge regression with spatial constraint (NLRRSC) parametric imaging algorithm was proposed in this study. For NLRRSC, a term that penalizes local spatial variation of parameters was added to the cost function of WNLS fitting. The initial estimates and spatial constraint were estimated by component representation model (CRM) with cluster analysis. A hierarchical cluster with average linkage method was used to extract components. The ridge parameter was determined by linear ridge regression theory at each iteration, and a modified Gauss-Newton algorithm was used for minimizing the cost function. Results from a computer simulation showed that the percent mean square error of estimates obtained by NLRRSC can be decreased by 60-80% compared to that of WNLS. The parametric images estimated by NLRRSC are significantly better than the ones generated by WNLS. A highly correlated linear relationship was found between the ROI values calculated from the microparametric images generated by NLRRSC and estimates from ROI kinetic fitting. NLRRSC provided a reliable estimate of glucose metabolite uptake rate with a comparable image quality compared to Patlak analysis. In conclusion, NLRRSC is a reliable and robust parametric imaging algorithm for dynamic PET studies.
能够同时呈现示踪剂动力学生理参数的空间分布和定量信息的参数图像,其价值早已得到认可。然而,动态正电子发射断层扫描(PET)像素动力学固有的高噪声水平,使其不适用于通过传统加权非线性最小二乘法(WNLS)拟合来生成示踪剂动力学模型微观参数的参数图像。基于应整合空间和时间信息以提高参数图像质量这一概念,本研究提出了一种具有空间约束的非线性岭回归(NLRRSC)参数成像算法。对于NLRRSC,在WNLS拟合的代价函数中添加了一个惩罚参数局部空间变化的项。初始估计值和空间约束通过带有聚类分析的成分表示模型(CRM)来估计。采用平均连锁法的层次聚类来提取成分。在每次迭代时,根据线性岭回归理论确定岭参数,并使用改进的高斯 - 牛顿算法来最小化代价函数。计算机模拟结果表明,与WNLS相比,NLRRSC获得的估计值的均方误差百分比可降低60 - 80%。由NLRRSC估计的参数图像明显优于由WNLS生成的图像。在从NLRRSC生成的微观参数图像计算得到的感兴趣区域(ROI)值与ROI动力学拟合估计值之间发现了高度相关的线性关系。与Patlak分析相比,NLRRSC在提供可靠的葡萄糖代谢物摄取率估计值的同时,具有相当的图像质量。总之,NLRRSC是一种用于动态PET研究的可靠且稳健的参数成像算法。