Joint reconstruction of image and motion in gated positron emission tomography

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.

Deformation-Compensated Learning for Image Reconstruction Without Ground Truth

Deep neural networks for medical image reconstruction are traditionally trained using high-quality ground-truth images as training targets. Recent work on Noise2Noise (N2N) has shown the potential of using multiple noisy measurements of the same object as an alternative to having a ground-truth. However, existing N2N-based methods are not suitable for learning from the measurements of an object undergoing nonrigid deformation. This paper addresses this issue by proposing the deformation-compensated learning (DeCoLearn) method for training deep reconstruction networks by compensating for object deformations. A key component of DeCoLearn is a deep registration module, which is jointly trained with the deep reconstruction network without any ground-truth supervision. We validate DeCoLearn on both simulated and experimentally collected magnetic resonance imaging (MRI) data and show that it significantly improves imaging quality.

Deep Learning Based Joint PET Image Reconstruction and Motion Estimation

Respiratory motion is one of the main sources of motion artifacts in positron emission tomography (PET) imaging. The emission image and patient motion can be estimated simultaneously from respiratory gated data through a joint estimation framework. However, conventional motion estimation methods based on registration of a pair of images are sensitive to noise. The goal of this study is to develop a robust joint estimation method that incorporates a deep learning (DL)-based image registration approach for motion estimation. We propose a joint estimation framework by incorporating a learned image registration network into a regularized PET image reconstruction. The joint estimation was formulated as a constrained optimization problem with moving gated images related to a fixed image via the deep neural network. The constrained optimization problem is solved by the alternating direction method of multipliers (ADMM) algorithm. The effectiveness of the algorithm was demonstrated using simulated and real data. We compared the proposed DL-ADMM joint estimation algorithm with a monotonic iterative joint estimation. Motion compensated reconstructions using pre-calculated deformation fields by DL-based (DL-MC recon) and iterative (iterative-MC recon) image registration were also included for comparison. Our simulation study shows that the proposed DL-ADMM joint estimation method reduces bias compared to the ungated image without increasing noise and outperforms the competing methods. In the real data study, our proposed method also generated higher lesion contrast and sharper liver boundaries compared to the ungated image and had lower noise than the reference gated image.

PET respiratory motion correction: quo vadis?

Positron emission tomography (PET) respiratory motion correction has been a subject of great interest for the last twenty years, prompted mainly by the development of multimodality imaging devices such as PET/computed tomography (CT) and PET/magnetic resonance imaging (MRI). PET respiratory motion correction involves a number of steps including acquisition synchronization, motion estimation and finally motion correction. The synchronization steps include the use of different external device systems or data driven approaches which have been gaining ground over the last few years. Patient specific or generic motion models using the respiratory synchronized datasets can be subsequently derived and used for correction either in the image space or within the image reconstruction process. Similar overall approaches can be considered and have been proposed for both PET/CT and PET/MRI devices. Certain variations in the case of PET/MRI include the use of MRI specific sequences for the registration of respiratory motion information. The proposed review includes a comprehensive coverage of all these areas of development in field of PET respiratory motion for different multimodality imaging devices and approaches in terms of synchronization, estimation and subsequent motion correction. Finally, a section on perspectives including the potential clinical usage of these approaches is included.

An optimized imaging protocol for [99mTc]Tc-DPD scintigraphy and SPECT/CT quantification in cardiac transthyretin (ATTR) amyloidosis

BACKGROUND

In [99mTc]Tc-DPD scintigraphy for myocardial ATTR amyloidosis, planar images 3 hour p.i. and SPECT/CT acquisition in L-mode are recommended. This study investigated if earlier planar images (1 hour p.i.) are beneficial and if SPECT/CT acquisition should be preferred in H-mode (180° detector angle) or L-mode (90°).

METHODS

In SPECT/CT phantom measurements (NaI cameras, N = 2; CZT, N = 1), peak contrast recovery (CRpeak) was derived from sphere inserts or myocardial insert (cardiac phantom; signal-to-background ratio [SBR], 10:1 or 5:1). In 25 positive and 38 negative patients (reference: endomyocardial biopsy or clinical diagnosis), Perugini scores and heart-to-contralateral (H/CL) count ratios were derived from planar images 1 hour and 3 hour p.i.

RESULTS

In phantom measurements, accuracy of myocardial CRpeak at SBR 10:1 (H-mode, 0.95-0.99) and reproducibility at 5:1 (H-mode, 1.02-1.14) was comparable for H-mode and L-mode. However, L-mode showed higher variability of background counts and sphere CRpeak throughout the field of view than H-mode. In patients, sensitivity/specificity were ≥ 95% for H/CL ratios at both time points and visual scoring 3 hour. At 1 hour, visual scores showed specificity of 89% and reduced reader's confidence.

CONCLUSIONS

Early DPD images provided no additional value for visual scoring or H/CL ratios. In SPECT/CT, H-mode is preferred over L-mode, especially if quantification is applied apart from the myocardium.

General simultaneous motion estimation and image reconstruction (G-SMEIR)

To achieve better performance for 4D multi-frame reconstruction with the parametric motion model (MF-PMM), a general simultaneous motion estimation and image reconstruction (G-SMEIR) method is proposed. In G-SMEIR, projection domain motion estimation and image domain motion estimation are performed alternatively to achieve better 4D reconstruction. This method can mitigate the local optimum trapping problem in either domain. To improve computational efficiency, the image domain motion estimation is accelerated by adapting fast convergent algorithms and graphics processing unit (GPU) computing. The proposed G-SMEIR method is tested using a cone-beam computed tomography (CBCT) simulation study of 4D XCAT phantom at different dose levels and compared with 3D total variation-based reconstruction (3D TV), 4D reconstruction with image domain motion estimation (IM4D), and SMEIR. G-SMEIR shows strong denoising capability and achieves similar performance at regular dose and half dose. The root mean squared error (RMSE) of G-SMEIR is the best among the four methods and improved about 12% over SMEIR for all respiratory phase images at full dose. G-SMEIR also achieved the best structural similarity index (SSIM) values among all methods. More importantly, G-SMEIR leads to more than 40% improvement of the mean deviation from the phantom tumor motion over SMEIR. A preliminary patient CBCT image reconstruction also shows better image quality of G-SMEIR than that of the frame-by-frame reconstruction (3D TV) and MF-PMM either using image domain motion estimation (IM4D) or using projection domain motion estimation (SMEIR) alone. G-SMEIR with a flexible combination of image domain and projection domain motion estimation provides an effective tool for 4D tomographic reconstruction.

Motion estimation and correction in SPECT, PET and CT

Patient motion impacts single photon emission computed tomography (SPECT), positron emission tomography (PET) and X-ray computed tomography (CT) by giving rise to projection data inconsistencies that can manifest as reconstruction artifacts, thereby degrading image quality and compromising accurate image interpretation and quantification. Methods to estimate and correct for patient motion in SPECT, PET and CT have attracted considerable research effort over several decades. The aims of this effort have been two-fold: to estimate relevant motion fields characterizing the various forms of voluntary and involuntary motion; and to apply these motion fields within a modified reconstruction framework to obtain motion-corrected images. The aims of this review are to outline the motion problem in medical imaging and to critically review published methods for estimating and correcting for the relevant motion fields in clinical and preclinical SPECT, PET and CT. Despite many similarities in how motion is handled between these modalities, utility and applications vary based on differences in temporal and spatial resolution. Technical feasibility has been demonstrated in each modality for both rigid and non-rigid motion, but clinical feasibility remains an important target. There is considerable scope for further developments in motion estimation and correction, and particularly in data-driven methods that will aid clinical utility. State-of-the-art machine learning methods may have a unique role to play in this context.

Simulation Study of a Frame-Based Motion Correction Algorithm for Positron Emission Imaging

Positron emission tomography (PET) is a functional non-invasive imaging modality that uses radioactive substances (radiotracers) to measure changes in metabolic processes. Advances in scanner technology and data acquisition in the last decade have led to the development of more sophisticated PET devices with good spatial resolution (1-3 mm of full width at half maximum (FWHM)). However, there are involuntary motions produced by the patient inside the scanner that lead to image degradation and potentially to a misdiagnosis. The adverse effect of the motion in the reconstructed image increases as the spatial resolution of the current scanners continues improving. In order to correct this effect, motion correction techniques are becoming increasingly popular and further studied. This work presents a simulation study of an image motion correction using a frame-based algorithm. The method is able to cut the acquired data from the scanner in frames, taking into account the size of the object of study. This approach allows working with low statistical information without losing image quality. The frames are later registered using spatio-temporal registration developed in a multi-level way. To validate these results, several performance tests are applied to a set of simulated moving phantoms. The results obtained show that the method minimizes the intra-frame motion, improves the signal intensity over the background in comparison with other literature methods, produces excellent values of similarity with the ground-truth (static) image and is able to find a limit in the patient-injected dose when some prior knowledge of the lesion is present.

Quantitative PET in the 2020s: a roadmap

Positron emission tomography (PET) plays an increasingly important role in research and clinical applications, catalysed by remarkable technical advances and a growing appreciation of the need for reliable, sensitive biomarkers of human function in health and disease. Over the last 30 years, a large amount of the physics and engineering effort in PET has been motivated by the dominant clinical application during that period, oncology. This has led to important developments such as PET/CT, whole-body PET, 3D PET, accelerated statistical image reconstruction, and time-of-flight PET. Despite impressive improvements in image quality as a result of these advances, the emphasis on static, semi-quantitative 'hot spot' imaging for oncologic applications has meant that the capability of PET to quantify biologically relevant parameters based on tracer kinetics has not been fully exploited. More recent advances, such as PET/MR and total-body PET, have opened up the ability to address a vast range of new research questions, from which a future expansion of applications and radiotracers appears highly likely. Many of these new applications and tracers will, at least initially, require quantitative analyses that more fully exploit the exquisite sensitivity of PET and the tracer principle on which it is based. It is also expected that they will require more sophisticated quantitative analysis methods than those that are currently available. At the same time, artificial intelligence is revolutionizing data analysis and impacting the relationship between the statistical quality of the acquired data and the information we can extract from the data. In this roadmap, leaders of the key sub-disciplines of the field identify the challenges and opportunities to be addressed over the next ten years that will enable PET to realise its full quantitative potential, initially in research laboratories and, ultimately, in clinical practice.

Variational multi-task MRI reconstruction: Joint reconstruction, registration and super-resolution

Motion degradation is a central problem in Magnetic Resonance Imaging (MRI). This work addresses the problem of how to obtain higher quality, super-resolved motion-free reconstructions from highly undersampled MRI data. In this work, we present for the first time a variational multi-task framework that allows joining three relevant tasks in MRI: reconstruction, registration and super-resolution. Our framework takes a set of multiple undersampled MR acquisitions corrupted by motion into a novel multi-task optimisation model, which is composed of an L² fidelity term that allows sharing representation between tasks, super-resolution foundations and hyperelastic deformations to model biological tissue behaviors. We demonstrate that this combination yields significant improvements over sequential models and other bi-task methods. Our results exhibit fine details and compensate for motion producing sharp and highly textured images compared to state of the art methods while keeping low CPU time. Our improvements are appraised on both clinical assessment and statistical analysis.

Joint compensation of motion and partial volume effects by iterative deconvolution incorporating wavelet-based denoising in oncologic PET/CT imaging

OBJECTIVES

We aim to develop and rigorously evaluate an image-based deconvolution method to jointly compensate respiratory motion and partial volume effects (PVEs) for quantitative oncologic PET imaging, including studying the impact of various reconstruction algorithms on quantification performance.

PROCEDURES

An image-based deconvolution method that incorporated wavelet-based denoising within the Lucy-Richardson algorithm was implemented and assessed. The method was evaluated using phantom studies with signal-to-background ratios (SBR) of 4 and 8, and clinical data of 10 patients with 42 lung lesions ≤30 mm in diameter. In each study, PET images were reconstructed using four different algorithms: OSEM-basic, PSF, TOF, and TOFPSF. The performance was quantified using contrast recovery (CR), coefficient of variation (COV) and contrast-to-noise-ratio (CNR) metrics. Further, in each study, variabilities arising due to the four different reconstruction algorithms were assessed.

RESULTS

In phantom studies, incorporation of wavelet-based denoising improved COV in all cases. Processing images using proposed method yielded significantly higher CR and CNR particularly in small spheres, for all reconstruction algorithms and all SBRs (P < 0.05). In patient studies, processing images using the proposed method yielded significantly higher CR and CNR (P < 0.05). The choice of the reconstruction algorithm impacted quantification performance for changes in motion amplitude, tumor size and SBRs.

CONCLUSIONS

Our results provide strong evidence that the proposed joint-compensation method can yield improved PET quantification. The choice of the reconstruction algorithm led to changes in quantitative accuracy, emphasizing the need to carefully select the right combination of reconstruction-image-based compensation methods.

Respiratory Motion Correction Using A Novel Positron Emission Particle Tracking Technique: A Framework Towards Individual Lesion-Based Motion Correction

Respiratory motion during PET/CT imaging is a matter of concern due to degraded image quality and reduced quantitative accuracy caused by motion artifacts. One class of motion correction methods relies on hardware-based respiratory motion tracking systems in order to use respiratory cycles for correcting motion artifacts. Another class of hardware-free methods extract motion information from the reconstructed images or sinograms. Hardware-based methods, however, are limited by calibration requirement, patient discomfort, lack of adaptability during scanning, presence of electronic drift during respiratory monitoring etc. Extracting motion information from reconstructed images is also limited by the fact that the original raw information requires significant processing before it can be used. Hence the motivation behind this work is to introduce a software-based approach that can be applied on raw 64-bit listmode data. The basic design of the proposed method is based on the fundamentals of Positron Emission Particle Tracking (PEPT) with additional incorporation of Time of Flight (TOF) information. Respiratory motion of patients has been extracted from the raw PET data by tracking a point source attached to the patient in areas on and near the chest. The key objective of this work is to describe a new process by which this particle tracking based motion correction system can eventually be lesion specific and correct the motion for a particular lesion within the patient. This work thus serves as a framework for lesion specific motion correction.

MRI-Driven PET Image Optimization for Neurological Applications

Positron emission tomography (PET) and magnetic resonance imaging (MRI) are established imaging modalities for the study of neurological disorders, such as epilepsy, dementia, psychiatric disorders and so on. Since these two available modalities vary in imaging principle and physical performance, each technique has its own advantages and disadvantages over the other. To acquire the mutual complementary information and reinforce each other, there is a need for the fusion of PET and MRI. This combined dual-modality (either sequential or simultaneous) could generate preferable soft tissue contrast of brain tissue, flexible acquisition parameters, and minimized exposure to radiation. The most unique superiority of PET/MRI is mainly manifested in MRI-based improvement for the inherent limitations of PET, such as motion artifacts, partial volume effect (PVE) and invasive procedure in quantitative analysis. Head motion during scanning significantly deteriorates the effective resolution of PET image, especially for the dynamic scan with lengthy time. Hybrid PET/MRI device can offer motion correction (MC) for PET data through MRI information acquired simultaneously. Regarding the PVE associated with limited spatial resolution, the process and reconstruction of PET data can be further optimized by using acquired MRI either sequentially or simultaneously. The quantitative analysis of dynamic PET data mainly relies upon an invasive arterial blood sampling procedure to acquire arterial input function (AIF). An image-derived input function (IDIF) method without the need of arterial cannulization, can serve as a potential alternative estimation of AIF. Compared with using PET data only, combining anatomical or functional information from MRI for improving the accuracy in IDIF approach has been demonstrated. Yet, due to the interference and inherent disparity between the two modalities, these methods for optimizing PET image based on MRI still have many technical challenges. This review discussed upon the most recent progress, current challenges and future directions of MRI-driven PET data optimization for neurological applications, with either sequential or simultaneous acquisition approach.

Respiratory- and cardiac motion-corrected simultaneous whole-heart PET and dual phase coronary MR angiography

PURPOSE

To develop a framework for efficient and simultaneous acquisition of motion-compensated whole-heart coronary MR angiography (CMRA) and left ventricular function by MR and myocardial integrity by PET on a 3T PET-MR system.

METHODS

An acquisition scheme based on a dual-phase CMRA sequence acquired simultaneously with cardiac PET data has been developed. The framework is integrated with a motion-corrected image reconstruction approach, so that non-rigid respiratory and cardiac deformation fields estimated from MR images are used to correct both the CMRA (respiratory motion correction for each cardiac phase) and the PET data (respiratory and cardiac motion correction). The proposed approach was tested in a cohort of 8 healthy subjects and 6 patients with coronary artery disease. Left ventricular (LV) function estimated from motion-corrected dual-phase CMRA was compared to the gold standard estimated from a stack of 2D CINE images for the healthy subjects. Relative increase of signal in motion-corrected PET images compared to uncorrected images was computed for standard 17-segment polar maps for each patient.

RESULTS

Motion-corrected dual-phase CMRA images allow for visualization of the coronary arteries in both systole and diastole for all healthy subjects and cardiac patients. LV functional indices from healthy subjects result in good agreement with the reference method, underestimating stroke volume by 3.07 ± 3.26 mL and ejection fraction by 0.30 ± 1.01%. Motion correction improved delineation of the myocardium in PET images, resulting in an increased 18 F-FDG signal of up to 28% in basal segments of the myocardial wall compared to uncorrected images.

CONCLUSION

The proposed motion-corrected dual-phase CMRA and cardiac PET produces co-registered good quality images in both modalities in a single efficient examination of ~13 min.

Multi-Tracer Guided PET Image Reconstruction

Multi-tracer positron emission tomography (PET) has the potential to enhance PET imaging by providing complementary information from different physiological processes. However, one or more of the images may present high levels of noise. Guided image reconstruction methods transfer information from a guide image into the PET image reconstruction to encourage edge-preserving noise reduction. In this work we aim to reduce noise in poorer quality PET datasets via guidance from higher quality ones by using a weighted quadratic penalty approach. In particular, we applied this methodology to [18F]fluorodeoxyglucose (FDG) and [11C]methionine imaging of gliomas. 3D simulation studies showed that guiding the reconstruction of methionine datasets using pre-existing FDG images reduced reconstruction errors across the whole-brain (-8%) and within a tumour (-36%) compared to maximum likelihood expectation-maximisation (MLEM). Furthermore, guided reconstruction outperformed a comparable non-local means filter, indicating that regularising during reconstruction is preferable to post-reconstruction approaches. Hyperparameters selected from the 3D simulation study were applied to real data, where it was observed that the proposed FDG-guided methionine reconstruction allows for better edge preservation and noise reduction than standard MLEM. Overall, the results in this work demonstrate that transferring information between datasets in multi-tracer PET studies improves image quality and quantification performance.

4D cone-beam computed tomography (CBCT) using a moving blocker for simultaneous radiation dose reduction and scatter correction

Four-dimensional (4D) x-ray cone-beam computed tomography (CBCT) is important for a precise radiation therapy for lung cancer. Due to the repeated use and 4D acquisition over a course of radiotherapy, the radiation dose becomes a concern. Meanwhile, the scatter contamination in CBCT deteriorates image quality for treatment tasks. In this work, we propose the use of a moving blocker (MB) during the 4D CBCT acquisition ('4D MB') and to combine motion-compensated reconstruction to address these two issues simultaneously. In 4D MB CBCT, the moving blocker reduces the x-ray flux passing through the patient and collects the scatter information in the blocked region at the same time. The scatter signal is estimated from the blocked region for correction. Even though the number of projection views and projection data in each view are not complete for conventional reconstruction, 4D reconstruction with a total-variation (TV) constraint and a motion-compensated temporal constraint can utilize both spatial gradient sparsity and temporal correlations among different phases to overcome the missing data problem. The feasibility simulation studies using the 4D NCAT phantom showed that 4D MB with motion-compensated reconstruction with 1/3 imaging dose reduction could produce satisfactory images and achieve 37% improvement on structural similarity (SSIM) index and 55% improvement on root mean square error (RMSE), compared to 4D reconstruction at the regular imaging dose without scatter correction. For the same 4D MB data, 4D reconstruction outperformed 3D TV reconstruction by 28% on SSIM and 34% on RMSE. A study of synthetic patient data also demonstrated the potential of 4D MB to reduce the radiation dose by 1/3 without compromising the image quality. This work paves the way for more comprehensive studies to investigate the dose reduction limit offered by this novel 4D MB method using physical phantom experiments and real patient data based on clinical relevant metrics.

Penalized maximum likelihood simultaneous longitudinal PET image reconstruction with difference-image priors

PURPOSE

Many clinical contexts require the acquisition of multiple positron emission tomography (PET) scans of a single subject, for example, to observe and quantitate changes in functional behaviour in tumors after treatment in oncology. Typically, the datasets from each of these scans are reconstructed individually, without exploiting the similarities between them. We have recently shown that sharing information between longitudinal PET datasets by penalizing voxel-wise differences during image reconstruction can improve reconstructed images by reducing background noise and increasing the contrast-to-noise ratio of high-activity lesions. Here, we present two additional novel longitudinal difference-image priors and evaluate their performance using two-dimesional (2D) simulation studies and a three-dimensional (3D) real dataset case study.

METHODS

We have previously proposed a simultaneous difference-image-based penalized maximum likelihood (PML) longitudinal image reconstruction method that encourages sparse difference images (DS-PML), and in this work we propose two further novel prior terms. The priors are designed to encourage longitudinal images with corresponding differences which have (a) low entropy (DE-PML), and (b) high sparsity in their spatial gradients (DTV-PML). These two new priors and the originally proposed longitudinal prior were applied to 2D-simulated treatment response [¹⁸ F]fluorodeoxyglucose (FDG) brain tumor datasets and compared to standard maximum likelihood expectation-maximization (MLEM) reconstructions. These 2D simulation studies explored the effects of penalty strengths, tumor behaviour, and interscan coupling on reconstructed images. Finally, a real two-scan longitudinal data series acquired from a head and neck cancer patient was reconstructed with the proposed methods and the results compared to standard reconstruction methods.

RESULTS

Using any of the three priors with an appropriate penalty strength produced images with noise levels equivalent to those seen when using standard reconstructions with increased counts levels. In tumor regions, each method produces subtly different results in terms of preservation of tumor quantitation and reconstruction root mean-squared error (RMSE). In particular, in the two-scan simulations, the DE-PML method produced tumor means in close agreement with MLEM reconstructions, while the DTV-PML method produced the lowest errors due to noise reduction within the tumor. Across a range of tumor responses and different numbers of scans, similar results were observed, with DTV-PML producing the lowest errors of the three priors and DE-PML producing the lowest bias. Similar improvements were observed in the reconstructions of the real longitudinal datasets, although imperfect alignment of the two PET images resulted in additional changes in the difference image that affected the performance of the proposed methods.

CONCLUSION

Reconstruction of longitudinal datasets by penalizing difference images between pairs of scans from a data series allows for noise reduction in all reconstructed images. An appropriate choice of penalty term and penalty strength allows for this noise reduction to be achieved while maintaining reconstruction performance in regions of change, either in terms of quantitation of mean intensity via DE-PML, or in terms of tumor RMSE via DTV-PML. Overall, improving the image quality of longitudinal datasets via simultaneous reconstruction has the potential to improve upon currently used methods, allow dose reduction, or reduce scan time while maintaining image quality at current levels.

Simultaneous maximum a posteriori longitudinal PET image reconstruction

Positron emission tomography (PET) is frequently used to monitor functional changes that occur over extended time scales, for example in longitudinal oncology PET protocols that include routine clinical follow-up scans to assess the efficacy of a course of treatment. In these contexts PET datasets are currently reconstructed into images using single-dataset reconstruction methods. Inspired by recently proposed joint PET-MR reconstruction methods, we propose to reconstruct longitudinal datasets simultaneously by using a joint penalty term in order to exploit the high degree of similarity between longitudinal images. We achieved this by penalising voxel-wise differences between pairs of longitudinal PET images in a one-step-late maximum a posteriori (MAP) fashion, resulting in the MAP simultaneous longitudinal reconstruction (SLR) method. The proposed method reduced reconstruction errors and visually improved images relative to standard maximum likelihood expectation-maximisation (ML-EM) in simulated 2D longitudinal brain tumour scans. In reconstructions of split real 3D data with inserted simulated tumours, noise across images reconstructed with MAP-SLR was reduced to levels equivalent to doubling the number of detected counts when using ML-EM. Furthermore, quantification of tumour activities was largely preserved over a variety of longitudinal tumour changes, including changes in size and activity, with larger changes inducing larger biases relative to standard ML-EM reconstructions. Similar improvements were observed for a range of counts levels, demonstrating the robustness of the method when used with a single penalty strength. The results suggest that longitudinal regularisation is a simple but effective method of improving reconstructed PET images without using resolution degrading priors.

PET Image Reconstruction and Deformable Motion Correction Using Unorganized Point Clouds

Quantitative positron emission tomography imaging often requires correcting the image data for deformable motion. With cyclic motion, this is traditionally achieved by separating the coincidence data into a relatively small number of gates, and incorporating the inter-gate image transformation matrices into the reconstruction algorithm. In the presence of non-cyclic deformable motion, this approach may be impractical due to a large number of required gates. In this paper, we propose an alternative approach to iterative image reconstruction with correction for deformable motion, wherein unorganized point clouds are used to model the imaged objects in the image space, and motion is corrected for explicitly by introducing a time-dependence into the point coordinates. The image function is represented using constant basis functions with finite support determined by the boundaries of the Voronoi cells in the point cloud. We validate the quantitative accuracy and stability of the proposed approach by reconstructing noise-free and noisy projection data from digital and physical phantoms. The point-cloud-based maximum likelihood expectation maximization (MLEM) and one-pass list-mode ordered-subset expectation maximization (OSEM) algorithms are validated. The results demonstrate that images reconstructed using the proposed method are quantitatively stable, with noise and convergence properties comparable to image reconstruction based on the use of rectangular and radially-symmetric basis functions.

List-mode image reconstruction for positron emission tomography using tetrahedral voxels

Image space decomposition based on tetrahedral voxels are interesting candidates for use in emission tomography. Tetrahedral voxels provide many of the advantages of point clouds with irregular spacing, such as being intrinsically multi-resolution, yet they also serve as a volumetric partition of the image space and so are comparable to more standard cubic voxels. Additionally, non-rigid displacement fields can be applied to the tetrahedral mesh in a straight-forward manner. So far studies incorporating tetrahedral decomposition of the image space have concentrated on pre-calculated, node-based, system matrix elements which reduces the flexibility of the tetrahedral approach and the capacity to accurately define regions of interest. Here, a list-mode on-the-fly calculation of the system matrix elements is described using a tetrahedral decomposition of the image space and volumetric elements-voxels. The algorithm is demonstrated in the context of awake animal PET which may require both rigid and non-rigid motion compensation, as well as quantification within small regions of the brain. This approach allows accurate, event based, motion compensation including non-rigid deformations.

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.

Motion correction options in PET/MRI

Subject motion is unavoidable in clinical and research imaging studies. Breathing is the most important source of motion in whole-body PET and MRI studies, affecting not only thoracic organs but also those in the upper and even lower abdomen. The motion related to the pumping action of the heart is obviously relevant in high-resolution cardiac studies. These two sources of motion are periodic and predictable, at least to a first approximation, which means certain techniques can be used to control the motion (eg, by acquiring the data when the organ of interest is relatively at rest). Additionally, nonperiodic and unpredictable motion can also occur during the scan. One obvious limitation of methods relying on external devices (eg, respiratory bellows or the electrocardiogram signal to monitor the respiratory or cardiac cycle, respectively) to trigger or gate the data acquisition is that the complex motion of internal organs cannot be fully characterized. However, detailed information can be obtained using either the PET or MRI data (or both) allowing the more complete characterization of the motion field so that a motion model can be built. Such a model and the information derived from simple external devices can be used to minimize the effects of motion on the collected data. In the ideal case, all the events recorded during the PET scan would be used to generate a motion-free or corrected PET image. The detailed motion field can be used for this purpose by applying it to the PET data before, during, or after the image reconstruction. Integrating all these methods for motion control, characterization, and correction into a workflow that can be used for routine clinical studies is challenging but could potentially be extremely valuable given the improvement in image quality and reduction of motion-related image artifacts.

The Use of Anatomical Information for Molecular Image Reconstruction Algorithms: Attenuation/Scatter Correction, Motion Compensation, and Noise Reduction

PET and SPECT are important tools for providing valuable molecular information about patients to clinicians. Advances in nuclear medicine hardware technologies and statistical image reconstruction algorithms enabled significantly improved image quality. Sequentially or simultaneously acquired anatomical images such as CT and MRI from hybrid scanners are also important ingredients for improving the image quality of PET or SPECT further. High-quality anatomical information has been used and investigated for attenuation and scatter corrections, motion compensation, and noise reduction via post-reconstruction filtering and regularization in inverse problems. In this article, we will review works using anatomical information for molecular image reconstruction algorithms for better image quality by describing mathematical models, discussing sources of anatomical information for different cases, and showing some examples.

Maximum-likelihood joint image reconstruction and motion estimation with misaligned attenuation in TOF-PET/CT

This work is an extension of our recent work on joint activity reconstruction/motion estimation (JRM) from positron emission tomography (PET) data. We performed JRM by maximization of the penalized log-likelihood in which the probabilistic model assumes that the same motion field affects both the activity distribution and the attenuation map. Our previous results showed that JRM can successfully reconstruct the activity distribution when the attenuation map is misaligned with the PET data, but converges slowly due to the significant cross-talk in the likelihood. In this paper, we utilize time-of-flight PET for JRM and demonstrate that the convergence speed is significantly improved compared to JRM with conventional PET data.

Maximum-Likelihood Joint Image Reconstruction/Motion Estimation in Attenuation-Corrected Respiratory Gated PET/CT Using a Single Attenuation Map

This work provides an insight into positron emission tomography (PET) joint image reconstruction/motion estimation (JRM) by maximization of the likelihood, where the probabilistic model accounts for warped attenuation. Our analysis shows that maximum-likelihood (ML) JRM returns the same reconstructed gates for any attenuation map (μ-map) that is a deformation of a given μ-map, regardless of its alignment with the PET gates. We derived a joint optimization algorithm accordingly, and applied it to simulated and patient gated PET data. We first evaluated the proposed algorithm on simulations of respiratory gated PET/CT data based on the XCAT phantom. Our results show that independently of which μ-map is used as input to JRM: (i) the warped μ-maps correspond to the gated μ-maps, (ii) JRM outperforms the traditional post-registration reconstruction and consolidation (PRRC) for hot lesion quantification and (iii) reconstructed gated PET images are similar to those obtained with gated μ-maps. This suggests that a breath-held μ-map can be used. We then applied JRM on patient data with a μ-map derived from a breath-held high resolution CT (HRCT), and compared the results with PRRC, where each reconstructed PET image was obtained with a corresponding cine-CT gated μ-map. Results show that JRM with breath-held HRCT achieves similar reconstruction to that using PRRC with cine-CT. This suggests a practical low-dose solution for implementation of motion-corrected respiratory gated PET/CT.

The development and initial evaluation of a realistic simulated SPECT dataset with simultaneous respiratory and cardiac motion for gated myocardial perfusion SPECT

We developed a realistic simulation dataset for simultaneous respiratory and cardiac (R&C) gated SPECT/CT using the 4D NURBS-based Cardiac-Torso (NCAT) Phantom and Monte Carlo simulation methods, and evaluated it for a sample application study. The 4D NCAT phantom included realistic respiratory motion and beating heart motion based on respiratory gated CT and cardiac tagged MRI data of normal human subjects. To model the respiratory motion, a set of 24 separate 3D NCAT phantoms excluding the heart was generated over a respiratory cycle. The beating heart motion was modeled separately with 48 frames per cardiac cycle for each of the 24 respiratory phases. The resultant set of 24  ×  48 3D NCAT phantoms provides a realistic model of a normal human subject at different phases of combined R&C motions. An almost noise-free SPECT projection dataset for each of the 1152 3D NCAT phantoms was generated using Monte Carlo simulation techniques and the radioactivity uptake distribution of (99m)Tc sestamibi in different organs. By grouping and summing the separate projection datasets, separate or simultaneous R&C gated acquired data with different gating schemes could be simulated. In the initial evaluation, we combined the projection datasets into ungated, 6 respiratory-gates only, 8 cardiac-gates only, and combined 6 respiratory-gates & 8 cardiac-gates projection datasets. Each dataset was reconstructed using 3D OS-EM without and with attenuation correction using the averaged and respiratory-gated attenuation maps, and the resulting reconstructed images were compared. These results were used to demonstrate the effects of R&C motions and the reduction of image artifact due to R&C motions by gating and attenuation corrections. We concluded that the realistic 4D NCAT phantom and Monte Carlo simulated SPECT projection datasets with R&C motions are powerful tools in the study of the effects of R&C motions, as well as in the development of R&C gating schemes and motion correction methods for improved SPECT/CT imaging.

The potential of positron emission tomography for intratreatment dynamic lung tumor tracking: a phantom study

PURPOSE

This study aims to evaluate the potential and feasibility of positron emission tomography for dynamic lung tumor tracking during radiation treatment. The authors propose a center of mass (CoM) tumor tracking algorithm using gated-PET images combined with a respiratory monitor and investigate the geometric accuracy of the proposed algorithm.

METHODS

The proposed PET dynamic lung tumor tracking algorithm estimated the target position information through the CoM of the segmented target volume on gated PET images reconstructed from accumulated coincidence events. The information was continuously updated throughout a scan based on the assumption that real-time processing was supported (actual processing time at each frame ≈ 10 s). External respiratory motion and list-mode PET data were acquired from a phantom programmed to move with measured respiratory traces (external respiratory motion and internal target motion) from human subjects, for which the ground truth target position was known as a function of time. The phantom was cylindrical with six hollow sphere targets (10, 13, 17, 22, 28, and 37 mm in diameter). The measured respiratory traces consisted of two sets: (1) 1D-measured motion from ten healthy volunteers and (2) 3D-measured motion from four lung cancer patients. The authors evaluated the geometric accuracy of the proposed algorithm by quantifying estimation errors (Euclidean distance) between the actual motion of targets (1D-motion and 3D-motion traces) and CoM trajectories estimated by the proposed algorithm as a function of time.

RESULTS

The time-averaged error of 1D-motion traces over all trajectories of all targets was 1.6 mm. The error trajectories decreased with time as coincidence events were accumulated. The overall error trajectory of 1D-motion traces converged to within 2 mm in approximately 90 s. As expected, more accurate results were obtained for larger targets. For example, for the 37 mm target, the average error over all 1D-motion traces was 1.1 mm; and for the 10 mm target, the average error over all 1D-motion traces was 2.8 mm. The overall time-averaged error of 3D-motion traces was 1.6 mm, which was comparable to that of the 1D-motion traces. There were small variations in the errors between the 3D-motion traces, although the motion trajectories were very different. The accuracy of the estimates was consistent for all targets except for the smallest.

CONCLUSIONS

The authors developed an algorithm for dynamic lung tumor tracking using list-mode PET data and a respiratory motion signal, and demonstrated proof-of-principle for PET-guided lung tumor tracking. The overall tracking error in phantom studies is less than 2 mm.

Motion freeze for respiration motion correction in PET/CT: a preliminary investigation with lung cancer patient data

PURPOSE

Respiratory motion presents significant challenges for accurate PET/CT. It often introduces apparent increase of lesion size, reduction of measured standardized uptake value (SUV), and the mismatch in PET/CT fusion images. In this study, we developed the motion freeze method to use 100% of the counts collected by recombining the counts acquired from all phases of gated PET data into a single 3D PET data, with correction of respiration by deformable image registration.

METHODS

Six patients with diagnosis of lung cancer confirmed by oncologists were recruited. PET/CT scans were performed with Discovery STE system. The 4D PET/CT with the Varian real-time position management for respiratory motion tracking was followed by a clinical 3D PET/CT scan procedure in the static mode. Motion freeze applies the deformation matrices calculated by optical flow method to generate a single 3D effective PET image using the data from all the 4D PET phases.

RESULTS

The increase in SUV and decrease in tumor size with motion freeze for all lesions compared to the results from 3D and 4D was observed in the preliminary data of lung cancer patients. In addition, motion freeze substantially reduced tumor mismatch between the CT image and the corresponding PET images.

CONCLUSION

Motion freeze integrating 100% of the PET counts has the potential to eliminate the influences induced by respiratory motion in PET data.

Management of respiratory motion in PET/computed tomography: the state of the art

Combined PET/computed tomography (CT) is of value in cancer diagnosis, follow-up, and treatment planning. For cancers located in the thorax or abdomen, the patient’s breathing causes artifacts and errors in PET and CT images. Many different approaches for artifact avoidance or correction have been developed; most are based on gated acquisition and synchronization between the respiratory signal and PET acquisition. The respiratory signal is usually produced by an external sensor that tracks a physiological characteristic related to the patient’s breathing. Respiratory gating is a compensation technique in which time or amplitude binning is used to exclude the motion in reconstructed PET images. Although this technique is performed in routine clinical practice, it fails to adequately correct for respiratory motion because each gate can mix several tissue positions. Researchers have suggested either selecting PET events from gated acquisitions or performing several PET acquisitions (corresponding to a breath-hold CT position). However, the PET acquisition time must be increased if adequate counting statistics are to be obtained in the different gates after binning. Hence, other researchers have assessed correction techniques that take account of all the counting statistics (without increasing the acquisition duration) and integrate motion information before, during, or after the reconstruction process. Here, we provide an overview of how motion is managed to overcome respiratory motion in PET/CT images.

Simultaneous motion estimation and image reconstruction (SMEIR) for 4D cone-beam CT

PURPOSE

Image reconstruction and motion model estimation in four-dimensional cone-beam CT (4D-CBCT) are conventionally handled as two sequential steps. Due to the limited number of projections at each phase, the image quality of 4D-CBCT is degraded by view aliasing artifacts, and the accuracy of subsequent motion modeling is decreased by the inferior 4D-CBCT. The objective of this work is to enhance both the image quality of 4D-CBCT and the accuracy of motion model estimation with a novel strategy enabling simultaneous motion estimation and image reconstruction (SMEIR).

METHODS

The proposed SMEIR algorithm consists of two alternating steps: (1) model-based iterative image reconstruction to obtain a motion-compensated primary CBCT (m-pCBCT) and (2) motion model estimation to obtain an optimal set of deformation vector fields (DVFs) between the m-pCBCT and other 4D-CBCT phases. The motion-compensated image reconstruction is based on the simultaneous algebraic reconstruction technique (SART) coupled with total variation minimization. During the forward- and backprojection of SART, measured projections from an entire set of 4D-CBCT are used for reconstruction of the m-pCBCT by utilizing the updated DVF. The DVF is estimated by matching the forward projection of the deformed m-pCBCT and measured projections of other phases of 4D-CBCT. The performance of the SMEIR algorithm is quantitatively evaluated on a 4D NCAT phantom. The quality of reconstructed 4D images and the accuracy of tumor motion trajectory are assessed by comparing with those resulting from conventional sequential 4D-CBCT reconstructions (FDK and total variation minimization) and motion estimation (demons algorithm). The performance of the SMEIR algorithm is further evaluated by reconstructing a lung cancer patient 4D-CBCT.

RESULTS

Image quality of 4D-CBCT is greatly improved by the SMEIR algorithm in both phantom and patient studies. When all projections are used to reconstruct a 3D-CBCT by FDK, motion-blurring artifacts are present, leading to a 24.4% relative reconstruction error in the NACT phantom. View aliasing artifacts are present in 4D-CBCT reconstructed by FDK from 20 projections, with a relative error of 32.1%. When total variation minimization is used to reconstruct 4D-CBCT, the relative error is 18.9%. Image quality of 4D-CBCT is substantially improved by using the SMEIR algorithm and relative error is reduced to 7.6%. The maximum error (MaxE) of tumor motion determined from the DVF obtained by demons registration on a FDK-reconstructed 4D-CBCT is 3.0, 2.3, and 7.1 mm along left-right (L-R), anterior-posterior (A-P), and superior-inferior (S-I) directions, respectively. From the DVF obtained by demons registration on 4D-CBCT reconstructed by total variation minimization, the MaxE of tumor motion is reduced to 1.5, 0.5, and 5.5 mm along L-R, A-P, and S-I directions. From the DVF estimated by SMEIR algorithm, the MaxE of tumor motion is further reduced to 0.8, 0.4, and 1.5 mm along L-R, A-P, and S-I directions, respectively.

CONCLUSIONS

The proposed SMEIR algorithm is able to estimate a motion model and reconstruct motion-compensated 4D-CBCT. The SMEIR algorithm improves image reconstruction accuracy of 4D-CBCT and tumor motion trajectory estimation accuracy as compared to conventional sequential 4D-CBCT reconstruction and motion estimation.

Towards coronary plaque imaging using simultaneous PET-MR: a simulation study

Coronary atherosclerotic plaque rupture is the main cause of myocardial infarction and the leading killer in the US. Inflammation is a known bio-marker of plaque vulnerability and can be assessed non-invasively using fluorodeoxyglucose-positron emission tomography imaging (FDG-PET). However, cardiac and respiratory motion of the heart makes PET detection of coronary plaque very challenging. Fat surrounding coronary arteries allows the use of MRI to track plaque motion during simultaneous PET-MR examination. In this study, we proposed and assessed the performance of a fat-MR based coronary motion correction technique for improved FDG-PET coronary plaque imaging in simultaneous PET-MR. The proposed methods were evaluated in a realistic four-dimensional PET-MR simulation study obtained by combining patient water-fat separated MRI and XCAT anthropomorphic phantom. Five small lesions were digitally inserted inside the patients coronary vessels to mimic coronary atherosclerotic plaques. The heart of the XCAT phantom was digitally replaced with the patient's heart. Motion-dependent activity distributions, attenuation maps, and fat-MR volumes of the heart, were generated using the XCAT cardiac and respiratory motion fields. A full Monte Carlo simulation using Siemens mMR's geometry was performed for each motion phase. Cardiac/respiratory motion fields were estimated using non-rigid registration of the transformed fat-MR volumes and incorporated directly into the system matrix of PET reconstruction along with motion-dependent attenuation maps. The proposed motion correction method was compared to conventional PET reconstruction techniques such as no motion correction, cardiac gating, and dual cardiac-respiratory gating. Compared to uncorrected reconstructions, fat-MR based motion compensation yielded an average improvement of plaque-to-background contrast of 29.6%, 43.7%, 57.2%, and 70.6% for true plaque-to-blood ratios of 10, 15, 20 and 25:1, respectively. Channelized Hotelling observer (CHO) signal-to-noise ratio (SNR) was used to quantify plaque detectability. CHO-SNR improvement ranged from 105% to 128% for fat-MR-based motion correction as compared to no motion correction. Likewise, CHO-SNR improvement ranged from 348% to 396% as compared to both cardiac and dual cardiac-respiratory gating approaches. Based on this study, our approach, a fat-MR based motion correction for coronary plaque PET imaging using simultaneous PET-MR, offers great potential for clinical practice. The ultimate performance and limitation of our approach, however, must be fully evaluated in patient studies.

Registration-based reconstruction of four-dimensional cone beam computed tomography

We present a new method for reconstruction of 4-D cone beam computed tomography from an undersampled set of X-ray projections. The novelty of the proposed method lies in utilizing optical flow based registration to facilitate that each temporal phase is reconstructed from the full set of acquired projections. The reconstruction of each phase thus exhibits limited aliasing despite significant intra-phase undersampling. The method is fully self-contained. Initially an approximate 4-D volume is reconstructed and an inter-phase registration based hereon. A subsequent reconstruction pass integrates the optical flow estimation in a cost function formulation in which the X-ray projections from all temporal phases are considered for the reconstruction of each individual phase. Quantitative and qualitative evaluations were performed through reconstruction of both a numerical phantom and a clinical dataset. The obtained reconstructions are compared to the state-of-the-art alternatives of total variation regularization and prior image constrained compressed sensing. Our studies show that the proposed method is the better overall "compromise" in the depiction of both moving and stationary anatomical structures.

Magnetic resonance-based motion correction for positron emission tomography imaging

Positron emission tomography (PET) image quality is limited by patient motion. Emission data are blurred owing to cardiac and/or respiratory motion. Although spatial resolution is 4 mm for standard clinical whole-body PET scanners, the effective resolution can be as low as 1 cm owing to motion. Additionally, the deformation of attenuation medium causes image artifacts. Previously, gating has been used to "freeze" the motion, but led to significantly increased noise level. Simultaneous PET/magnetic resonance (MR) modality offers a new way to perform PET motion correction. MR can be used to measure 3-dimensional motion fields, which can then be incorporated into the iterative PET reconstruction to obtain motion-corrected PET images. In this report, we present MR imaging techniques to acquire dynamic images, a nonrigid image registration algorithm to extract motion fields from acquired MR images, and a PET reconstruction algorithm with motion correction. We also present results from both phantom and in vivo animal PET/MR studies. We demonstrate that MR-based PET motion correction using simultaneous PET/MR improves image quality and lesion detectability compared with gating and no motion correction.

Online tracking of deformable objects under occlusion using dominant points

This paper deals with tracking of deformable objects in the presence of occlusion using dominant point representation of the boundary contour. A novel nonintegral time propagation model for propagating the dominant points is proposed. It uses an initial guess generated from a linear operation and an analytical conjugate gradient approach for online robust learning of the shape deformation and motion model. A scheme is presented to automatically detect and correct the region of large local deformation. In order to deal with occlusion, admissible restrictions on deformation and motion of the object are automatically determined. The proposed method overcomes the need of offline learning and learns the deformation and motion model of the object using very few initial frames of the input video. The performance of the method is demonstrated using varieties of videos of different objects.

List-mode reconstruction for the Biograph mCT with physics modeling and event-by-event motion correction

Whole-body PET/CT scanners are important clinical and research tools to study tracer distribution throughout the body. In whole-body studies, respiratory motion results in image artifacts. We have previously demonstrated for brain imaging that, when provided with accurate motion data, event-by-event correction has better accuracy than frame-based methods. Therefore, the goal of this work was to develop a list-mode reconstruction with novel physics modeling for the Siemens Biograph mCT with event-by-event motion correction, based on the MOLAR platform (Motion-compensation OSEM List-mode Algorithm for Resolution-Recovery Reconstruction). Application of MOLAR for the mCT required two algorithmic developments. First, in routine studies, the mCT collects list-mode data in 32 bit packets, where averaging of lines-of-response (LORs) by axial span and angular mashing reduced the number of LORs so that 32 bits are sufficient to address all sinogram bins. This degrades spatial resolution. In this work, we proposed a probabilistic LOR (pLOR) position technique that addresses axial and transaxial LOR grouping in 32 bit data. Second, two simplified approaches for 3D time-of-flight (TOF) scatter estimation were developed to accelerate the computationally intensive calculation without compromising accuracy. The proposed list-mode reconstruction algorithm was compared to the manufacturer's point spread function + TOF (PSF+TOF) algorithm. Phantom, animal, and human studies demonstrated that MOLAR with pLOR gives slightly faster contrast recovery than the PSF+TOF algorithm that uses the average 32 bit LOR sinogram positioning. Moving phantom and a whole-body human study suggested that event-by-event motion correction reduces image blurring caused by respiratory motion. We conclude that list-mode reconstruction with pLOR positioning provides a platform to generate high quality images for the mCT, and to recover fine structures in whole-body PET scans through event-by-event motion correction.

Building motion models of lung tumours from cone-beam CT for radiotherapy applications

A method is presented to build a surrogate-driven motion model of a lung tumour from a cone-beam CT scan, which does not require markers. By monitoring an external surrogate in real time, it is envisaged that the motion model be used to drive gated or tracked treatments. The motion model would be built immediately before each fraction of treatment and can account for inter-fraction variation. The method could also provide a better assessment of tumour shape and motion prior to delivery of each fraction of stereotactic ablative radiotherapy. The two-step method involves enhancing the tumour region in the projections, and then fitting the surrogate-driven motion model. On simulated data, the mean absolute error was reduced to 1 mm. For patient data, errors were determined by comparing estimated and clinically identified tumour positions in the projections, scaled to mm at the isocentre. Averaged over all used scans, the mean absolute error was under 2.5 mm in superior-inferior and transverse directions.

Cardiac motion compensation and resolution modeling in simultaneous PET-MR: a cardiac lesion detection study

Cardiac motion and partial volume effects (PVE) are two of the main causes of image degradation in cardiac PET. Motion generates artifacts and blurring while PVE lead to erroneous myocardial activity measurements. Newly available simultaneous PET-MR scanners offer new possibilities in cardiac imaging as MRI can assess wall contractility while collecting PET perfusion data. In this perspective, we develop a list-mode iterative reconstruction framework incorporating both tagged-MR derived non-rigid myocardial wall motion and position dependent detector point spread function (PSF) directly into the PET system matrix. In this manner, our algorithm performs both motion 'deblurring' and PSF deconvolution while reconstructing images with all available PET counts. The proposed methods are evaluated in a beating non-rigid cardiac phantom whose hot myocardial compartment contains small transmural and non-transmural cold defects. In order to accelerate imaging time, we investigate collecting full and half k-space tagged MR data to obtain tagged volumes that are registered using non-rigid B-spline registration to yield wall motion information. Our experimental results show that tagged-MR based motion correction yielded an improvement in defect/myocardium contrast recovery of 34-206% as compared to motion uncorrected studies. Likewise, lesion detectability improved by respectively 115-136% and 62-235% with MR-based motion compensation as compared to gating and no motion correction and made it possible to distinguish non-transmural from transmural defects, which has clinical significance given the inherent limitations of current single modality imaging in identifying the amount of residual ischemia. The incorporation of PSF modeling within the framework of MR-based motion compensation significantly improved defect/myocardium contrast recovery (5.1-8.5%, p < 0.01) and defect detectability (39-56%, p < 0.01). No statistical difference was found in PET contrast and lesion detectability based on motion fields obtained with half and full k-space tagged data.

Noise properties of motion-compensated tomographic image reconstruction methods

Motion-compensated image reconstruction (MCIR) methods incorporate motion models to improve image quality in the presence of motion. MCIR methods differ in terms of how they use motion information and they have been well studied separately. However, there have been less theoretical comparisions of different MCIR methods. This paper compares the theoretical noise properties of three popular MCIR methods assuming known nonrigid motion. We show the relationship among three MCIR methods-motion-compensated temporal regularization (MTR), the parametric motion model (PMM), and post-reconstruction motion correction (PMC)-for penalized weighted least square cases. These analyses show that PMM and MTR are matrix-weighted sums of all registered image frames, while PMC is a scalar-weighted sum. We further investigate the noise properties of MCIR methods with Poisson models and quadratic regularizers by deriving accurate and fast variance prediction formulas using an "analytical approach." These theoretical noise analyses show that the variances of PMM and MTR are lower than or comparable to the variance of PMC due to the statistical weighting. These analyses also facilitate comparisons of the noise properties of different MCIR methods, including the effects of different quadratic regularizers, the influence of the motion through its Jacobian determinant, and the effect of assuming that total activity is preserved. Two-dimensional positron emission tomography simulations demonstrate the theoretical results.

Joint image and motion reconstruction for PET using a B-spline motion model

We present a novel joint image and motion reconstruction method for PET. The method is based on gated data and reconstructs an image together with a motion function. The motion function can be used to transform the reconstructed image to any of the input gates. All available events (from all gates) are used in the reconstruction. The presented method uses a B-spline motion model, together with a novel motion regularization procedure that does not need a regularization parameter (which is usually extremely difficult to adjust). Several image and motion grid levels are used in order to reduce the reconstruction time. In a simulation study, the presented method is compared to a recently proposed joint reconstruction method. While the presented method provides comparable reconstruction quality, it is much easier to use since no regularization parameter has to be chosen. Furthermore, since the B-spline discretization of the motion function depends on fewer parameters than a displacement field, the presented method is considerably faster and consumes less memory than its counterpart. The method is also applied to clinical data, for which a novel purely data-driven gating approach is presented.

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

Cardiac and respiratory movements pose significant challenges to image quality and quantitative accuracy in PET imaging. Cardiac and/or respiratory gating attempt to address this issue, but instead lead to enhanced noise levels. Direct four-dimensional (4D) PET image reconstruction incorporating motion compensation has the potential to minimize noise amplification while removing considerable motion blur. A wide-ranging choice of such techniques is reviewed in this work. Future opportunities and the challenges facing the adoption of 4D PET reconstruction and its role in basic and clinical research are also discussed.