Abstract
Purpose:
Image segmentation can be time-consuming and lacks consistency between different oncologists, which is essential in conformal radiotherapy techniques. We aimed to evaluate automatic delineation results generated by convolutional neural networks (CNNs) from geometry and dosimetry perspectives and explore the reliability of these segmentation tools in rectal cancer.
Methods:
Forty-seven rectal cancer cases treated from February 2018 to April 2019 were randomly collected retrospectively in our cancer center. The oncologists delineated regions of interest (ROIs) on planning CT images as the ground truth, including clinical target volume (CTV), bladder, small intestine, and femoral heads. The corresponding automatic segmentation results were generated by DeepLabv3+ and ResUNet, and we also used Atlas-Based Autosegmentation (ABAS) software for comparison. The geometry evaluation was carried out using the volumetric Dice similarity coefficient (DSC) and surface DSC, and critical dose parameters were assessed based on replanning optimized by clinically approved or automatically generated CTVs and organs at risk (OARs), i.e., the Planref and Plantest. Pearson test was used to explore the correlation between geometric metrics and dose parameters.
Results:
In geometric evaluation, DeepLabv3+ performed better in DCS metrics for the CTV (volumetric DSC, mean = 0.96, P< 0.01; surface DSC, mean = 0.78, P< 0.01) and small intestine (volumetric DSC, mean = 0.91, P< 0.01; surface DSC, mean = 0.62, P< 0.01), ResUNet had advantages in volumetric DSC of the bladder (mean = 0.97, P< 0.05). For critical dose parameters analysis between Planref and Plantest, there was a significant difference for target volumes (P< 0.01), and no significant difference was found for the ResUNet-generated small intestine (P > 0.05). For the correlation test, a negative correlation was found between DSC metrics (volumetric, surface DSC) and dosimetric parameters (δD95, δD95, HI, CI) for target volumes (P< 0.05), and no significant correlation was found for most tests of OARs (P > 0.05).
Conclusions:
CNNs show remarkable repeatability and time-saving in automatic segmentation, and their accuracy also has a certain potential in clinical practice. Meanwhile, clinical aspects, such as dose distribution, may need to be considered when comparing the performance of auto-segmentation methods.
1 Introduction
Preoperative radiotherapy is currently considered to be the standard treatment for locally advanced rectal cancer and has been proven to reduce local recurrence (–). With the development of radiotherapy technology, such as intensity modulated radiotherapy (IMRT) and volumetric modulated arc therapy (VMAT), the target volume can receive a highly conformal dose distribution (). In addition, it has been proven that IMRT and VMAT are dosimetrically superior to other conformal techniques in protecting organs at risk (OARs) in rectal cancer (). Thus, the accurate delineation of the clinical target volume (CTV) and OARs is crucial for treatment planning in rectal cancer.
Interobserver differences occur during manual delineation, which depend on oncologists’ clinical experience, resulting in significant changes in dose distributions (). Multiple studies have applied deep learning methods to automatic segmentation to solve the problem of time consumption and the lack of consistency in manual contouring (–). Based on the planning computed tomography (pCT) images, the oncologists’ delineated regions of interest (ROIs) as a training set. These structures are imported into deep learning models with CT images, and their corresponding features are extracted to train models according to the framework characteristics of different models.
The accuracy of automatic segmentation requires clinical evaluation. Objective evaluation metrics such as the volumetric Dice similarity coefficient (volumetric DSC) and Hausdorff distance (HD) are widely used, and some studies have carried out dosimetry assessments (–). However, clinical evaluation of the quality of deep learning delineation has limitations (). Considering the different accuracy requirements of CTVs and OARs in clinical practice, it is necessary to combine their clinical importance and tolerant errors and carry out a comprehensive evaluation from the perspectives of geometry and dosimetry.
We carried out a retrospective study of radiotherapy patients with rectal cancer. CTV and OARs were segmented manually as the ground truth (GT), two convolutional neural networks we have trained—DeepLabv3+ and ResUNet—were used for automatic delineation (), and a common method Atlas-based Auto segmentation (ABAS) was used as a comparison. Our research aimed to explore the clinical impact of auto-segmentation results from a dosimetric perspective.
2 Materials and methods
2.1 Patient data
The retrospective study was approved by the ethics committee of West China Hospital in 2019, with no extra health risks and no need for patient consent. Rectal cancer patients from February 2018 to April 2019 at West China Hospital were chosen randomly and metastases in advanced patients were ignored. Each patient was immobilized in a supine position with arms over the head using a radiotherapy thermoplastic mold, and this position was applied during simulation and treatment. The contrast-enhanced CT images (tube voltage, 120 kVp; matrix size, 512 × 512; voxel resolution, 0.9 × 0.9 × 3.0 mm in left-right, antero-posterior and cranio-caudal directions) were acquired as patient pCT on the same CT scanner (SOMATOM Definition AS, Siemens Healthcare).
Based on the pCT, a radiation oncologist manually segmented the CTV and OARs of rectal cancer by referring to Radiation Therapy Oncology Group consensus guidelines. Then, the structures were modified and approved by a senior expert physician and labelled as ground truth, including CTVGT, bladderGT, small intestineGT, left femoral headGT and right femoral headGT.
2.2 Automatic segmentation
DeepLabv3+ and ResUNet, two typical CNNs, were used for automating delineation. DeepLabv3+ employs an atrous spatial pyramid pooling module and concat aggregation for the extraction and integration of high-level features, and ResUNet has shortcut connections for each level of features (). The pCT images of all patients were imported into the models to obtain the mask of each structure on every CT slice. The two-dimensional masks were then converted to a three-dimensional structure in DICOM format, imported to the treatment planning system (TPS), and labelled as ROIDeepLabv3+ and ROIResUNet, respectively. The network models were uploaded onto github (https://github.com/hujunjiescu/DeepRadiology_rectum), and the architecture diagram was shown in Figure 1.
CNN models had the same training settings. The contouring tasks were worked out based on the Pytorch deep learning framework using Python. The experiments were performed on a Linux operating system workstation with the CPU Intel Xeon E5-2620 v3@ 2.4GHz, GPU NVIDIA Tesla K40 Xp, and 64 GB RAM. The loss function for the optimization was the weighted cross-entropy, which was defined as:
where N, C, w, y, and a denoted the batch size, number of classes, weight factor, ground truth sets, and prediction sets, respectively. The batch size N was set at 10, the weight factor w at 2, and the total training epoch T at 100. The stochastic gradient descent method was used to optimize the network with the initial learning rate set as 0.01, which was multiplied by for the epoch t. The segmentation results were rewritten into DICOM RT structure (RS) files based on their original spatial resolutions.
The ABAS worked on CT datasets using a multi-patient atlas. We randomly selected 5 atlas patients from the CNN training set, then their pCT images and manual contoured structures were imported to ABAS software (Version 2.01.00, Elekta CMS, Inc.). The Simultaneous Truth And Performance Level Evaluation (STAPLE) algorithm was used to fuse the multiple single-subject atlas auto-segmentation sets into one multi-subject auto-segmentation set ().
2.3 Treatment plans
To evaluate the clinical dosimetry value of the two automatic delineation methods, a two-round optimization protocol was performed using TPS (Raystation, version 4.7.5, Raystation Laboratories, Stockholm, Sweden): 1) The corresponding PTV was obtained based on CTV expanded with a three-dimensional margin of 5 mm; 2) The dose prescription was set to 50.4Gy/28 fraction to the PTV; 3) Two full arcs, one from 181 to 180° clockwise and the other from 180 to 181° counterclockwise, were designed using the VMAT technique and 6 MV photons; 4) The initial optimization parameters applied to the first round VMAT planning were shown in Table 1, and in the second round, the weight of Parameter4 was set to 100, and the weight of Max EUD objectives was set to 0.01.
Table 1
| Parameter (Pi) | ROI | Description | Weight |
|---|---|---|---|
| P1 | PTV | Min Dose 50.4Gy | 60 |
| P2 | PTV | Max Dose 52.42Gy | 90 |
| P3 | PTV | Min DVH 50.9Gy to 95% volume | 100 |
| P4 | PTV | Uniform Dose 51.21Gy | 0.2 |
| P5 | Bladder Avoid | Max DVH 40Gy to 52% volume | 20 |
| P6 | Bladder Avoid | Max EUD 28Gy, Parameter α 1 | 1 |
| P7 | Small intestine Avoid | Max DVH 30Gy to 30% volume | 20 |
| P8 | Small intestine Avoid | Max EUD 20Gy, Parameter α 1 | 1 |
| P9 | Femoral Head Right | Max DVH 40Gy to 5% volume | 20 |
| P10 | Femoral Head Right | Max EUD 1500, Parameter α 1 | 1 |
| P11 | Femoral Head Left | Max DVH 40Gy to 5% volume | 20 |
| P12 | Femoral Head Left | Max EUD, Parameter α 1 | 1 |
| P13 | External | Dose Fall-Off [H]50.4Gy [L]25.2Gy, [D]2.8 cm | 15 |
Initial planning objective set for plan optimization.
Taking the difficulty of CTV delineation into consideration, we divided the results of the three auto-segmentation methods into CTV and OAR groups and introduced them as optimization objectives separately to obtain Plantest: 1) the plan optimized using CTVDeepLabv3+ and OARGT labelled Plan1; 2) the plan optimized using CTVResUNet and OARGT labelled Plan2; 3) the plan optimized using CTVABAS and OARGT labelled Plan3; 4) the plan optimized using CTVGT and OARDeepLabv3+ labelled Plan4; 5) the plan optimized using CTVGT and OARResUNet labelled Plan5; and 6) the plan optimized using CTVGT and OARABAS labelled Plan6. In addition, we obtained a Planref optimized using ROIGT, calculated characteristic dose parameters of ROIGT in all plans, and compared parameters extracted from Plantest with those from the Planref respectively.
Figure 1
2.4 Evaluation metrics and statistical analysis
In terms of geometry, two DSCs were used to evaluate quantitatively, which were calculated on the overlap of the ROI structures. The ROIs were converted from RS files to thresholding masks, and the masks were divided into slices corresponding to the CT images. The mathematical operations were carried out based on all mask slices of a certain structure and the average was obtained as DSC value.
The volume similarity was usually evaluated by volumetric DSC, calculated using:
where V1 were the ROIs of ground truth set and V2 were the corresponding auto-segmentation structures. Volumetric DSC varies between 0 (no overlap) and 1 (complete overlap), which indicates the degree of overlap between ROIGT and auto-segmentation results.
To characterize the proportion of the contour edges that need to be redrawn, the surface DSC was applied to assess the agreement between just the surface of two contours (
Figure 2

Calculation method of surface DSC. (A) acceptable tolerance τ value, (B) surface DSC formula, (C) the calculation process taken CTV as an example, in which the red lines in and were the part that exceeds the tolerance.
In addition, the clinical practicability of contours delineated automatically is evaluated by the accuracy of the dose distribution in plan design. The characteristic dosimetry parameters of ROIGT in every plan were extracted for comparison. D2 (Dn representing the absolute dose of n% volume) and D98 were extracted to signify hot spots and cold spots for all structures, respectively (
The collected data were analyzed using SPSS Statistics software (version 26.0, SPSS Inc., Chicago, IL, United States). Normality tests were performed on all datasets of geometric and dosimetric parameters. Paired samples t tests or Wilcoxon signed rank tests were chosen for group comparison with statistical significance set at P value< 0.05 (2-tailed). To make a more intuitive comparison, we calculated the absolute difference between the dose parameters extracted from the Planref and the Plantest, denoted as DAbs, and carried out a statistical description. In particular, Vn, HI, and CI were relative values and directly subtracted, while Dn were absolute values and converted to normalized dose difference () (
3 Results
Forty-seven rectal cancer patients were included in the study. The median age was 54 years, with a interquartile range (IQR) of 13.97, and other characteristics are shown in Table 2. For patients diagnosed with stage IV disease, the study ignored metastases in the training and evaluation. In 5 cases, the structures were not successfully generated from the ABAS software.
Table 2
| Characteristic | Value |
|---|---|
| Sex | |
| Male | 23 (49%) |
| Female | 24 (51%) |
| Age | |
| Median (range) | 54 (28-83) |
| Cancer classification | |
| I | 4 (9%) |
| II | 5 (11%) |
| III | 32 (68%) |
| IV | 6 (13%) |
Characteristics of 47 patients.
The statistical analysis results of volumetric and surface DSC are shown in Figure 3. In general, the volumetric and surface DSC of the three automatic segmentation results were significantly different, except for the surface DSC of the bladder structure delineated by DeepLabv3+ and ResUNet (P = 0.78). For CTV, DeepLabv3+ showed the best performance on volumetric DSC (mean = 0.96, P< 0.01) and surface DSC (mean = 0.96, P< 0.01). The DSCs of CNNs automatic contouring bladder were significantly higher than those of ABAS (P< 0.01), although there was no significant difference in surface DSC between the two CNNs, the mean value of ResUNet was slightly higher than that of DeepLabv3+ (BladderDeepLabv3+vs. BladderResUNet, 0.82 vs. 0.85, P = 0.78). In the delineation of the small intestine, DeepLabv3+ showed significant advantages, whose mean DSCs (volumetric DSC, 0.91; surface DSC, 0.62) were greater than those of the other two groups (P< 0.01), as well as a lower standard deviation (volumetric DSC, 0.05; surface DSC, 0.10). For the segmentation of the right and left femoral head, ABAS achieved the best performance, then the ResUNet, and DeepLabv3+ ranked the last based on the volumetric DSC (mean for the right femoral head, ABAS vs. ResUNet vs. DeepLabv3+, 0.94 vs. 0.85 vs. 0.84, P< 0.01) and surface DSC (mean for the right femoral head, ABAS vs. ResUNet vs. DeepLabv3+, 0.84 vs. 0.70 vs. 0.67, P< 0.01), despite more outliers in volumetric DSC of ABAS. The ground truth and the automatic delineation of a random case were shown in Figure 4.
Figure 3

Statistical analysis results of volumetric and surface DSC. The paired- sample tests were performed between the three auto-segmentation results at a significance level of 0.05(2-tailed), and the missing data in the ABAS dataset (n=42) were replaced with the mean.
Figure 4

A case of structures comparison. The red line represents ground truth, the yellow line represents automatic segmentation of DeepLabv3+, the blue line represents automatic segmentation of ResUNet.
The DAbs of the dose parameters between Plantest and Planref were shown in Table 3, and it contained descriptive statistics and results of statistical tests. We observed a statistical difference in dose distribution of real CTV between the Planref and Plan 1-3 (P< 0.01), which used automatically delineated CTV as an inverse optimization parameter. The difference, however, was that some dose parameters of real OAR in Plan1 were not significantly different from the Planref (P > 0.05), which showed a similar trend to the outperformance of DeepLabv3+ in geometric evaluation. When we introduced the auto-segmentation OAR groups into the inverse plan, we found that only critical dose parameters of the small intestine between Plan5 and the Planref had no statistical difference (P > 0.05), but the small intestine generated by ResUNet was not optimal in the geometric assessment (volumetric DSC, mean = 0.84; surface DSC, mean = 0.52).
Table 3
| ROI | Dose Parameter | Plan1 (CTVDeepLabv3+ and OARGT) | Plan2 (CTVResUNet and OARGT) | Plan3 (CTVABASt and OARGT) | Plan4 (OARDeepLabv3+ and CTVGT) | Plan5 (OARResUNet and CTVGT) | Plan6 (OARABAS and CTVGT) | ||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Pa | DAbsb | P | DAbs | P | DAbs | P | DAbs | P | DAbs | P | DAbs | ||
| CTVGT | V50.4 | <0.01 | 3.74(4.34) | <0.01* | 4.14(5.27) | <0.01 | 6.39(9.31) | 0.04 | 1.87(2.80) | 0.24 | 1.61(2.23) | <0.01 | 2.17(3.95) |
| PTVGT | D2 | <0.01* | 0.40(0.30) | <0.01 | 0.34(0.35) | <0.01 | 0.61(0.42) | <0.01 | 0.23(0.19) | 0.02 | 0.19(0.35) | 0.03 | 0.32(0.42) |
| D95 | <0.01* | 4.15(4.20) | <0.01* | 6.82(4.66) | <0.01 | 20.36(80.90) | <0.01 | 0.15(0.19) | 0.12 | 0.13(0.17) | 0.04 | 0.32(0.46) | |
| D98 | <0.01* | 10.43(7.60) | <0.01 | 12.71(7.93) | <0.01 | 41.28(77.17) | <0.01 | 0.19(0.25) | 0.02 | 0.18(0.22) | 0.16 | 0.44(0.68) | |
| CI | <0.01 | 0.11(0.04) | <0.01 | 0.13(0.05) | <0.01 | 0.20(0.08) | 0.03 | 0.01(0.02) | 0.23 | 0.01(0.02) | 0.08 | 0.03(0.05) | |
| HI | <0.01* | 0.10(0.08) | <0.01* | 0.12(0.08) | <0.01 | 0.40(0.76) | <0.01 | 0.003(0.004) | <0.01 | 0.003(0.003) | <0.01 | 0.01(0.01) | |
| BladderGT | D2 | <0.01 | 0.40(0.61) | 0.01 | 0.40(0.53) | <0.01 | 0.44(0.61) | 0.11* | 0.19(0.39) | 0.14 | 0.29(0.40) | 0.83 | 0.36(0.48) |
| V30 | 0.72 | 5.34(5.29) | 0.11 | 5.78(6.55) | <0.01 | 18.49(10.37) | 0.01* | 2.21(3.07) | 0.01* | 2.56(4.14) | <0.01 | 23.89(9.01) | |
| V40 | 0.28 | 4.58(7.34) | 0.01 | 6.01(8.22) | <0.01* | 16.71(13.28) | 0.36* | 1.59(2.07) | 0.35* | 1.53(1.91) | <0.01 | 22.87(12.91) | |
| V50 | 0.12 | 5.67(6.65) | 0.01 | 5.84(5.65) | 0.09* | 8.04(7.19) | <0.01* | 0.90(1.68) | <0.01* | 1.11(1.31) | 0.06 | 1.11(2.39) | |
| Small intestineGT | D2 | 0.16* | 0.32(0.37) | 0.01* | 0.29(0.42) | <0.01 | 0.47(0.76) | 0.82 | 0.11(0.19) | 0.49 | 0.18(0.17) | 0.14 | 0.26(0.49) |
| V15 | <0.01* | 3.10(3.04) | <0.01* | 3.44(3.60) | 0.04 | 8.63(17.07) | 0.20* | 1.04(1.84) | 0.56* | 1.08(1.82) | <0.01 | 3.83(7.82) | |
| V45 | 0.02 | 2.06(3.42) | <0.01 | 2.64(3.68) | 0.97* | 3.30(5.35) | 0.02* | 0.26(0.66) | 0.13* | 0.24(0.60) | <0.01* | 4.09(4.33) | |
| V50 | 0.04 | 1.92(2.57) | <0.01 | 1.92(3.44) | <0.01* | 2.05(3.84) | 0.12* | 0.24(0.30) | 0.54* | 0.18(0.37) | 0.95 | 0.26(0.67) | |
| Right Femoral HeadGT | D2 | 0.17 | 4.75(7.34) | 0.06 | 4.14(5.28) | <0.01* | 5.29(4.68) | <0.01 | 9.40(8.06) | <0.01 | 9.12(10.53) | 0.81* | 4.42(5.65) |
| V40 | 0.12* | 0.38(1.15) | 0.01* | 0.21(0.46) | <0.01 | 3.43(2.41) | <0.01* | 0.99(2.02) | <0.01* | 1.16(2.95) | <0.01 | 5.80(5.44) | |
| V45 | 0.10* | 0.03(0.25) | 0.03* | 0.002(0.05) | 0.05 | 0.03(0.19) | <0.01* | 0.05(0.39) | <0.01* | 0.08(0.35) | <0.01 | 0.32(1.18) | |
| Left Femoral HeadGT | D2 | 0.71 | 4.95(6.22) | <0.01 | 5.00(5.82) | <0.01* | 10.07(8.78) | <0.01 | 8.14(9.27) | <0.01 | 8.11(7.12) | 0.20* | 5.19(7.16) |
| V40 | 0.97* | 0.63(1.17) | <0.01 | 0.85(1.23) | <0.01 | 2.63(2.56) | <0.01* | 1.94(2.89) | <0.01 | 1.86(2.58) | <0.01* | 7.17(7.15) | |
| V45 | 0.98* | 0.13(0.54) | <0.01* | 0.09(0.39) | 0.13 | 0.14(0.40) | <0.01* | 0.36(0.73) | <0.01* | 0.44(0.74) | <0.01 | 1.11(1.72) | |
Statistical analysis of dose parameter difference between test and reference plans.
a. P values (2-tailed) marked with * indicated results of paired samples Wilcoxon signed rank test for original dose parameters, while the unmarked P values were the results of paired samples t tests.
b. DAbs were the absolute difference between specific dose parameters extracted from the test plans and the reference plan, described by the median (IQR). In particular, the difference of Dn was normalized according to the δDn formula.
P values with no statistical significance were bold.
Although the volumetric and surface DSCs of all structures were numerically different, there was a correlation between them (P< 0.01). The correlation analysis results of geometric metrics and dose parameters were shown in Table 4. The volumetric and surface DSCs of CTV generated by ResUNet were correlated with all dose parameters in target volume (P< 0.05) in Plan2, on the other hand, the volumetric DSC of the three CTV groups were correlated with δD95, δD98, HI, and CI (P< 0.05) respectively in Plan 1-3 . For the bladder, the volumetric DSC of ResUNet results was correlated with all dose parameters of the bladder in Plan5, and δD2 in Plan 4-6 was correlated with both DSC metrics. There were few correlation indexes in the small intestine, only volumetric DSC vs. V15 and surface DSC vs. V45/V50. In the results of bilateral femoral heads, both DSC metrics of two CNNs were correlated with the corresponding δD2 in Plan 4-5.
Table 4
| Volumetric DSC | Surface DSC | ||||||
|---|---|---|---|---|---|---|---|
| DeepLabv3+ | ResUNet | ABAS | DeepLabv3+ | ResUNet | ABAS | ||
| CTV | δD2 | -0.31* | -0.34* | 0.11 | -0.28 | -0.38** | 0.01 |
| δD95 | -0.54** | -0.79** | -0.79** | -0.41** | -0.66** | -0.15 | |
| δD98 | -0.43** | -0.68** | -0.77** | -0.31* | -0.57** | -0.19 | |
| HI | -0.39** | -0.67** | -0.77** | -0.28 | -0.57** | -0.19 | |
| CI | -0.31* | -0.57** | -0.45** | -0.30* | -0.60** | -0.56** | |
| Bladder | V30 | -0.27 | -0.30* | -0.41** | -0.24 | -0.28 | -0.55** |
| V40 | -0.28 | -0.34* | -0.25 | -0.26 | -0.31* | -0.38* | |
| V50 | -0.23 | -0.34* | 0.25 | -0.20 | -0.27 | 0.22 | |
| δD2 | -0.37* | -0.46** | -0.31* | -0.36* | -0.41** | -0.42** | |
| Small Intestine | V15 | 0.20 | 0.17 | 0.50** | 0.01 | <0.01 | -0.08 |
| V45 | -0.05 | -0.11 | 0.14 | -0.16 | -0.20 | -0.37* | |
| V50 | -0.06 | -0.09 | 0.11 | -0.15 | -0.16 | -0.39* | |
| δD2 | 0.23 | <-0.01 | -0.06 | 0.19 | -0.11 | -0.02 | |
| Right Femoral Head | V40 | -0.21 | -0.25 | 0.21 | -0.21 | -0.26 | 0.18 |
| V45 | -0.14 | -0.14 | -0.02 | -0.14 | -0.13 | -0.02 | |
| δD2 | -0.45** | -0.51** | -0.42** | -0.46** | -0.52** | -0.30 | |
| Left Femoral Head | V40 | -0.28 | -0.20 | 0.39** | -0.31* | -0.22 | 0.44** |
| V45 | -0.17 | -0.11 | 0.16 | -0.19 | -0.12 | 0.28 | |
| δD2 | -0.42** | -0.38** | 0.07 | -0.43** | -0.40** | 0.08 | |
Pearson correlation analysis between geometric parameters and corresponding multiple dose parameters of the same structure.
The values in the table were Pearson correlation coefficient “r”, r values marked with ** indicated a significant correlation at test level 0.01 (2-tailed), r values marked with * indicated a significant correlation at test level 0.05 (2-tailed), and unmarked r values indicated no significant correlation. PTV was generated by CTV plus a uniform margin, so we tested the dose parameters of PTV with the geometric parameters of CTV.
4 Discussion
In our retrospective study, geometric and dosimetric evaluations of CTV and OARs for rectal cancer were carried out using manually segmented structures as the ground truth, while commercial software ABAS generated structures as reference. The results showed that the CNN models had a remarkable performance in accuracy and repeatability of automatic segmentation, but their performance in geometric metrics and dose parameters was not completely consistent.
The effect of automatic segmentation requires objective metrics for evaluation. Volumetric DSC is the most commonly used metric and describes the degree of overlap between two structures; however, it weights all misplaced segmentations equally and cannot characterize the distance of the ROI surface. For example, a structure with more proportion needs to be modified slightly and takes a long time may obtain a high volumetric DSC, while a structure requiring a large amount of modification locally and a short time-consuming may have a low volumetric DSC. For the description of surface distance, a commonly used metric is the HD, which represents the maximum of the shortest distances from any delineated point to the other contour (
Besides geometric accuracy, the effect of automatically generated ROIs on treatment planning should also be considered. Since the dose distribution is affected by mechanical and physical factors and cannot fully fit the edge of the structure, parts of the automatic contouring that are not perfectly consistent with clinical ground truth may be covered by the isodose lines, which perhaps can be considered as the “robustness” of the structure. Therefore, the evaluation of automatic delineation results should be combined with dosimetry results rather than a single geometric evaluation. There are many methods for dose evaluation; the simplest is to transplant a reference plan into different ROI sets and calculate the dose parameters for comparison (
At present, our institution has established an integrated platform for automatic delineation (including head, chest, abdomen, pelvis, etc.), which is connected to TPS and CT workstations through the hospital’s internal network. The platform can realize CT image transmission, conversion between RS files and masks, and continuous input of abnormal cases (such as recognizing the skull as the femoral head) to improve network performance. This is also a key step in the real application of artificial intelligence to the clinic. The auto-segmentation assessment should integrate subjective and objective methods, but the subjective assessment will introduce inter-observer variability, so it needs multi-center external validation (
5 Conclusion
In this study, we evaluated the automatic segmentation results from the perspectives of geometry and dosimetry. The results showed the advantages of speed and repeatability of deep learning in ROI delineation, which is of great help to the routine workflow of radiotherapy. The auto-segmentation function of CNNs is a stability tool for VMAT and IMRT treatment plan design, and it may have further potential in adaptive radiotherapy, which requires repeat CT scans and CTV delineation before each treatment fraction (
The characteristics of convolutional neural networks are different, and the segmentation effect on the ROIs of rectal cancer also differs. We can integrate the two networks or classify them according to the advantageous structure of each network; however, whether further exploration can bring better results requires a comprehensive clinical evaluation.
Data availability statment
The original contributions presented in the study are included in the article. Further inquiries can be directed to the corresponding author.
Statements
Ethics statement
The studies involving humans were approved by Ethics Committee on Biomedical Research, West China Hospital of Sichuan University. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Author contributions
JL, YS, and SB contributed to conception and design of the study. SB and GL contributed to administrative support and provision of study materials. YS organized the database. JL, YW, and LL collected the data. JL performed the statistical analysis. All authors contributed to the manuscript and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China, grant numbers 81972848; Cancer Precision Radiotherapy Spark Program of China International Medical Foundation, grant number 2019-N-11-04; Sichuan Province Science and Technology Support Program, grant number 2021YFS0143; Sichuan University Innovation Research Project, grant number 2022SCUH0021.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AbrahaIAristeiCPalumboILupattelliMTrastulliSCirocchiRet al. Preoperative radiotherapy and curative surgery for the management of localised rectal carcinoma. Cochrane Database Syst Rev (2018) 10(10):CD002102. doi: 10.1002/14651858.CD002102.pub3
2
QiaoliWYongpingHWeiXGuoqiangXYunheJQiuyanLet al. Preoperative short-course radiotherapy (5 × 5 Gy) with delayed surgery versus preoperative long-course radiotherapy for locally resectable rectal cancer: a meta-analysis. Int J Colorectal Dis (2019) 34(12):2171–83. doi: 10.1007/s00384-019-03433-9
3
RahbariNNElbersHAskoxylakisVMotschallEBorkUBüchlerMWet al. Neoadjuvant radiotherapy for rectal cancer: meta-analysis of randomized controlled trials. Ann Surg Oncol (2013) 20(13):4169–82. doi: 10.1245/s10434-013-3198-9
4
LobefaloFBignardiMReggioriGTozziATomatisSAlongiFet al. Dosimetric impact of inter-observer variability for 3D conformal radiotherapy and volumetric modulated arc therapy: the rectal tumor target definition case. Radiat Oncol (2013) 8:176. doi: 10.1186/1748-717X-8-176
5
ClivioAClivioAFogliataAFranzetti-PellandaANicoliniGVanettiEet al. Volumetric-modulated arc radiotherapy for carcinomas of the anal canal: A treatment planning comparison with fixed field IMRT. Radiother Oncol (2009) 92(1):118–24. doi: 10.1016/j.radonc.2008.12.020
6
TaoCJYiJLChenNYRenWChengJTungSet al. Multi-subject atlas-based auto-segmentation reduces interobserver variation and improves dosimetric parameter consistency for organs at risk in nasopharyngeal carcinoma: A multi-institution clinical study. Radiother Oncol (2015) 115(3):407–11. doi: 10.1016/j.radonc.2015.05.012
7
KosminMLedsamJRomera-ParedesBMendesRMoinuddinSde SouzaDet al. Rapid advances in auto-segmentation of organs at risk and target volumes in head and neck cancer. Radiother Oncol (2019) 135:130–40. doi: 10.1016/j.radonc.2019.03.004
8
McQuinlanYBrouwerCLLinZGanYSung KimJvan ElmptWet al. An investigation into the risk of population bias in deep learning autocontouring. Radiother Oncol (2023) 186:109747. doi: 10.1016/j.radonc.2023.109747
9
LiaoWLuoXHeYDongYLiCLiKet al. Comprehensive evaluation of a deep learning model for automatic organs at risk segmentation on heterogeneous computed tomography images for abdominal radiation therapy. Int J Radiat Oncol Biol Phys (2023) S0360-3016(23):00520–5. doi: 10.1016/j.ijrobp.2023.05.034
10
WangTLeiYTianSJiangXZhouJLiuTet al. Learning-based automatic segmentation of arteriovenous malformations on contrast CT images in brain stereotactic radiosurgery. Med Phys (2019) 46(7):3133–41. doi: 10.1002/mp.13560
11
BurbachJPKleijnenJPReerinkOSeravalliEPhilippensMESchakelTet al. Inter-observer agreement of MRI-based tumor delineation for preoperative radiotherapy boost in locally advanced rectal cancer. Radiother Oncol (2016) 118(2):399–407. doi: 10.1016/j.radonc.2015.10.030
12
van DijkLVVan den BoschLAljabarPPeressuttiDBothSSteenbakkers RJHMet al. Improving automatic delineation for head and neck organs at risk by Deep Learning Contouring. Radiother Oncol (2020) 142:115–23. doi: 10.1016/j.radonc.2019.09.022
13
SongYHuJWuQXuFNieSZhaoYet al. Automatic delineation of the clinical target volume and organs at risk by deep learning for rectal cancer postoperative radiotherapy. Radiother Oncol (2020) 145:186–92. doi: 10.1016/j.radonc.2020.01.020
14
NikolovSBlackwellSZverovitchAMendesRLivneMDe FauwJet al. Clinically applicable segmentation of head and neck anatomy for radiotherapy: deep learning algorithm development and validation study. J Med Internet Res (2021) 23(7):e26151. doi: 10.2196/26151
15
StanleyJDunscombePLauHBurnsPLimGLiuHWet al. The effect of contouring variability on dosimetric parameters for brain metastases treated with stereotactic radiosurgery. Int J Radiat Oncol Biol Phys (2013) 87(5):924–31. doi: 10.1016/j.ijrobp.2013.09.013
16
GrégoireVMackieTR. State of the art on dose prescription, reporting and recording in Intensity-Modulated Radiation Therapy (ICRU report No. 83) Cancer Radiother (2011) 15(6-7):555–9. doi: 10.1016/j.canrad.2011.04.003
17
KiliçDÇatliSÜlgerŞKapucuLÖ.Is there any impact of PET/CT on radiotherapy planning in rectal cancer patients undergoing preoperative IMRT? Turk J Med Sci (2015) 45(1):129–35. doi: 10.3906/sag-1312-149
18
RobertsonJMLockmanDYanDWallaceM. The dose-volume relationship of small bowel irradiation and acute grade 3 diarrhea during chemoradiotherapy for rectal cancer. Int J Radiat Oncol Biol Phys (2008) 70(2):413–8. doi: 10.1016/j.ijrobp.2007.06.066
19
GallagherMJBreretonHDRostockRAZeroJMZekoskiDAPoyssLFet al. A prospective study of treatment techniques to minimize the volume of pelvic small bowel with reduction of acute and late effects associated with pelvic irradiation. Int J Radiat Oncol Biol Phys (1986) 12(9):1565–73. doi: 10.1016/0360-3016(86)90279-8
20
LetschertJGLebesqueJVde BoerRWHartAABartelinkH. Dose-volume correlation in radiation-related late small-bowel complications: a clinical study. Radiother Oncol (1990) 18(4):307–20. doi: 10.1016/0167-8140(90)90111-9
21
ShiraishiSMooreKL. Knowledge-based prediction of three-dimensional dose distributions for external beam radiotherapy. Med Phys (2016) 43(1):378. doi: 10.1118/1.4938583
22
SongYHuJLiuYHuHHuangYBaiSet al. Dose prediction using a deep neural network for accelerated planning of rectal cancer radiotherapy. Radiother Oncol (2020) 149:111–6. doi: 10.1016/j.radonc.2020.05.005
23
ShaXWangHShaHXieLZhouQZhangWet al. Clinical target volume and organs at risk segmentation for rectal cancer radiotherapy using the Flex U-Net network. Front Oncol (2023) 13:1172424. doi: 10.3389/fonc.2023.1172424
24
QiuWZhangWMaXKongYShiPFuMet al. Auto-segmentation of important centers of growth in the pediatric skeleton to consider during radiation therapy based on deep learning. Med Phys (2023) 50(1):284–96. doi: 10.1002/mp.15919
25
JansenEPNijkampJGubanskiMLindPAVerheijM. Interobserver variation of clinical target volume delineation in gastric cancer. Int J Radiat Oncol Biol Phys (2010) 77(4):1166–70. doi: 10.1016/j.ijrobp.2009.06.023
26
XiaXWangJLiYPengJFanJZhangJet al. An artificial intelligence-based full-process solution for radiotherapy: A proof of concept study on rectal cancer. Front Oncol (2021) 10:616721. doi: 10.3389/fonc.2020.616721
27
DefeudisAMazzettiSPanicJMicilottaMVassalloLGiannettoGet al. MRI-based radiomics to predict response in locally advanced rectal cancer: comparison of manual and automatic segmentation on external validation in a multicenter study. Eur Radiol Exp (2022) 6(1):19. doi: 10.1186/s41747-022-00272-2
28
NijkampJMarijnenCvan HerkMvan TriestBSonkeJJ. Adaptive radiotherapy for long course neo-adjuvant treatment of rectal cancer. Radiother Oncol (2012) 103(3):353–9. doi: 10.1016/j.radonc.2012.02.013
29
LiLXuBZhuangZLiJHuYYangHet al. Accurate tumor segmentation and treatment outcome prediction with DeepTOP. Radiother Oncol (2023) 183:109550. doi: 10.1016/j.radonc.2023.109550
30
BretoALSpielerBZavala-RomeroOAlhusseiniMPatelNVAsherDAet al. Deep learning for per-fraction automatic segmentation of gross tumor volume (GTV) and organs at risk (OARs) in adaptive radiotherapy of cervical cancer. Front Oncol (2022) 12:854349. doi: 10.3389/fonc.2022.854349
Summary
Keywords
automatic segmentation, deep learning, rectal neoplasms, radiotherapy, CNN
Citation
Li J, Song Y, Wu Y, Liang L, Li G and Bai S (2023) Clinical evaluation on automatic segmentation results of convolutional neural networks in rectal cancer radiotherapy. Front. Oncol. 13:1158315. doi: 10.3389/fonc.2023.1158315
Received
03 February 2023
Accepted
11 August 2023
Published
05 September 2023
Volume
13 - 2023
Edited by
Gozde Akar, Middle East Technical University, Türkiye
Reviewed by
Qiu Tang, School of Medicine, Zhejiang University, China; Jiazhou Wang, Fudan University, China
Updates

Check for updates
Copyright
© 2023 Li, Song, Wu, Liang, Li and Bai.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sen Bai, baisen@scu.edu.cn
†These authors have contributed equally to this work and shared first authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.