Optimisation of the air fraction correction for lung PET/CT: addressing resolution mismatch

Abstract

BackgroundIncreased pulmonary 18\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{18}{}$$\\end{document}F-FDG metabolism in patients with idiopathic pulmonary fibrosis, and other forms of diffuse parenchymal lung disease, can predict measurements of health and lung physiology. To improve PET quantification, voxel-wise air fractions (AF) determined from CT can be used to correct for variable air content in lung PET/CT. However, resolution mismatches between PET and CT can cause artefacts in the AF-corrected image.MethodsThree methodologies for determining the optimal kernel to smooth the CT are compared with noiseless simulations and non-TOF MLEM reconstructions of a patient-realistic digital phantom: (i) the point source insertion-and-subtraction method, hpts\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$h_{pts}$$\\end{document}; (ii) AF-correcting with varyingly smoothed CT to achieve the lowest RMSE with respect to the ground truth (GT) AF-corrected volume of interest (VOI), hAFC\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$h_{AFC}$$\\end{document}; iii) smoothing the GT image to match the reconstruction within the VOI, hPVC\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$h_{PVC}$$\\end{document}. The methods were evaluated both using VOI-specific kernels, and a single global kernel optimised for the six VOIs combined. Furthermore, hPVC\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$h_{PVC}$$\\end{document} was implemented on thorax phantom data measured on two clinical PET/CT scanners with various reconstruction protocols.ResultsThe simulations demonstrated that at <200\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$<200$$\\end{document} iterations (200 i), the kernel width was dependent on iteration number and VOI position in the lung. The hpts\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$h_{pts}$$\\end{document} method estimated a lower, more uniform, kernel width in all parts of the lung investigated. However, all three methods resulted in approximately equivalent AF-corrected VOI RMSEs (<10%) at ≥\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\ge$$\\end{document}200i. The insensitivity of AF-corrected quantification to kernel width suggests that a single global kernel could be used. For all three methodologies, the computed global kernel resulted in an AF-corrected lung RMSE <10% at ≥\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\ge$$\\end{document}200i, while larger lung RMSEs were observed for the VOI–specific kernels. The global kernel approach was then employed with the hPVC\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$h_{PVC}$$\\end{document} method on measured data. The optimally smoothed GT emission matched the reconstructed image well, both within the VOI and the lung background. VOI RMSE was <10%, pre-AFC, for all reconstructions investigated.ConclusionsSimulations for non-TOF PET indicated that around 200i were needed to approach image resolution stability in the lung. In addition, at this iteration number, a single global kernel, determined from several VOIs, for AFC, performed well over the whole lung. The hPVC\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$h_{PVC}$$\\end{document} method has the potential to be used to determine the kernel for AFC from scans of phantoms on clinical scanners.