Investigation of practical initial attenuation image estimates in TOF-MLAA reconstruction for PET/MR

J. C. Cheng, A. Salomon, M. Yaqub, R. Boellaard

Research output: Contribution to journalArticleAcademicpeer-review

Abstract

Purpose: Time-of-flight joint attenuation and activity positron emission tomography reconstruction requires additional calibration (scale factors) or constraints during or post-reconstruction to produce a quantitative mu-map. In this work, the impact of various initializations of the joint reconstruction was investigated, and the initial average mu-value (IAM) method was introduced such that the forward-projection of the initial mu-map is already very close to that of the reference mu-map, thus reducing/minimizing the offset (scale factor) during the early iterations of the joint reconstruction. Consequently, the accuracy and efficiency of unconstrained joint reconstruction such as time-of-flight maximum likelihood estimation of attenuation and activity (TOF-MLAA) can be improved by the proposed IAM method. Methods: 2D simulations of brain and chest were used to evaluate TOF-MLAA with various initial estimates which include the object filled with water uniformly (conventional initial estimate), bone uniformly, the average mu-value uniformly (IAM magnitude initialization method), and the perfect spatial mu-distribution but with a wrong magnitude (initialization in terms of distribution). 3D GATE simulation was also performed for the chest phantom under a typical clinical scanning condition, and the simulated data were reconstructed with a fully corrected list-mode TOF-MLAA algorithm with various initial estimates. The accuracy of the average mu-values within the brain, chest, and abdomen regions obtained from the MR derived mu-maps was also evaluated using computed tomography mu-maps as the gold-standard. Results: The estimated mu-map with the initialization in terms of magnitude (i.e., average mu-value) was observed to reach the reference more quickly and naturally as compared to all other cases. Both 2D and 3D GATE simulations produced similar results, and it was observed that the proposed IAM approach can produce quantitative mu-map/emission when the corrections for physical effects such as scatter and randoms were included. The average mu-value obtained from MR derived mu-map was accurate within 5% with corrections for bone, fat, and uniform lungs. Conclusions: The proposed IAM-TOF-MLAA can produce quantitative mu-map without any calibration provided that there are sufficient counts in the measured data. For low count data, noise reduction and additional regularization/rescaling techniques need to be applied and investigated. The average mu-value within the object is prior information which can be extracted from MR and patient database, and it is feasible to obtain accurate average mu-value using MR derived mu-map with corrections as demonstrated in this work. (C) 2016 American Association of Physicists in Medicine.
Original languageEnglish
Pages (from-to)4163-4173
Number of pages11
JournalMedical Physics
Volume43
DOIs
Publication statusPublished - 2016

Cite this

@article{3e0db865beeb498bbe6064a259c53dd5,
title = "Investigation of practical initial attenuation image estimates in TOF-MLAA reconstruction for PET/MR",
abstract = "Purpose: Time-of-flight joint attenuation and activity positron emission tomography reconstruction requires additional calibration (scale factors) or constraints during or post-reconstruction to produce a quantitative mu-map. In this work, the impact of various initializations of the joint reconstruction was investigated, and the initial average mu-value (IAM) method was introduced such that the forward-projection of the initial mu-map is already very close to that of the reference mu-map, thus reducing/minimizing the offset (scale factor) during the early iterations of the joint reconstruction. Consequently, the accuracy and efficiency of unconstrained joint reconstruction such as time-of-flight maximum likelihood estimation of attenuation and activity (TOF-MLAA) can be improved by the proposed IAM method. Methods: 2D simulations of brain and chest were used to evaluate TOF-MLAA with various initial estimates which include the object filled with water uniformly (conventional initial estimate), bone uniformly, the average mu-value uniformly (IAM magnitude initialization method), and the perfect spatial mu-distribution but with a wrong magnitude (initialization in terms of distribution). 3D GATE simulation was also performed for the chest phantom under a typical clinical scanning condition, and the simulated data were reconstructed with a fully corrected list-mode TOF-MLAA algorithm with various initial estimates. The accuracy of the average mu-values within the brain, chest, and abdomen regions obtained from the MR derived mu-maps was also evaluated using computed tomography mu-maps as the gold-standard. Results: The estimated mu-map with the initialization in terms of magnitude (i.e., average mu-value) was observed to reach the reference more quickly and naturally as compared to all other cases. Both 2D and 3D GATE simulations produced similar results, and it was observed that the proposed IAM approach can produce quantitative mu-map/emission when the corrections for physical effects such as scatter and randoms were included. The average mu-value obtained from MR derived mu-map was accurate within 5{\%} with corrections for bone, fat, and uniform lungs. Conclusions: The proposed IAM-TOF-MLAA can produce quantitative mu-map without any calibration provided that there are sufficient counts in the measured data. For low count data, noise reduction and additional regularization/rescaling techniques need to be applied and investigated. The average mu-value within the object is prior information which can be extracted from MR and patient database, and it is feasible to obtain accurate average mu-value using MR derived mu-map with corrections as demonstrated in this work. (C) 2016 American Association of Physicists in Medicine.",
author = "Cheng, {J. C.} and A. Salomon and M. Yaqub and R. Boellaard",
note = "M1 - 7 ISI Document Delivery No.: DQ4KE Times Cited: 0 Cited Reference Count: 23 Cheng, Ju-Chieh (Kevin) Salomon, Andre Yaqub, Maqsood Boellaard, Ronald Boellaard, Ronald/0000-0002-0313-5686 Netherlands Organization for Health Research and Development IMDI Grant [104002006] This work was supported by the Netherlands Organization for Health Research and Development IMDI Grant No. 104002006 and was presented in part at the annual Society of Nuclear Medicine and Molecular Imaging Meeting, Baltimore, MD, June 2015 and at the IEEE Medical Imaging Conference, San Diego, CA, November 2015. 0 1 6 AMER ASSOC PHYSICISTS MEDICINE AMER INST PHYSICS MELVILLE MED PHYS",
year = "2016",
doi = "10.1118/1.4953634",
language = "English",
volume = "43",
pages = "4163--4173",
journal = "Medical Physics",
issn = "0094-2405",
publisher = "AAPM - American Association of Physicists in Medicine",

}

Investigation of practical initial attenuation image estimates in TOF-MLAA reconstruction for PET/MR. / Cheng, J. C.; Salomon, A.; Yaqub, M.; Boellaard, R.

In: Medical Physics, Vol. 43, 2016, p. 4163-4173.

Research output: Contribution to journalArticleAcademicpeer-review

TY - JOUR

T1 - Investigation of practical initial attenuation image estimates in TOF-MLAA reconstruction for PET/MR

AU - Cheng, J. C.

AU - Salomon, A.

AU - Yaqub, M.

AU - Boellaard, R.

N1 - M1 - 7 ISI Document Delivery No.: DQ4KE Times Cited: 0 Cited Reference Count: 23 Cheng, Ju-Chieh (Kevin) Salomon, Andre Yaqub, Maqsood Boellaard, Ronald Boellaard, Ronald/0000-0002-0313-5686 Netherlands Organization for Health Research and Development IMDI Grant [104002006] This work was supported by the Netherlands Organization for Health Research and Development IMDI Grant No. 104002006 and was presented in part at the annual Society of Nuclear Medicine and Molecular Imaging Meeting, Baltimore, MD, June 2015 and at the IEEE Medical Imaging Conference, San Diego, CA, November 2015. 0 1 6 AMER ASSOC PHYSICISTS MEDICINE AMER INST PHYSICS MELVILLE MED PHYS

PY - 2016

Y1 - 2016

N2 - Purpose: Time-of-flight joint attenuation and activity positron emission tomography reconstruction requires additional calibration (scale factors) or constraints during or post-reconstruction to produce a quantitative mu-map. In this work, the impact of various initializations of the joint reconstruction was investigated, and the initial average mu-value (IAM) method was introduced such that the forward-projection of the initial mu-map is already very close to that of the reference mu-map, thus reducing/minimizing the offset (scale factor) during the early iterations of the joint reconstruction. Consequently, the accuracy and efficiency of unconstrained joint reconstruction such as time-of-flight maximum likelihood estimation of attenuation and activity (TOF-MLAA) can be improved by the proposed IAM method. Methods: 2D simulations of brain and chest were used to evaluate TOF-MLAA with various initial estimates which include the object filled with water uniformly (conventional initial estimate), bone uniformly, the average mu-value uniformly (IAM magnitude initialization method), and the perfect spatial mu-distribution but with a wrong magnitude (initialization in terms of distribution). 3D GATE simulation was also performed for the chest phantom under a typical clinical scanning condition, and the simulated data were reconstructed with a fully corrected list-mode TOF-MLAA algorithm with various initial estimates. The accuracy of the average mu-values within the brain, chest, and abdomen regions obtained from the MR derived mu-maps was also evaluated using computed tomography mu-maps as the gold-standard. Results: The estimated mu-map with the initialization in terms of magnitude (i.e., average mu-value) was observed to reach the reference more quickly and naturally as compared to all other cases. Both 2D and 3D GATE simulations produced similar results, and it was observed that the proposed IAM approach can produce quantitative mu-map/emission when the corrections for physical effects such as scatter and randoms were included. The average mu-value obtained from MR derived mu-map was accurate within 5% with corrections for bone, fat, and uniform lungs. Conclusions: The proposed IAM-TOF-MLAA can produce quantitative mu-map without any calibration provided that there are sufficient counts in the measured data. For low count data, noise reduction and additional regularization/rescaling techniques need to be applied and investigated. The average mu-value within the object is prior information which can be extracted from MR and patient database, and it is feasible to obtain accurate average mu-value using MR derived mu-map with corrections as demonstrated in this work. (C) 2016 American Association of Physicists in Medicine.

AB - Purpose: Time-of-flight joint attenuation and activity positron emission tomography reconstruction requires additional calibration (scale factors) or constraints during or post-reconstruction to produce a quantitative mu-map. In this work, the impact of various initializations of the joint reconstruction was investigated, and the initial average mu-value (IAM) method was introduced such that the forward-projection of the initial mu-map is already very close to that of the reference mu-map, thus reducing/minimizing the offset (scale factor) during the early iterations of the joint reconstruction. Consequently, the accuracy and efficiency of unconstrained joint reconstruction such as time-of-flight maximum likelihood estimation of attenuation and activity (TOF-MLAA) can be improved by the proposed IAM method. Methods: 2D simulations of brain and chest were used to evaluate TOF-MLAA with various initial estimates which include the object filled with water uniformly (conventional initial estimate), bone uniformly, the average mu-value uniformly (IAM magnitude initialization method), and the perfect spatial mu-distribution but with a wrong magnitude (initialization in terms of distribution). 3D GATE simulation was also performed for the chest phantom under a typical clinical scanning condition, and the simulated data were reconstructed with a fully corrected list-mode TOF-MLAA algorithm with various initial estimates. The accuracy of the average mu-values within the brain, chest, and abdomen regions obtained from the MR derived mu-maps was also evaluated using computed tomography mu-maps as the gold-standard. Results: The estimated mu-map with the initialization in terms of magnitude (i.e., average mu-value) was observed to reach the reference more quickly and naturally as compared to all other cases. Both 2D and 3D GATE simulations produced similar results, and it was observed that the proposed IAM approach can produce quantitative mu-map/emission when the corrections for physical effects such as scatter and randoms were included. The average mu-value obtained from MR derived mu-map was accurate within 5% with corrections for bone, fat, and uniform lungs. Conclusions: The proposed IAM-TOF-MLAA can produce quantitative mu-map without any calibration provided that there are sufficient counts in the measured data. For low count data, noise reduction and additional regularization/rescaling techniques need to be applied and investigated. The average mu-value within the object is prior information which can be extracted from MR and patient database, and it is feasible to obtain accurate average mu-value using MR derived mu-map with corrections as demonstrated in this work. (C) 2016 American Association of Physicists in Medicine.

U2 - 10.1118/1.4953634

DO - 10.1118/1.4953634

M3 - Article

VL - 43

SP - 4163

EP - 4173

JO - Medical Physics

JF - Medical Physics

SN - 0094-2405

ER -