Abstract
Background/Aim: To evaluate the utility of high spatial resolution digital positron emission tomography images with the source-to-background ratio (SBR) algorithm for gross tumour volume (GTV) delineation. Materials and Methods: The bowl and spheres (10-37 mm) were filled with fluoro-2-deoxy-D-glucose to achieve 4-16 times background radioactivity. The images were reconstructed using three isotropic voxel sizes. The SBR and percentage threshold (TH) to SUVmax were calculated. The plots between SBR and TH were fitted using a regression equation. The contoured volumes (CVs) of the spheres were calculated by applying TH. Results: TH was 38.6+75.0/SBR for 4 mm voxel size; 39.6+37.0/SBR for 2 mm; and 38.8+35.2/SBR for 1 mm. The mean relative errors between CV and true volume for 4, 2, and 1 mm voxel sizes were 15%, 7%, and 7%, respectively. Conclusion: The present technique is useful for GTV delineation with reduced contouring error.
- Positron emission tomography
- radiotherapy planning
- gross tumour volume
- auto-delineation
- digital positron emission tomography
Functional imaging using [18F]-fluoro-2-deoxy-D-glucose (FDG)-positron emission tomography (PET), which provides information pertaining to glucose metabolism, has been widely used for radiotherapy (RT) planning (1-3). Auto-delineation of gross tumour volume (GTV) using threshold based on the standardised uptake value (SUV) is a useful method for RT planning (4-6). The optimal percentage threshold (TH) to the maximum SUV (SUVmax) of the tumour reportedly ranges from 40% to 50% (4, 6, 7). However, some studies have shown that tumour size and background affect the TH; i.e., smaller lesions and a higher background activity necessitate a higher TH (8, 9). The source-to-background ratio (SBR) algorithm, which is based on the inverse function (y=a+b×1/x) between the SBR and TH, is widely accepted as a sophisticated and practical approach to address the effect of tumour size and background on TH (10, 11). Previous studies have demonstrated the clinical value of the SBR algorithm for RT planning (10, 11). Van Baardwijk et al. have compared the SBR-based auto-delineation with pathological findings in 23 patients with non–small-cell lung cancer (11). They found a strong correlation between the maximal tumour diameter of the SBR-based auto-delineation and the macroscopic diameter of primary tumours (Pearson correlation coefficient: 0.90). They have also evaluated the interobserver variability by calculating the concordance index (CI) for tumour delineation among the five observers. They found that the integration of SBR-based auto-delineation into manual delineation based on computed tomography (CT) images significantly increased the CI from 0.70 to 0.74 for primary tumours (p<0.001). Daisne et al. have estimated the volume of SBR-based auto-delineation and manual delineation based on the CT and magnetic resonance (MR) images and compared the volumes with surgical specimen obtained from nine patients with laryngeal cancer (10). The average mismatch of the volume between the imaging modalities and the surgical specimen was 46% for the SBR-based auto-delineation, 81% for manual delineation based on the CT images, and 107% for manual delineation based on the MR images. They concluded that SBR-based auto-delineation is the most accurate method for delineation of GTV among these three methods.
Positron emission tomography images with three radioactivity levels and three voxel sizes.
Conventional PET images have a limited spatial resolution (typical voxel size: 4-5 mm), and are influenced by lesion size and background on the TH (8). State-of-the-art digital PET provides higher spatial resolution images with a minimal voxel size of 1 mm (12-14). To the best of our knowledge, the introduction of the high spatial resolution digital PET images in radiotherapy planning has not been discussed. We conducted a phantom study to evaluate the utility of high spatial resolution digital PET images with SBR algorithm for auto-delineation of GTV in radiotherapy planning.
Materials and Methods
Phantoms and FDG-PET/CT imaging. We used the National Electrical Manufacturers Association International Electrotechnical Commission body phantom set (Data Spectrum Corporation, Hillsborough, NC, USA), which consists of an elliptical bowl (volume: 10.3 l) and six spheres (diameter range=10-37 mm). The bowl was filled with 2.65 kBq/ml FDG (Nihon Medi-Physics, Tokyo, Japan), and the spheres were filled with FDG to achieve 4, 8, and 16 times the background radioactivity to obtain the various SBRs (15). PET imaging was conducted using a whole-body 3D digital PET/CT scanner (Vereos; Philips Medical Systems, Cleveland, OH, USA). We acquired CT images [414 mA, 120 kV, 512×512 matrix, 600 mm field of view (FOV), 0.8 mm slice thickness] and then performed a PET scan using one bed position (FOV 576 mm; scan duration: 5 min). Attenuation correction was conducted with CT transmission data, and emission images were reconstructed using list-mode ordered subset expectation maximisation time-of-flight PET imaging, with three isotropic voxel sizes of 4, 2, and 1 mm and without point spread function. For each voxel size, we used a fixed number of iterations and subsets based on the previous report and our preliminary experiment (three iterations and 13 subsets for 4 and 2mm; three iterations and six subsets for 1 mm) (13). Post-reconstruction smoothing filter was not applied.
Image analysis. PET image files were transferred to the Velocity oncology imaging informatics system version 3.0.2 (Varian Medical Systems, Inc, Palo Alto, CA, USA) (Figure 1). We recorded the SUVmax of the spheres and background activity, which is calculated from the average activity in volume of interest of approximately 100 ml delineated in the phantom outside the sphere (16). The SUVmax, SBR, and TH were estimated for each sphere in the different reconstructed images with the three voxel sizes. The plots between the SBR and TH were fitted using a regression equation (y=a+b×1/x). To verify the validity of the relationships between SBR and TH established above, the contoured volume (CV) was calculated by applying the TH determined using regression equations. To compare CV with the true volume of sphere (TV), we estimated the mean relative error for each voxel size. Relative error was calculated using the following equation:
Statistical analysis. All statistical analyses were performed using SPSS version 25.0 (IBM, Armonk, NY, USA). Friedman test was performed to compare the SUVmax for six sphere sizes and reconstructed images with the three voxel sizes. Dunn-Bonferroni test was performed to compare between images with different voxel sizes as a post hoc analysis if the Friedman test result was significant (8). Regression analyses were performed to evaluate the relationship between the SBR and TH. p-Values <0.05 were considered significant.
Results
Eighteen groups of the SUVmax, SBR and TH were acquired from the images of each voxel size. Differences in the SUVmax were significant between six sphere sizes in the images with 4 mm voxel size (p=0.010), 2 mm (p=0.018), and 1 mm (p=0.012) (Table I, Figure 2). Differences in the SUVmax were significant between the three voxel sizes (p<0.001). The SUVmax of the images with 4 mm voxel size were significantly lower than the SUVmax of the images with 2 mm (p=0.003) and 1 mm (p<0.001). There was no significant difference in the SUVmax between the images with 2 mm and 1 mm voxel size (p=0.29). The relationship between the SBR and TH was described as follows: TH=38.6+75.0/SBR for 4 mm voxel size (R2=0.69, p<0.001); TH=39.6+37.0/SBR for 2 mm voxel size (R2=0.58, p<0.001); and TH=38.8+35.2/SBR for 1 mm voxel size (R2=0.60, p<0.001) (Figure 3). The mean relative error for 4, 2, and 1 mm voxel size was 15%, 7%, and 7%, respectively (Figure 4). Figure 5 shows the waterfall plots of relative error of CVs based on the above formulae for three voxel sizes.
Discussion
One of the limitations of conventional PET, which has a photomultiplier tube, is the relatively low spatial resolution; this results in a partial volume effect that affects the images both visually and quantitatively and encompasses a low signal-to-noise ratio (8, 9). Digital PET, which has a solid-state digital photon counting PET detector, helps overcome this limitation (12-14, 17, 18). However, our results suggested that the partial volume effect was unavoidable even in the high spatial resolution images with 2 mm and 1 mm voxel sizes. Our results also suggested that the degree of the partial volume effect depended on the spatial resolution of the images. Particularly, the partial volume effect was significantly greater with 4 mm voxel size than with 2 mm or smaller.
The SUVmax of six spheres for three source-to-background ratios and three voxel sizes.
Our results demonstrated that the SBR algorithm is useful not only for conventional resolution images with voxel size of 4 mm but also for high spatial resolution images with voxel size of 2 and 1 mm. However, the SBR algorithm is strongly affected by the spatial resolution of PET images: the relationship between the SBR and TH for 4 mm voxel size differed greatly from that for 2 and 1 mm. This difference in the SBR algorithm may be due to the difference in the partial volume effect between the images of the voxel sizes. Thus, the use of SBR algorithm based on the conventional PET images for auto-delineation of GTV based on high spatial resolution images may result in inappropriate GTV delineation, and the SBR algorithm for auto-delineation of GTV should be optimised on the basis of the spatial resolution of the PET images. Our results also suggested that the use of high spatial resolution images with isotropic voxel sizes of 2 or 1 mm reduces the contouring error compared with conventional resolution images with voxel size of 4 mm. Introduction of the combination of high spatial resolution PET images and dedicated SBR algorithm for auto-delineation facilitate precise definition of GTV.
There are some limitations of our study. First, this was an experimental phantom study and could not replicate the clinical setting with respect to complex tumour shape and heterogeneous background activity. Second, we did not evaluate the influence of reconstruction algorithm on the SBR algorithm. This factor may slightly influence SUVmax and the relationship between the SBR and TH (16). Further investigations are underway to address these issues.
The SUVmax of the six spheres and three voxel sizes for radioactivities of 4 (A), 8 (B), and 16 (C) times the background activity.
In conclusion, high spatial resolution digital PET images with SBR algorithm are useful for auto-delineation of GTV and for reducing contouring error compared with using conventional resolution PET images. A dedicated SBR algorithm for the high-resolution digital PET images is required for precise tumour delineation.
Relationship between the source-to-background ratio and the threshold with voxel sizes of 4 mm (A), 2 mm (B), and 1 mm (C).
Relationship between the true volume and the contoured volume of the spheres based on the source-to-background algorithm with voxel sizes of 4 mm (A), 2 mm (B), and 1 mm (C).
Relative error of contoured volume of the spheres based on the source-to-background algorithm with voxel sizes of 4 mm (A), 2 mm (B), and 1 mm (C).
Acknowledgements
This work was partially supported by the JSPS KAKENHI [grant number 18K07758].
Footnotes
Authors' Contributions
TW developed the study design; collected, analyzed, and interpreted the data; performed statistical analysis; and wrote the manuscript. RT developed the study design, analyzed and interpreted the data, and wrote the manuscript. TS, TM, and SS developed the study design and analyzed and interpreted the data. YK and YS collected the data. NO developed the study design and interpreted the data. All Authors have read and approved the final manuscript.
This article is freely accessible online.
Presentation at a conference: The 61st Annual Meeting of the American Society for Radiation Oncology, September 15–18, 2019, Chicago, IL, USA.
Conflicts of Interest
The Authors declare that they have no conflicts of interest regarding this study.
- Received March 11, 2020.
- Revision received March 28, 2020.
- Accepted March 29, 2020.
- Copyright© 2020, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved