Variance prediction in SPECT reconstruction based on the Fisher information using a novel angular blurring algorithm for computation of the system matrix

In this paper we propose a Gaussian decomposition approach to compute the SPECT imaging system matrix. The flexibility of our method allows the implementation of a wide range of imaging systems. With the system matrix, we predict the variance in reconstructed images using the Fisher information matrix and local block circulant approximation. We present results of noise prediction for 3 multi-pinhole systems and 3 multi-slit slit-slat systems, each of them being designed to be inserted in an MRI system. Results show that for a particular phantom 2x2 multi-pinhole and 2 slits slit-slat systems achieve lowest variance. 2013 IEEE

Concurrent respiratory motion correction of abdominal PET and dynamic contrast-enhanced-MRI using a compressed sensing approach

We present an approach for concurrent reconstruction of respiratory motion-compensated abdominal dynamic contrast-enhanced (DCE)-MRI and PET data in an integrated PET/MR scanner. The MR and PET reconstructions share the same motion vector fields derived from radial MR data; the approach is robust to changes in respiratory pattern and does not increase the total acquisition time. Methods: PET and DCE-MRI data of 12 oncologic patients were simultaneously acquired for 6 min on an integrated PET/MR system after administration of 18F-FDG and gadoterate meglumine. Goldenangle radial MR data were continuously acquired simultaneously with PET data and sorted into multiple motion phases on the basis of a respiratory signal derived directly from the radial MR data. The resulting multidimensional dataset was reconstructed using a compressed sensing approach that exploits sparsity among respiratory phases. Motion vector fields obtained using the full 6-min (MC6-min) and only the last 1 min (MC1-min) of data were incorporated into the PET reconstruction to obtain motion-corrected PET images and in an MR iterative reconstruction algorithm to produce a series of motion-corrected DCE-MR images (moco-GRASP). The motioncorrection methods (MC6-min and MC1-min) were evaluated by qualitative analysis of the MR images and quantitative analysis of SUVmax and SUVmean, contrast, signal-to-noise ratio (SNR), and lesion volume in the PET images. Results: Motion-corrected MC6-min PET images demonstrated 30%, 23%, 34%, and 18% increases in average SUVmax, SUVmean, contrast, and SNR and an average 40% reduction in lesion volume with respect to the non-motion-corrected PET images. The changes in these figures of merit were smaller but still substantial for the MCMC1-min protocol: 19%, 10%, 15%, and 9% increases in average SUVmax, SUVmean, contrast, and SNR; and a 28% reduction in lesion volume. Moco-GRASP images were deemed of acceptable or better diagnostic image quality with respect to conventional breath-hold Cartesian volumetric interpolated breath-hold examination acquisitions. Conclusion: We presented a method that allows the simultaneous acquisition of respiratory motion-corrected diagnostic quality DCE-MRI and quantitatively accurate PET data in an integrated PET/MR scanner with negligible prolongation in acquisition time compared with routine PET/DCE-MRI protocols

PET/MRI in the presence of metal implants: Completion of the attenuation map from PET emission data

We present a novel technique for accurate whole-body attenuation correction in the presence of metallic endoprosthesis, on integrated non-time-of-flight (non-TOF) PET/MRI scanners. The proposed implant PET-based attenuationmap completion (IPAC) method performs a joint reconstruction of radioactivity and attenuation from the emission data to determine the position, shape, and linear attenuation coefficient (LAC) of metallic implants. Methods: The initial estimate of the attenuation map was obtained using the MR Dixon method currently available on the Siemens Biograph mMR scanner. The attenuation coefficients in the area of the MR image subjected to metal susceptibility artifacts are then reconstructed fromthe PET emission data using the IPAC algorithm. The method was tested on 11 subjects presenting 13 different metallic implants, who underwent CT and PET/MR scans. Relative mean LACs and Dice similarity coefficients were calculated to determine the accuracy of the reconstructed attenuation values and the shape of the metal implant, respectively. The reconstructed PET images were compared with those obtained using the reference CT-based approach and the Dixon-based method. Absolute relative change (aRC) images were generated in each case, and voxel-based analyses were performed. Results: The error in implant LAC estimation, using the proposed IPAC algorithm, was 15.7%± 7.8%, which was significantly smaller than the Dixon- (100%) and CT-(39%) derived values. A mean Dice similarity coefficient of 73% ± 9% was obtained when comparing the IPAC- with the CT-derived implant shape. The voxel-based analysis of the reconstructed PET images revealed quantification errors (aRC) of 13.2% ± 22.1% for the IPACwith respect to CT-corrected images. The Dixon-based method performed substantially worse, with a mean aRC of 23.1% ± 38.4%. Conclusion: We have presented a non-TOF emission-based approach for estimating the attenuation map in the presence of metallic implants, to be used for whole-body attenuation correction in integrated PET/MR scanners. The Graphics Processing Unit implementation of the algorithm will be included in the open-source reconstruction toolbox Occiput.io