Abstract
18F-fluoride PET/CT offers the opportunity for accurate skeletal metastasis staging, compared with conventional imaging methods. 18F-fluoride is a bone-specific tracer whose uptake depends on osteoblastic activity. Because of the resulting increase in bone mineralization and sclerosis, the osteoblastic process can also be detected morphologically in CT images. Although CT is characterized by high resolution, the potential of PET is limited by its lower spatial resolution and the resulting partial-volume effect. In this context, the synergy between PET and CT presents an opportunity to resolve this limitation using a novel multimodal approach called synergistic functional–structural resolution recovery (SFS-RR). Its performance is benchmarked against current resolution recovery technology using the point-spread function (PSF) of the scanner in the reconstruction procedure. Methods: The SFS-RR technique takes advantage of the multiresolution property of the wavelet transform applied to both functional and structural images to create a high-resolution PET image that exploits the structural information of CT. Although the method was originally conceived for PET/MR imaging of brain data, an ad hoc version for whole-body PET/CT is proposed here. Three phantom experiments and 2 datasets of metastatic bone 18F-fluoride PET/CT images from primary prostate and breast cancer were used to test the algorithm performances. The SFS-RR images were compared with the manufacturer’s PSF-based reconstruction using the standardized uptake value (SUV) and the metabolic volume as metrics for quantification. Results: When compared with standard PET images, the phantom experiments showed a bias reduction of 14% in activity and 1.3 cm3 in volume estimates for PSF images and up to 20% and 2.5 cm3 for the SFS-RR images. The SFS-RR images were characterized by a higher recovery coefficient (up to 60%) whereas noise levels remained comparable to those of standard PET. The clinical data showed an increase in the SUV estimates for SFS-RR images up to 34% for peak SUV and 50% for maximum SUV and mean SUV. Images were also characterized by sharper lesion contours and better lesion detectability. Conclusion: The proposed methodology generates PET images with improved quantitative and qualitative properties. Compared with standard methods, SFS-RR provides superior lesion segmentation and quantification, which may result in more accurate tumor characterization.
Up to 70% of patients with prostate and breast cancer will develop bone metastases (1,2). 18F-fluoride has shown efficacy in both diagnosis and treatment response assessment (3–6), and recent studies on skeletal metastases report improved diagnostic sensitivity and specificity when morphologic evaluation from CT scans are combined with functional evaluation of 18F-fluoride PET (7,8). This radiotracer accumulates at skeletal metastatic sites as a result of increased blood flow, osteoblastic activity, and bone mineralization (9–11). In prostate cancer, 18F-fluoride accumulation corresponds to sites of osteosclerosis and increased bone density that are usually visible on CT (12).
Evaluation criteria for tumor staging and response assessment include visual or quantitative evaluation of the extent, intensity, and changes in 18F-fluoride uptake in bone lesions (13). In this perspective, the influence of the partial-volume effect (PVE) is of importance (14) when comparing activity and morphologic changes before and after therapy (15,16) considering the poor image resolution and quantification bias resulting from activity spill-over. PVE in PET has been addressed with several image-based partial-volume correction (PVC) methods (17) that can be classified as either voxel-based methods, such as partition-based (18), multiresolution (19,20), or region-of-interest (ROI)–based techniques (21,22), which are limited by assumptions of radiotracer distribution homogeneity in the area of interest. A distinct approach consists of the incorporation of a model for the system point-spread function (PSF) within the image reconstruction algorithm to account for resolution degradation (23,24).
In this study, our aim was to correct for the PVE in whole-body 18F-fluoride PET/CT to allow a more robust lesion classification in terms of activity quantification and volume definition. The methodology was developed from previous work by Shidahara et al. (20) and exploits the local functional–structural relationship of PET/CT in a synergistic fashion for a realistic noise-controlled resolution recovery of PET images, hence the name synergistic functional–structural resolution recovery (SFS-RR). Here, the SFS-RR algorithm is optimized for 18F-fluoride PET given the correlation between functional (fluoride uptake) and morphologic (sclerosis) signals on PET and CT images. For benchmarking, the resulting images were compared with standard reconstructed PET images and images reconstructed with the inclusion of the PSF model (25).
MATERIALS AND METHODS
Image Resolution Recovery
The SFS-RR algorithm was first introduced by Shidahara et al. (20) for PVC of brain PET/MR data. The structural information was exploited by segmenting a T1-MR image through a probability atlas (26) defining 83 anatomic regions. Hence, the resolution recovery is ROI-based and relies on good coregistration between PET and MR images as well as between MR imaging and the probability atlas. The idea stems from previous concepts on wavelet-based resolution recovery (19) and denoising (27).
In this work, we developed SFS-RR further to fit a novel clinical requirement, specifically 18F-fluoride PET/CT for detecting and monitoring bone metastases. The choice of the application is not fortuitous; in the first instance, PET/CT images provide synergistic information (i.e., both modalities show high image intensity in correspondence of lesions) and, second, they do not require additional coregistration as for 2 separate PET and MR acquisitions. Furthermore, all the structural information of interest is contained in the CT and can be automatically segmented for each subject with no need for a universal atlas.
The algorithm decomposes both functional (PET) and anatomic (CT) images into several resolution elements by means of a wavelet transform. The high-resolution components of both modalities are then combined together via a statistical model with appropriate scaling, resolution correction, and weighting, to create a high-resolution PET image that exploits the structural information, when present, but preserves PET data when matching structural data are not present.
Anatomic Image Segmentation
In their original work, Shidahara et al. (20) proposed the use of an anatomic brain atlas to obtain suitable anatomic images. In brain studies, this is a reasonable procedure given the possibility of normalizing to a common space (e.g., probabilistic atlas).
In whole-body PET/CT, the atlas-based approach is not feasible. In 18F-fluoride PET/CT acquisitions, the good spatial correlation between morphologic and functional information is such that the CT images of each subject can be processed individually to highlight the structures of interest supplying the required structural base.
The first step consists of an initial coarse segmentation of the CT images based on thresholding the Hounsfield unit (HU) values as follows: bone (100 ≤ HU < 1,400), soft tissue (0 ≤ HU < 100), and fat (−150 ≤ HU < 0) (28).
Bone is further segmented into 100 bins after image histogram equalization. New intensity values are assigned to all the segmented regions. These intensity values are calculated from the average of each corresponding region in the original PET image to obtain the subject-specific structural reference image used as anatomic information for the SFS-RR algorithm (Supplemental Fig. 1; supplemental materials are available at http://jnm.snmjournals.org) (20,29).
Image segmentation and the SFS-RR algorithm implementation were both performed in Matlab R2011b (The MathWorks Inc.). The whole procedure is schematically described in Figure 1; for a more detailed mathematic exposition the reader is referred to the supplemental materials (Synergistic-Functional-Structural Resolution Recovery Algorithm section).
Phantom Data
For the evaluation of the SFS-RR method, we used the National Electrical Manufacturers Association International Electrotechnical Commission body phantom and an insert with 6 spheres of different volumes: 26.52 cm3 (S1), 11.49 cm3 (S2), 5.57 cm3 (S3), 2.57 cm3 (S4), 1.15 cm3 (S5), and 0.52 cm3 (S6). Compartments were filled with both iodinated contrast medium (CM) Omnipaque300 (GE Healthcare) (organic iodine [300 mg/mL]) and radioactive tracer 18F-fluoride. We aimed to reproduce contrast levels between different structures in both the PET and the CT images as observed in clinical 18F-fluoride bone scans. Specifically, we reproduced PET and CT contrasts as observed in normal soft tissue, normal bone, and metastatic bone. We performed 3 different experiments, changing the layout of CT and PET contrasts. This procedure aimed to account for possible mismatches between functional and anatomic images (i.e., whereas a lesion would be detectable in only 1 imaging modality) resulting in a more robust method validation. A summary of CM and radiotracer concentrations used in each experiment is reported in Table 1. For a detailed description of the experimental procedure, we refer to Grecchi et al. (30).
Images were acquired on a Discovery 710 PET/CT scanner (GE Healthcare). CT scans were obtained with a routine clinical protocol (115 mA, 140 kVp, and 0.5-s gantry rotation speed) followed by a fully 3-dimensional PET time-of-flight acquisition. PET data were reconstructed using our routine clinical protocol, a standard time-of-flight ordered-subsets expectation maximization algorithm (24 subsets, 2 iterations) (Q.Core VuePoint FX [GE Healthcare], henceforth called standard PET). The resulting images were then processed with the SFS-RR algorithm. In addition, we reconstructed the same data with the inclusion of a PSF model into the standard time-of-flight ordered-subsets expectation maximization algorithm (Q.Core VuePoint FX-S [GE Healthcare], henceforth called PET-PSF).
All the images were finally smoothed with a gaussian 3-dimensional filter (full width at half maximum, 6.4 mm); this is required by PSF-reconstructed images as well when used in clinical routine given that increased noise levels hamper visual assessment.
Clinical Dataset
The impact of the proposed resolution recovery technique was tested with 2 different sets of oncologic patient data, both characterized by the presence of bone metastases. The institutional review board approved this study, and all subjects signed a written informed consent form.
The first dataset was a prospective observational study of patients with bone-predominant metastatic prostate cancer, at first diagnosis or at progression of disease, who were embarking on docetaxel chemotherapy. The second dataset was a prospective observational study of patients with bone-predominant metastatic breast cancer, at first diagnosis or at progression of disease, who were embarking on a new line of endocrine treatment in combination with bone-targeted therapy. In total, 7 patients with active skeletal metastases were included in the analysis, 4 with prostate cancer and 3 with breast cancer.
In each patient, a whole-body 18F-fluoride PET/CT scan was acquired with a total of 8 bed positions, from the base of the skull to upper thighs, 60 min after injection of approximately 250 MBq.
The image reconstruction protocol for both datasets was the same as for the phantom experiment, with the exclusion of the PSF reconstruction, which was not performed as it was not included in the clinical protocol.
Data Analysis
In-house software was used to perform quantitative analysis on both phantom and patient data. An ROI was manually drawn on the outer border of each lesion to completely contain the whole lesion volume (or sphere for the phantom data) and then automatically segmented with a threshold of 40% of the maximum value of the ROI. The mean standardized uptake value (SUVmean), maximum SUV (SUVmax), peak SUV (SUVpeak), and lesion metabolic active tumor volume (MATV) were then computed for the automatically segmented ROI. SUVpeak is computed here as the SUVmean measured over a fixed small circular volume of about 1 cm3, in the hottest area of the tumor (more active region). It is considered more reproducible because it involves the mean value of a few voxels involving and surrounding the hottest tumor area.
Solely for the phantom experiments, knowing the ground truth, we used the root mean square error (RMSE) and contrast-to-noise ratio (CNR) as additional metrics for image-quality assessment. Specifically, we evaluated the quantification accuracy and the trade-off between contrast improvement and image noise. The 2 metrics are defined as follow represents the mean activity estimated inside a sphere and the corresponding ground truth whereas is the mean activity estimated in the phantom background and its SD.
RESULTS
Phantom Data
Figure 2 shows representative transaxial views of structural (CT) and functional (PET) images for the 3 phantom experiments, 1 for each line. PET images from PSF reconstruction (PET-PSF) and from the application of the resolution recovery algorithm (PET-SFS) are also displayed. Supplemental Figure 2 shows representative line profiles for the 3 experiments for spheres 4–5.
Improved qualitative resolution for the smaller structures (red marker in Fig. 2) is noted when the SFS-RR algorithm is applied. Even though larger spheres (green and blue markers in Fig. 2) are easily detectable in the images from all modalities, it is possible to appreciate a reduction in the blurring surrounding the structure when the resolution recovery algorithm is implemented.
Furthermore, it is possible to appreciate the robustness of the anatomy-based resolution recovery algorithm to unexpected mismatches between anatomy and functional acquisition by studying images from experiment 1. Even though spheres 4–6 cannot be detected on the CT image, they are not lost in the new functional image returned by the SFS-RR algorithm.
The quantitative evaluation of functional images obtained with the 3 different methods is reported in Table 2. The table reports for each phantom compartment in all the experiments the estimates of SUVmean, SUVmax, SUVpeak, and MATV together with the corresponding ground-truth values.
A summary of method performances is summarized in the supplemental materials (Supplemental Fig. 3) as the average among the 3 experiments. The general trend shows that the smaller the sphere, the bigger the bias in the activity estimation, regardless of the method used. However, with the SFS-RR application the bias decreases with an average range of 1%–5% in the PET-PSF images and 5%–19% in PET-SFS images. The same trend applies to lesion size estimation in which the bias decreases in a range of 0.46–0.95 cm for PET-PSF data and 0.56–1.09 cm for PET-SFS data.
The better performances of the SFS-RR algorithm are upheld by the RMSE and CNR comparison in Figure 3. Images resulting from the application of the SFS algorithm show lower RMSE on average (up to 15%, compared with standard PET for the smallest sphere) while being consistent with the trend of the RMSE increasing for smaller spheres. The improved image resolution does not come with reduced image quality because noise levels are contained with CNR either higher than or comparable with the standard PET image.
Clinical Dataset
The influence of using the SFS-RR algorithm on real patient data can be appreciated qualitatively in Figure 4 (and Supplemental Figs. 4–6) in which the maximum-intensity projection and 2 different transaxial views are displayed for a representative subject.
There is a clear increase in lesion sharpness after the application of a resolution recovery technique. The quantitative characterization of all lesions of this specific subject is reported in Supplemental Fig. 7 and Supplemental Table 1. The transaxial views of Figure 4 are a good example of the effect in lesion definition and characterization using the SFS-RR algorithm. Sharper contours and the activity recovery in the PET-SFS images for the rib lesion (Fig. 4, transaxial view, red marker) and also in the spine (Fig. 4, transaxial view, blue marker) are evident, and in the spine lesion it is easier to appreciate that the activity is in the periphery of the lesion in which there is greatest osteoblastic activity compared with the relatively photopenic center. In terms of quantitative characterization, there is an increase in the SUVmean estimates of 60% (rib lesion, SUV PET = 30.7 and SUV SFS = 49.1) and 43% (spine lesion, SUV PET = 23.2 and SUV SFS = 33.1) from the standard PET to the PET-SFS. In contrast, the automatic segmented MATV has a relative reduction of 25% (rib lesion, MATV PET = 2.7 cm3 and MATV SFS = 2.1 cm3) and 31% (spine lesion, MATV PET = 7.5 cm3 and MATV SFS = 5.2 cm3).
A comprehensive comparison of the quantitative differences due to application of the SFS-RR algorithm is reported in Figure 5. As previously performed for patient 01 (Supplemental Fig. 7; Supplemental Table 1), we segmented all lesions of the remaining patients and collected the corresponding values of SUVmean, SUVmax, SUVpeak, and MATV. Figure 5 reports the relative differences (SUVs and MATV) between lesions segmented in standard PET and PET-SFS images for the entire datasets.
There is a general increment in activity estimates for PET-SFS compared with standard PET. The average increments for different indices are as follow: ΔSUVmean = 49%, ΔSUVmax = 47%, and ΔSUVpeak = 34%. For low-activity lesions, the range of differences in the estimates is larger than for more active lesions, indicating that the algorithm efficiency is dependent on signal-to-noise ratio. For the MATV, there is an average reduction of 1.4 cm3 when segmentation is performed on PET-SFS images. When lesion size increased, the difference reached values of 4–5 cm3, which might be relevant if patient classification were based on characterization of the larger detectable lesions.
DISCUSSION
In this work, we evaluated the influence of a multimodal PVC technique on the quantification and assessment of metastatic bone lesions from primary prostate and breast cancers. This work targets specifically to 18F-fluoride PET/CT bone scanning because SFS-RR obviously produces better results when the correspondence between functional and structural signals is stronger. Our results showed an average 50% increase in SUVmax and SUVmean and a 30% increase in the SUVpeak for partial-volume-corrected images when compared with the standard PET, depending on lesion size (lesion volume range, 0.5–25 cm3). Our results are in agreement with findings in similar experimental settings from previous studies (31,32). Although SUVmax estimates depend on image noise, the 50% increase is not a consequence of noise bursts given the comparable RMSE between the SFS-RR–corrected images and standard PET images.
The higher activity recovery and the good noise control from the phantom analysis indicate a better image quality when the SFS-RR algorithm is applied. Indeed, patient images show lesions with sharper and better-defined contours, which result in improved lesion conspicuity and segmentation even for smaller volumes. The CNR depends on the absolute activity value, thus explaining why one set of spheres shows higher CNR than the other on Figure 3.
The results from the phantom experiments showed that the SFS-RR images outperform both standard PET and PSF images in terms of image quality and quantification accuracy. PSF-based image reconstruction is known to contribute to the appearance of artifacts (33) and is computationally cumbersome hence is not performed in routine clinical studies in our unit; for this reason, standard PET was used as a reference for SFS-RR images when it came to patient image analysis. In this regard, no artifacts have been generally observed in this and previous applications of SFS-RR, which is now a mature-enough technology worthy of further testing in the clinical setting.
It is worth highlighting the robustness of the methodology regarding possible mismatch between PET and CT images. Phantom acquisitions showed that even if some structures were visible only on the functional images, they were preserved after the application of the algorithm. This is of importance because lesions that might be lost in the CT segmentation, for example, due to their small size or the fact that the metastasis does not show sufficient bone mineralization to appear sclerotic, will still be visible in the final enhanced PET images.
Although the SFS-RR algorithm showed qualitatively and quantitatively better images than standard PET, further analyses are necessary to quantify the influence of the improved image quality on the assessment of patient skeletal staging and therapy response. Additional analyses may allow better definition and quantification of lesions after therapy or allow greater detectability and segmentation of metastatic spread at staging. Of additional interest would be the evaluation of whether lesion heterogeneity is affected by higher resolution and evaluation of the consequent impact on textural analysis, given the increasing oncologic applications of textural analysis (34,35).
CONCLUSION
We have proposed and tested on a set of phantom studies and demonstrated on clinical data a multimodal methodology for quantitative resolution recovery for whole-body PET/CT, here specifically designed for 18F-fluoride PET imaging of bone metastases. The technique allows rapid and straightforward application and produces images of significantly improved visual quality and quantitative properties.
DISCLOSURE
The costs of publication of this article were defrayed in part by the payment of page charges. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734. The project is supported by the EPSRC (scholarship EP/K502868/1 and grant WT 088641/Z/09/Z); the NIHR Biomedical Research Centre at Guy's and St. Thomas' NHS Foundation Trust and King’s College London; UCL Comprehensive Cancer Imaging Centre, funded by the CRUK and EPSRC in association with the MRC and DoH (England); Centre of Excellence in Medical Engineering; Wellcome Trust; and MRC PET Methodology Program grant G1100809/1. The views expressed are those of the authors and not necessarily those of the NHS, the NIHR, or the Department of Health. No other potential conflict of interest relevant to this article was reported.
Acknowledgments
We wish to acknowledge Dr. Benjamin Taylor for providing us anonymized patient data.
Footnotes
Published online Jul. 16, 2015.
- © 2015 by the Society of Nuclear Medicine and Molecular Imaging, Inc.
REFERENCES
- Received for publication May 8, 2015.
- Accepted for publication July 8, 2015.