Four-Dimensional Image Reconstruction Strategies in Cardiac-Gated and Respiratory-Gated PET Imaging

Spatio-Temporal Positron Emission Tomography Reconstruction with Attenuation and Motion Correction

The detection of cancer lesions of a comparable size to that of the typical system resolution of modern scanners is a long-standing problem in Positron Emission Tomography. In this paper, the effect of composing an image-registering convolutional neural network with the modeling of the static data acquisition (i.e., the forward model) is investigated. Two algorithms for Positron Emission Tomography reconstruction with motion and attenuation correction are proposed and their performance is evaluated in the detectability of small pulmonary lesions. The evaluation is performed on synthetic data with respect to chosen figures of merit, visual inspection, and an ideal observer. The commonly used figures of merit-Peak Signal-to-Noise Ratio, Recovery Coefficient, and Signal Difference-to-Noise Ration-give inconclusive responses, whereas visual inspection and the Channelised Hotelling Observer suggest that the proposed algorithms outperform current clinical practice.

Monte Carlo Simulation and Reconstruction: Assessment of Myocardial Perfusion Imaging of Tracer Dynamics With Cardiac Motion Due to Deformation and Respiration Using Gamma Camera With Continuous Acquisition

Purpose

Myocardial perfusion imaging (MPI) with single photon emission computed tomography (SPECT) is routinely used for stress testing in nuclear medicine. Recently, our group extended its potential going from 3D visual qualitative image analysis to 4D spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration and the estimated myocardial blood flow (MBF) and coronary flow reserve (CFR). However, the quality of reconstructed image is compromised due to cardiac deformation and respiration. The work presented here develops an algorithm that reconstructs the dynamic sequence of separate respiratory and cardiac phases and evaluates the algorithm with data simulated with a Monte Carlo simulation for the continuous image acquisition and processing with a slowly rotating SPECT camera.

Methods

A clinically realistic Monte Carlo (MC) simulation is developed using the 4D Extended Cardiac Torso (XCAT) digital phantom with respiratory and cardiac motion to model continuous data acquisition of dynamic cardiac SPECT with slowly rotating gamma cameras by incorporating deformation and displacement of the myocardium due to cardiac and respiratory motion. We extended our previously developed 4D maximum-likelihood expectation-maximization (MLEM) reconstruction algorithm for a data set binned from a continuous list mode (LM) simulation with cardiac and respiratory information. Our spatiotemporal image reconstruction uses splines to explicitly model the temporal change of the tracer for each cardiac and respiratory gate that delineates the myocardial spatial position as the tracer washes in and out. Unlike in a fully list-mode data acquisition and reconstruction the accumulated photons are binned over a specific but very short time interval corresponding to each cardiac and respiratory gate. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it continuously deforms. These results are then compared with the conventional 4D spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. Mean Stabilized Activity (MSA), signal to noise ratio (SNR) and Bias for the myocardium activities for three different target-to-background ratios (TBRs) are evaluated. Dynamic quantitative indices such as wash-in (K₁) and wash-out (k₂) rates at each gate were also estimated.

Results

The MSA and SNR are higher with higher TBRs while biases were improved with higher TBRs to less than 10%. The correlation between exhalation-inhalation sequence with the ground truth during respiratory cycle was excellent. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Estimated values of K₁ and k₂ were also consistent with the ground truth.

Conclusion

The continuous image acquisition for dynamic scan using conventional two-head gamma cameras can provide valuable information for MPI. Our study demonstrated the viability of using a continuous image acquisition method on a widely used clinical two-head SPECT system. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Precise implementation of reconstruction algorithms, better segmentation techniques by generating images of different tissue types and background activity would improve the feasibility of the method in real clinical environment.

Motion in nuclear cardiology imaging: types, artifacts, detection and correction techniques

In this paper, the authors review the field of motion detection and correction in nuclear cardiology with single photon emission computed tomography (SPECT) and positron emission tomography (PET) imaging systems. We start with a brief overview of nuclear cardiology applications and description of SPECT and PET imaging systems, then explaining the different types of motion and their related artefacts. Moreover, we classify and describe various techniques for motion detection and correction, discussing their potential advantages including reference to metrics and tasks, particularly towards improvements in image quality and diagnostic performance. In addition, we emphasize limitations encountered in different motion detection and correction methods that may challenge routine clinical applications and diagnostic performance.

Motion correction for PET data using subspace-based real-time MR imaging in simultaneous PET/MR

Image quality of positron emission tomography (PET) reconstructions is degraded by subject motion occurring during the acquisition. Magnetic resonance (MR)-based motion correction approaches have been studied for PET/MR scanners and have been successful at capturing regular motion patterns, when used in conjunction with surrogate signals (e.g. navigators) to detect motion. However, handling irregular respiratory motion and bulk motion remains challenging. In this work, we propose an MR-based motion correction method relying on subspace-based real-time MR imaging to estimate motion fields used to correct PET reconstructions. We take advantage of the low-rank characteristics of dynamic MR images to reconstruct high-resolution MR images at high frame rates from highly undersampled k-space data. Reconstructed dynamic MR images are used to determine motion phases for PET reconstruction and estimate phase-to-phase nonrigid motion fields able to capture complex motion patterns such as irregular respiratory and bulk motion. MR-derived binning and motion fields are used for PET reconstruction to generate motion-corrected PET images. The proposed method was evaluated on in vivo data with irregular motion patterns. MR reconstructions accurately captured motion, outperforming state-of-the-art dynamic MR reconstruction techniques. Evaluation of PET reconstructions demonstrated the benefits of the proposed method in terms of motion artifacts reduction, improving the contrast-to-noise ratio by up to a factor 3 and achieveing a target-to-background ratio up to 90% superior compared to standard/uncorrected methods. The proposed method can improve the image quality of motion-corrected PET reconstructions in clinical applications.

Partial volume and motion correction in cardiac PET: First results from an in vs ex vivo comparison using animal datasets

BACKGROUND

In a previous study on ex vivo, static cardiac datasets, we investigated the benefits of performing partial volume correction (PVC) in cardiac 18F-Fluorodeoxyglucose(FDG) PET datasets. In the present study, we extend the analysis to in vivo cardiac datasets, with the aim of defining which reconstruction technique maximizes quantitative accuracy and, ultimately, makes PET a better diagnostic tool for cardiac pathologies.

METHODS

In vivo sheep datasets were acquired and reconstructed with/without motion correction and using several reconstruction algorithms (with/without resolution modeling, with/without non-anatomical priors). Corresponding ex vivo scans of the excised sheep hearts were performed on a small-animal PET scanner (Siemens Focus 220, microPET) to provide high-resolution reference data unaffected by respiratory and cardiac motion. A comparison between the in vivo cardiac reconstructions and the corresponding ex vivo ground truth was performed.

RESULTS

The use of an edge-preserving prior (Total Variation (TV) prior in this work) in combination with motion correction reduces the bias in absolute quantification when compared to the standard clinical reconstructions (- 0.83 vs - 3.74 SUV units), when the end-systolic gate is considered. At end-diastole, motion correction improves absolute quantification but the PVC with priors does not improve the similarity to the ground truth more than a regular iterative reconstruction with motion correction and without priors. Relative quantification was not influenced much by the chosen reconstruction algorithm.

CONCLUSIONS

The relative ranking of the algorithms suggests superiority of the PVC reconstructions with dual gating in terms of overall absolute quantification and noise properties. A well-tuned edge-preserving prior, such as TV, enhances the noise properties of the resulting images of the heart. The end-systolic gate yields the most accurate quantification of cardiac datasets.

MR-based cardiac and respiratory motion correction of PET: application to static and dynamic cardiac 18F-FDG imaging

Motion of the myocardium deteriorates the quality and quantitative accuracy of cardiac PET images. We present a method for MR-based cardiac and respiratory motion correction of cardiac PET data and evaluate its impact on estimation of activity and kinetic parameters in human subjects. Three healthy subjects underwent simultaneous dynamic 18F-FDG PET and MRI on a hybrid PET/MR scanner. A cardiorespiratory motion field was determined for each subject using navigator, tagging and golden-angle radial MR acquisitions. Acquired coincidence events were binned into cardiac and respiratory phases using electrocardiogram and list mode-driven signals, respectively. Dynamic PET images were reconstructed with MR-based motion correction (MC) and without motion correction (NMC). Parametric images of 18F-FDG consumption rates (Ki) were estimated using Patlak's method for both MC and NMC images. MC alleviated motion artifacts in PET images, resulting in improved spatial resolution, improved recovery of activity in the myocardium wall and reduced spillover from the myocardium to the left ventricle cavity. Significantly higher myocardium contrast-to-noise ratio and lower apparent wall thickness were obtained in MC versus NMC images. Likewise, parametric images of Ki calculated with MC data had improved spatial resolution as compared to those obtained with NMC. Consistent with an increase in reconstructed activity concentration in the frames used during kinetic analyses, MC led to the estimation of higher Ki values almost everywhere in the myocardium, with up to 18% increase (mean across subjects) in the septum as compared to NMC. This study shows that MR-based motion correction of cardiac PET results in improved image quality that can benefit both static and dynamic studies.

Cardiac contraction motion compensation in gated myocardial perfusion SPECT: A comparative study

INTRODUCTION

Cardiac contraction significantly degrades quality and quantitative accuracy of gated myocardial perfusion SPECT (MPS) images. In this study, we aimed to explore different techniques in motion-compensated temporal processing of MPS images and their impact on image quality and quantitative accuracy.

MATERIAL AND METHOD

50 patients without known heart condition underwent gated MPS. 3D motion compensation methods using Motion Freezing by Cedars Sinai (MF), Log-domain Diffeomorphic Demons (LDD) and Free-Form Deformation (FFD) were applied to warp all image phases to fit the end-diastolic (ED) phase. Afterwards, myocardial wall thickness, myocardial to blood pool contrast, and image contrast-to noise ratio (CNR) were measured in summed images with no motion compensation (NoMC) and compensated images (MF, LDD and FFD). Total Perfusion Defect (TPD) was derived from Cedars-Sinai software, on the basis of sex-specific normal limits.

RESULT

Left ventricle (LV) lateral wall thickness was reduced after applying motion compensation (p < 0.05). Myocardial to blood pool contrast and CNR in compensated images were greater than NoMC (p < 0.05). TPD_LDD was in good agreement with the corresponding TPD_MF (p = 0.13).

CONCLUSION

All methods have improved image quality and quantitative performance relative to NoMC. LDD and FFD are fully automatic and do not require any manual intervention, while MF is dependent on contour definition. In terms of diagnostic parameters LDD is in good agreement with MF which is a clinically accepted method. Further investigation along with diagnostic reference standards, in order to specify diagnostic value of each technique is recommended.

Impact of motion compensation and partial volume correction for 18F-NaF PET/CT imaging of coronary plaque

Recent studies have suggested that ¹⁸F-NaF-PET enables visualization and quantification of plaque micro-calcification in the coronary tree. However, PET imaging of plaque calcification in the coronary arteries is challenging because of the respiratory and cardiac motion as well as partial volume effects. The objective of this work is to implement an image reconstruction framework, which incorporates compensation for respiratory as well as cardiac motion (MoCo) and partial volume correction (PVC), for cardiac ¹⁸F-NaF PET imaging in PET/CT. We evaluated the effect of MoCo and PVC on the quantification of vulnerable plaques in the coronary arteries. Realistic simulations (Biograph TPTV, Biograph mCT) and phantom acquisitions (Biograph mCT) were used for these evaluations. Different uptake values in the calcified plaques were evaluated in the simulations, while three 'plaque-type' lesions of 36, 31 and 18 mm³ were included in the phantom experiments. After validation, the MoCo and PVC methods were applied in four pilot NaF-PET patient studies. In all cases, the MoCo-based image reconstruction was performed using the STIR software. The PVC was obtained from a local projection (LP) method, previously evaluated in preclinical and clinical PET. The results obtained show a significant increase of the measured lesion-to-background ratios (LBR) in the MoCo  +  PVC images. These ratios were further enhanced when using directly the tissue-activities from the LP method, making this approach more suitable for the quantitative evaluation of coronary plaques. When using the LP method on the MoCo images, LBR increased between 200% and 1119% in the simulated data, between 212% and 614% in the phantom experiments and between 46% and 373% in the plaques with positive uptake observed in the pilot patients. In conclusion, we have built and validated a STIR framework incorporating MoCo and PVC for ¹⁸F-NaF PET imaging of coronary plaques. First results indicate an improved quantification of plaque-type lesions.

Enhancing ejection fraction measurement through 4D respiratory motion compensation in cardiac PET imaging

ECG gated cardiac PET imaging measures functional parameters such as left ventricle (LV) ejection fraction (EF), providing diagnostic and prognostic information for management of patients with coronary artery disease (CAD). Respiratory motion degrades spatial resolution and affects the accuracy in measuring the LV volumes for EF calculation. The goal of this study is to systematically investigate the effect of respiratory motion correction on the estimation of end-diastolic volume (EDV), end-systolic volume (ESV), and EF, especially on the separation of normal and abnormal EFs. We developed a respiratory motion incorporated 4D PET image reconstruction technique which uses all gated-frame data to acquire a motion-suppressed image. Using the standard XCAT phantom and two individual-specific volunteer XCAT phantoms, we simulated dual-gated myocardial perfusion imaging data for normally and abnormally beating hearts. With and without respiratory motion correction, we measured the EDV, ESV, and EF from the cardiac-gated reconstructed images. For all the phantoms, the estimated volumes increased and the biases significantly reduced with motion correction compared with those without. Furthermore, the improvement of ESV measurement in the abnormally beating heart led to better separation of normal and abnormal EFs. The simulation study demonstrated the significant effect of respiratory motion correction on cardiac imaging data with motion amplitude as small as 0.7 cm. The larger the motion amplitude the more improvement respiratory motion correction brought about on the EF measurement. Using data-driven respiratory gating, we also demonstrated the effect of respiratory motion correction on estimating the above functional parameters from list mode patient data. Respiratory motion correction has been shown to improve the accuracy of EF measurement in clinical cardiac PET imaging.

Attenuation correction in 4D-PET using a single-phase attenuation map and rigidity-adaptive deformable registration

PURPOSE

Four-dimensional positron emission tomography (4D-PET) imaging is a potential solution to the respiratory motion effect in the thoracic region. Computed tomography (CT)-based attenuation correction (AC) is an essential step toward quantitative imaging for PET. However, due to the temporal difference between 4D-PET and a single attenuation map from CT, typically available in routine clinical scanning, motion artifacts are observed in the attenuation-corrected PET images, leading to errors in tumor shape and uptake. We introduced a practical method to align single-phase CT with all other 4D-PET phases for AC.

METHODS

A penalized non-rigid Demons registration between individual 4D-PET frames without AC provides the motion vectors to be used for warping single-phase attenuation map. The non-rigid Demons registration was used to derive deformation vector fields (DVFs) between PET matched with the CT phase and other 4D-PET images. While attenuated PET images provide useful data for organ borders such as those of the lung and the liver, tumors cannot be distinguished from the background due to loss of contrast. To preserve the tumor shape in different phases, an ROI-covering tumor was excluded from nonrigid transformation. Instead the mean DVF of the central region of the tumor was assigned to all voxels in the ROI. This process mimics a rigid transformation of the tumor along with a nonrigid transformation of other organs. A 4D-XCAT phantom with spherical lung tumors, with diameters ranging from 10 to 40 mm, was used to evaluate the algorithm. The performance of the proposed hybrid method for attenuation map estimation was compared to (a) the Demons nonrigid registration only and (b) a single attenuation map based on quantitative parameters in individual PET frames.

RESULTS

Motion-related artifacts were significantly reduced in the attenuation-corrected 4D-PET images. When a single attenuation map was used for all individual PET frames, the normalized root-mean-square error (NRMSE) values in tumor region were 49.3% (STD: 8.3%), 50.5% (STD: 9.3%), 51.8% (STD: 10.8%) and 51.5% (STD: 12.1%) for 10-mm, 20-mm, 30-mm, and 40-mm tumors, respectively. These errors were reduced to 11.9% (STD: 2.9%), 13.6% (STD: 3.9%), 13.8% (STD: 4.8%), and 16.7% (STD: 9.3%) by our proposed method for deforming the attenuation map. The relative errors in total lesion glycolysis (TLG) values were -0.25% (STD: 2.87%) and 3.19% (STD: 2.35%) for 30-mm and 40-mm tumors, respectively, in proposed method. The corresponding values for Demons method were 25.22% (STD: 14.79%) and 18.42% (STD: 7.06%). Our proposed hybrid method outperforms the Demons method especially for larger tumors. For tumors smaller than 20 mm, nonrigid transformation could also provide quantitative results.

CONCLUSION

Although non-AC 4D-PET frames include insignificant anatomical information, they are still useful to estimate the DVFs to align the attenuation map for accurate AC. The proposed hybrid method can recover the AC-related artifacts and provide quantitative AC-PET images.

Effect of Bayesian-penalized likelihood reconstruction on [13N]-NH3 rest perfusion quantification

ABSTACT

OBJECTIVES: Myocardial blood flow (MBF) imaging is used in patients with suspected cardiac sarcoidosis, and also in stress/rest studies. The accuracy of MBF is dependent on imaging parameters such as new reconstruction methodologies. In this work, we aim to assess the impact of a novel PET reconstruction algorithm (Bayesian-penalized likelihood-BPL) on the values determined from the calculation of [13N]-NH3 MBF values.

METHODS

Data from 21 patients undergoing rest MBF evaluation [13N]-NH3 as part of sarcoidosis imaging were retrospectively analyzed. Each scan was reconstructed with a range of BPL coefficients (1-500), and standard clinical FBP and OSEM reconstructions. MBF values were calculated via an automated software routine for all datasets.

RESULTS

Reconstruction of [13N]-NH3 dynamic data using the BPL, OSEM, or FBP reconstruction showed no quantitative differences for the calculation of territorial or global MBF (P = .97). Image noise was lower using OSEM or BPL reconstructions than FBP and noise from BPL reached levels seen in OSEM images between B = 300 and B = 400. Intrasubject differences between all reconstructions over all patients in respect of all cardiac territories showed a maximum coefficient of variation of 9.74%.

CONCLUSION

Quantitation of MBF via kinetic modeling of cardiac rest MBF by [13N]-NH3 is minimally affected by the use of a BPL reconstruction technique, with BPL images presenting with less noise.

Respiratory motion correction in 4D-PET by simultaneous motion estimation and image reconstruction (SMEIR)

In conventional 4D positron emission tomography (4D-PET), images from different frames are reconstructed individually and aligned by registration methods. Two issues that arise with this approach are as follows: (1) the reconstruction algorithms do not make full use of projection statistics; and (2) the registration between noisy images can result in poor alignment. In this study, we investigated the use of simultaneous motion estimation and image reconstruction (SMEIR) methods for motion estimation/correction in 4D-PET. A modified ordered-subset expectation maximization algorithm coupled with total variation minimization (OSEM-TV) was used to obtain a primary motion-compensated PET (pmc-PET) from all projection data, using Demons derived deformation vector fields (DVFs) as initial motion vectors. A motion model update was performed to obtain an optimal set of DVFs in the pmc-PET and other phases, by matching the forward projection of the deformed pmc-PET with measured projections from other phases. The OSEM-TV image reconstruction was repeated using updated DVFs, and new DVFs were estimated based on updated images. A 4D-XCAT phantom with typical FDG biodistribution was generated to evaluate the performance of the SMEIR algorithm in lung and liver tumors with different contrasts and different diameters (10-40 mm). The image quality of the 4D-PET was greatly improved by the SMEIR algorithm. When all projections were used to reconstruct 3D-PET without motion compensation, motion blurring artifacts were present, leading up to 150% tumor size overestimation and significant quantitative errors, including 50% underestimation of tumor contrast and 59% underestimation of tumor uptake. Errors were reduced to less than 10% in most images by using the SMEIR algorithm, showing its potential in motion estimation/correction in 4D-PET.

Influence of respiratory motion correction on quantification of myocardial perfusion SPECT

BACKGROUND

Respiratory-related cardiac motion could have considerable effects on myocardial perfusion imaging, leading to misinterpretation of the images. In this study, we examined the influence of respiratory correction on ECG-gated myocardial perfusion SPECT (RC-GSPECT) concerning regional myocardial perfusion and function.

MATERIALS AND METHODS

Using the NCAT phantom, a typical torso phantom was generated. SimSET, a Monte Carlo simulator, was used to image the photon emerging from the phantom. Twenty-six patients underwent a 2-day stress-rest ECG-gated myocardial perfusion SPECT (GSPECT) imaging. A separate study was also performed by simultaneous respiratory and cardiac triggering with the real-time position management (RPM) for respiratory correction (RC).

RESULTS

In simulation study, count density in the inferior and inferoseptal walls increased in the lower bin of the respiratory cycle. On the other hand, there was a higher correlation between RC-GSPECT and echocardiography for left ventricular ejection fraction (LVEF) (r = 0.95, P < .01 vs r = 0.88, P < .01 for GSPECT).

CONCLUSION

We proposed a new approach for respiratory and cardiac-gated SPECT to eliminate respiratory motion artifacts. RC-GSPECT is a feasible method in MPI studies and may play an important role to improve the quality of MPI images, particularly in the inferior wall.

Quantitative myocardial perfusion PET parametric imaging at the voxel-level

Quantitative myocardial perfusion (MP) PET has the potential to enhance detection of early stages of atherosclerosis or microvascular dysfunction, characterization of flow-limiting effects of coronary artery disease (CAD), and identification of balanced reduction of flow due to multivessel stenosis. We aim to enable quantitative MP-PET at the individual voxel level, which has the potential to allow enhanced visualization and quantification of myocardial blood flow (MBF) and flow reserve (MFR) as computed from uptake parametric images. This framework is especially challenging for the (82)Rb radiotracer. The short half-life enables fast serial imaging and high patient throughput; yet, the acquired dynamic PET images suffer from high noise-levels introducing large variability in uptake parametric images and, therefore, in the estimates of MBF and MFR. Robust estimation requires substantial post-smoothing of noisy data, degrading valuable functional information of physiological and pathological importance. We present a feasible and robust approach to generate parametric images at the voxel-level that substantially reduces noise without significant loss of spatial resolution. The proposed methodology, denoted physiological clustering, makes use of the functional similarity of voxels to penalize deviation of voxel kinetics from physiological partners. The results were validated using extensive simulations (with transmural and non-transmural perfusion defects) and clinical studies. Compared to post-smoothing, physiological clustering depicted enhanced quantitative noise versus bias performance as well as superior recovery of perfusion defects (as quantified by CNR) with minimal increase in bias. Overall, parametric images obtained from the proposed methodology were robust in the presence of high-noise levels as manifested in the voxel time-activity-curves.

Patient-dependent count-rate adaptive normalization for PET detector efficiency with delayed-window coincidence events

Quantitative PET imaging is widely used in clinical diagnosis in oncology and neuroimaging. Accurate normalization correction for the efficiency of each line-of- response is essential for accurate quantitative PET image reconstruction. In this paper, we propose a normalization calibration method by using the delayed-window coincidence events from the scanning phantom or patient. The proposed method could dramatically reduce the 'ring' artifacts caused by mismatched system count-rates between the calibration and phantom/patient datasets. Moreover, a modified algorithm for mean detector efficiency estimation is proposed, which could generate crystal efficiency maps with more uniform variance. Both phantom and real patient datasets are used for evaluation. The results show that the proposed method could lead to better uniformity in reconstructed images by removing ring artifacts, and more uniform axial variance profiles, especially around the axial edge slices of the scanner. The proposed method also has the potential benefit to simplify the normalization calibration procedure, since the calibration can be performed using the on-the-fly acquired delayed-window dataset.

Challenges of radiotherapy: report on the 4D treatment planning workshop 2013

This report, compiled by experts on the treatment of mobile targets with advanced radiotherapy, summarizes the main conclusions and innovations achieved during the 4D treatment planning workshop 2013. This annual workshop focuses on research aiming to advance 4D radiotherapy treatments, including all critical aspects of time resolved delivery, such as in-room imaging, motion detection, motion managing, beam application, and quality assurance techniques. The report aims to revise achievements in the field and to discuss remaining challenges and potential solutions. As main achievements advances in the development of a standardized 4D phantom and in the area of 4D-treatment plan optimization were identified. Furthermore, it was noticed that MR imaging gains importance and high interest for sequential 4DCT/MR data sets was expressed, which represents a general trend of the field towards data covering a longer time period of motion. A new point of attention was work related to dose reconstructions, which may play a major role in verification of 4D treatment deliveries. The experimental validation of results achieved by 4D treatment planning and the systematic evaluation of different deformable image registration methods especially for inter-modality fusions were identified as major remaining challenges. A challenge that was also suggested as focus for future 4D workshops was the adaptation of image guidance approaches from conventional radiotherapy into particle therapy. Besides summarizing the last workshop, the authors also want to point out new evolving demands and give an outlook on the focus of the next workshop.

Self-gated radial MRI for respiratory motion compensation on hybrid PET/MR systems

Accurate localization and uptake quantification of lesions in the chest and abdomen using PET imaging is challenging due to the respiratory motion during the exam. The advent of hybrid PET/MR systems offers new ways to compensate for respiratory motion without exposing the patient to additional radiation. The use of self-gated reconstructions of a 3D radial stack-of-stars GRE acquisition is proposed to derive a high-resolution MRI motion model. The self-gating signal is used to perform respiratory binning of the simultaneously acquired PET raw data. Matching mu-maps are generated for every bin, and post-reconstruction registration is performed in order to obtain a motion-compensated PET volume from the individual gates. The proposed method is demonstrated in-vivo for three clinical patients. Motion-corrected reconstructions are compared against ungated and gated PET reconstructions. In all cases, motion-induced blurring of lesions in the liver and lung was substantially reduced, without compromising SNR as it is the case for gated reconstructions.

Dynamic whole-body PET parametric imaging: I. Concept, acquisition protocol optimization and clinical application

Static whole-body PET/CT, employing the standardized uptake value (SUV), is considered the standard clinical approach to diagnosis and treatment response monitoring for a wide range of oncologic malignancies. Alternative PET protocols involving dynamic acquisition of temporal images have been implemented in the research setting, allowing quantification of tracer dynamics, an important capability for tumor characterization and treatment response monitoring. Nonetheless, dynamic protocols have been confined to single-bed-coverage limiting the axial field-of-view to ~15-20 cm, and have not been translated to the routine clinical context of whole-body PET imaging for the inspection of disseminated disease. Here, we pursue a transition to dynamic whole-body PET parametric imaging, by presenting, within a unified framework, clinically feasible multi-bed dynamic PET acquisition protocols and parametric imaging methods. We investigate solutions to address the challenges of: (i) long acquisitions, (ii) small number of dynamic frames per bed, and (iii) non-invasive quantification of kinetics in the plasma. In the present study, a novel dynamic (4D) whole-body PET acquisition protocol of ~45 min total length is presented, composed of (i) an initial 6 min dynamic PET scan (24 frames) over the heart, followed by (ii) a sequence of multi-pass multi-bed PET scans (six passes × seven bed positions, each scanned for 45 s). Standard Patlak linear graphical analysis modeling was employed, coupled with image-derived plasma input function measurements. Ordinary least squares Patlak estimation was used as the baseline regression method to quantify the physiological parameters of tracer uptake rate Ki and total blood distribution volume V on an individual voxel basis. Extensive Monte Carlo simulation studies, using a wide set of published kinetic FDG parameters and GATE and XCAT platforms, were conducted to optimize the acquisition protocol from a range of ten different clinically acceptable sampling schedules examined. The framework was also applied to six FDG PET patient studies, demonstrating clinical feasibility. Both simulated and clinical results indicated enhanced contrast-to-noise ratios (CNRs) for Ki images in tumor regions with notable background FDG concentration, such as the liver, where SUV performed relatively poorly. Overall, the proposed framework enables enhanced quantification of physiological parameters across the whole body. In addition, the total acquisition length can be reduced from 45 to ~35 min and still achieve improved or equivalent CNR compared to SUV, provided the true Ki contrast is sufficiently high. In the follow-up companion paper, a set of advanced linear regression schemes is presented to particularly address the presence of noise, and attempt to achieve a better trade-off between the mean-squared error and the CNR metrics, resulting in enhanced task-based imaging.

The effect of regularization in motion compensated PET image reconstruction: a realistic numerical 4D simulation study

Following continuous improvement in PET spatial resolution, respiratory motion correction has become an important task. Two of the most common approaches that utilize all detected PET events to motion-correct PET data are the reconstruct-transform-average method (RTA) and motion-compensated image reconstruction (MCIR). In RTA, separate images are reconstructed for each respiratory frame, subsequently transformed to one reference frame and finally averaged to produce a motion-corrected image. In MCIR, the projection data from all frames are reconstructed by including motion information in the system matrix so that a motion-corrected image is reconstructed directly. Previous theoretical analyses have explained why MCIR is expected to outperform RTA. It has been suggested that MCIR creates less noise than RTA because the images for each separate respiratory frame will be severely affected by noise. However, recent investigations have shown that in the unregularized case RTA images can have fewer noise artefacts, while MCIR images are more quantitatively accurate but have the common salt-and-pepper noise. In this paper, we perform a realistic numerical 4D simulation study to compare the advantages gained by including regularization within reconstruction for RTA and MCIR, in particular using the median-root-prior incorporated in the ordered subsets maximum a posteriori one-step-late algorithm. In this investigation we have demonstrated that MCIR with proper regularization parameters reconstructs lesions with less bias and root mean square error and similar CNR and standard deviation to regularized RTA. This finding is reproducible for a variety of noise levels (25, 50, 100 million counts), lesion sizes (8 mm, 14 mm diameter) and iterations. Nevertheless, regularized RTA can also be a practical solution for motion compensation as a proper level of regularization reduces both bias and mean square error.