A major advantage of hybrid PET/MR systems is the radiation-free high spatial and temporal resolution of cardiac MR imaging that can be used to estimate respiratory and cardiac motion present during PET data acquisition. This information can be incorporated into algorithms to correct for the effects of motion in the PET data to reduce blurring and increase target-to-background ratios of PET hotspots. In this work, we demonstrate a method for respiratory and cardiac motion correction in patients with cardiac sarcoidosis.
A patient with proven extra-cardiac sarcoidosis and suspected cardiac involvement underwent hybrid PET/MR (Biograph mMR, Siemens) to test the feasibility of MR-based motion correction of the PET data. 18F-FDG (10 mCi) was administered. After patient set-up and initial image scouts, list-mode PET data acquisition began after 30 minutes and lasted for 60 minutes.
Motion Estimation: Respiratory motion was estimated using a free-breathing 3D golden-angle radial VIBE sequence (based on Siemens WIP-793, [2]). A whole-body coronal slab with 3x3x3mm-resolution and 1600 spokes was acquired over 6-7min. The amplitude of the center of k-space was extracted from the raw data and plotted over time. This quantity was taken to reflect the volume of liver in the field of view and subsequently to provide an estimate of the respiratory phase during the acquisition. Spokes were then divided into 4 respiratory frames based on signal amplitude and reconstructed in Matlab using NUFFT algorithms [3]. Motion vector fields (MVFs) between respiratory frames were then estimated using freely-available non-rigid registration algorithms [4].
Cardiac motion was estimated separately using a similar acquisition but with higher resolution of 1.4x1.4x1.4mm and a coronal slab just covering the heart. In addition, imaging proceeded during infusion of 0.2mmol/kg gadolinium-based contrast agent (Multihance, Bracco) to provide contrast between the blood-pool, coronary vessels and myocardium. Respiratory phase was estimated in the same manner before finding the actual head-to-foot displacement from the images. k-Space data for each respiratory phase were phase-shifted to correct for the head-to-foot displacement. Corrected k-space data were then sorted into 3 cardiac frames based on recorded ECG trigger timing (0-300, 300-600, 600-inf ms after trigger) before offline reconstruction. MVFs between cardiac frames were then estimated.
Motion-Phase Monitoring: To employ PET data acquired during the entire 60 minutes, respiratory phase was estimated from the variation of total PET counts over time, similarly to the method used for MR data. ECG trigger times were then used to double-sort all PET data into a matrix of 4 respiratory and 3 cardiac frames.
Motion Corrected PET Reconstruction: Motion correction of PET data employed the reconstruct-transform-average (RTA) approach. Double-gated list-mode PET data were reconstructed offline (e7tools, Siemens) using an iterative algorithm (OP-OSEM, 3 iterations, 21 subsets). Attenuation correction maps [5] of the body were transformed using the estimated respiratory MVFs. Finally, RTA was performed in two steps. All cardiac frames were transformed to the end expiration position using corresponding MVFs. Then each respiratory-motion-corrected cardiac frame was transformed to the diastolic position using the cardiac MVFs before averaging all frames. A flow diagram of the steps involved is given in Fig. 1.
Image analysis: Non-motion corrected (non-MC), and double respiratory and cardiac motion corrected (MC) 18F-FDG-PET images were evaluated qualitatively for image blurring and overall quality and by finding target-to-background ratio (TBR) in cardiac hotspots identified and compared to non-corrected images (blood in the left atrium was used to estimate background).