List-mode Positron Emission Tomography Data Research Articles

High-resolution motion compensation for brain PET imaging using real-time electromagnetic motion tracking.

Substantial improvements in spatial resolution in brain positron emission tomography (PET) scanners have greatly reduced partial volume effect, making head movement the main source of image blur. To achieve high-resolution PET neuroimaging, precise real-time estimation of both head position and orientation is essential for accurate motion compensation. A high-resolution electromagnetic motion tracking (EMMT) system with an event-by-event motion correction is developed for PET-CTscanners. EMMT is comprised of a source, an array of sensors, and a readout electronic unit (REU). The source acts as a transmitter and emits an EM dipole field. It is placed in close proximity to the sensor array and detects changes in EM flux density due to sensor movement. The REU digitizes signals from each sensor and captures precise rotational and translational movements in real time. Tracked motion in the EMMT coordinate system is synchronized with the PET list-mode data and transformed into the scanner coordinate system by locating paired positions in both systems. The optimal rigid motion is estimated using singular value decomposition. The rigid motion and depth-of-interaction (DOI) parallax effect are corrected by event-by-event rebinning of mispositioned lines-of-response (LORs). We integrated the EMMT with our recently developed ultra-high resolution Prism-PET prototype brain scanner and a commercial Siemens Biograph mCT PET-CT scanner. We assessed the imaging performance of the Prism-PET/EMMT system using multi-frame motion of point sources and phantoms. The mCT/EMMT system was validated using a set of point sources attached to both a mannequin head and a human volunteer, for simulating multiframe and continuous motions, respectively. Additionally, a human subject for [18F]MK6240 PET imaging wasincluded. The tracking accuracy of the Prism-PET/EMMT system was quantified as a root-mean-square (RMS) error of 0.49 for 100 axial rotations, and an RMS error of 0.15 mm for 100 mm translations.The percent difference (%diff) in average full width at half maximum (FWHM) of point source between motion-corrected and static images, within a motion range of and 10 mm from the center of the scanner's field-of-view (FOV), was 3.9%. The measured recovery coefficients of the 2.5-mm diameter sphere in the activity-filled partial volume correction phantom were 23.9%, 70.8%, and 74.0% for the phantom with multi-frame motion, with motion and motion compensation, and without motion, respectively. In the mCT/EMMT system, the %diff in average FWHM of point sources between motion-corrected and static images, within a motion range of and 10 mm from the center of the FOV, was 14%. Applying motion correction to the [18F]MK6240 PET imaging reduced the motion-induced spill-in artifact in the lateral ventricle region, lowering its standardized uptake value ratio (SUVR) from 0.70 to 0.34. The proposed EMMT system is a cost-effective, high frame-rate, and none-line-of-sight alternative to infrared camera-based tracking systems and is capable of achieving high rotational and translational tracking accuracies for mitigating motion-induced blur in high-resolution brain dedicated PETscanners.

Refining penalty parameter selection in whole-body PET image reconstruction for lung cancer patients using the cross-validation log-likelihood method

Objective.Penalty parameters in penalized likelihood positron emission tomography (PET) reconstruction are typically determined empirically. The cross-validation log-likelihood (CVLL) method has been introduced to optimize these parameters by maximizing a CVLL function, which assesses the likelihood of reconstructed images using one subset of a list-mode dataset based on another subset. This study aims to validate the efficacy of the CVLL method in whole-body imaging for cancer patients using a conventional clinical PET scanner.Approach.Fifteen lung cancer patients were injected with 243.7 ± 23.8 MBq of [18F]FDG and underwent a 22 min PET scan on a Biograph mCT PET/CT scanner, starting at 60 ± 5 min post-injection. The PET list-mode data were partitioned by subsampling without replacement, with 20 minutes of data for image reconstruction using an in-house ordered subset expectation maximization algorithm and the remaining 2 minutes of data for cross-validation. Two penalty parameters, penalty strengthβand Fair penalty function parameterδ, were subjected to optimization. Whole-body images were reconstructed, and CVLL values were computed across various penalty parameter combinations. The optimal image corresponding to the maximum CVLL value was selected by a grid search for each patient.Main results.Theδvalue required to maximize the CVLL value was notably small (⩽10-6in this study). The influences of voxel size and scan duration on image optimization were investigated. A correlation analysis revealed a significant inverse relationship between optimalβand scan count level, with a correlation coefficient of -0.68 (p-value = 3.5 × 10-5). The optimal images selected by the CVLL method were compared with those chosen by two radiologists based on their diagnostic preferences. Differences were observed in the selection of optimal images.Significance.This study demonstrates the feasibility of incorporating the CVLL method into routine imaging protocols, potentially allowing for a wide range of combinations of injected radioactivity amounts and scan durations in modern PET imaging.

Unsupervised deep learning framework for data-driven gating in positron emission tomography.

Physiological motion, such as respiratory motion, has become a limiting factor in the spatial resolution of positron emission tomography (PET) imaging as the resolution of PET detectors continue to improve. Motion-induced misregistration between PET and CT images can also cause attenuation correction artifacts. Respiratory gating can be used to freeze the motion and to reduce motion induced artifacts. In this study, we propose a robust data-driven approach using an unsupervised deep clustering network that employs an autoencoder (AE) to extract latent features for respiratory gating. We first divide list-mode PET data into short-time frames. The short-time frame images are reconstructed without attenuation, scatter, or randoms correction to avoid attenuation mismatch artifacts and to reduce image reconstruction time. The deep AE is then trained using reconstructed short-time frame images to extract latent features for respiratory gating. No additional data are required for the AE training. K-means clustering is subsequently used to perform respiratory gating based on the latent features extracted by the deep AE. The effectiveness of our proposed Deep Clustering method was evaluated using physical phantom and real patient datasets. The performance was compared against phase gating based on an external signal (External) and image based principal component analysis (PCA) with K-means clustering (Image PCA). The proposed method produced gated images with higher contrast and sharper myocardium boundaries than those obtained using the External gating method and Image PCA. Quantitatively, the gated images generated by the proposed Deep Clustering method showed larger center of mass (COM) displacement and higher lesion contrast than those obtained using the other two methods. The effectiveness of our proposed method was validated using physical phantom and real patient data. The results showed our proposed framework could provide superior gating than the conventional External method and Image PCA.

Deep learning based registration for head motion correction in positron emission tomography as a strategy for improved image quantification

Objectives: Positron emission tomography (PET) is affected by various kinds of patient movement during a scan. Frame-by-frame image registration is one of the most practiced motion correction techniques. In recent years, deep learning has shown a remarkable ability to quickly and accurately register images once trained. This paper studies the feasibility of using a deep learning framework to correct 3D positron emission tomography image volumes for head motion in routine positron emission tomography imaging to improve quantification in motion impacted data.Materials and Methods: A neural network was trained with 3D positron emission tomography image volumes in an unsupervised manner to predict transformation parameters required to perform image registration. A multi-step convolutional neural network (CNN) was combined with a spatial transform layer. Pairs of target and source images were used as input to the network. To prepare the training dataset, a previously published TOF-PEPT algorithm was applied to automatically detect static frames where the patient remained in a relatively steady position and transitional frames where they underwent abrupt motion. A single image volume was reconstructed for each static frame. The image reconstructed from the first static frame served as the target image with images from subsequent static frames being used as source images. The trained neural network predicted transformation parameters that could be used to perform frame-by-frame image-based motion correction but also enabled raw listmode positron emission tomography data correction where individual line-of-responses were repositioned. Line profiles and ROIs were drawn across the reconstructed image volumes to compare performance and quantitative results between standard registration tools and the deep learning technique. Corrected volumes were further compared to motion free images quantitatively using Dice indices.Results: In total, one hundred 3D positron emission tomography image volumes were used to train the network. Cross-validation was carried out using a 4:1 ratio for the training and test data. A conventional algorithm for affine registration from the Advanced Normalization Tools (ANTs) software package served as a baseline. To evaluate the correction performance, the mean Dice index and standardized uptake value (SUV) were used. Application of the algorithm to clinical data showed good performance with respect to registration accuracy as well as processing time. The neural network yielded a mean Dice index of ∼0.87 which was similar to the advanced Normalization Tools algorithm and did so ∼3x faster using a multi-core CPU and ∼20x faster with a GPU. Standardized uptake value analysis showed that quantitative results were 30%–60% higher in the motion-corrected images, and the neural network performed better than or close to the advanced Normalization Tools.Conclusion: The aim of this work was to study the quantitative impact of using a data-driven deep learning motion correction technique for positron emission tomography data and assess its performance. The results showed the technique is capable of producing high quality registrations that compensate for patient motion that occurs during a scan and improve quantitative accuracy.

Cross-Attention for Improved Motion Correction in Brain PET.

Head movement during long scan sessions degrades the quality of reconstruction in positron emission tomography (PET) and introduces artifacts, which limits clinical diagnosis and treatment. Recent deep learning-based motion correction work utilized raw PET list-mode data and hardware motion tracking (HMT) to learn head motion in a supervised manner. However, motion prediction results were not robust to testing subjects outside the training data domain. In this paper, we integrate a cross-attention mechanism into the supervised deep learning network to improve motion correction across test subjects. Specifically, cross-attention learns the spatial correspondence between the reference images and moving images to explicitly focus the model on the most correlative inherent information - the head region the motion correction. We validate our approach on brain PET data from two different scanners: HRRT without time of flight (ToF) and mCT with ToF. Compared with traditional and deep learning benchmarks, our network improved the performance of motion correction by 58% and 26% in translation and rotation, respectively, in multi-subject testing in HRRT studies. In mCT studies, our approach improved performance by 66% and 64% for translation and rotation, respectively. Our results demonstrate that cross-attention has the potential to improve the quality of brain PET image reconstruction without the dependence on HMT. All code will be released on GitHub: https://github.com/OnofreyLab/dl_hmc_attention_mlcn2023.

A back-projection-and-filtering-like (BPF-like) reconstruction method with the deep learning filtration from listmode data in TOF-PET.

The time-of-flight (TOF) information improves signal-to-noise ratio (SNR) for positron emission tomography (PET) imaging. Existing analytical algorithms for TOF PET usually follow a filtered back-projection process on reconstructing images from the sinogram data. This work aims to develop a back-projection-and-filtering-like (BPF-like) algorithm that reconstructs the TOF PET image directly from listmode data rapidly. We extended the 2D conventional non-TOF PET projection model to a TOF case, where projection data are represented as line integrals weighted by the one-dimensional TOF kernel along the projection direction. After deriving the central slice theorem and the TOF back-projection of listmode data, we designed a deep learning network with a modified U-net architecture to perform the spatial filtration (reconstruction filter). The proposed BP-Net method was validated via Monte Carlo simulations of TOF PET listmode data with three different time resolutions for two types of activity phantoms. The network was only trained on the simulated full-dose XCAT dataset and then evaluated on XCAT and Jaszczak data with different time resolutions and dose levels. Reconstructed images show that when compared with the conventional BPF algorithm and the MLEM algorithm proposed for TOF PET, the proposed BP-Net method obtains better image quality in terms of peak signal-to-noise ratio, relative mean square error, and structure similarity index; besides, the reconstruction speed of the BP-Net is 1.75 times faster than BPF and 29.05 times faster than MLEM using 15 iterations. The results also indicate that the performance of the BP-Net degrades with worse time resolutions and lower tracer doses, but degrades less than BPF or MLEM reconstructions. In this work, we developed an analytical-like reconstruction in the form of BPF with the reconstruction filtering operation performed via a deep network. The method runs even faster than the conventional BPF algorithm and provides accurate reconstructions from listmode data in TOF-PET, free of rebinning data to a sinogram.

Fast Energy Dependent Scatter Correction for List-Mode PET Data.

Improvements in energy resolution of modern positron emission tomography (PET) detectors have created opportunities to implement energy-based scatter correction algorithms. Here, we use the energy information of auxiliary windows to estimate the scatter component. Our method is directly implemented in an iterative reconstruction algorithm, generating a scatter-corrected image without the need for sinograms. The purpose was to implement a fast energy-based scatter correction method on list-mode PET data, when it was not possible to use an attenuation map as a practical approach for the scatter degradation. The proposed method was evaluated using Monte Carlo simulations of various digital phantoms. It accurately estimated the scatter fraction distribution, and improved the image contrast in the simulated studied cases. We conclude that the proposed scatter correction method could effectively correct the scattered events, including multiple scatters and those originated in sources outside the field of view.

Feasibility of equivalent performance of 3D TOF [18F]-FDG PET/CT with reduced acquisition time using clinical and semiquantitative parameters

BackgroundHigh-performance time-of-flight (TOF) positron emission tomography (PET) systems have the capability for rapid data acquisition while preserving diagnostic image quality. However, determining a reliable and clinically applicable cut-off of the acquisition time plays an important role in routine practice. This study aimed to assess the diagnostic equivalence of short acquisition time of 57 with routine 75 seconds per bed position (s/BP) of [18F]-fluoro-deoxy-glucose (FDG) PET.Phantom studies applying EARL criteria suggested the feasibility of shortened acquisition time in routine clinical imaging by 3D TOF PET/CT scanners. Ninety-six patients with melanoma, lung or head and neck cancer underwent a standard whole-body, skull base-to-thigh or vertex-to-thigh [18F]-FDG PET/CT examination using the 3D TOF Ingenuity TF PET/CT system (Philips, Cleveland, OH). The [18F]-FDG activity applied was equal to 4MBq per kg body weight. Retrospectively, PET list-mode data were used to calculate a second PET study per patient with a reduced acquisition time of 57 s instead of routine 75 s/BP. PET/CT data were reconstructed using a 3D OSEM TOF algorithm. Blinded patient data were analysed by two nuclear medicine physicians. The number of [18F]-FDG-avid lesions per body region (head&neck, thorax, abdomen, bone, extremity) and image quality (grade 1–5) were evaluated. Semiquantitative analyses were performed by standardized uptake value (SUV) measurements using 3D volume of interests (VOI). The visual and semiquantitative diagnostic equivalence of 214 [18F]-FDG-avid lesions were analysed in the routine standard (75 s/BP) as well as the calculated PET/CT studies with short acquisition time. Statistical analyses were performed by equivalence testing and Bland–Altman plots.ResultsLesion detection rate per patient’s body region agreed in > 98% comparing 57 s/BP and 75 s/BP datasets. Overall image quality was determined as equal or superior to 75 s in 80% and 69%, respectively. In the semiquantitative lesion-based analyses, a significant equivalence was found between the 75 s/BP and 57 s/BP PET/CT images both for SUVmax (p = 0.004) and SUVmean (p = 0.003).ConclusionThe results of this study demonstrate significant clinical and semiquantitative equivalence between short acquisition time of 57 s/BP and standard 75 s/BP 3D TOF [18F]-FDG PET/CT scanning, which may improve the patient’s workflow in routine practice.

A data-driven respiratory motion estimation approach for PET based on time-of-flight weighted positron emission particle tracking.

PurposeRespiratory motion of patients during positron emission tomography (PET)/computed tomography (CT) imaging affects both image quality and quantitative accuracy. Hardware‐based motion estimation, which is the current clinical standard, requires initial setup, maintenance, and calibration of the equipment, and can be associated with patient discomfort. Data‐driven techniques are an active area of research with limited exploration into lesion‐specific motion estimation. This paper introduces a time‐of‐flight (TOF)‐weighted positron emission particle tracking (PEPT) algorithm that facilitates lesion‐specific respiratory motion estimation from raw listmode PET data.MethodsThe TOF‐PEPT algorithm was implemented and investigated under different scenarios: (a) a phantom study with a point source and an Anzai band for respiratory motion tracking; (b) a phantom study with a point source only, no Anzai band; (c) two clinical studies with point sources and the Anzai band; (d) two clinical studies with point sources only, no Anzai band; and (e) two clinical studies using lesions/internal regions instead of point sources and no Anzai band. For studies with radioactive point sources, they were placed on patients during PET/CT imaging. The motion tracking was performed using a preselected region of interest (ROI), manually drawn around point sources or lesions on reconstructed images. The extracted motion signals were compared with the Anzai band when applicable. For the purposes of additional comparison, a center‐of‐mass (COM) algorithm was implemented both with and without the use of TOF information. Using the motion estimate from each method, amplitude‐based gating was applied, and gated images were reconstructed.ResultsThe TOF‐PEPT algorithm is shown to successfully determine the respiratory motion for both phantom and clinical studies. The derived motion signals correlated well with the Anzai band; correlation coefficients of 0.99 and 0.94‐0.97 were obtained for the phantom study and the clinical studies, respectively. TOF‐PEPT was found to be 13–38% better correlated with the Anzai results than the COM methods. Maximum Standardized Uptake Values (SUVs) were used to quantitatively compare the reconstructed‐gated images. In comparison with the ungated image, a 14–39% increase in the max SUV across several lesion areas and an 8.7% increase in the max SUV on the tracked lesion area were observed in the gated images based on TOF‐PEPT. The distinct presence of lesions with reduced blurring effect and generally sharper images were readily apparent in all clinical studies. In addition, max SUVs were found to be 4–10% higher in the TOF‐PEPT‐based gated images than in those based on Anzai and COM methods.ConclusionA PEPT‐ based algorithm has been presented for determining movement due to respiratory motion during PET/CT imaging. Gating based on the motion estimate is shown to quantifiably improve the image quality in both a controlled point source phantom study and in clinical data patient studies. The algorithm has the potential to facilitate true motion correction where the reconstruction algorithm can use all data available.

Data-driven respiratory gating based on localized diaphragm sensing in TOF PET

It is important to measure the respiratory cycle in positron emission tomography (PET) to enhance the contrast of the tumor as well as the accuracy of its localization in organs such as the lung and liver. Several types of data-driven respiratory gating methods, such as center of mass and principal component analysis, have been developed to directly measure the breathing cycle from PET images and listmode data. However, the breathing cycle is still hard to detect in low signal-to-noise ratio (SNR) data, particularly in low dose PET/CT scans. To address this issue, a time-of-flight (TOF) PET is currently utilized for the data-driven respiratory gating because of its higher SNR and better localization of the region of interest. To further improve the accuracy of respiratory gating with TOF information, we propose an accurate data-driven respiratory gating method, which retrospectively derives the respiratory signal using a localized sensing method based on a diaphragm mask in TOF PET data. To assess the accuracy of the proposed method, the performance is evaluated with three patient datasets, and a pressure-belt signal as the ground truth is compared. In our experiments, we validate that the respiratory signal using the proposed data-driven gating method is well matched to the pressure-belt respiratory signal with less than 5% peak time errors and over 80% trace correlations. Based on gated signals, the respiratory-gated image of the proposed method provides more clear edges of organs compared to images using conventional non-TOF methods. Therefore, we demonstrate that the proposed method can achieve improvements for the accuracy of gating signals and image quality.

Body motion detection and correction in cardiac PET: Phantom and human studies.

Patient body motion during a cardiac positron emission tomography (PET) scan can severely degrade image quality. We propose and evaluate a novel method to detect, estimate, and correct body motion in cardiac PET. Our method consists of three key components: motion detection, motion estimation, and motion-compensated image reconstruction. For motion detection, we first divide PET list-mode data into 1-s bins and compute the center of mass (COM) of the coincidences' distribution in each bin. We then compute the covariance matrix within a 25-s sliding window over the COM signals inside the window. The sum of the eigenvalues of the covariance matrix is used to separate the list-mode data into "static" (i.e., body motion free) and "moving" (i.e. contaminated by body motion) frames. Each moving frame is further divided into a number of evenly spaced sub-frames (referred to as "sub-moving" frames), in which motion is assumed to be negligible. For motion estimation, we first reconstruct the data in each static and sub-moving frame using a rapid back-projection technique. We then select the longest static frame as the reference frame and estimate elastic motion transformations to the reference frame from all other static and sub-moving frames using nonrigid registration. For motion-compensated image reconstruction, we reconstruct all the list-mode data into a single image volume in the reference frame by incorporating the estimated motion transformations in the PET system matrix. We evaluated the performance of our approach in both phantom and human studies. Visually, the motion-corrected (MC) PET images obtained using the proposed method have better quality and fewer motion artifacts than the images reconstructed without motion correction (NMC). Quantitative analysis indicates that MC yields higher myocardium to blood pool concentration ratios. MC also yields sharper myocardium than NMC. The proposed body motion correction method improves image quality of cardiac PET.

Ultra-low-dose PET reconstruction using generative adversarial network with feature matching and task-specific perceptual loss.

Our goal was to use a generative adversarial network (GAN) with feature matching and task-specific perceptual loss to synthesize standard-dose amyloid Positron emission tomography (PET) images of high quality and including accurate pathological features from ultra-low-dose PET images only. Forty PET datasets from 39 participants were acquired with a simultaneous PET/MRI scanner following injection of 330±30MBq of the amyloid radiotracer 18F-florbetaben. The raw list-mode PET data were reconstructed as the standard-dose ground truth and were randomly undersampled by a factor of 100 to reconstruct 1% low-dose PET scans. A 2D encoder-decoder network was implemented as the generator to synthesize a standard-dose image and a discriminator was used to evaluate them. The two networks contested with each other to achieve high-visual quality PET from the ultra-low-dose PET. Multi-slice inputs were used to reduce noise by providing the network with 2.5D information. Feature matching was applied to reduce hallucinated structures. Task-specific perceptual loss was designed to maintain the correct pathological features. The image quality was evaluated by peak signal-to-noise ratio (PSNR), structural similarity (SSIM), and root mean square error (RMSE) metrics with and without each of these modules. Two expert radiologists were asked to score image quality on a 5-point scale and identified the amyloid status (positive or negative). With only low-dose PET as input, the proposed method significantly outperformed Chen et al.'s method (Chen et al. Radiology. 2018;290:649-656) (which shows the best performance in this task) with the same input (PET-only model) by 1.87dB in PSNR, 2.04% in SSIM, and 24.75% in RMSE. It also achieved comparable results to Chen et al.'s method which used additional magnetic resonance imaging (MRI) inputs (PET-MR model). Experts' reading results showed that the proposed method could achieve better overall image quality and maintain better pathological features indicating amyloid status than both PET-only and PET-MR models proposed by Chen et al. CONCLUSION: Standard-dose amyloid PET images can be synthesized from ultra-low-dose images using GAN. Applying adversarial learning, feature matching, and task-specific perceptual loss are essential to ensure image quality and the preservation of pathological features.

Self-Gating: An Adaptive Center-of-Mass Approach for Respiratory Gating in PET.

The goal is to develop an adaptive center-of-mass (COM)-based approach for device-less respiratory gating of list-mode positron emission tomography (PET) data. Our method contains two steps. The first is to automatically extract an optimized respiratory motion signal from the list-mode data during acquisition. The respiratory motion signal was calculated by tracking the location of COM within a volume of interest (VOI). The signal prominence (SP) was calculated based on Fourier analysis of the signal. The VOI was adaptively optimized to maximize SP. The second step is to automatically correct signal-flipping effects. The sign of the signal was determined based on the assumption that the average patient spends more time during expiration than inspiration. To validate our methods, thirty-one 18F-FDG patient scans were included in this paper. An external device-based signal was used as the gold standard, and the correlation coefficient of the data-driven signal with the device-based signal was measured. Our method successfully extracted respiratory signal from 30 out of 31 datasets. The failure case was due to lack of uptake in the field of view. Moreover, our sign determination method obtained correct results for all scans excluding the failure case. Quantitatively, the proposed signal extraction approach achieved a median correlation of 0.85 with the device-based signal. Gated images using optimized data-driven signal showed improved lesion contrast over static image and were comparable to those using device-based signal. We presented a new data-driven method to automatically extract respiratory motion signal from list-mode PET data by optimizing VOI for COM calculation, as well as determine motion direction from signal asymmetry. Successful application of the proposed method on most clinical datasets and comparison with device-based signal suggests its potential of serving as an alternative to external respiratory monitors.

Quantitative and Visual Assessments toward Potential Sub-mSv or Ultrafast FDG PET Using High-Sensitivity TOF PET in PET/MRI.

Newer high-performance time-of-flight (TOF) positron emission tomography (PET) systems have the capability to preserve diagnostic image quality with low count density, while maintaining a high raw photon detection sensitivity that would allow for a reduction in injected dose or rapid data acquisition. To assess this, we performed quantitative and visual assessments of the PET images acquired using a highly sensitive (23.3cps/kBq) large field of view (25-cm axial) silicon photomultiplier (SiPM)-based TOF PET (400-ps timing resolution) integrated with 3T-MRI in comparison to PET images acquired on non-TOF PET/x-ray computed tomography (CT) systems. Whole-body 2-deoxy-2-[18F]fluoro-D-glucose ([18F]FDG) PET/CT was acquired for 15 patients followed by whole body PET/magnetic resonance imaging (MRI) with an average injected dose of 325±84MBq. The PET list mode data from PET/MRI were reconstructed using full datasets (4min/bed) and reduced datasets (2, 1, 0.5, and 0.25min/bed). Qualitative assessment between PET/CT and PET/MR images were made. A Likert-type scale between 1 and 5, 1 for non-diagnostic, 3 equivalent to PET/CT, and 5 superior quality, was used. Maximum and mean standardized uptake values (SUVmax and SUVmean) of normal tissues and lesions detected were measured and compared. Mean visual assessment scores were 3.54±0.32, 3.62±0.38, and 3.69±0.35 for the brain and 3.05±0.49, 3.71±0.45, and 4.14±0.44 for the whole-body maximum intensity projections (MIPs) for 1, 2, and 4min/bed PET/MR images, respectively. The SUVmean values for normal tissues were lower and statistically significant for images acquired at 4, 2, 1, 0.5, and 0.25min/bed on the PET/MR, with values of -18±28% (p < 0.001), - 16±29% (p = 0.001), - 16±31% (p = 0.002), - 14±35% (p < 0.001), and - 13±34% (p = 0.002), respectively. SUVmax and SUVpeak values of all lesions were higher and statistically significant (p < 0.05) for 4, 2, 1, 0.50, and 0.25min/bed PET/MR datasets. High-sensitivity TOF PET showed comparable but still better visual image quality even at a much reduced activity in comparison to lower-sensitivity non-TOF PET. Our data translates to a seven times reduction in either injection dose for the same time or total scan time for the same injected dose. This "ultra-sensitivity" PET system provides a path to clinically acceptable extremely low-dose FDG PET studies (e.g., sub 1mCi injection or sub-mSv effective dose) or PET studies as short as 1min/bed (e.g., 6min of total scan time) to cover whole body without compromising diagnostic performance.

Regularization parameter selection for penalized-likelihood list-mode image reconstruction in PET

Penalized likelihood (PL) reconstruction has demonstrated potential to improve image quality of positron emission tomography (PET) over unregularized ordered-subsets expectation-maximization (OSEM) algorithm. However, selecting proper regularization parameters in PL reconstruction has been challenging due to the lack of ground truth and variation of penalty functions. Here we present a method to choose regularization parameters using a cross-validation log-likelihood (CVLL) function. This new method does not require any knowledge of the true image and is directly applicable to list-mode PET data. We performed statistical analysis of the mean and variance of the CVLL. The results show that the CVLL provides an unbiased estimate of the log-likelihood function calculated using the noise free data. The predicted variance can be used to verify the statistical significance of the difference between CVLL values. The proposed method was validated using simulation studies and also applied to real patient data. The reconstructed images using optimum parameters selected by the proposed method show good image quality visually.

Super-Resolution Reconstruction of 3D PET Images Using Two Respiratory-Phase Low-Dose CT Images

Positron emission tomography (PET) images usually suffer from incorrect quantification of the radioactive uptake of small lesions due to low spatial resolution. To improve the spatial resolution, we previously proposed super-resolution (SR) algorithms based on wobble scanning. The proposed algorithms, however, require mechanical motion of the patient bed or a system gantry for wobble scanning. In this paper, we propose a framework for SR reconstruction of 3D PET images based on the use of respiratory motion rather than wobble motion. As in the conventional protocol of PET imaging, gated list-mode PET data are acquired in a free breathing condition. In addition, we acquire two low-dose CT images in a breath-hold manner at exhale and inhale phases, without increasing the radiation burden to a patient. Using the two low-dose CT images, we estimate the 4D motion vector field (MVF) and correspondingly generate a virtual 4D CT image that are matched to the 4D PET image. The 3D CT images have much better spatial resolution than PET images and therefore the corresponding estimated 3D MVFs can be considered reliable for PET SR reconstruction. We then estimate space-variant point spread functions (PSFs) in the imaging field of view using a minimum number of PSFs obtained through Monte-Carlo simulations. Finally, SR reconstruction is performed by incorporating the estimated 3D MVFs and space-variant PSFs. In the SR reconstruction, to avoid the resolution degradation in the evenly spaced parallel-beam rebinning and to reduce the computational time on the graphics processing unit, we introduce a parallel-friendly spanned line of response reconstruction technique based on fan-beam reordering. The proposed framework is evaluated via Monte-Carlo simulations with the digital XCAT phantom and via experiments with several patient datasets including moving lung lesions. Both the simulation and experiment results show that the proposed framework provides much clearer organ boundaries as well as more accurate quantitative lesion information than the conventional methods, with a considerable reduction of computational time.

Repurposing the Microsoft Kinect for Windows v2 for external head motion tracking for brain PET

Medical imaging systems such as those used in positron emission tomography (PET) are capable of spatial resolutions that enable the imaging of small, functionally important brain structures. However, the quality of data from PET brain studies is often limited by subject motion during acquisition. This is particularly challenging for patients with neurological disorders or with dynamic research studies that can last 90 min or more. Restraining head movement during the scan does not eliminate motion entirely and can be unpleasant for the subject. Head motion can be detected and measured using a variety of techniques that either use the PET data itself or an external tracking system. Advances in computer vision arising from the video gaming industry could offer significant benefits when re-purposed for medical applications. A method for measuring rigid body type head motion using the Microsoft Kinect v2 is described with results presenting ⩽0.5 mm spatial accuracy. Motion data is measured in real-time at 30 Hz using the KinectFusion algorithm. Non-rigid motion is detected using the residual alignment energy data of the KinectFusion algorithm allowing for unreliable motion to be discarded. Motion data is aligned to PET listmode data using injected pulse sequences into the PET/CT gantry allowing for correction of rigid body motion. Pilot data from a clinical dynamic PET/CT examination is shown.

Sensitivity estimation in time-of-flight list-mode positron emission tomography.

An accurate quantification of the images in positron emission tomography (PET) requires knowing the actual sensitivity at each voxel, which represents the probability that a positron emitted in that voxel is finally detected as a coincidence of two gamma rays in a pair of detectors in the PET scanner. This sensitivity depends on the characteristics of the acquisition, as it is affected by the attenuation of the annihilation gamma rays in the body, and possible variations of the sensitivity of the scanner detectors. In this work, the authors propose a new approach to handle time-of-flight (TOF) list-mode PET data, which allows performing either or both, a self-attenuation correction, and self-normalization correction based on emission data only. The authors derive the theory using a fully Bayesian statistical model of complete data. The authors perform an initial evaluation of algorithms derived from that theory and proposed in this work using numerical 2D list-mode simulations with different TOF resolutions and total number of detected coincidences. Effects of randoms and scatter are not simulated. The authors found that proposed algorithms successfully correct for unknown attenuation and scanner normalization for simulated 2D list-mode TOF-PET data. A new method is presented that can be used for corrections for attenuation and normalization (sensitivity) using TOF list-mode data.

Automatic Cardiac Self-Gating of Small-Animal PET Data.

The cardiac gating signal (CGS) in positron emission tomography (PET) studies is usually obtained from an electrocardiography (ECG) monitor. In this work, we propose a method to obtain the CGS in small-animal PET using the acquired list-mode data without using any hardware or end-user input. The CGS was obtained from the number of coincidences over time acquired in the lines-of-response connected with the cardiac region. This region is identified in the image as its value changes with frequencies in the range of 3 to 12 Hz. The procedure was tested in a study with 29 Wistar rats and 6 mice injected with 2-deoxy-2-[(18)F]fluoro-D-glucose, which underwent a 45-min single-bed list-mode PET scan of the heart syncronized with an ECG. The estimated signals and the reconstructed images using eight-gated frames were compared with the ones obtained using the ECG signal from the monitor. The differences of the PET-based CGS with respect to the ECG relative to the duration of the heartbeats were 5.6 % in rats and 11.0 % in mice. The reconstructed gated images obtained from the proposed method do not differ qualitatively with respect to the ones obtained with the ECG. The quantitative analysis of both set of images were performed measuring the volume of the left ventricle (LV) of the rats in the end-of-systole and end-of-diastole phase. The differences found in these parameters between both methods were below 12.1 % in the ESV and 9.3 % in the EDV with a 95 % confidence interval, which are comparable to the accuracy (7 %) of the method used for segmenting the LV. The proposed method is able to provide a valid and accurate CGS in small-animal PET list-mode data.

Single-cell tracking with PET using a novel trajectory reconstruction algorithm.

Virtually all biomedical applications of positron emission tomography (PET) use images to represent the distribution of a radiotracer. However, PET is increasingly used in cell tracking applications, for which the "imaging" paradigm may not be optimal. Here, we investigate an alternative approach, which consists in reconstructing the time-varying position of individual radiolabeled cells directly from PET measurements. As a proof of concept, we formulate a new algorithm for reconstructing the trajectory of one single moving cell directly from list-mode PET data. We model the trajectory as a 3-D B-spline function of the temporal variable and use nonlinear optimization to minimize the mean-square distance between the trajectory and the recorded list-mode coincidence events. Using Monte Carlo simulations (GATE), we show that this new algorithm can track a single source moving within a small-animal PET system with 3 mm accuracy provided that the activity of the cell [Bq] is greater than four times its velocity [mm/s]. The algorithm outperforms conventional ML-EM as well as the "minimum distance" method used for positron emission particle tracking (PEPT). The new method was also successfully validated using experimentally acquired PET data. In conclusion, we demonstrated the feasibility of a new method for tracking a single moving cell directly from PET list-mode data, at the whole-body level, for physiologically relevant activities and velocities.