- Open Access
Motion Correction of Whole-Body PET Data with a Joint PET-MRI Registration Functional
© Fieseler et al.; licensee BioMed Central Ltd. 2014
- Published: 28 February 2014
Respiratory motion is known to degrade image quality in PET imaging. The necessary acquisition time of several minutes per bed position will inevitably lead to a blurring effect due to organ motion. A lot of research has been done with regards to motion correction of PET data. As full-body PET-MRI became available recently, the anatomical data provided by MRI is a promising source of motion information. Current PET-MRI-based motion correction approaches, however, do not take into account the available information provided by PET data. PET data, though, may add valuable additional information to increase motion estimation robustness and precision.
In this work we propose a registration functional that is capable of performing motion detection in gated data of two modalities simultaneously. Evaluation is performed using phantom data. We demonstrate that performing a joint registration of both modalities does improve registration accuracy and PET image quality.
- Motion Correction
- Image Registration
- Joint Registration
- Multi-Modal Data
- Positron Emission Tomography
- Magnetic Resonance Imaging
Respiratory motion is known to impair image quality as well as quantification in positron emission tomography (PET) . As the acquisition of PET takes several minutes per bed position, organ motion due to respiration cannot be avoided and will thus result in blurred images. By using gating methods, the acquired PET data can be divided into different motion phases. Gating reduces the amount of motion contained within each gate to a large extent, yet at the expense of reduced statistics and thereby image quality [2, 3]. To alleviate this, motion between gates can be estimated and subsequently be used to correct PET data for motion, resulting in a single image volume with reduced motion artifacts and full statistics. Various approaches for motion estimation of gated PET data have been studied, including optical flow , B-spline based methods , and registration methods including mass-preservation .
The use of 4D-CT data for motion correction has been proposed . The advantage of this approach lies in the usage of anatomical data, which is independent of tracer uptake. Acquisition of 4D-CT data, however, increases the radiation burden for the patient.
Whole body PET-MRI is promising regarding PET motion correction. High resolution MR data may allow for a precise motion estimation independent of tracer uptake and without additional radiation burden for the patient. The feasibility of PET-MRI-based motion correction of PET data has been demonstrated already, e.g., using hardware phantoms , animals , and simulation data . A commonality of current approaches to MR-based motion correction is that a 4D MR dataset is acquired from which motion is estimated and subsequently used in the correction of PET data. Strategies for the generation of 4D MR include acquisition of 2D slices with subsequent reordering , fast, consecutive acquisitions of 3D volumes , and sorting of k-space data during or after acquisition [13, 14]. Fayad et al. propose an approach where motion and image data are estimated simultaneously from the acquired MR data .
The motion information contained in PET data, however, remains unused in these approaches. In clinical routine, time for the acquisition of MR data needed for motion correction may be limited, as clinical protocols may demand for further, diagnostic sequences . Since time is proportional to image quality in MR, limited acquisition time may not allow to exploit the full potential of MR, resulting in poorer image quality than technically possible. Accordingly, all available information for motion estimation, including PET data, should be used.
Further, both modalities may contribute valuable information to the motion detection process. In MR, e.g., the lungs give relatively little signal due to their low proton density . Integrating information from PET data, if, e.g., active lesions are present in the lungs, may benefit motion estimation. Using information from both modalities should result in more reliable registration results. In the present work we propose an approach that uses information from both modalities by combining two registration functionals into a joint functional.
For motion correction approaches in PET-MRI, proper synchronisation of MR and PET is mandatory. The motion determined from MR has to be related to the PET data with respect to time. For the approach described in the following, we assume a gated PET dataset as well as a series of MR datasets. We assume that each PET gate corresponds to one MR gate with respect to its motion phase. The feasibility of creating corresponding phases has been demonstrated, e.g., in .
In the following we describe the proposed registration functional, followed by a description of the phantom data used in this work.
where is a distance functional, is the reference volume, is the template volume to be registered, and is a regularizer penalizing unfavourable transformations. The scalar value α weights the influence of the regularizer. For non-rigid transformations, y is chosen as a non-parametric transformation (one vector per voxel).
Here, and denote two reference volumes and and the template volumes. The scalar value β allows to weight the influence of the data term for PET.
In the registration functional in Equation (2), the deformation is represented by a common grid y for both modalities. Since the input data will not necessarily share the same resolution, resampling to a common grid is performed. Here, we chose a common grid of 2 mm3 voxel size.
where yref is a reference grid, given by the identity transformation in our case. The three summands control changes in length, surface area, and volume. Parameters for the hyperelastic regularizer were chosen empirically as α l = 1, α a = 0.1, α v = 1. Throughout all experiments in this paper, we keep α l , α a , α v fixed and vary the regularization strength by changing α in Eq. (2). For further details regarding the regularizer we refer to .
As the distance functional , we choose the sum of squared differences (SSD) for both PET and MR. The functional was implemented using the Matlab-based FAIR toolbox . For all registration experiments, linear interpolation was used and a multi-level approach using a downscaling factor of 0.5 was applied. Optimiziation was performed using a Gauss-Newton scheme with a preconditioned conjugate gradient solver.
In the present work we use data generated using a software phantom for evaluation. This allows us to compare motion estimates against ground-truth motion data. The XCAT phantom  is widely used as the basis for the simulation of imaging modalities and the evaluation of correction methods, e.g., in [6, 20, 21]. Using the XCAT phantom, we created an artificial PET-MRI dataset as described in the following.
The XCAT phantom was set to a maximum diaphragm motion of 2 cm. We selected eight frames representing the full range from inspiration to expiration .
For the creation of MR data a labelled XCAT dataset of 1 mm3 resolution was created. This dataset and known tissue values for T 1, T 2, and proton density [23, 24] were used as an input to the freely available MR simulation software SIMRI . We simulated an MR acquisition of stacked 2D slices covering the thorax. For respiratory motion, the largest amount of motion can be expected in the cranio-caudal and anterior-posterior directions. Thus, a sagittal slice orientation was chosen to capture these directions in-plane. For the phantom dataset used here, the maximum extents of motion are 2.5 mm (left-right), 10.9 mm (anterior-posterior), and 25.2 mm (head-feet). The following MR parameters were used: gradient echo, TE/TR 10 ms/30 ms, 12° flip angle, 256 × 256 pixels, 2 mm pixel size in-plane, slice spacing 1 cm, slice thickness 1 cm.
Using the phantom data described above allows for evaluation based on ground-truth motion as well as ground-truth activity data. First, we set β = 0, thereby performing an MR-based registration and use varying values of α to determine the best result for an MR-based registration.
Using the determined registration parameters, we compute registrations for increasing values of β, thus adding increasing amounts of PET information to the registration.
where Ω is the image domain, y i the i-th component of vector y, and yGT the ground-truth vector. Here, we examine averaged values of the AEE for all gates.
Further, the computed motion estimates are used to perform a motion correction of the dataset. The PET gates are warped using the computed motion and averaged. We evaluate correlation values of the corrected images as well as the recovered activities in the three lesion regions. For all registrations, attenuated PET data are used. For evaluation of correlation coefficients and activity recovery, attenuation corrected PET data are used.
Recovered Activity in Lesions
Using the joint registration approach, local improvements are observable. For the heart, the lungs, and the lesions (lungs and liver), improvements in terms of registration error (AEE) are achieved. In particular the lesions in the lungs show a large reduction in registration error when PET data are added to the registration. For the liver region, an increase of the AEE is observable. For the area of the lesion added to the liver, though, the AEE is decreased. With the exception of the liver region in toto, those regions that exhibit tracer uptake in PET seem to contribute to a better motion estimation result. Correlation of PET data is improved if PET data is added to the registration, indicating a benefit. Additionally, the lower registration error leads to a slightly better recovery of the activity in the lesions. Globally, a slight increase of the AEE is observable with increasing values of β.
The results presented here indicate that a benefit of using motion information from both modalities, MR and PET, is achievable. Certainly, the extent of the benefit will depend on many factors, with one major factor being the image quality of the two image modalities. The amount of remaining motion-induced blurring within the individual MR and PET gates will certainly limit the precision to which motion can be estimated. The extent to which PET can contribute will as well depend on the image quality, determined by factors like, e.g., type of injected tracer, injected dose, tracer uptake, and acquisition time. The phantom data used here does not contain artifacts other than noise. Particularly, we did not simulate motion artifacts which are likely to occur during MR acquisitions.
We have presented a joint registration functional that makes use of motion information derived from PET and MR data simultaneously. In this approach, motion information from both modalities is used. We demonstrated that the proposed method leads to a lower local registration error and better recovery of lesion activity, thus using information from both modalities simultaneously is beneficial regarding motion correction. As a result, clinical scenarios involving lesion quantification might in particular benefit from the proposed method.
In future work, we will evaluate our approach for a broader range of data, including cardiac motion. This will include the addition of mass-preservation . Certainly, the phantom data used in this work does only approximate reality. Thus, we will evaluate the proposed approach on real phantom data.
This work was partly funded by the German Research Foundation (Deutsche Forschungsgemeinschaft), Sonderforschungsbereich SFB 656 MoBil (project B3), and a research grant from Siemens Healthcare, Erlangen, Germany. We acknowledge support by Deutsche Forschungsgemeinschaft and Open Access Publication Fund of University of Muenster.
This article has been published as part of BioMedical Engineering OnLine Volume 13 Supplement 1, 2014: Selected articles from the 35th Annual International Conference of the IEEE Engineering in Medicine and Biology Society: Workshop on Current Challenging Image Analysis and Information Processing in Life Sciences. The full contents of the supplement are available online at http://www.biomedical-engineering-online.com/supplements/13/S1
- Nehmeh S, Erdi Y, Ling C, Rosenzweig K, Squire O, Braban L, Ford E, Sidhu K, Mageras G, Larson S, et al.: Effect of respiratory gating on reducing lung motion artifacts in pet imaging of lung cancer. Medical physics 2002, 29: 366. 10.1118/1.1448824View ArticleGoogle Scholar
- Büther F, Ernst I, Dawood M, Kraxner P, Schäfers M, Schober O, Schäfers KP: Detection of respiratory tumour motion using intrinsic list mode-driven gating in positron emission tomography. European Journal of Nuclear Medicine and Molecular Imaging 2010, 37: 2315–2327. 10.1007/s00259-010-1533-yView ArticleGoogle Scholar
- Dawood M, Buether F, Lang N, Schober O, Schaefers KP: Respiratory gating in positron emission tomography: A quantitative comparison of different gating schemes. Medical Physics 2007,34(7):3067. 10.1118/1.2748104View ArticleGoogle Scholar
- Dawood M, Büther F, Jiang X, Schäfers KP: Respiratory motion correction in 3-D PET data with advanced optical flow algorithms. IEEE Transactions on Medical Imaging 2008,27(8):1164–1175.View ArticleGoogle Scholar
- Bai W, Brady M: Regularized b-spline deformable registration for respiratory motion correction in PET images. Physics in Medicine and Biology 2009,54(9):2719. 10.1088/0031-9155/54/9/008View ArticleGoogle Scholar
- Gigengack F, Ruthotto L, Burger M, Wolters CH, Jiang X, Schäfers KP: Motion correction in dual gated cardiac PET using mass-preserving image registration. IEEE Transactions on Medical Imaging 2012,31(3):698–712.View ArticleGoogle Scholar
- Lamare F, Carbayo MJL, Cresson T, Kontaxakis G, Santos A, Rest CCL, Reader AJ, Visvikis D: List-mode-based reconstruction for respiratory motion correction in PET using non-rigid body transformations. Physics in Medicine and Biology 2007,52(17):5187. 10.1088/0031-9155/52/17/006View ArticleGoogle Scholar
- Tsoumpas C, Mackewn J, Halsted P, King A, Buerger C, Totman J, Schaeffter T, Marsden P: Simultaneous PET-MR acquisition and MR-derived motion fields for correction of non-rigid motion in PET. Annals of Nuclear Medicine 2010,24(10):745–750. 10.1007/s12149-010-0418-2View ArticleGoogle Scholar
- Chun SY, Reese TG, Ouyang J, Guerin B, Catana C, Zhu X, Alpert NM, El Fakhri G: MRI-based nonrigid motion correction in simultaneous PET/MRI. Journal of Nuclear Medicine 2012,53(8):1284–1291. 10.2967/jnumed.111.092353View ArticleGoogle Scholar
- Dikaios N, Izquierdo-Garcia D, Graves M, Mani V, Fayad Z, Fryer T: MRI-based motion correction of thoracic PET: initial comparison of acquisition protocols and correction strategies suitable for simultaneous PET/MRI systems. European Radiology 2012, 22: 439–446. 10.1007/s00330-011-2274-4View ArticleGoogle Scholar
- Würslin C, Schmidt H, Martirosian P, Boss A, Schwenzer NF, Stegger L: Respiratory motion correction in oncologic PET using t1-weighted MR imaging on a simultaneous whole-body PET/MR system. Journal of Nuclear Medicine 2013,54(3):464–471. 10.2967/jnumed.112.105296View ArticleGoogle Scholar
- Tsoumpas C, Buerger C, King AP, Mollet P, Keereman V, Vandenberghe S, Schulz V, Schleyer P, Schaeffter T, Marsden PK: Fast generation of 4D PET-MR data from real dynamic MR acquisitions. Physics in Medicine and Biology 2011,56(20):6597. 10.1088/0031-9155/56/20/005View ArticleGoogle Scholar
- Ouyang J, Li Q, Fakhri GE: Magnetic resonance-based motion correction for positron emission tomography imaging. Seminars in Nuclear Medicine 2013,43(1):60–67. 10.1053/j.semnuclmed.2012.08.007View ArticleGoogle Scholar
- Dutta J, El Fakhri G, Huang C, Petibon Y, Reese TG, Li Q: Respiratory motion compensation in simultaneous PET/MR using a maximum a posteriori approach. In Biomedical Imaging (ISBI), 2013 IEEE 10th International Symposium On. IEEE; 2013:800–803.View ArticleGoogle Scholar
- Fayad H, Odille F, Felblinger J, Visvikis D: A generic PET/MRI respiratory motion correction using a generalized reconstruction by inversion of coupled systems (GRICS) approach. In Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC), 2012 IEEE. IEEE; 2012:2813–2816.View ArticleGoogle Scholar
- Bailey DL, Beyer T, Boellaard R, Gückel B, Hellwig D, Herzog H, Pichler BJ, Quick HH, Sabri O, et al.: Summary report of the first international workshop on PET/MR imaging, march 19–23, 2012, tübingen, germany. Molecular Imaging and Biology 2013, 1–11.Google Scholar
- Burger M, Modersitzki J, Ruthotto L: A hyperelastic regularization energy for image registration. SIAM Journal on Scientific Computing 2012.Google Scholar
- Modersitzki J: Fair: Flexible Algorithms for Image Registration (Fundamentals of Algorithms). Society for Industrial and Applied Mathematics, Philadelphia 2009.Google Scholar
- Segars WP, Sturgeon G, Mendonca S, Grimes J, Tsui BM: 4D XCAT phantom for multimodality imaging research. Medical physics 2010,37(9):4902–4915. 10.1118/1.3480985View ArticleGoogle Scholar
- Xia T, Alessio AM, De Man B, Manjeshwar R, Asma E, Kinahan PE: Ultra-low dose CT attenuation correction for PET/CT. Physics in Medicine and Biology 2012,57(2):309. 10.1088/0031-9155/57/2/309View ArticleGoogle Scholar
- Tabary J, Marache-Francisco S, Valette S, Segars WP, Lartizien C: Realistic X-ray CT simulation of the XCAT phantom with SINDBAD. Nuclear Science Symposium Conference Record (NSS/MIC), 2009 IEEE 2009, 3980–3983.View ArticleGoogle Scholar
- Dawood M, Büther F, Stegger L, Jiang X, Schober O, Schäfers M, Schäfers KP: Optimal number of respiratory gates in positron emission tomography: A cardiac patient study. Medical physics 2009, 36: 1775. 10.1118/1.3112422View ArticleGoogle Scholar
- Stanisz GJ, Odrobina EE, Pun J, Escaravage M, Graham SJ, Bronskill MJ, Henkelman RM: T1, T2 relaxation and magnetization transfer in tissue at 3T. Magnetic Resonance in Medicine 2005,54(3):507–512. 10.1002/mrm.20605View ArticleGoogle Scholar
- Tobon-Gomez C, Sukno F, Bijnens BH, Huguet M, Frangi A: Realistic simulation of local image appearance of cardiac magnetic resonance imaging using a virtual phantom population. CI2BM09 - MICCAI Workshop on Cardiovascular Interventional Imaging and Biophysical Modelling, London, Royaume-Uni 2009, 132–140.Google Scholar
- Benoit-Cattin H, Collewet G, Belaroussi B, Saint-Jalmes H, Odet C: The SIMRI project: a versatile and interactive MRI simulator. Journal of Magnetic Resonance 2005,173(1):97–115. 10.1016/j.jmr.2004.09.027View ArticleGoogle Scholar
- Kösters T, Schäfers KP, Wübbeling F: EMrecon: An expectation maximization based image reconstruction framework for emission tomography data. Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC), 2011 IEEE 2011, 4365–4368.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.