MR-based attenuation correction for PET/MRI neurological studies with continuous-valued attenuation coefficients for bone through a conversion from R2 to CT-Hounsfield units

• Published in: NeuroImage (Orlando, Fla.) Vol. 112; pp. 160 - 168
• Main Authors: Juttukonda, Meher R., Mersereau, Bryant G., Chen, Yasheng, Su, Yi, Rubin, Brian G., Benzinger, Tammie L.S., Lalush, David S., An, Hongyu
• Format: Journal Article
• Language: English
• Published: United States Elsevier Inc 15.05.2015; Elsevier Limited
• Subjects: Accuracy; Adipose Tissue - anatomy & histology; Aged; Air; Algorithms; Aniline Compounds; Bone and Bones - anatomy & histology; Bone and Bones - diagnostic imaging; Brain - anatomy & histology; Brain - diagnostic imaging; Ethylene Glycols; Female; Humans; Image Processing, Computer-Assisted; Magnetic Resonance Imaging - methods; Magnetic Resonance Imaging - statistics & numerical data; Male; Methods; Middle Aged; Neuroimaging - statistics & numerical data; Positron-Emission Tomography - methods; Positron-Emission Tomography - statistics & numerical data; Radiopharmaceuticals; Skull - anatomy & histology; Skull - diagnostic imaging; Studies; Ute
• ISSN: 1053-8119; 1095-9572; 1095-9572
• DOI: 10.1016/j.neuroimage.2015.03.009

MR-based correction for photon attenuation in PET/MRI remains challenging, particularly for neurological applications requiring quantitation of data. Existing methods are either not sufficiently accurate or are limited by the computation time required. The goal of this study was to develop an MR-based attenuation correction method that accurately separates bone tissue from air and provides continuous-valued attenuation coefficients for bone. PET/MRI and CT datasets were obtained from 98 subjects (mean age [±SD]: 66yrs [±9.8], 57 females) using an IRB-approved protocol and with informed consent. Subjects were injected with 352±29MBq of 18F-Florbetapir tracer, and PET acquisitions were begun either immediately or 50min after injection. CT images of the head were acquired separately using a PET/CT system. Dual echo ultrashort echo-time (UTE) images and two-point Dixon images were acquired. Regions of air were segmented via a threshold of the voxel-wise multiplicative inverse of the UTE echo 1 image. Regions of bone were segmented via a threshold of the R2* image computed from the UTE echo 1 and UTE echo 2 images. Regions of fat and soft tissue were segmented using fat and water images decomposed from the Dixon images. Air, fat, and soft tissue were assigned linear attenuation coefficients (LACs) of 0, 0.092, and 0.1cm−1, respectively. LACs for bone were derived from a regression analysis between corresponding R2* and CT values. PET images were reconstructed using the gold standard CT method and the proposed CAR-RiDR method. The RiDR segmentation method produces mean Dice coefficient±SD across subjects of 0.75±0.05 for bone and 0.60±0.08 for air. The CAR model for bone LACs greatly improves accuracy in estimating CT values (28.2%±3.0 mean error) compared to the use of a constant CT value (46.9%±5.8, p<10−6). Finally, the CAR-RiDR method provides a low whole-brain mean absolute percent-error (MAPE±SD) in PET reconstructions across subjects of 2.55%±0.86. Regional PET errors were also low and ranged from 0.88% to 3.79% in 24 brain ROIs. We propose an MR-based attenuation correction method (CAR-RiDR) for quantitative PET neurological imaging. The proposed method employs UTE and Dixon images and consists of two novel components: 1) accurate segmentation of air and bone using the inverse of the UTE1 image and the R2* image, respectively and 2) estimation of continuous LAC values for bone using a regression between R2* and CT-Hounsfield units. From our analysis, we conclude that the proposed method closely approaches (<3% error) the gold standard CT-scaled method in PET reconstruction accuracy. •Good bone/air segmentation can be achieved using UTE and Dixon images.•A sigmoid model is used to estimate continuous CT-HU values using MR R2* values.•PET reconstructions with a mean error of 2.6% in whole-brain are produced.•Attenuation map computation time is less than 15s per map.