Motion estimation and correction for simultaneous PET/MR using SIRF and CIL

SIRF is a powerful PET/MR image reconstruction research tool for processing data and developing new algorithms. In this research, new developments to SIRF are presented, with focus on motion estimation and correction. SIRF's recent inclusion of the adjoint of the resampling operator allows gradient propagation through resampling, enabling the MCIR technique. Another enhancement enabled registering and resampling of complex images, suitable for MRI. Furthermore, SIRF's integration with the optimization library CIL enables the use of novel algorithms. Finally, SPM is now supported, in addition to NiftyReg, for registration. Results of MR and PET MCIR reconstructions are presented, using FISTA and PDHG, respectively. These demonstrate the advantages of incorporating motion correction and variational and structural priors. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 2'

Motion estimation and correction for simultaneous PET/MR using SIRF and CIL

SIRF is a powerful PET/MR image reconstruction research tool for processing data and developing new algorithms. In this research, new developments to SIRF are presented, with focus on motion estimation and correction. SIRF's recent inclusion of the adjoint of the resampling operator allows gradient propagation through resampling, enabling the MCIR technique. Another enhancement enabled registering and resampling of complex images, suitable for MRI. Furthermore, SIRF's integration with the optimization library CIL enables the use of novel algorithms. Finally, SPM is now supported, in addition to NiftyReg, for registration. Results of MR and PET MCIR reconstructions are presented, using FISTA and PDHG, respectively. These demonstrate the advantages of incorporating motion correction and variational and structural priors. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 2'

Developments in PET-MRI for Radiotherapy Planning Applications

The hybridization of magnetic resonance imaging (MRI) and positron emission tomography (PET) provides the benefit of soft-tissue contrast and specific molecular information in a simultaneous acquisition. The applications of PET-MRI in radiotherapy are only starting to be realised. However, quantitative accuracy of PET relies on accurate attenuation correction (AC) of, not only the patient anatomy but also MRI hardware and current methods, which are prone to artefacts caused by dense materials. Quantitative accuracy of PET also relies on full characterization of patient motion during the scan. The simultaneity of PET-MRI makes it especially suited for motion correction. However, quality assurance (QA) procedures for such corrections are lacking. Therefore, a dynamic phantom that is PET and MR compatible is required. Additionally, respiratory motion characterization is needed for conformal radiotherapy of lung. 4D-CT can provide 3D motion characterization but suffers from poor soft-tissue contrast. In this thesis, I examine these problems, and present solutions in the form of improved MR-hardware AC techniques, a PET/MRI/CT-compatible tumour respiratory motion phantom for QA measurements, and a retrospective 4D-PET-MRI technique to characterise respiratory motion. Chapter 2 presents two techniques to improve upon current AC methods that use a standard helical CT scan for MRI hardware in PET-MRI. One technique uses a dual-energy computed tomography (DECT) scan to construct virtual monoenergetic image volumes and the other uses a tomotherapy linear accelerator to create CT images at megavoltage energies (1.0 MV) of the RF coil. The DECT-based technique reduced artefacts in the images translating to improved μ-maps. The MVCT-based technique provided further improvements in artefact reduction, resulting in artefact free μ-maps. This led to more AC of the breast coil. In chapter 3, I present a PET-MR-CT motion phantom for QA of motion-correction protocols. This phantom is used to evaluate a clinically available real-time dynamic MR images and a respiratory-triggered PET-MRI protocol. The results show the protocol to perform well under motion conditions. Additionally, the phantom provided a good model for performing QA of respiratory-triggered PET-MRI. Chapter 4 presents a 4D-PET/MRI technique, using MR sequences and PET acquisition methods currently available on hybrid PET/MRI systems. This technique is validated using the motion phantom presented in chapter 3 with three motion profiles. I conclude that our 4D-PET-MRI technique provides information to characterise tumour respiratory motion while using a clinically available pulse sequence and PET acquisition method

Uncertainty analysis of MR-PET image registration for precision neuro-PET imaging

Accurate regional brain quantitative PET measurements, particularly when using partial volume correction, rely on robust image registration between PET and MR images. We argue here that the precision, and hence the uncertainty, of MR-PET image registration is mainly driven by the registration implementation and the quality of PET images due to their lower resolution and higher noise compared to the structural MR images. We propose a dedicated uncertainty analysis for quantifying the precision of MR-PET registration, centred around the bootstrap resampling of PET list-mode events to generate multiple PET image realisations with different noise (count) levels. The effects of PET image reconstruction parameters, such as the use of attenuation and scatter corrections and different number of iterations, on the precision and accuracy of MR-PET registration were investigated. In addition, the performance of four software packages with their default settings for rigid inter-modality image registration were considered: NiftyReg, Vinci, FSL and SPM. Four distinct PET image distributions made of two early time frames (similar to cortical FDG) and two late frames using two amyloid PET dynamic acquisitions of one amyloid positive and one amyloid negative participants were investigated. For the investigated four PET frames, the biggest impact on the uncertainty was observed between registration software packages (up to 10-fold difference in precision) followed by the reconstruction parameters. On average, the lowest uncertainty for different PET frames and brain regions was observed with SPM and two iterations of fully quantitative image reconstruction. The observed uncertainty for the varying PET count-level (from 5% to 60%) was slightly lower than for the reconstruction parameters. We also observed that the registration uncertainty in quantitative PET analysis depends on amyloid status of the considered PET frames, with increased uncertainty (up to three times) when using post-reconstruction partial volume correction. This analysis is applicable for PET data obtained from either PET/MR or PET/CT scanners

Penalized PET/CT Reconstruction Algorithms With Automatic Realignment for Anatomical Priors

Two algorithms for solving misalignment issues in penalized PET/CT reconstruction using anatomical priors are proposed. Both approaches are based on a recently published joint motion estimation and image reconstruction method. The first approach deforms the anatomical image to align it with the functional one while the second approach deforms both images to align them with the measured data. Our current implementation alternates between alignment estimation and image reconstruction. We have chosen parallel level sets (PLSs) as a representative anatomical penalty, incorporating a spatially variant penalty strength. The performance was evaluated using simulated nontime-of-flight data generated with an XCAT phantom in the thorax region. We used the attenuation map in the anatomical prior. The results demonstrated that both methods can estimate the misalignment and deform the anatomical image accordingly. However, the performance of the first approach depends highly on the workflow of the alternating process. The second approach shows a faster convergence rate to the correct alignment and is less sensitive to the workflow. The presence of anatomical information improves the convergence rate of misalignment estimation for the second approach but slow it down for the first approach. Both approaches show improved performance in misalignment estimation as the data noise level decreases

Approaches for the optimization of MR protocols in clinical hybrid PET/MRI studies

Magnetic resonance imaging (MRI) is the examination method of choice for the diagnosis of a variety of diseases. MRI allows us to obtain not only anatomical information but also identification of physiological and functional parameters such as networks in the brain and tumor cellularity, which plays an increasing role in oncologic imaging, as well as blood flow and tissue perfusion. However, in many cases such as in epilepsy, degenerative neurological diseases and oncological processes, additional metabolic and molecular information obtained by PET can provide essential complementary information for better diagnosis. The combined information obtained from MRI and PET acquired in a single imaging session allows a more accurate localization of pathological findings and better assessment of the underlying physiopathology, thus providing a more powerful diagnostic tool. Two hundred and twenty-one patients were scanned from April 2011 to January 2012 on a Philips Ingenuity TF PET/MRI system. The purpose of this review article is to provide an overview of the techniques used for the optimization of different protocols performed in our hospital by specialists in the following fields: neuroradiology, head and neck, breast, and prostate imaging. This paper also discusses the different problems encountered, such as the length of studies, motion artifacts, and accuracy of image fusion including physical and technical aspects, and the proposed solution

MRI-Based Attenuation Correction in Emission Computed Tomography

The hybridization of magnetic resonance imaging (MRI) with positron emission tomography (PET) or single photon emission computed tomography (SPECT) enables the collection of an assortment of biological data in spatial and temporal register. However, both PET and SPECT are subject to photon attenuation, a process that degrades image quality and precludes quantification. To correct for the effects of attenuation, the spatial distribution of linear attenuation coefficients (μ-coefficients) within and about the patient must be available. Unfortunately, extracting μ-coefficients from MRI is non-trivial. In this thesis, I explore the problem of MRI-based attenuation correction (AC) in emission tomography. In particular, I began by asking whether MRI-based AC would be more reliable in PET or in SPECT. To this end, I implemented an MRI-based AC algorithm relying on image segmentation and applied it to phantom and canine emission data. The subsequent analysis revealed that MRI-based AC performed better in SPECT than PET, which is interesting since AC is more challenging in SPECT than PET. Given this result, I endeavoured to improve MRI-based AC in PET. One problem that required addressing was that the lungs yield very little signal in MRI, making it difficult to infer their μ-coefficients. By using a pulse sequence capable of visualizing lung parenchyma, I established a linear relationship between MRI signal and the lungs’ μ-coefficients. I showed that applying this mapping on a voxel-by-voxel basis improved quantification in PET reconstructions compared to conventional MRI-based AC techniques. Finally, I envisaged that a framework for MRI-based AC methods would potentiate further improvements. Accordingly, I identified three ways an MRI can be converted to μ-coefficients: 1) segmentation, wherein the MRI is divided into tissue types and each is assigned an μ-coefficient, 2) registration, wherein a template of μ-coefficients is aligned with the MRI, and 3) mapping, wherein a function maps MRI voxels to μ-coefficients. I constructed an algorithm for each method and catalogued their strengths and weaknesses. I concluded that a combination of approaches is desirable for MRI-based AC. Specifically, segmentation is appropriate for air, fat, and water, mapping is appropriate for lung, and registration is appropriate for bone

Improved UTE-based attenuation correction for cranial PET-MR using dynamic magnetic field monitoring.

Purpose: Ultrashort echo time (UTE) MRI has been proposed as a way to produce segmented attenuation maps for PET, as it provides contrast between bone, air, and soft tissue. However, UTE sequences require samples to be acquired during rapidly changing gradient fields, which makes the resulting images prone to eddy current artifacts. In this work it is demonstrated that this can lead to misclassification of tissues in segmented attenuation maps (AC maps) and that these effects can be corrected for by measuring the true k-space trajectories using a magnetic field camera. Methods: The k-space trajectories during a dual echo UTE sequence were measured using a dynamic magnetic field camera. UTE images were reconstructed using nominal trajectories and again using the measured trajectories. A numerical phantom was used to demonstrate the effect of reconstructing with incorrect trajectories. Images of an ovine leg phantom were reconstructed and segmented and the resulting attenuation maps were compared to a segmented map derived from a CT scan of the same phantom, using the Dice similarity measure. The feasibility of the proposed method was demonstrated in in vivo cranial imaging in five healthy volunteers. Simulated PET data were generated for one volunteer to show the impact of misclassifications on the PET reconstruction. Results: Images of the numerical phantom exhibited blurring and edge artifacts on the bone-tissue and air-tissue interfaces when nominal k-space trajectories were used, leading to misclassification of soft tissue as bone and misclassification of bone as air. Images of the tissue phantom and the in vivo cranial images exhibited the same artifacts. The artifacts were greatly reduced when the measured trajectories were used. For the tissue phantom, the Dice coefficient for bone in MR relative to CT was 0.616 using the nominal trajectories and 0.814 using the measured trajectories. The Dice coefficients for soft tissue were 0.933 and 0.934 for the nominal and measured cases, respectively. For air the corresponding figures were 0.991 and 0.993. Compared to an unattenuated reference image, the mean error in simulated PET uptake in the brain was 9.16% when AC maps derived from nominal trajectories was used, with errors in the SUV max for simulated lesions in the range of 7.17%-12.19%. Corresponding figures when AC maps derived from measured trajectories were used were 0.34% (mean error) and -0.21% to +1.81% (lesions). CONCLUSIONS: Eddy current artifacts in UTE imaging can be corrected for by measuring the true k-space trajectories during a calibration scan and using them in subsequent image reconstructions. This improves the accuracy of segmented PET attenuation maps derived from UTE sequences and subsequent PET reconstruction

Positron Emission Tomography: Current Challenges and Opportunities for Technological Advances in Clinical and Preclinical Imaging Systems

Positron emission tomography (PET) imaging is based on detecting two time-coincident high-energy photons from the emission of a positronemitting radioisotope. The physics of the emission, and the detection of the coincident photons, give PET imaging unique capabilities for both very high sensitivity and accurate estimation of the in vivo concentration of the radiotracer. PET imaging has been widely adopted as an important clinical modality for oncological, cardiovascular, and neurological applications. PET imaging has also become an important tool in preclinical studies, particularly for investigating murine models of disease and other small-animal models. However, there are several challenges to using PET imaging systems. These include the fundamental trade-offs between resolution and noise, the quantitative accuracy of the measurements, and integration with X-ray computed tomography and magnetic resonance imaging. In this article, we review how researchers and industry are addressing these challenges.This work was supported in part by National Institutes of Health grants R01-CA042593, U01-CA148131, R01CA160253, R01CA169072, and R01CA164371; by Human Frontier Science Program grant RGP0004/2013; and by the Innovative Medicines Initiative under grant agreement 115337, which comprises financial contributions from the European Union’s Seventh Framework Program (FP7/2007–2013