Patch Based Convolutional Neural Network for Whole Slide Tissue Image Classification

From statwiki
Revision as of 21:14, 15 November 2021 by Alivochk (talk | contribs) (Previous Work)
Jump to: navigation, search

Presented by

Cassandra Wong, Anastasiia Livochka, Maryam Yalsavar, David Evans

Introduction

Despite the fact that CNN are well-known for their success in image classification, it is computationally impossible to use them for cancer classification. This problem is due to high-resolution images that cancer classification is dealing with. As a result, this paper argues that using a patch level CNN can outperform an image level based one and considers two main challenges in patch level classification – aggregation of patch-level classification results and existence of non-discriminative patches. For dealing with these challenges, training a decision fusion model and an Expectation-Maximization (EM) based method for locating the discriminative patches are suggested respectively. At the end the authors proved their claims and findings by testing their model to the classification of glioma and non-small-cell lung carcinoma cases.

Previous Work

The proposed patch-level CNN and training a decision fusion model as a two-level model was made apparent by the various breakthroughs and results noted below:

  • Majority of Whole Slide Tissue Images classification methods fixate on classifying or obtaining features on patches [1, 2, 3]. These methods excel when an abundance of patch labels are provided [1, 2], allowing patch-level supervised classifiers to learn the assortment of cancer subtypes. However, labeling patches requires specialized annotators; an excessive task at a large scale.
  • Multiple Instance Learning (MIL) based classification [4, 5] utilizes unlabeled patches to predict a label of a bag. For a binary classification problem, the main assumption (Standard Multi-Instance assumption, SMI) states that a bag is positive if and only if there exists at least one positive instance in the bag. Some authors combine MIL with Neural Networks[43, 57] and model SMI by max-pooling. This approach is inefficient due to only one instance with a maximum score (because of max-pooling) being trained in one training iteration on the entire bag.
  • Other works sometimes apply average pooling (voting). However, it has been shown that many decision fusion models can outperform simple voting[42, 45]. The choice of the decision fusion function would depend heavily on the domain.

EM-based Method with CNN

Figure 2. Top: A CNN is trained on patches and EM-based method iteratively eliminates non-discriminative patches. Bottom: An image-level decision fusion model is trained on histograms of patch-level predictions to predict the image-level label

The high-resolution image is modelled as a bag, and patches extracted from it are instances that form a specific bag. The ground truth labels are provided for the bag only, so we model the labels of an instance (discriminative or not) as a hidden binary variable. Hidden binary variables are estimated by the Expectation-Maximization algorithm. A summary of the proposed approach can be found in Fig.2. Please note that this approach will work for any discriminative model.

In this paper [math]X = \{X_1, \dots, X_N\}[/math] denotes dataset containing [math]N[/math] bags. A bag [math]X_i= \{X_{i,1}, X_{i,2}, \dots, X_{i, N_i}\}[/math] consists of [math]N_i[/math] pathes (instances) and [math]X_{i,j} = \lt x_{i,j}, y_j\gt [/math] denotes j-th instance and it’s label in i-th bag. We assume bags are i.i.d. (independent identically distributed), [math]X[/math] and associated hidden labels [math]H[/math] are generated by the following model: $$P(X, H) = \prod_{i = 1}^N P(X_{i,1}, \dots , X_{i,N_i}| H_i)P(H_i) \quad \quad \quad \quad (1) $$ [math]Hi = {H_{i, 1}, \dots, H_{i, Ni}}[/math] denotes the set of hidden variables for instances in the bag [math]X_i[/math] and [math]H_{i, j}[/math] indicates whether the patch [math]X_{i,j}[/math] is discriminative for [math]y_i[/math] (it is discriminative if estimated label of the instance coincides with the label of the whole bag). Authors assume that [math]X_{i, j}[/math] is independent from hidden labels of all other instances in the i-th bag, therefore [math](1)[/math] can be simplified as: $$P(X, H) = \prod_{i = 1}^{N} \prod_{j=1}^{N_i} P(X_{i, j}| H_{i, j})P(H_{i, j}) \quad \quad (2)$$ Authors propose to estimate the hidden labels of the individual patches [math]H[/math] by maximizing the data likelihood [math]P(X)[/math] using Expectation Maximization. In one iteration of EM ​​we alternate between performing E step (Expectation) where we estimate hidden variables [math]H_{i, j}[/math] and M step (Maximization) where we update the parameters of the model [math](2)[/math] such that data likelihood [math]P(X)[/math] is maximized. Let's denote [math]D[/math] the set of discriminative instances. We start by assuming all instances are in [math]D[/math] (all [math]H_{i, j}=1[/math]).

Discriminative Patch Selection

References

[1] A. Cruz-Roa, A. Basavanhally, F. Gonzalez, H. Gilmore, M. Feldman, S. Ganesan, N. Shih, J. Tomaszewski, and A. Madabhushi. Automatic detection of invasive ductal carcinoma in whole slide images with convolutional neural networks. In Medical Imaging, 2014. 2, 3
[2] H. S. Mousavi, V. Monga, G. Rao, and A. U. Rao. Automated discrimination of lower and higher grade gliomas based on histopathological image analysis. JPI, 2015. 2, 6
[3] Y. Xu, Z. Jia, Y. Ai, F. Zhang, M. Lai, E. I. Chang, et al. Deep convolutional activation features for large scale brain tumor histopathology image classification and segmentation. In ICASSP, 2015. 2, 5, 6
[4] E. Cosatto, P.-F. Laquerre, C. Malon, H.-P. Graf, A. Saito, T. Kiyuna, A. Marugame, and K. Kamijo. Automated gastric cancer diagnosis on h&e-stained sections; ltraining a classifier on a large scale with multiple instance machine learning. In Medical Imaging, 2013. 2
[5] Y. Xu, T. Mo, Q. Feng, P. Zhong, M. Lai, E. I. Chang, et al. Deep learning of feature representation with multiple instance learning for medical image analysis. In ICASSP, 2014. 2
[6]