Improved Regional Activity Quantitation in Nuclear Medicine Using a New Approach to Correct for Tissue Partial Volume and Spillover Effects

• Published in: IEEE transactions on medical imaging Vol. 31; no. 2; pp. 405 - 416
• Main Authors: Moore, S. C., Southekal, S., Mi-Ae Park, McQuaid, S. J., Kijewski, M. F., Muller, S. P.
• Format: Journal Article
• Language: English
• Published: United States IEEE 01.02.2012; The Institute of Electrical and Electronics Engineers, Inc. (IEEE)
• Subjects: Algorithms; Animals; Beads; Bias; Computer simulation; Estimates; Image Enhancement - methods; Image Interpretation, Computer-Assisted - methods; Image reconstruction; Image resolution; Image segmentation; Imaging; Liver; Lungs; Mathematical model; Mice; Neoplasms - diagnostic imaging; Ordered subsets expectation-maximization (OS-EM); partial-volume effect; Pattern Recognition, Automated - methods; Phantoms, Imaging; positron emission tomography (PET); Projection; Reproducibility of Results; Sensitivity and Specificity; single photon emission computed tomography (SPECT); standardized uptake value (SUV); Studies; Tomography, Emission-Computed, Single-Photon - instrumentation; Tomography, Emission-Computed, Single-Photon - methods; Tumors
• ISSN: 0278-0062; 1558-254X
• DOI: 10.1109/TMI.2011.2169981

We have developed a neσw method of compensating for effects of partial volume and spillover in dual-modality imaging. The approach requires segmentation of just a few tissue types within a small volume-of-interest (VOI) surrounding a lesion; the algorithm estimates simultaneously, from projection data, the activity concentration within each segmented tissue inside the VOI. Measured emission projections were fitted to the sum of resolution-blurred projections of each such tissue, scaled by its unknown activity concentration, plus a global background contribution obtained by reprojection through the reconstructed image volume outside the VOI. The method was evaluated using multiple-pinhole μSPECT data simulated for the MOBY mouse phantom containing two spherical lung tumors and one liver tumor, as well as using multiple-bead phantom data acquired on μSPECT and μCT scanners. Each VOI in the simulation study was 4.8 mm (12 voxels) cubed and, depending on location, contained up to four tissues (tumor, liver, heart, lung) with different values of relative 99m Tc concentration. All tumor activity estimates achieved <; 3% bias after ~ 15 ordered-subsets expectation maximization (OSEM) iterations (×10 subsets), with better than 8% precision (≤ 25% greater than the Cramer-Rao lowσer bound). The projection-based fitting approach also outperformed three standardized uptake value (SUV)-like metrics, one of which was corrected for count spillover. In the bead phantom experiment, the mean ± standard deviation of the bias of VOI estimates of bead concentration were 0.9±9.5%, comparable to those of a perturbation geometric transfer matrix (pGTM) approach (-5.4±8.6%); however, VOI estimates were more stable with increasing iteration number than pGTM estimates, even in the presence of substantial axial misalignment between μCT and μSPECT image volumes.