Sm: enhanced localization in Multiple Instance Learning for medical imaging classification

Multiple Instance Learning (MIL) is widely used in medical imaging classification to reduce the labeling effort. While only bag labels are available for training, one typically seeks predictions at both bag and instance levels (classification and localization tasks, respectively). Early MIL methods treated the instances in a bag independently. Recent methods account for global and local dependencies among instances. Although they have yielded excellent results in classification, their performance in terms of localization is comparatively limited. We argue that these models have been designed to target the classification task, while implications at the instance level have not been deeply investigated. Motivated by a simple observation -- that neighboring instances are likely to have the same label -- we propose a novel, principled, and flexible mechanism to model local dependencies. It can be used alone or combined with any mechanism to model global dependencies (e.g., transformers). A thorough empirical validation shows that our module leads to state-of-the-art performance in localization while being competitive or superior in classification. Our code is at https://github.com/Franblueee/SmMIL.

Paper

Similar papers

Peer review

Reviewer uFQt5/10 · confidence 4/52024-07-09

Summary

The paper is well-written and proposes a new modular mechanism that can (selectively) combine local and global interactions in a model. The authors discuss the theory behind their proposed smoothness operator and pooling mechanism and provide detailed proof for their claims. They have also evaluated their model against other methods on three different datasets (1 CT and 2 H&E) with reporting AUC and F1 score as their evaluation metrics. Their method is superior to other methods and the proper ablation study supports their claims.

Strengths

The writing is concise and the contributions are clear. The author proposes a new smoothness mechanism, generally an interesting idea, to integrate local and global information. Besides the quantitative experiments, there are some qualitative experiments such as attention maps that show the model's relative importance on different regions in a slide.

Weaknesses

There are three major concerns with this work: 1. The choice of encoder. The authors have used Resnet18 for RSNA and Camelyon16 and Resnet50 for PANDA. I could not find any justification for why different encoders have been used! Also, with the trend toward foundation models, transformer-based backbones are now of interest to the community. Therefore, for the sake of consistency, the author should use the same encoder for different datasets or report both resent50 and resnet18. And, for the sake of the generality of their work, they should add one encoder such as ViT or Swin as well (w/ ImageNet weights) to support the generality of their claims. 2. The body of research has been founded on local-to-global interactions. There are quite a few standard graph-based methodologies in the literature that the authors need to compare their work against. Currently, there are no representatives from those families of the MIL method in the benchmark. Two examples of such methods are (1) and (2). 3. There is a family of MIL methods in the literature that try to pseudo-label the instances during training, which is essentially equivalent to localization in this work. For instance, two of the most recent such methods are (3) and (4). The author should compare their method against these as they are essentially from the same family. (1) Chen, Richard J., et al. "Whole slide images are 2d point clouds: Context-aware survival prediction using patch-based graph convolutional networks." Medical Image Computing and Computer Assisted Intervention–MICCAI 2021: 24th International Conference, Strasbourg, France, September 27–October 1, 2021, Proceedings, Part VIII 24. Springer International Publishing, 2021. (2) Li, R., Yao, J., Zhu, X., Li, Y., Huang, J. (2018). Graph CNN for Survival Analysis on Whole Slide Pathological Images. In: Frangi, A., Schnabel, J., Davatzikos, C., Alberola-López, C., Fichtinger, G. (eds) Medical Image Computing and Computer Assisted Intervention – MICCAI 2018. MICCAI 2018. Lecture Notes in Computer Science(), vol 11071. Springer, Cham. https://doi.org/10.1007/978-3-030-00934-2_20 (3) Z. Shao et al., "LNPL-MIL: Learning from Noisy Pseudo Labels for Promoting Multiple Instance Learning in Whole Slide Image," 2023 IEEE/CVF International Conference on Computer Vision (ICCV), Paris, France, 2023, pp. 21438-21438, doi: 10.1109/ICCV51070.2023.01965. (4) Ren, Q. et al. (2023). IIB-MIL: Integrated Instance-Level and Bag-Level Multiple Instances Learning with Label Disambiguation for Pathological Image Analysis. In: Greenspan, H., et al. Medical Image Computing and Computer Assisted Intervention – MICCAI 2023. MICCAI 2023. Lecture Notes in Computer Science, vol 14225. Springer, Cham. https://doi.org/10.1007/978-3-031-43987-2_54

Questions

I am curious if the author expects the same performance boost obtained on top of resent, to be obtained on top of a histopathology-specific feature extractor.

Rating

5

Confidence

4

Soundness

3

Presentation

3

Contribution

2

Limitations

It has been justified :)

Reviewer uFQt2024-08-08

Questions regarding the new experiments

I would like to thank the authors for extending their experiments and including new insights. I would appreciate the authors if they could provide more details regarding the following questions. It has been mentioned in line 45 that "This is a flexible module that can be used alone on top of classical MIL approaches, or in combination with transformers to also account for global dependencies." If I have understood it right, the main contribution is to have sm as a local-global module that helps to improve performance by capturing both local and global. Judging only from Table 1 and Resnet50 with ImageNet weights, how do authors justify the models without global interactions that are superior to their proposed smTAP with global interactions? If I am not missing anything, this is a bit contradictory to the claims of the paper. Also, comparing Resnet50+BT and Resnet50 only for the Calmelyon16 dataset, it seems that their proposed model benefits heavily from pre-trained features (improvement comes mostly from the encoder rather than the aggregation proposed). I am wondering if the authors can justify this observation. Does this come from any biological underlying meaning or can the Sm method be intuitively linked to anything?

Area Chair vXTG2024-08-08

Please read the rebuttal to check if the authors addressed your concerns

Dear Reviewer uFQt, Can you have a look at the rebuttal and see if your concerns have been addressed? Best regards Your AC.

Authorsrebuttal2024-08-10

On the new questions

We would like to thank Reviewer uFQt for their positive comments on the new experiments and the new insights we have provided. We address the new questions below. > **Q. It has been mentioned in line 45 that "This is a flexible module that can be used alone on top of classical MIL approaches, or in combination with transformers to also account for global dependencies." If I have understood it right, the main contribution is to have sm as a local-global module that helps to improve performance by capturing both local and global. Judging only from Table 1 and Resnet50 with ImageNet weights, how do authors justify the models without global interactions that are superior to their proposed smTAP with global interactions? If I am not missing anything, this is a bit contradictory to the claims of the paper.** As we indicate in the paper, the proposed Sm operator introduces local interactions in a principled manner (through Dirichlet energy minimization with the presented theoretical guarantees). This operator can be naturally combined with self-attention layers (transformer layers) to account for both local and global interactions. As the reviewer points out, in Table 1 of the rebuttal PDF, when looking at ResNet50 with Imagenet weights, the models without global interactions (i.e., no transformer layers) do better at the instance level (localization task), with the proposed SmAP achieving the highest rank (see the table in the rebuttal to Reviewer uFQt). However, when looking at the classification task (bag level, Table 2), the models with global interactions outperform to those without global interactions, with the proposed SmTAP achieving the highest rank (see the table in the rebuttal to Reviewer uFQt). These results suggest that models without global interactions (SmAP) perform better at the instance level (localization task, Table 1), while models with global interactions perform better at the bag level (classification task, Table 2). This behaviour was also pointed out by Rev. YAE1. As we argue in the rebuttal to Reviewer YAE1, this is explained by the fact that self-attention layers of transformers obtain a richer bag representation (which leads to better bag level results). However, the signal that identifies each token in the sequence becomes weaker after each self-attention layer (the tokens become more similar). This leads to worse instance level results. Note that this phenomenon is known as over-smoothing and has been observed in (1), (2), (3) (see references below). We believe that the results in Table 1 are not contradictory with the claims in the paper, but rather reinforce them. The proposed Sm improves instance level results (localization) and remains competitive at the bag level (classification). This holds true whether it is combined with transformer layers (SmTAP) or not (SmAP). > **Q. Also, comparing Resnet50+BT and Resnet50 only for the Calmelyon16 dataset, it seems that their proposed model benefits heavily from pre-trained features (improvement comes mostly from the encoder rather than the aggregation proposed). I am wondering if the authors can justify this observation. Does this come from any biological underlying meaning or does the Sm method can be intuitively linked to anything?** As the reviewer points out, the proposed Sm benefits from SSL pre-trained features. Our intuition is that if an SSL pre-trained encoder is used, instances that are similar from a biological point of view will tend to cluster together in the feature space. In particular, instances with the same label will be clustered together, and far away from instances with a different label. Since neighboring instances are naturally (a priori) expected to have similar biological properties (e.g., the same label), the smoothness introduced by the proposed Sm will be particularly beneficial. This interesting observation is related to that in CAMIL [13], where it was shown that a SSL pre-trained encoder is crucial for the optimal performance of their neighbor-constrained attention mask. ---- (1) Shi, Han, et al. "Revisiting Over-smoothing in BERT from the Perspective of Graph." International Conference on Learning Representations (ICLR). 2022. (2) Wang, Peihao, et al. "Anti-Oversmoothing in Deep Vision Transformers via the Fourier Domain Analysis: From Theory to Practice." International Conference on Learning Representations (ICLR). 2022. (3) Dong, Yihe, Jean-Baptiste Cordonnier, and Andreas Loukas. "Attention is not all you need: Pure attention loses rank doubly exponentially with depth." International Conference on Machine Learning (ICML). PMLR, 2021.

Authorsrebuttal2024-08-12

Concern whether our reply was satisfactory

In case it is necessary, we want to summarize our response to the new questions by the reviewer. > **Q. It has been mentioned in line 45 that "This is a flexible module that can be used alone on top of classical MIL approaches, or in combination with transformers to also account for global dependencies." If I have understood it right, the main contribution is to have sm as a local-global module that helps to improve performance by capturing both local and global. Judging only from Table 1 and Resnet50 with ImageNet weights, how do authors justify the models without global interactions that are superior to their proposed smTAP with global interactions? If I am not missing anything, this is a bit contradictory to the claims of the paper.** Instance level results worsen with transformers. This is not a consequence of the proposed Sm, but of the oversmoothing phenomenon observed in transformers. This has been studied in the works referenced in our previous response. The proposed method without a transformer outperforms methods that do not a transformer; and the proposed method with a transformer outperforms methods that use a transformer. > **Q. Also, comparing Resnet50+BT and Resnet50 only for the Calmelyon16 dataset, it seems that their proposed model benefits heavily from pre-trained features (improvement comes mostly from the encoder rather than the aggregation proposed). I am wondering if the authors can justify this observation. Does this come from any biological underlying meaning or does the Sm method can be intuitively linked to anything?** Note that every method has access to the same features, and every method benefits from using them (see bag level results in Table 2). That our method benefits at the instance level from using, we interpret it as a sign that local dependencies indeed play an important role in MIL localization. Please let us know if you have any questions so that we can respond in a timely manner. We look forward to your response, Authors

Reviewer uFQt2024-08-12

I appreciate the authors' effort in providing further details and explanations for my questions. The response provided addressed my concern partly, yet I agree with Reviewer YAE1 regarding the improvement being marginal. However, given the efforts in providing newer benchmarks and also considering that the content needs to be updated with consistent backbones for different datasets (as provided in the previous reply). Therefore, I would raise my rating to Borderline Accept).

Reviewer jUey6/10 · confidence 4/52024-07-13

Summary

The authors propose a technique to improve the localization capabilities of the current MIL, especially for CAD models that perform CT and WSI analysis. The method is based on seeing the attention attributed to each patch as a graph and minimizing its Dirichlet energy, promoting smooth transitions on the attention values of neighbor patches, which is in accordance with the assumption that neighboring patches share the same labels.

Strengths

The work is well-developed and well written, with a methodology explained through concise equations and theorems that guide the reader through the development of the idea. The results validate what is proposed, but minor adjustments and additional information will make the work even more valuable.

Weaknesses

The method introduces the smooth operator on top of the existing ABMIL model. The proposed approach improves the performance of ABMIL on the CAMELYON dataset, both on localization and classification. Nevertheless, when analyzing the attention maps provided in the appendix, the differences between ABMIL and smAP are minimal. Although we have a numerical evaluation of the Dirichlet energy, we cannot tell if the drop in energy is significant only by these numbers, so this visual empiric evaluation is important for a better insight regarding the developed work. Observing the appendix figures, ABMIL already performs well in matching its attention to the tiles annotations. I suggest the authors apply the smooth operator on top of other models that did not perform well on CAMELYON for this task, such as CLAM, CAMIL, and SETMIL, to evaluate how minimizing the Dirichlet energy of the attention improves the localization task quantitatively and qualitatively. I suggest the authors perform an ablation study to evaluate how minimizing the Dirichlet energy compares to other smoothing strategies. For instance, summing the variance of the attentions to the model’s loss may provide similar results as minimizing the Dirichlet energy. The analysis of the alpha parameter is great, but again, a qualitative approach is valuable for this case. Seeing how the attention map changes as the parameter alpha increases is highly valuable for the reader. Lines 127 to 131: trying to debate equations not present in the paper is not good for the reader. I suggest writing the debated equations in the article or explaining the arguments without mentioning them.

Questions

The formulation seems to assume bags of equal size, which is not the case in MIL (in general). The proposed method (loss) may have the drawback of over-smoothing the transition between tiles, making the localization regions blurry. Is this true? Could this be overcome?

Rating

6

Confidence

4

Soundness

3

Presentation

3

Contribution

3

Limitations

yes

Area Chair vXTG2024-08-08

Please read the rebuttal to check if the authors addressed your concerns

Dear Reviewer jUey, Can you have a look at the rebuttal and see if your concerns have been addressed? Best regards Your AC.

Area Chair vXTG2024-08-08

Please read the rebuttal to check if the authors addressed your concerns

Dear Reviewer YAE1, Can you have a look at the rebuttal and see if your concerns have been addressed? Best regards Your AC.

Reviewer YAE12024-08-09

reviewer feedback

Thanks for the authors' effort to present additional results and make clarifications. I raised the score despite the concern that this still seems to be a marginal improvement.

Program Chairsdecision2024-09-25

Decision

Accept (poster)

© 2026 NYSGPT2525 LLC