The Effect of MR-based Motion Correction on PET Kinetic Parameters Estimation
Rong Guo1, Yoann Petibon2, Yixin Ma1, Kui Ying1, and Jinsong Ouyang2

1Engineering Physics, Tsinghua University, Beijing, China, People's Republic of, 2Center for Advanced Radiological Sciences, Massachusetts General Hospital, Boston, MA, United States

Synopsis

Bias may be introduced in the estimation of the PET myocardial kinetic parameters by both cardiac and respiratory motion. Simultaneous PET-MR makes it possible to perform MR-based PET motion correction. We have investigated the performance of MR-based motion correction on the estimation of myocardial PET kinetic parametermsat

TARGET AUDIENCE

People who are interested in cardiac PET-MR.

PURPOSE

Both cardiac and respiratory motion may introduce bias in the estimation of the PET myocardial kinetic parameters. Simultaneous PET-MR makes it possible to perform MR-based PET motion correction. We have investigated the performance of MR-based motion correction on the estimation of myocardial PET kinetic parameters.

METHOD

Fig.1 shows the flowchart of PET-MR simulation, image reconstruction, and estimation of PET kinetic parameters. We performed simulation studies using an XCAT torso phantom [1], which includes heart, lungs, liver, and soft-tissue compartments and generates both cardiac and respiratory motion fields. The time activity curves of all the compartments were simulated according to a one-tissue compartmental model with realistic kinetic parameters and arterial input function from previously reported human ammonia perfusion studies [2]. The simulation data were equivalent to a 9-min dynamic PET scan with 0.5 mCi injection dose and the framing scheme of 8×5 sec, 4×10 sec, 2×20 sec, 1×40 sec, 1×2 min and 1×4 min. We simulated a defect in the left ventricle myocardium by lowering K1 and k2 values by 60% and 28%, respectively, of the original values. PET sinograms for each time frame were generated using a forward projection model, which incorporates attenuation, point spread function, and Poisson noise. Both scatter and random events were not accounted in the simulations. MRI simulation was performed using MRILAB [3] at 3T with published the T1, T2, and spin density values and standard GRE sequence. Three different types of motion were investigated: cardiac motion only, respiratory motion only, and both cardiac and respiratory motion. For each type of motion, a reference motion phase was selected. The motion fields transforming from a given motion phase to the reference motion phase were obtained by applying non-rigid Demon registration algorithm to the reconstructed MR images. For each time frame, PET image was reconstructed using Filtered Back Projection (FBP) for each motion phase. The resulting reconstructed PET images were then transformed to the reference motion phase using the motion fields. The motion corrected PET image for the reference motion phase in the time frame was then obtained by summing up all the transformed PET images. Finally, voxel-wise kinetic parameters were estimated using the motion corrected PET images by curve-fitting.

RESULTS and DISCUSSION

For each type of motion, Fig.2 shows the reference (i.e. gating, only the PET events within the reference motion phase were used), motion-uncorrected, and motion- corrected static PET images. Fig.3 shows the estimated K1 maps for the reference motion phase and each type of motion. Tab.1 shows the estimated kinetic parameters in different parts of the myocardium. Our results show that MR-based motion correction reduces motion blurring effects as well as the bias of the estimated kinetic parameters. Moreover, respiratory motion appears to contribute more to the blurring effect and the bias of the estimated kinetic parameters than cardiac motion if no motion correction is applied.

CONCLUSION

We have shown that MR-based motion correction reduces the bias of estimated kinetic parameters in dynamic PET.

Acknowledgements

No acknowledgement found.

References

[1] W. P. Segarsa, etc. 4D XCAT phantom for multimodality imaging research. Med. Phys. 37 (9), Sep 2010.

[2] Otto Muzik .etc, Validation of Nitrogen- 13-Ammonia Tracer Kinetic Model for Quantification of Myocardial Blood Flow Using PET. J NucIMed 1993; 34:83-91.

[3]D.Kroon. http://www.mathworks.com/matlabcentral/fileexchange/21451-multimodality-non-rigid-demon-algorithm-image-registration

Figures

Fig.1 Flowchart of the study process including PET-MRI simulation, image reconstruction, motion correction and curve fitting

Fig.2 Reconstructed static PET images with gating (i.e., for the reference motion phase), motion-uncorrected, and motion corrected methods applied to three different types of motion, which include cardiac motion only, respiratory motion only, and both cardiac and respiratory motion.

Fig.3 Estimated K1 maps with gating (i.e., for the reference motion phase), motion-uncorrected, and motion corrected methods applied to three different types of motion, which include cardiac motion only, respiratory motion only, and both cardiac and respiratory motion.

Tab.1 Estimated kinetic parameters. Card: cardiac motion only; Resp: respiratory motion only; Both: both cardiac and respiratory motion; Bott: bottom of the myocardium; Left: left ventricle myocardium; Defe: defect. Kinetic parameters used in the simulation: K1=0.8, k2=0.14 for myocardium; K1=0.316, k2=0.1 for defect.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
1868