One model to use them all: training a segmentation model with complementary datasets

Purpose Understanding surgical scenes is crucial for computer-assisted surgery systems to provide intelligent assistance functionality. One way of achieving this is via scene segmentation using machine learning (ML). However, such ML models require large amounts of annotated training data, containing examples of all relevant object classes, which are rarely available. In this work, we propose a method to combine multiple partially annotated datasets, providing complementary annotations, into one model, enabling better scene segmentation and the use of multiple readily available datasets. Methods Our method aims to combine available data with complementary labels by leveraging mutual exclusive properties to maximize information. Specifically, we propose to use positive annotations of other classes as negative samples and to exclude background pixels of these binary annotations, as we cannot tell if a positive prediction by the model is correct. Results We evaluate our method by training a DeepLabV3 model on the publicly available Dresden Surgical Anatomy Dataset, which provides multiple subsets of binary segmented anatomical structures. Our approach successfully combines 6 classes into one model, significantly increasing the overall Dice Score by 4.4% compared to an ensemble of models trained on the classes individually. By including information on multiple classes, we were able to reduce the confusion between classes, e.g. a 24% drop for stomach and colon. Conclusion By leveraging multiple datasets and applying mutual exclusion constraints, we developed a method that improves surgical scene segmentation performance without the need for fully annotated datasets. Our results demonstrate the feasibility of training a model on multiple complementary datasets. This paves the way for future work further alleviating the need for one specialized large, fully segmented dataset but instead the use of already existing datasets.


Introduction
The understanding of the visible surgical scene is key for computer-assisted surgery (CAS) systems to understand the current situation and provide adapted assistance functions.One approach of achieving this, is via the use of full scene semantic segmentation models, which are able to classify every visible part of the scene.These models provide the basis for recognizing the current situation or actions, and enable useful assistance functions.
In recent work, progress has been made in improving surgical semantic scene segmentation by the use of temporal context [1], stereo vision [2], simulated data [3], or weak label annotations [4].Further steps towards solving this task in the surgical setting have been taken by the Robotic Scene Segmentation Challenge [5], and the HeiSurF Challenge [6].Both EndoVis1 sub-challenges provided annotations of 11 and 21 classes, respectively, including surgical tools and human anatomy, and challenged participants to semantically segment all of them.
Nevertheless, a major bottleneck for clinical translation of surgical data science (SDS) applications remains the availability of such datasets [7], due to the high amount of time required for experts, to create such segmentation annotations.This is even more challenging for full scene semantic segmentation, which requires a pixel-wise annotation of multiple classes for the complete frame.
This issue was sidestepped in the recently published Dresden Surgical Anatomy Dataset (DSAD) [8] by simply providing binary segmentations.The dataset contains 11 classes of anatomical structures, split into multiple subsets providing one class each.This way, the authors were able to publish over 13000 expert-approved semantic segmentations.
At present there are a handful of datasets providing different annotations relevant to the field of SDS [7].These datasets can differ in the granularity of classes and the classes annotated in general, depending on the protocol used during annotation.Recent works have presented different methods that rely only on partial annotations for training a CT segmentation model, improving the usage of existing knowledge [9][10][11][12].The used approaches range from dataset specific backbones and pseudo label generation [12], over merging unlabeled classes with the background class and adding a mutual exclusion constraint [9], to simply masking unlabeled classes during loss calculation [10].
We assume that learning to segment multiple classes causes a single model to develop a better understanding of the scene, leading to better segmentation performance.We therefore propose the usage of information from mutual-exclusivity, which can be applied on top of any state-of-the-art approach.In this work we apply it in addition to the masking during loss calcation [10].We demonstrate the feasibility of training a model for multi-class organ segmentation of laparoscopic surgery images on multiple complementary dataset, thereby overcoming the data bottleneck challenge in SDS.The code and models are publicly available on https://gitlab.com/nct tso public/dsad-segmentation/.

Methods
In this section, we introduce our proposed method and define an upper and lower baseline to compare the performance of the proposed method to.A visualization of the architectures is shown in Figure 1.Further, we provide a definition of how the average dice scores are calculated in this work.

Baselines
A naive approach to combine multiple datasets is to train one model per dataset, and subsequently building an ensemble to merge their predictions into one final multi-class prediction.In case of one binary segmenting model per class, this kind of ensemble prediction can be achieved by applying the argmax over the different sigmoid outputs per pixel.This is followed by a threshold to determine if the most likely class is predicted positive, otherwise the background class is assigned.In this work this approach will be used as the lower baseline and is referred to as ensemble (EN).
Alternatively a single model can be trained on a single fully labeled dataset.As this requires a large fully labeled dataset, which is hard to obtain but provides the maximum amount of information, this approach will be used as the upper baseline in this work and is referred to as fully supervised (FS).

Implication based labeling
To overcome the ambiguity inherent to merging multiple models and lost information between classes, while sidestepping the need for a single large fully annotated dataset, we propose to combine the classes provided by multiple datasets into one model.This reduces computational costs during inference as only one output from one model is required and lowers the requirements we pose against the dataset, as not all classes need to be included in every dataset.In addition, the model could benefit from the shared knowledge among the different classes.
The model takes single images of the datasets as input and outputs a class probability for every pixel and every class.For every pixel, the final class is selected based on the highest predicted class probability.The outputs are normalized by a sigmoid and therefore each value is independent of the outputs of other classes.On one hand, this setup does not prevent the model from predicting multiple classes per pixel, instead this has to be enforced by the loss during training.On the other hand, this setup allows us to access the predicted probability of every class separately, which is required for dealing with the problem of training from incomplete knowledge in the datasets, as each dataset provides only information on its contained class and not the classes introduced by other datasets.
Assuming every pixel in the target semantic scene segmentation problem is exactly part of one class (mutually exclusivity of classes), this can be used to maximize the information provided by a partially annotated dataset by applying the following rules: 1.A positive annotation of one class implies a negative annotation of all other classes.2. A negative annotation of a class (in a binary annotation the background) provides no information if this region contains other classes.As no implication to other classes can be made, the annotation stays unknown.
The application of these rules is visualized in the blue box of Figure 1.
In this work we use these implication rules to obtain additional negative samples from datasets not containing the examined class.Cases in which the annotation stays unknown are excluded from the loss calculation by masking, as no decision on the correctness of the prediction can be made.This method is called implicit labeling (IL) in the following.
The loss per class c can be formulated as:

Metrics
In this work the model performance is evaluated using the dice score [13].For the average dice score per class, the dice is calculated per image and averaged subsequently over all images.If a class does not occur in either the target nor the prediction, the F1 is not defined, we therefore set the score to one in those cases as the model performed as expected.The average dice score over all classes is calculated by averaging all average dice scores per class.
Statistical significance is calculated using a two-sided Wilcoxon signed-rank test on the image-wise and class-wise dice scores of our approach against the lower baseline.The significance is calculated per class.The significance of the mean dice score is calculated by evaluationg all image-wise and class-wise dice scores over all classes.

Dataset
We evaluate our proposed approach using the publicly available DSAD dataset [8].This dataset consists of 13,195 laparoscopic images split into 11 subsets with a minimum of 1,000 frames from a minimum of 20 surgeries each.In each subset, binary segmentations for one of the 11 classes (abdominal wall, colon, inferior mesenteric artery, intestinal veins, liver, pancreas, small intestine, spleen, stomach, ureter, vesicular glands) are provided.For the stomach subset, additional masks are available, annotating six of the remaining ten classes (abdominal wall, colon, liver, pancreas, small intestine, spleen) visible in this subset, resulting in one multi-class subset.This work uses the proposed split [14] into training, validation, and test set.To be able to compare to the multi-class subset, only the binary subsets of the contained classes are used, with the exception of the spleen.The spleen class was excluded in this work due to the lack of positive examples in the validation and test split in the multi-class subset.To examine the ability of IL and EN to join datasets, they were trained by splitting the multi-class subset into six binary sets, one for each class.The methods trained on these subsets are comparable to the FS approach trained on the combined annotations of the multi-class subset.As the binary subsets have no overlapping classes, we interpret them as separate datasets in this work, which we aim to join into one model.All classes fulfill the required mutual exclusiveness.

Trials
To validate our approach, we conducted five trials to compare our proposed implicit labeling approach to the ensemble approach and the fully supervised approach.In all experiments, the DeepLabV3 architecture with a ResNet50 backbone [15] is used.Models are initialized using the default PyTorch pretraining on COCO [16].All models were trained using PyTorch [17] v1.12 on Ubuntu 18.04 and Nvidia Tesla V100 GPUs, evaluation was done on Ubuntu 20.04 with Nvidia RTX A5000 GPUs.The images were down scaled to a size of 640x512 pixels for memory and time reasons.For all trials an initial learning rate of 3 × 10 −4 was used with a scheduler reducing the learn rate by a gamma of 0.9 every 10 epochs.Further all model were trained for 100 epochs using an AdamW optimizer with weight decay of 0.1 and an cross-entropy loss with a positive weight factor for balancing positive and negative pixels.The best model per trial was selected according to the dice score on the validation set.
The ensemble approach evaluates an ensemble of six models, each trained on one of the classes, serving as the lower baseline.The prediction of the ensemble was obtained by selecting the class with the highest value over all models per pixel.The positive weight for the background was set to 1 and the remaining classes to the negative to positive pixel ratio per class.This approach was used in two trials, once trained on the binary subset, and once on the binary sets extracted from the multi-class subset.The fully supervised approach examines the fully supervised method, serving as the upper baseline.For this a single DeepLabV3 was trained on the multi-class subset.The positive weight was calculated by the share of positive pixels of the class, negating the share, and normalizing all negated shares with the softmax function.As this approach requires a fully annotated dataset, this method is only trained once on the multi-class subset.In the implicit labeling approach, our proposed method was used to train a single DeepLabV3 on the six binary subsets, which represent the same classes that are available in the multi-class subset.The positive weight was set to the negative to positive pixel ratio, while including the positive pixels of other classes as negatives, as described before.The output of the loss function was masked before aggregation to ignore pixels as required by our approach.This approach was used in two trials, once trained on the binary subset, and once on the binary sets extracted from the multi-class subset.

Results
The trials were evaluated twice on the up-to-now unseen test split of the dataset, once using the binary subset for each class and once using the multi-class subset.As the multi-class subset is based on the stomach subset, the frames and, therefore, the results for this class are identical.The results on the binary testset are shown in Table 1.The overall best performance is reached by our proposed IL approach on the binary trainset with an average dice score of 72% over all classes outperforming EN by 4% On the binary trainset all scores of the IL approach are either significantly higher or not significantly lower than the EN approach.For the trials trained on the multi-class subset the upper baseline, FS, reaches the highest score of 46%.For the highest score per class distributes among the FS and IL approach.On the multi-class trainset all scores of the IL approach are significantly above the EN approach, except for the stomach.The performance of all methods is significantly lower for Pancreas and Liver on the multi-class trainset compared to the binary trainset, while IL maintains the highest score.The results on the multi-class testset are shown in Table 2.The overall best performance is reached by the FS approach with an average dice score of 78%.On the binary trainset IL outperforms EN by 7%.Especially in the classes Colon, Pancreas, Small Intestine and Stomach IL performs significantly better than EN, while EN is better for the Liver.On the multi-class testset FS has the hightest score on all classes.The second highest class is distributed among IL and EN.The average performance of IL outperforms EN by 4%.
Figure 2 shows the pixelwise confusion of the classes for the trials tested on the multi-class testset.For trials trained on the multi-class trainset IL shows lower confusion with the background for the pancreas and the small intestine than EN.This means pixels of both classes are less often missed to detect.Compared to FS, IL less often confuses the small intestine with the colon.The liver is more often confused with the abdominal wall by IL compared to FS and EN.For the binary trials, IL lowers the confusion of the stomach with the colon and small intestine.The confusion of the liver with the abdominal wall and background slightly increases.
Figure 3 shows four examples of segmentation results.In the rows (a) and (b), FS and IL are able to segment a single structure with one class while the EN approach is mixing multiple classes.In rows (c) andd (d) all three approaches produce patched results per organ and detect classes not present.
Table 3 shows the results of out ablation studies trained and evaluated on the binary subset.The first row shows the dice score of the separate models in the ensemble before applying the argmax.All classes except the abdominal wall reach higher scores compared to the merged ensemble.IL is outperformed on all classes except the stomach.The second row shows the results of our IL approach if no additional negative samples are inferred due to mutual exclusion and only the loss masking is applied.While this model is better than our proposed IL approach for the classes abdominal wall and liver the remaining classes and the average performance are below IL.Table 3 Dice score of each class calculated on the testsets of the binary subsets, µ being the mean over all classes.The first row shows the performance of each separate models in the ensemble before applying the argmax.The second row shows a variant of the IL approach if no additional negative samples are inferred.The classes Colon, Liver, Pancreas, Small Intestine and Stomach are abbreviated by their first three letters.The trials were evaluated with respect to the inference time by inferring 1000 random inputs on an Nvidia RTX A5000.The time needed from loading the image to the GPU to downloading the prediction back to the CPU was averaged.All models were already prepared on the GPU.The ensemble approach required 136 ± 5.8 ms per frame, the fully supervised and implicit labeling approaches both required 23 ± 0.1 ms per frame.Resulting in 7fps and 43fps, respectively, and therefore a 6x speed up.Memory consumption on the GPU stayed below 3.2 GB for the ensemble and 2.3 GB for the fully supervised and implicit labeling approaches.

Discussion
In summary, this work demonstrates that the implicit labeling approach is able to leverage multiple complementary datasets into one model.We find that models benefit from a better scene understanding through more learned classes and reach better performance.Further, we point out limitations due to the lack of data diversity we came across in the different data subsets.
As Table 3 shows the performance on EN merged by argmax drops compared to the performance of its separate models.This again proves the importance of being able to train a single model being able to understand and segment all classes.Best results are reached by fully annotated datasets, as shown by FS in Tables 1 and 2. But as these are rarely available, our proposed method is a valid alternative in cases where multiple datasets need to be merged.
Tables 1 and 2 clearly show, that IL significantly outperforms EN in all four combinations of test and trainset of the binary and multi-class subsets.On the binary train and binary test combination our proposed approach is either within 1% difference on the dice score or significantly better than the lower baseline.This supports our initial assumption that the models benefit from better scene understanding due to complementary information.As shown in Figure 3 (a) and (b) the knowledge of multiple classes, as in the FS and IL approach, helps the model to select one class.This missing information in the EN approach leads to segmentations being patched together of multiple classes.IL is able to solve this issue without the need for fully segmented annotations.This is demonstrated further in Figure 2 by comparing the class confusion of IL to EN.For all classes except the liver the confusion drops.Especially the confusion of colon, small intestine and stomach on the binary trainset, where those classes do look similar in many cases, are reduced by the knowledge of all classes.Interestingly the fully supervised model is not able to prevent confusion of the small intestine and the colon, which might be caused by rare occurrence of those classes due to the stomach centric multi-class subset.Interestingly, the large increase of negative samples in the IL approach compared to the FS approach does not lead to the model to prefer the background class rather the opposite is the case.This might be due to the higher weighting factors of positive pixels used in the loss.This proves our assumption that models benefit from better scene understanding through more learned classes.
As shown in Table 1 and 2, the models perform worse when tested on the subset they are not trained on.Especially on the Pancreas and Small Intestine the FS approach drops in performance.This is explained by the fact that the multi-class subset is based on the stomach.Therefore, all frames contain this structure, limiting the possible viewing angle and distances of other annotated classes.IL is able to maintain the highest cross-domain performance in both directions.This cross-domain capabilities can further be seen in the first row of Table 2 where IL is trained on the multi-class trainset and tested on the binary subset without any performance drop on the average dice score compared to the third row which shows the performance of IL tested on the multi-class testset.EN, in rows two and four, is dropping by 5% due to the different appearance in the subsets.
The second row of Table 3 shows the importance of the implicit labels, as the models performance drops without them.This ablation is comparable to the approach presented by Ulrich et al. [10], which further proves the benefit of our method, as by simply adding the implicit knowledge we outperform related work.
When looking at the situation given by the dataset we are limited to a small fully annotated dataset and have more binary data on hand.While it is possible to train the FS approach on the small dataset and achieve a good performance when testing on similar data, as shown in Table 2, we do see the performance drop when the model is applied to different data, as shown in Table 1.The classes Pancreas and Small Intestine drop by 51% and 79%, respectively, resulting in dice scores below 20%.Considering that the binary dataset is focused around the organ it is tested for, these misclasifications can be considered critical, as the central organ is not recognized.When training our IL approach, on the other hand, not only more data is available, but also the performance drop in changing data is less severe.While we do see dropping performance on classes, if they are no longer in the focus of the frame, no class drops below a dice score of 25%, as shown in Table 2.This proves our IL approach is more applicable in a realistic setting with limited data availability, producing more robust results without the need for complicated preprocessing or much higher computational costs.
Finally, as the inference time results show, the implicit labeling approach combines fast inference with the ability to learn on not fully annotated datasets.The ensemble is 6 times slower, as it needs to infer 6 models instead of one, setting a limit to the scalability.The runtime of the ensemble increases linearly with the number of classes, while our implicit labeling approach allows adding more segmentation data with a negligible effect on runtime.

Conclusion
In this paper, we presented, to the best of our knowledge, the first approach for laparoscopic organ segmentation that combines multiple datasets into one model.We accomplished this by applying a combination of masking during loss and mutual exclusion constraints.We were able to show that segmentation models benefit from a better understanding of the scene in the sense of knowing more classes, improving the overall dice score, reducing confusion between classes and improving generalization to changes in the appearance of the classes.Further, we were able to show that we do not require all classes to be annotated in a single dataset but rather can combine complementary ones.The resulting model was able to achieve real-time capable inference speeds of 43fps on an Nvidia RTX A5000 GPU, and additional classes can be added without linearly increasing the runtime.
Even though the already good results, we see potential to further improve the approach.For example, the model might include weak labels during training in the form of binary presence of classes.Also, the effects of the use of datasets with overlapping classes should be investigated, as well as the possibility to combine noncomplementary datasets, like different levels of detail in the classes.For example, some datasets only segment instruments in one class, and others split them up into different types.Further, the use of active learning to selectively annotate missing data might be promising, as it matches very well with the structure of complementary datasets.Overall, we see great potential in applying our method in multiple settings of semantic segmentation of surgical data science.

Declarations
Competing interests The authors declare no conflict of interest.Ethical approval For this type of study, formal consent is not required.Informed consent This article contains patient data from publicly available datasets.Funding This work is funded by the German Federal Ministry of Health (BMG), on the basis of a decision by the German Bundestag, within the "Surgomics" project (Grant Number BMG 2520DAT82), the German Cancer Research Center (CoBot 2.0), the German Research Foundation (Deutsche Forschungsgemeinschaft, DFG) as part of Germany's Excellence Strategy (EXC 2050/1, Project ID 390696704) within the Cluster of Excellence "Centre for Tactile Internet with Human-in-the-Loop" (CeTI) of the Dresden University of Technology and by the European Union through NEAR-

Fig. 1 A
Fig. 1 A flowchart describing the flow of annotation data through the three architectures used in this work.The proposed implicit label generation is shown in the blue box.
where B is the number of images, b ∈ [1, B], in the batch, and P is the number of Pixels, p ∈ P , where λ = 1, ŷ(c) b,p is the prediction of class c for pixel p of image b and y (c) b,p the respective ground truth.A pixel p is included in the loss calculation if either the image b is annotated for class c or if y (ĉ) b,p = 1 for any other class ĉ ∈ C \c, otherwise the pixel is excluded by λ b,p,c being set to zero.In the second case the ground truth for class c is false, y (c) b,p = 0, due to the mutually exclusivity of classes.

Fig. 2
Fig. 2 Pixelwise confusion matrix per approach and trainset tested on the multi-class testset.

Fig. 3
Fig. 3 Segmentation results for 3 example images.The columns show the original frame, the ground truth, and the predicted segmentations of the fully supervised (FS), implicit labeling (IL), and ensemble (EN) approaches.Rows (a) and (b) show examples of good IL performance, (c) and (d) show more difficult cases.

Table 1
Dice score of each class calculated on the testsets of the binary subsets, µ being the mean over all classes.Every row represents one of the trials, fully supervised (FS), ensemble (EN), and implicit labeling (IL), trained on the binary and multi-class sets.The highest value per column and trainset is highlighted in bold.The classes Abdominal Wall, Colon, Liver, Pancreas, Small Intestine and Stomach are abreviated by their first three letters.Cells are highlited in green if the difference between IL and EN is strongly significant (p<.01) or yellow if significant (p<.05).

Table 2
Dice score of each class calculated on the testset of the multi-class subset, µ being the mean over all classes.Every row represents one of the trials, fully supervised (FS), ensemble (EN), and implicit labeling (IL), trained on the binary and multi-class sets.The highest value per column and trainset is highlighted in bold.The classes Background, Abdominal Wall, Colon, Liver, Pancreas, Small Intestine and Stomach are abreviated by their first three letters.Cells are highlited in green if the difference between IL and EN is strongly significant (p<.01) or yellow if significant (p<.05).