Characterization and reduction of noise in dynamic PET data using masked volumewise principal component analysis

Quantitation of dynamic total-body PET imaging: recent developments and future perspectives

BACKGROUND

Positron emission tomography (PET) scanning is an important diagnostic imaging technique used in disease diagnosis, therapy planning, treatment monitoring, and medical research. The standardized uptake value (SUV) obtained at a single time frame has been widely employed in clinical practice. Well beyond this simple static measure, more detailed metabolic information can be recovered from dynamic PET scans, followed by the recovery of arterial input function and application of appropriate tracer kinetic models. Many efforts have been devoted to the development of quantitative techniques over the last couple of decades.

CHALLENGES

The advent of new-generation total-body PET scanners characterized by ultra-high sensitivity and long axial field of view, i.e., uEXPLORER (United Imaging Healthcare), PennPET Explorer (University of Pennsylvania), and Biograph Vision Quadra (Siemens Healthineers), further stimulates valuable inspiration to derive kinetics for multiple organs simultaneously. But some emerging issues also need to be addressed, e.g., the large-scale data size and organ-specific physiology. The direct implementation of classical methods for total-body PET imaging without proper validation may lead to less accurate results.

CONCLUSIONS

In this contribution, the published dynamic total-body PET datasets are outlined, and several challenges/opportunities for quantitation of such types of studies are presented. An overview of the basic equation, calculation of input function (based on blood sampling, image, population or mathematical model), and kinetic analysis encompassing parametric (compartmental model, graphical plot and spectral analysis) and non-parametric (B-spline and piece-wise basis elements) approaches is provided. The discussion mainly focuses on the feasibilities, recent developments, and future perspectives of these methodologies for a diverse-tissue environment.

Spectral-resolved cone-beam X-ray luminescence computed tomography with principle component analysis

Cone-beam X-ray luminescence computed tomography (CB-XLCT) has become a promising technique for its higher utilization of X-ray and shorter scanning time compared to the narrow-beam XLCT, but it suffers from the low-spatial resolution that results in the insufficiency to resolve the adjacent multiple probes. In multispectral CB-XLCT, multiple probes show different emission behaviors in the dimension of the spectrum. In this work, a spectral-resolved CB-XLCT method combining multispectral CB-XLCT with principle component analysis (PCA) was proposed to improve the imaging resolution. Results of digital simulation and the phantom experiment illustrated that the proposed method was capable of resolving adjacent multiple probes accurately and had better performance than the common multispectral CB-XLCT with spectrum information priori.

Improved accuracy of quantitative parameter estimates in dynamic contrast-enhanced CT study with low temporal resolution

PURPOSE

A previously proposed method to reduce radiation dose to patient in dynamic contrast-enhanced (DCE) CT is enhanced by principal component analysis (PCA) filtering which improves the signal-to-noise ratio (SNR) of time-concentration curves in the DCE-CT study. The efficacy of the combined method to maintain the accuracy of kinetic parameter estimates at low temporal resolution is investigated with pixel-by-pixel kinetic analysis of DCE-CT data.

METHODS

The method is based on DCE-CT scanning performed with low temporal resolution to reduce the radiation dose to the patient. The arterial input function (AIF) with high temporal resolution can be generated with a coarsely sampled AIF through a previously published method of AIF estimation. To increase the SNR of time-concentration curves (tissue curves), first, a region-of-interest is segmented into squares composed of 3 × 3 pixels in size. Subsequently, the PCA filtering combined with a fraction of residual information criterion is applied to all the segmented squares for further improvement of their SNRs. The proposed method was applied to each DCE-CT data set of a cohort of 14 patients at varying levels of down-sampling. The kinetic analyses using the modified Tofts' model and singular value decomposition method, then, were carried out for each of the down-sampling schemes between the intervals from 2 to 15 s. The results were compared with analyses done with the measured data in high temporal resolution (i.e., original scanning frequency) as the reference.

RESULTS

The patients' AIFs were estimated to high accuracy based on the 11 orthonormal bases of arterial impulse responses established in the previous paper. In addition, noise in the images was effectively reduced by using five principal components of the tissue curves for filtering. Kinetic analyses using the proposed method showed superior results compared to those with down-sampling alone; they were able to maintain the accuracy in the quantitative histogram parameters of volume transfer constant [standard deviation (SD), 98th percentile, and range], rate constant (SD), blood volume fraction (mean, SD, 98th percentile, and range), and blood flow (mean, SD, median, 98th percentile, and range) for sampling intervals between 10 and 15 s.

CONCLUSIONS

The proposed method of PCA filtering combined with the AIF estimation technique allows low frequency scanning for DCE-CT study to reduce patient radiation dose. The results indicate that the method is useful in pixel-by-pixel kinetic analysis of DCE-CT data for patients with cervical cancer.

The 2D Hotelling filter - a quantitative noise-reducing principal-component filter for dynamic PET data, with applications in patient dose reduction

Background

In this paper we apply the principal-component analysis filter (Hotelling filter) to reduce noise from dynamic positron-emission tomography (PET) patient data, for a number of different radio-tracer molecules. We furthermore show how preprocessing images with this filter improves parametric images created from such dynamic sequence.

We use zero-mean unit variance normalization, prior to performing a Hotelling filter on the slices of a dynamic time-series. The Scree-plot technique was used to determine which principal components to be rejected in the filter process. This filter was applied to [¹¹C]-acetate on heart and head-neck tumors, [¹⁸F]-FDG on liver tumors and brain, and [¹¹C]-Raclopride on brain. Simulations of blood and tissue regions with noise properties matched to real PET data, was used to analyze how quantitation and resolution is affected by the Hotelling filter. Summing varying parts of a 90-frame [¹⁸F]-FDG brain scan, we created 9-frame dynamic scans with image statistics comparable to 20 MBq, 60 MBq and 200 MBq injected activity. Hotelling filter performed on slices (2D) and on volumes (3D) were compared.

Results

The 2D Hotelling filter reduces noise in the tissue uptake drastically, so that it becomes simple to manually pick out regions-of-interest from noisy data. 2D Hotelling filter introduces less bias than 3D Hotelling filter in focal Raclopride uptake. Simulations show that the Hotelling filter is sensitive to typical blood peak in PET prior to tissue uptake have commenced, introducing a negative bias in early tissue uptake. Quantitation on real dynamic data is reliable. Two examples clearly show that pre-filtering the dynamic sequence with the Hotelling filter prior to Patlak-slope calculations gives clearly improved parametric image quality. We also show that a dramatic dose reduction can be achieved for Patlak slope images without changing image quality or quantitation.

Conclusions

The 2D Hotelling-filtering of dynamic PET data is a computer-efficient method that gives visually improved differentiation of different tissues, which we have observed improve manual or automated region-of-interest delineation of dynamic data. Parametric Patlak images on Hotelling-filtered data display improved clarity, compared to non-filtered Patlak slope images without measurable loss of quantitation, and allow a dramatic decrease in patient injected dose.

Separating structures of different fluorophore concentrations by principal component analysis on multispectral excitation-resolved fluorescence tomography images

Multispectral excitation-resolved fluorescence tomography (MEFT) uses excitation light of different wavelengths to illuminate the fluorophores and obtains the reconstruction image frame which is fluorescence yield at each corresponding wavelength. For structures containing fluorophores of different concentrations, fluorescence yields show different variation trends with the excitation spectrum. In this study, principal component analysis (PCA) is used to analyze the MEFT reconstructed image frames. By taking advantage of the different variation trends of fluorescence yields, PCA can provide a set of principal components (PCs) in which structures containing different concentrations of fluorophores are shown separately. Simulations and experiments are both performed to test the performance of the proposed algorithm. The results suggest that the location and structure of fluorophores with different concentrations can be obtained and the contrast of fluorophores can be improved further by using this algorithm.