1 October 2018 Evaluation of deep learning methods for parotid gland segmentation from CT images
Author Affiliations +
Abstract
The segmentation of organs at risk is a crucial and time-consuming step in radiotherapy planning. Good automatic methods can significantly reduce the time clinicians have to spend on this task. Due to its variability in shape and low contrast to surrounding structures, segmenting the parotid gland is challenging. Motivated by the recent success of deep learning, we study the use of two-dimensional (2-D), 2-D ensemble, and three-dimensional (3-D) U-Nets for segmentation. The mean Dice similarity to ground truth is ∼0.83 for all three models. A patch-based approach for class balancing seems promising for false-positive reduction. The 2-D ensemble and 3-D U-Net are applied to the test data of the 2015 MICCAI challenge on head and neck autosegmentation. Both deep learning methods generalize well onto independent data (Dice 0.865 and 0.88) and are superior to a selection of model- and atlas-based methods with respect to the Dice coefficient. Since appropriate reference annotations are essential for training but often difficult and expensive to obtain, it is important to know how many samples are needed for training. We evaluate the performance after training with different-sized training sets and observe no significant increase in the Dice coefficient for more than 250 training cases.

1.

Introduction

High radiation dose on organs at risk (OAR) during radiotherapy (RT) treatment can have severe side effects. An important set of OAR during treatment of head and neck cancer are the parotid glands, which are a type of salivary glands that are very sensitive to radiation. The most common radiation-induced side effect on the parotid glands is xerostomia (dry mouth), which can significantly decrease life quality.1 Therefore, modern radiotherapy planning requires accurate segmentation of target structures and OAR for precise and highly localized dose planning.2 As manual segmentation of the planning images is very time-consuming and user-dependent, radiotherapy planning could highly benefit from automatic methods for contouring. However, fully automatic segmentation of the parotid glands from computed tomography (CT) images is difficult due to their high variability in shape and appearance and often low soft-tissue contrast to surrounding structures. Because of their anatomical location, they are also prone to be affected by dental metal artifacts.

Various methods other than deep learning have been proposed for automatic segmentation of OAR in the head and neck region. The approaches include, among others, (multi) atlas-based methods,34.5.6 model-based methods,78.9.10 or their combinations.11,12 Some of the methods6,89.10.11 have been evaluated in the 2015 MICCAI challenge on head and neck autosegmentation, where the best mean Dice score on the parotid glands was 0.84.13

Since about 2012, deep learning methods have been extensively used for medical image processing problems, with a remarkable proportion of published papers on segmentation applications.14,15 Some of the most popular deep neural network architectures for segmentation include the U-Net,16,17 V-Net,18 and multiresolution architectures such as DeepMedic19 and F-Net.20 Some approaches have also been made with recurrent neural networks,21 generative adversarial networks,22 and other neural network architectures proposed by the computer vision community. Moreover, deep neural networks have proven to be powerful in several recent segmentation challenges such as brain tumor segmentation,23 liver tumor segmentation,24 and ischemic stroke lesion segmentation.25

Deep learning methods have also already been applied to various steps in the RT workflow, including but not limited to automatic contouring.26,27 In first clinical validations, deep learning segmentation methods have shown to decrease the total time required for OAR contouring in comparison to manual and atlas-based contouring.28 Head and neck segmentation from CT images has also been previously addressed using deep learning.29,30 However, while Ibragimov and Xing30 could observe superior performance on many OAR to previously reported results using their neural network, the performance on parotid gland segmentation only was comparable to previously reported Dice coefficients in the literature. Fritscher et al.29 also reported on the promising results of using deep learning methods for head and neck segmentation. All of these recent successes of deep learning in various medical domains motivate us to further explore deep learning for parotid gland segmentation from CT images.

In this contribution, we study parotid gland segmentation using deep learning, focusing on the U-Net architecture16,17 that we apply in two-dimensional (2-D), three-dimensional (3-D), and in a 2-D ensemble mode. We also evaluate the performance of the trained neural networks on the publicly available test dataset of the 2015 MICCAI challenge on head and neck autosegmentation and compare to the model and atlas-based methods that competed in the challenge. We show that the deep learning approaches yield superior results based on the obtained Dice coefficients, which are statistically significant for all but one compared method. Another important aspect of deep learning is the availability of a large, annotated training dataset that ideally captures most of the anatomical variability. In practice, medical datasets are often rather small compared to the typical datasets used in the computer vision community such as the ImageNet dataset.31 We investigate the influence of the number of training samples on the resulting Dice coefficients on a validation set.

2.

Methods

2.1.

Data

The image data used in the study consisted of 254 head and neck CT scans from two different clinical sites. For each CT scan, clinical routine level, uncurated reference segmentations of the left and right parotid glands, created by a single physician per image, were available. As the contours come from clinical routine, several doctors were involved in the contour creation at both clinical sites. The axial in-plane resolution was either 0.977 or 1.172 mm for all but two images, for which it was 1.219 mm. The slice spacing was either 2 or 3 mm. Due to the clinical uncurated nature of the contours, inconsistencies in the contouring of different scans are to be expected.

To simplify the segmentation problem, we decided to focus on a binary segmentation task (see Sec. 2.3). Therefore, all neural networks were trained to segment the left parotid gland only. The image data and reference segmentations of the right parotid gland were used for data augmentation by mirroring it and thus included into the training and validation. In total, 507 reference segmentations (253 left side, 254 right side mirrored, 1 left parotid gland was resected) were used and divided into 467 examples for training and 40 examples for validation.

In addition, for testing on an independent and publicly available test set, the 10 off-site and 5 on-site test cases of the 2015 MICCAI challenge on head and neck autosegmentation13 were used. The test data, including image data and carefully drawn segmentations based on best practices and scientific literature, are available at the Public Domain Database for Computational Anatomy.32

2.2.

Preprocessing

The axial in-plane resolution of the training and validation data was close to 1×1  mm2 for all images but the slice spacing was either 2 mm or 3 mm. Therefore, in order to preserve details in the axial plane but unify the resolution along the transverse axis, the data were resampled along the transverse axis only to 2 mm slice spacing using a Lanczos kernel with support size 3. For inference, the independent MICCAI challenge test data were resampled to 1×1×2  mm3 to match the resolution used during training of the neural networks. Furthermore, the data were preprocessed by automatically removing the treatment couch via masking of the patient volume.

2.3.

Neural Networks

We chose the U-Net architecture16,17 as the basis of all experiments. Three different U-Net models were implemented: a 2-D U-Net trained on axial slices, a 3-D U-Net, and a 2-D U-Net ensemble, as schematically shown in Fig. 1. The 2-D U-Net ensemble consisted of three 2-D U-Nets trained on axial, coronal, and sagittal reformatted slices, respectively, whose predictions were combined via a majority voting. This ensemble approach for combining orthogonal view directions is an alternative to 2.5-D networks, where multiple slices from orthogonal image planes are simultaneously fed into one neural network.33 All neural networks were trained to solve the binary segmentation task of segmenting the left parotid gland only. As a consequence, they had two output channels for foreground and background, followed by a softmax layer. All neural networks were implemented in-house in accordance with the architecture descriptions by Ronneberger et al.16 and Ciçek et al.,17 using the deep learning framework Lasagne.34 However, in contrast to the U-Net originally described in Ref. 16, batch normalization35 was also used for the 2-D U-Nets and only three of the U-Net’s internal resolution levels were implemented. The resulting U-Net has a maximum receptive field size of 44×44 voxels in the 2-D case and 44×44×44 voxels in the 3-D case, which contribute to the classification of a single voxel. For an image resolution of 1×1×2  mm3, this corresponds to a receptive field size of 44×44×88  mm3, as visualized in Fig. 2. The receptive field depends on the order and number of convolutional and pooling layers and their kernel sizes and should not be confused with the patch size used for training.36 From Fig. 2, the receptive field size seems sufficient as the parotid gland is almost completely inside the receptive field for a voxel chosen at the gland’s center. Moreover, the resulting U-Nets have fewer parameters than a U-Net with more resolution levels, hence they are easier to train on a GPU with limited memory.

Fig. 1

Schematic drawing of the different neural network architectures used: (1) 2-D U-Net working on 2-D patches, (2) 3-D U-Net working on 3-D patches, and (3) ensemble of three 2-D U-Nets working on axial, coronal, and sagittal patches, respectively. The individual predictions of the ensemble are combined to a single prediction via a majority vote.

JMI_6_1_011005_f001.png

Fig. 2

Receptive field visualization of a 3-D U-Net with three resolution levels: for the classification of the voxel marked with a cross hair in (a), the context that the neural network can use for classification is a cuboid of size 44×44×44 voxels or 44×44×88  mm3. The central planes through the cuboid in axial, sagittal and coronal view are visualized in (b).

JMI_6_1_011005_f002.png

2.4.

Training

All training of neural networks was performed on a desktop computer with an NVIDIA GeForce GTX 1080 graphics card with 8 GB graphical memory. Batch sizes, number of epochs, and learning rates were adapted to the memory consumption and training convergence of each individual neural network.

As the structure to be segmented is small compared to the full volume, a strategy for increasing the foreground percentage during training is necessary. Otherwise, the network will mostly see background and not learn the segmentation task. The problem is addressed to some extent by using a soft Dice loss for training of all U-Nets, which accounts for low foreground percentage during the loss calculation.18 Furthermore, two different strategies to actually increase the foreground percentage during training were applied for the 2-D and 3-D U-Nets, respectively.

2.4.1.

Training in a region of interest

For the 2-D U-Nets, a region of interest (ROI)-based training strategy was applied. During training, slices to be presented to the neural network were drawn only from within the ROI. It was defined by any slices (axial, sagittal, or coronal reformatted depending on the view direction of the neural network) containing the left or mirrored right parotid gland. Additionally, the U-Net should learn not to segment anything on slices directly neighboring the organ. Therefore, an additional margin of five slices in view direction was included into the ROI. As the scan range of head and neck scans can be large, a margin of 25 voxels above and below the reference segmentation was additionally applied to sagittal and coronal reformatted slices in order to further increase the foreground percentage.

2.4.2.

Patch-based training with controlled batch composition

For the 3-D U-Net, the ROI training strategy could not be applied as the 3-D ROI did not fit into GPU memory together with the neural network. Therefore, a patch-based training with controlled batch composition was implemented. Patches of size 72×72×56 voxels were extracted from the full volume instead of an ROI. From these patches, the center 32×32×16 voxels were classified and the remainder was padding, which was needed for the U-Net’s specific architecture with unpadded convolutions, where the image size shrinks with each convolutional layer. The padding also serves to provide context for the output neurons toward the borders of the unpadded patch. Class balancing of foreground and background was performed by composing each minibatch during training so that 50% of the patches overlapped with the parotid gland.

In addition to training on the full training dataset with 467 cases, the 2-D and 3-D U-Nets were trained on subsets of 50, 150, 250, 350, and 450 training cases in order to investigate the impact of training data quantity on the segmentation results.

2.5.

Postprocessing

The U-Nets’ predictions were binarized by thresholding at 0.5. A connected component (CC) analysis was performed and the largest component was taken as the final segmentation result with the aim to eliminate small false-positive findings. In our application, this basic postprocessing seems to be sufficient, as we expect to find a single left parotid gland inside each image volume as the largest segmented structure. Furthermore, the use of the Dice loss for training produces a raw neural network output, which is already close to binary with values either very close to 0 or 1. Therefore, no complex algorithm for converting soft predictions to hard labels is necessary in this case. In more complex detection tasks such as brain lesion segmentation, where the number of structures to be detected is unknown and prediction maps may be soft, additional advanced postprocessing using conditional random fields or other machine learning algorithms can be necessary.19

Finally, all segmentation results were resampled to the original image resolution using nearest-neighbor interpolation, in order to be able to compare to the original reference contours.

2.6.

Evaluation

The quantitative evaluation was mainly based on the Dice coefficient.37 For two sets of voxels X and Y, which represent segmentation result and reference segmentation, it is given by

(1)

DSC(X,Y)=2|XY||X|+|Y|.
To evaluate the overall over- versus underestimation in Sec. 3.2, we further compute an average signed surface distance (ASSD) integrated into the development platform MeVisLab.38 For sets of voxels X and Y, it is given as

(2)

ASSD(X,Y)=1|X|+|Y|[xXd(x,Y)yYd(y,X)],
where |X| denotes the border voxels of X and

(3)

d(x,Y)={infyYxy2,if  xYCinfyYxy2,if  xY,
is the signed Euclidean distance of a voxel x to the set Y. The distance measure is derived from the average or mean surface distance often used for segmentation evaluation and challenges.13,39 However, instead of calculating absolute differences, signed distances are used. This way, the distance measure is positive outside of the reference mask and negative inside it, therefore it can be used to distinguish over- and underestimation. For statistical evaluation, we used the Wilcoxon-signed rank test at significance level 0.05 to compare the performance of different neural networks.

3.

Results

The results are presented in four parts. First, the 2-D U-Net, 2-D U-Net ensemble, and 3-D U-Net are evaluated on the validation data and the effect of the training strategies is described. Second, the performance of the 2-D ensemble and 3-D U-Net on validation versus MICCAI test data is examined. Then, the deep learning results are compared to the quantitative results of the 2015 MICCAI challenge. Finally, the impact of the number of training samples on the 2-D ensemble and 3-D U-Net performance is presented.

3.1.

Comparison of 2-D, 2-D Ensemble, and 3-D U-Net Performance and Training Strategies

Figure 3 shows the impact of the ROI-based training strategy onto the axial 2-D U-Net’s segmentation performance. The Dice coefficient was calculated on all 40 validation cases once on an ROI around the target structure (extracted the same way as during training) and once on the full volume. Additionally, it was computed before and after the postprocessing step of automatically selecting the largest CC. Within the ROI, the 2-D U-Net achieves a median Dice coefficient of 0.830 before and 0.831 after CC selection. The difference between the two median Dice scores is very small, as there are very few false positives next to the segmented parotid gland within the ROI. However, when evaluated on the full volume, the median Dice coefficient drops significantly to 0.542 before and 0.725 after CC selection, with 18 cases with a Dice coefficient of 0 for the latter. This means that within the ROI, the 2-D U-Net can on average accurately segment the parotid gland, but there are many outliers in distant body regions that are also large in volume. This could have been expected, as with the ROI-based training strategy, the U-Net never saw patches from outside the ROI.

Fig. 3

Dice coefficients of contours generated on the validation data by the axial 2-D U-Net when evaluated on an ROI versus the full volume. Dice scores were calculated once before and once after automatic selection of the largest CC in the binarized result mask.

JMI_6_1_011005_f003.png

In Fig. 4, the 2-D U-Net performance is compared to the 2-D ensemble and the 3-D U-Net. With the idea to provide an automatic ROI detection in the future, the 2-D U-Net performance is shown on the ROI only as discussed before, whereas the other two neural networks models are evaluated on the full volume of all validation cases. All Dice coefficients in the plot were computed after the CC selection step. The median Dice coefficients of the 2-D ensemble (0.835) and the 3-D U-Net (0.830) are comparable to the 2-D U-Net (0.831) and the differences shown in Fig. 4 are not significant (p>0.05). The main observation is that both 2-D ensemble and 3-D U-Net segmentations contain only small false positives, which can be eliminated via CC selection so that the Dice coefficient on the full volume is comparable to that of the 2-D U-Net on the ROI. In the case of the 2-D ensemble, the combination of three individual predictions in orthogonal view directions can eliminate false positives of each single prediction. In the 3-D case, the patch-based training strategy helps to eliminate large false positives outside the ROI, as during the training, background patches from the whole training volumes are presented to the neural network.

Fig. 4

Dice coefficients of results generated by the axial 2-D U-Net (on an ROI only), the 2-D U-Net ensemble, and the 3-D U-Net on the validation data.

JMI_6_1_011005_f004.png

Figure 5 contains reference segmentations and deep learning results on selected exemplary validation cases that highlight typical observations. First of all, the segmentations by the three different neural networks are in general very similar, as also reflected in the Dice coefficients in Fig. 4. Most differences can be observed at the elongated lateral part of the parotid gland (first row), which is challenging for all neural networks. One common inconsistency between reference and autogenerated contour that can be observed in the second row is the inclusion of vessels at the medial part, which is consistent among all neural networks. This might be due to inconsistent annotations in the reference set. All trained neural networks are robust to dental metal artifacts (third row). In the exemplary case, they generate contours that even seem superior to the reference contours especially at the lateral anterior part. The fourth row shows contours of the third case in a sagittal plane. All neural networks, also the pure 2-D axial network, produce smooth contours in contrast to the reference segmentation, which often has an anatomically inconsistent shape in sagittal and coronal reformatted views. This might partially be due to resampling of the contours to the original image grid, which takes the segmentation from neighboring slices into account.

Fig. 5

Reference contours and segmentation results by the axial 2-D U-Net, the 2-D U-Net ensemble, and the 3-D U-Net for a selection of validation cases (one case per row).

JMI_6_1_011005_f005.png

3.2.

Comparison of Results on Validation and Independent Test Data

Both the 2-D ensemble and the 3-D U-Net were evaluated on the 15 test cases of the 2015 MICCAI challenge on head and neck autosegmentation. The axial 2-D U-Net was omitted due to the inability to compute the correct segmentation on the full volume after CC selection as seen in Sec. 3.1. Figure 6(a) shows the resulting Dice coefficients on the MICCAI test data compared to the results on the validation data from Sec. 3.1. The median Dice coefficient on the MICCAI data (0.865/0.880 for 2-D ensemble/3-D U-Net) is higher than on the validation data (0.835/0.830) and the variance is lower. This might be because uncurated data from clinical routine was used for validation (and training) that may contain errors that lead to lower Dice coefficients even for correct segmentations by the deep neural network. In Fig. 6(b), the average signed surface distance of the segmentation results to the reference is plotted. On the validation data, the median is close to zero for both 2-D ensemble (0.172 mm) and 3-D U-Net (0.239 mm), which means there is little bias toward under- or overestimation. The 2-D ensemble underestimates (negative ASSD) all but one MICCAI test case, the 3-D ensemble underestimates 21 of the 30 test cases. One common mistake of both neural network models is the underestimation of the medial part of the parotid gland that can be observed in Fig. 7. This might be due to the use of different contouring guidelines for the creation of the training versus MICCAI test contours, which leads to a systematic error during inference and Dice computation.

Fig. 6

(a) Dice coefficients and (b) ASSD of segmentation results generated by the 2-D ensemble and the 3-D U-Net on the validation data and on the independent MICCAI test data.

JMI_6_1_011005_f006.png

Fig. 7

Reference contours and segmentation results by the 2-D U-Net ensemble and the 3-D U-Net for three MICCAI test cases (one case per row).

JMI_6_1_011005_f007.png

3.3.

Comparison to Results of the 2015 MICCAI Challenge on Head and Neck Autosegmentation

In Fig. 8, the Dice coefficients computed from the segmentations by the 2-D ensemble and the 3-D U-Net are plotted together with the results of the MICCAI challenge on head and neck autosegmentation held in 2015. The original challenge results were published by Raudaschl et al.13 and kindly provided to us for comparison to our deep learning results. Results are plotted for the challenge’s on-site, off-site, and combined (on- and off-site) test cases separately, as some of the methods differ significantly for the on-site and off-site test cases. The two neural networks, however, yield similar quantitative results on each subset (see Table 1), which implies that the test cases in the two subsets are of similar difficulty. On the combined test set, the 3-D U-Net segmentations result in significantly (p<0.05) higher Dice coefficients than all methods in the challenge. For the 2-D ensemble, all differences are significant except for the comparison to Mannion-Haworth et al. (p=0.786). One has to keep in mind that in contrast to the methods in the challenge, our neural network models were trained on a completely independent training dataset. These results are promising as they demonstrate that deep learning methods can outperform other methods and are robust even when applied to independent data.

Fig. 8

Dice coefficients of contours produced by the 2-D ensemble and the 3-D U-Net on the 2015 MICCAI challenge test data in comparison to the methods that competed in the challenge. For the challenge, the test data were split into off-site and on-site test cases, which is reflected in the plot.

JMI_6_1_011005_f008.png

Table 1

Median Dice coefficients of the teams in the 2015 MICCAI challenge and the U-Net approaches for on-site, off-site, and combined test dataset.

Segmentation method/test casesOff-siteOn-siteAll
Albrecht et al.110.8330.8050.827
Chen et al.60.8250.8230.825
Jung et al.80.6990.6870.694
Mannion-Haworth et al.90.8540.8460.852
Orbes-Arteaga et al.100.8030.7460.795
2-D U-Net ensemble0.8670.8630.865
3-D U-Net0.8800.8800.880

3.4.

Impact of Training Data Quantity

Figure 9 shows the results of the training of the 2-D ensemble and 3-D U-Net on increasing amounts of training data. Dice scores are reported on the validation data only. Even though the differences in the Dice coefficients on the validation data are small, the median Dice coefficient increases significantly (p<0.05) from 50 to 150 to 250 training cases and for both models. When adding more training cases, the median Dice coefficient for the 3-D U-Net still increases, but the differences are not significant for either neural network.

Fig. 9

Dice coefficients of contours produced by the 2-D ensemble and the 3-D U-Net on the validation data when trained using an increasing number of training samples.

JMI_6_1_011005_f009.png

4.

Discussion

While the results for all investigated neural network models were similar throughout our experiments, the 3-D U-Net seems most promising with respect to contour quality and also model simplicity compared to an ensemble approach with three distinct models. Moreover, the ensemble sometimes produces artifacts at the contour borders due to the majority voting, which will require further postprocessing such as smoothing. An important aspect to keep in mind when comparing the different approaches are the different numbers of parameters in the 2-D and 3-D U-Nets. With more parameters, the 3-D U-Net has a higher learning capacity than the 2-D U-Net and the ensemble, which may also partially explain the observed differences in performance. An option would be to use a deeper 2-D U-Net with one more resolution level, so that more parameters are available to the 2-D neural network. However, this would automatically increase the receptive field size of the 2-D U-Net compared to the 3-D U-Net (with respect to the axial plane), so that again no fair comparison could be granted.

In general, fair comparison of different neural network architectures with different numbers of parameters is difficult, in the sense that changing one parameter such as network depth directly impacts other network properties (such as receptive field). Hence, while comparing to other popular architectures different from the U-Net may yield interesting results, it would not be evident why one method outperforms the other and which conclusions to draw. Therefore, comparison of different architectures and training strategies as we have done in this work may give first directions toward which network architecture or sampling strategy should be investigated in more detail. The chosen parameters should then be optimized in a hyperparameter search, in which only minimal changes are made to the training. That way, the impact of single parameters can be analyzed in a fair comparison. For example, future work could focus on binary versus nonbinary segmentation with the 3-D U-Net, such as simultaneous segmentation of right and left parotid gland or inclusion of further OAR in the head and neck region. Such a neural network would have only slightly more parameters due to the increased number of output channels but produce multiple segmentations in only slightly increased amount of time. It should then be investigated how the generalization capacity of a neural network changes, when it has to segment several distinct anatomical structures instead of a single one.

Furthermore, as the results produced by different neural networks can be similar as observed in Fig. 5, it would be interesting to compute local, application-specific distance measures instead of global general ones such as Dice score or surface distances. For example, one could then focus on the lateral elongated part of the parotid gland and assess model robustness with respect to this specific anatomical feature.

A limitation of quantitative results on the validation data as seen in Secs. 3.1 and 3.4 is the fact that they are not completely reliable as uncurated reference contours from clinical routine were used. The effect of this could be seen in Fig. 6, where the variance of the Dice scores on the uncurated validation data was much higher than on the curated test data. Assuming that the validation cases are not significantly harder to segment than the test cases, the increased variance may likely stem from the uncertainty in the reference masks. For future work, having a subset of the data for validation with well-curated reference contours would be highly desirable for quantitative evaluation. Still, the difference in results between the neural networks indicate that a patch-based training strategy with controlled batch composition is to be preferred over an ROI sampling strategy with respect to false-positive reduction. An automatic ROI detection, e.g., using a smaller neural network working on a coarse resolution would, however, still be desirable in order to speed up inference.

The use of clinical quality reference contours in Sec. 3.4 also poses the question whether the amount of available data is sufficient for the problem at hand as one might read from the plots. Another interpretation is that with the currently available clinical reference quality, no better results can be achieved no matter how much more data of similar quality is added to the training. We have conducted first experiments on training with curated versus clinical quality reference contours.40 Our results suggest that the difference in segmentation performance between training on a large dataset of clinical quality references versus a smaller dataset of curated quality references is small. After this rather technical study, an important further step will be to do a clinical qualitative evaluation of the clinical acceptance of the contours generated using deep learning. It needs to be investigated whether contours derived from clinical uncurated quality data via deep learning can fulfill clinical contouring guidelines. Another train of thought is that deeper networks with more parameters may also be able to integrate more information and further improve the segmentation performance with more samples. However, a fair comparison to a deeper neural network would be difficult due to the problems already discussed.

Finally, the results of Sec. 3.3 demonstrate the generalization potential of deep learning methods for autocontouring of CT data from different sites. We suppose that this is primarily due to the availability of the image data in well-calibrated Hounsfield units. Autosegmentation of MR planning images or cone-beam CT images for adaptive radiotherapy could be more challenging with respect to homogeneous datasets and generalization. In that case, training of neural networks for a specific site and scanner may be necessary. Moreover, the deep learning methods show a higher robustness or lower variance in performance compared to the atlas- and model-based methods in the segmentation challenge. This may be because the learned features can represent a large anatomical variety without prior assumptions and have additionally been trained on a larger dataset than the methods in the challenge.

5.

Conclusion

We have presented results on parotid gland segmentation using deep learning, including an analysis of the amount of training data. Using the U-Net architecture and using uncurated reference segmentations from clinical routine for training, we can achieve results on a publicly available test set that are significantly better than those by several model- and atlas-based methods. This demonstrates the high potential of deep learning-based autosegmentation methods for radiotherapy planning, where manual or semiautomatic contouring is still a bottleneck in the workflow.

Disclosures

Tobias Gass, Tomasz Morgas, and Benjamin Haas are employed by Varian Medical Systems. Varian Medical Systems partly funded the research conducted by employees of Fraunhofer MEVIS.

Acknowledgments

The work presented was funded by Varian Medical Systems and the Fraunhofer project QuantMed. We thank Patrik F. Raudaschl for providing quantitative result data of the 2015 MICCAI challenge on head and neck autosegmentation,13 which we used for comparison to our results in Sec. 3.3. An earlier version of parts of this manuscript has been published in the conference proceedings of SPIE Medical Imaging 2018.41

References

1. 

Grundmann O., Mitchell G. C., Limesand K. H., “Sensitivity of salivary glands to radiation: from animal models to therapies,” J. Dent. Res., 88 (10), 894 –903 (2009). https://doi.org/10.1177/0022034509343143 JDREAF 0022-0345 Google Scholar

2. 

Ng J., Shuryak I., “Minimizing second cancer risk following radiotherapy: current perspectives,” Cancer Manage. Res., 7 1 –11 (2015). https://doi.org/10.2147/CMAR.S47220 Google Scholar

3. 

Han X. et al., “Atlas-based auto-segmentation of head and neck CT images,” Lect. Notes Comput. Sci., 5242 434 –441 (2008). https://doi.org/10.1007/978-3-540-85990-1 LNCSD9 0302-9743 Google Scholar

4. 

Fortunati V. et al., “Tissue segmentation of head and neck CT images for treatment planning: a multiatlas approach combined with intensity modeling,” Med. Phys., 40 (7), 071905 (2013). https://doi.org/10.1118/1.4810971 MPHYA6 0094-2405 Google Scholar

5. 

Wachinger C. et al., “Contour-driven atlas-based segmentation,” IEEE Trans. Med. Imaging, 34 2492 –2505 (2015). https://doi.org/10.1109/TMI.2015.2442753 ITMID4 0278-0062 Google Scholar

6. 

Chen A., Dawant B., “A multi-atlas approach for the automatic segmentation of multiple structures in head and neck CT images,” MIDAS J., (2015) http://hdl.handle.net/10380/3540 Google Scholar

7. 

Qazi A. A. et al., “Auto-segmentation of normal and target structures in head and neck CT images: a feature-driven model-based approach,” Med. Phys., 38 (11), 6160 –6170 (2011). https://doi.org/10.1118/1.3654160 MPHYA6 0094-2405 Google Scholar

8. 

Jung F., Knapp O., Wesarg S., “Automatic segmentation of structures in CT head and neck images using a coupled shape model,” MIDAS J., (2015) http://hdl.handle.net/10380/3543 Google Scholar

9. 

Mannion-Haworth R. et al., “Fully automatic segmentation of head and neck organs using active appearance models,” MIDAS J., (2015) http://hdl.handle.net/10380/3537 Google Scholar

10. 

Orbes-Arteaga M., Pea D. C., Dominguez G. C., “Head and neck auto segmentation challenge based on non-local generative models,” MIDAS J., (2015) http://hdl.handle.net/10380/3539 Google Scholar

11. 

Albrecht T. et al., “Multi atlas segmentation with active shape model refinement for multi-organ segmentation in head and neck cancer radiotherapy planning,” MIDAS J., (2015) http://hdl.handle.net/10380/3531 Google Scholar

12. 

Fritscher K. D. et al., “Automatic segmentation of head and neck CT images for radiotherapy treatment planning using multiple atlases, statistical appearance models, and geodesic active contours,” Med. Phys., 41 (5), 051910 (2014). https://doi.org/10.1118/1.4871623 MPHYA6 0094-2405 Google Scholar

13. 

Raudaschl P. F. et al., “Evaluation of segmentation methods on head and neck CT: auto-segmentation challenge 2015,” Med. Phys., 44 (5), 2020 –2036 (2017). https://doi.org/10.1002/mp.12197 MPHYA6 0094-2405 Google Scholar

14. 

Litjens G. et al., “A survey on deep learning in medical image analysis,” Med. Image Anal., 42 60 –88 (2017). https://doi.org/10.1016/j.media.2017.07.005 Google Scholar

15. 

Shen D., Wu G., Suk H.-I., “Deep learning in medical image analysis,” Annu. Rev. Biomed. Eng., 19 221 –248 (2017). https://doi.org/10.1146/annurev-bioeng-071516-044442 ARBEF7 1523-9829 Google Scholar

16. 

Ronneberger O., Fischer P., Brox T., “U-Net: convolutional networks for biomedical image segmentation,” Lect. Notes Comput. Sci., 9351 234 –241 (2015). https://doi.org/10.1007/978-3-319-24574-4 LNCSD9 0302-9743 Google Scholar

17. 

Ciçek O. et al., “3D U-Net: learning dense volumetric segmentation from sparse annotation,” Lect. Notes Comput. Sci., 9901 424 –432 (2016). https://doi.org/10.1007/978-3-319-46723-8 LNCSD9 0302-9743 Google Scholar

18. 

Milletari F., Navab N., Ahmadi S.-A., “V-Net: fully convolutional neural networks for volumetric medical image segmentation,” in IEEE Fourth Int. Conf. on 3D Vision, 565 –571 (2016). https://doi.org/10.1109/3DV.2016.79 Google Scholar

19. 

Kamnitsas K. et al., “Efficient multi-scale 3D CNN with fully connected CRF for accurate brain lesion segmentation,” Med. Image Anal., 36 61 –78 (2017). https://doi.org/10.1016/j.media.2016.10.004 Google Scholar

20. 

Brosch T., Saalbach A., “Foveal fully convolutional nets for multi-organ segmentation,” Proc. SPIE, 10574 105740U (2018). https://doi.org/10.1117/12.2293528 PSISDG 0277-786X Google Scholar

21. 

Andermatt S., Pezold S., Cattin P., “Multi-dimensional gated recurrent units for the segmentation of biomedical 3D-data,” Lect. Notes Comput. Sci., 10008 142 –151 (2016). https://doi.org/10.1007/978-3-319-46976-8 LNCSD9 0302-9743 Google Scholar

22. 

Xue Y. et al., “SegAN: adversarial network with multi-scale L1 loss for medical image segmentation,” Neuroinformatics, 16 383 –392 (2018). https://doi.org/10.1007/s12021-018-9377-x 1539-2791 Google Scholar

23. 

Brain tumor segmentation, “Multimodal brain tumor segmentation (BraTS), Challenge 2017,” (2018) https://www.med.upenn.edu/sbia/brats2017.html July ). 2018). Google Scholar

24. 

Liver tumor segmentation, “Liver tumor segmentation (LiTS), Challenge 2017,” (2018) https://competitions.codalab.org/competitions/15595 July ). 2018). Google Scholar

25. 

Ischemic stroke lesion segmentation, “Ischemic stroke lesion segmentation (ISLES), Challenge 2017,” (2018) http://www.isles-challenge.org/ISLES2017/ July ). 2018). Google Scholar

26. 

Meyer P. et al., “Survey on deep learning for radiotherapy,” Comput. Biol. Med., 98 126 –146 (2018). https://doi.org/10.1016/j.compbiomed.2018.05.018 CBMDAW 0010-4825 Google Scholar

27. 

Thompson R. F. et al., “Artificial intelligence in radiation oncology: a specialty-wide disruptive transformation?,” Radiother. Oncol., (2018). https://doi.org/10.1016/j.radonc.2018.05.030 RAONDT 0167-8140 Google Scholar

28. 

Lustberg T. et al., “Clinical evaluation of atlas and deep learning based automatic contouring for lung cancer,” Radiother. Oncol., 126 (2), 312 –317 (2018). https://doi.org/10.1016/j.radonc.2017.11.012 RAONDT 0167-8140 Google Scholar

29. 

Fritscher K. et al., “Deep neural networks for fast segmentation of 3D medical images,” Lect. Notes Comput. Sci., 9901 158 –165 (2016). https://doi.org/10.1007/978-3-319-46723-8 LNCSD9 0302-9743 Google Scholar

30. 

Ibragimov B., Xing L., “Segmentation of organs-at-risks in head and neck CT images using convolutional neural networks,” Med. Phys., 44 (2), 547 –557 (2017). https://doi.org/10.1002/mp.12045 MPHYA6 0094-2405 Google Scholar

31. 

Deng J. et al., “ImageNet: a large-scale hierarchical image database,” in Proc. IEEE Conf. on Computer Vision and Pattern Recognition, 248 –255 (2009). https://doi.org/10.1109/CVPR.2009.5206848 Google Scholar

32. 

Sharp G. C. et al., “PDDCA: a public domain database for computational anatomy,” (2016) https://www.imagenglab.com/newsite/pddca/ November 2017). Google Scholar

33. 

Prasoon A. et al., “Deep feature learning for knee cartilage segmentation using a triplanar convolutional neural network,” Lect. Notes Comput. Sci., 8150 246 –253 (2013). https://doi.org/10.1007/978-3-642-40763-5 LNCSD9 0302-9743 Google Scholar

34. 

Dieleman S. et al., “Lasagne: first release,” (2015) https://dx.doi.org/10.5281/zenodo.27878 Google Scholar

35. 

Ioffe S., Szegedy C., “Batch normalization: accelerating deep network training by reducing internal covariate shift,” in Proc. of Machine Learning Research, 448 –456 (2015). Google Scholar

36. 

Luo W. et al., “Understanding the effective receptive field in deep convolutional neural networks,” in Proc. of the 30th Int. Conf. on Neural Information Processing Systems, 4898 –4906 (2016). Google Scholar

37. 

Dice L. R., “Measures of the amount of ecologic association between species,” Ecology, 26 (3), 297 –302 (1945). https://doi.org/10.2307/1932409 ECGYAQ 0094-6621 Google Scholar

38. 

Ritter F. et al., “Medical image analysis: a visual approach,” IEEE Pulse, 2 (6), 60 –70 (2011). https://doi.org/10.1109/MPUL.2011.942929 Google Scholar

39. 

Heimann T. et al., “Comparison and evaluation of methods for liver segmentation from CT datasets,” IEEE Trans. Med. Imaging, 28 (8), 1251 –1265 (2009). https://doi.org/10.1109/TMI.2009.2013851 ITMID4 0278-0062 Google Scholar

40. 

Hänsch A. et al., “Parotid gland segmentation with deep learning using clinical vs. curated training data,” Radiother. Oncol., 127 (Suppl. 1), S281 –S282 (2018). https://doi.org/10.1016/S0167-8140(18)30840-5 RAONDT 0167-8140 Google Scholar

41. 

Hänsch A. et al., “Comparison of different deep learning approaches for parotid gland segmentation from CT images,” Proc. SPIE, 10575 1057519 (2018). https://doi.org/10.1117/12.2292962 PSISDG 0277-786X Google Scholar

Biography

Annika Hänsch is a PhD student at Fraunhofer MEVIS in Bremen, Germany. She received her BS degree in mathematics from the University of Oldenburg in 2013 and her MS degree in computational life science from the University of Lübeck in 2016. Her current main research interest is medical image analysis, especially deep learning applied to segmentation problems.

Biographies for the other authors are not available.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Annika Hänsch, Annika Hänsch, Michael Schwier, Michael Schwier, Tobias Gass, Tobias Gass, Tomasz Morgas, Tomasz Morgas, Benjamin Haas, Benjamin Haas, Volker Dicken, Volker Dicken, Hans Meine, Hans Meine, Jan Klein, Jan Klein, Horst K. Hahn, Horst K. Hahn, "Evaluation of deep learning methods for parotid gland segmentation from CT images," Journal of Medical Imaging 6(1), 011005 (1 October 2018). https://doi.org/10.1117/1.JMI.6.1.011005 . Submission: Received: 31 May 2018; Accepted: 31 August 2018
Received: 31 May 2018; Accepted: 31 August 2018; Published: 1 October 2018
JOURNAL ARTICLE
9 PAGES


SHARE
Back to Top