Quantile Representation for Indirect Immunofluorescence Image Classification

02/06/2014 ∙ by David M. J. Tax, et al. ∙ 0

In the diagnosis of autoimmune diseases, an important task is to classify images of slides containing several HEp-2 cells. All cells from one slide share the same label, and by classifying cells from one slide independently, some information on the global image quality and intensity is lost. Considering one whole slide as a collection (a bag) of feature vectors, however, poses the problem of how to handle this bag. A simple, and surprisingly effective, approach is to summarize the bag of feature vectors by a few quantile values per feature. This characterizes the full distribution of all instances, thereby assuming that all instances in a bag are informative. This representation is particularly useful when each bag contains many feature vectors, which is the case in the classification of the immunofluorescence images. Experiments on the classification of indirect immunofluorescence images show the usefulness of this approach.

READ FULL TEXT VIEW PDF
POST COMMENT

Comments

There are no comments yet.

Authors

page 1

page 2

page 3

page 4

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

Anti-nuclear antibodies (ANAs) are autoantibodies that bind to cell nuclei. In healthy individuals, the immune system produces antibodies against foreign proteins or antigens, but not against human proteins (autoantigens). However, in individuals suffering from autoimmune disorders, antibodies for autoantigens are produced as well. Therefore, ANAs are important for diagnostic purposes and different types of ANAs are indicative of different illnesses [9].

A common test for detecting ANAs is indirect immunofluorescence (IIF) with slides of HEp-2 cells. ANAs bind to the cell nuclei, and the type of ANAs (and therefore type of illness) can be distinguished by examining the staining pattern of the cells. Problems such as the presence of noise, low level of standardization and errors in interpretation all lead to uncertainty in the diagnoses. Due to these challenges, computer aided diagnosis has been suggested as an alternative to examining the staining patterns [12].

Recently, an annotated dataset of six types of staining patterns from 28 subjects has been provided [7]

to facilitate the development of classifiers for the HEp-2 staining patterns. The dataset consists of 28 images of slides (one per subject), where each image contains several cells of the same class. A competition organized at the International Conference of Pattern Recognition 2012 revealed that there is large variability between cells of the same class, but across different subjects, leading to poor generalization. Our goal is to address this variability using a multiple instance learning 

[11] approach.

In multiple instance learning an object is represented by a collection of feature vectors [5]. Typically, this set is called a bag, and the feature vectors inside a bag are called the instances. This is an extension to the standard pattern recognition approach where objects are represented by a single feature vector only. By using a set of feature vectors, the representational capacity is enriched, which potentially increases the discriminability between classes. The obvious drawback is that for the classification of an object, the full set of feature vectors has to be taken into account. Depending on assumptions on how to deal with the set of feature vectors, several multiple instance learning approaches have been proposed. One branch of classifiers focuses on finding the single most informative feature vector (approaches such as miSVM [1], Diverse Density [11]

, and Expectation-Maximization 

[17] fall into this category). Another branch focuses on describing the full distribution of the set of feature vectors (with approaches such as Citation-NN [16], MILES [3], kernels [8] and bag dissimilarities [15, 4]).

When classifiying HEp-2 cells, the slides are the bags and the individual cells are the instances. Because cells in the same slide have the same label, selecting a single most informative cell per slide does not seem appropriate. Therefore we focus on approaches that describe the full distribution of the bags. We propose representing a bag by the quantiles of its instance distribution, an approach that is a generalization of the minimum-maximum representation in [8]

. In the minimum-maximum representation each bag is represented by the minimum and maximum feature value that appears in the bag. When the number of instances in a bag is large, these minima and maxima may be heavily influenced by outliers. By selecting appropriate quantiles, noisy feature values are avoided.

In Section 2 a short explanation of learning on sets of feature vectors is given. Next, in Section 3 the dataset and its preprocessing is explained, followed by the proposed approach for the classification of HEp-2 cells in Section 4. Finally, in Sections 5 and 6 experiments are shown and conclusions are drawn.

2 Learning with sets

In multiple instance learning (MIL) [11] and group-based classification (GBC) [13], objects are represented as sets of feature vectors, rather than single feature vectors only. In MIL, labels are available only at a coarse level, for a set (bag) of feature vectors (instances). Both training and test objects are bags. In GBC, only the test objects are bags, with the added information that all instances inside one bag are from the same class.

In the case of HEp-2 cell classification, we can consider slides as bags, and individual cells as the instances. There are two advantages to doing so.

In the training phase, a clear advantage of MIL is that learning is possible even if only coarse bag labels are availble. Although in the cell classification problem, the instance labels are available as well, there is an another advantage to using MIL in the training phase. In the case that the cells within one image are not independent, considering these cells jointly can be more informative than considering all cells individually.

In the testing phase, considering several cells jointly can also help to improve performance. This is illustrated in the following example. Consider the 1-dimensional binary classification problem in Fig. 1, and assume that we have found the Bayes optimal classifier. The shaded circles are the test set, and their true labels are . If we were to classify these instances independently, the error would be equal to , because the leftmost object will be misclassified to class . However, with the added information that these instances belong to a bag of objects from the same class, we could apply a majority voting combination to classify the bag as , and propagate the label to all the individual instances, reducing the error to 0.

)

)

Figure 1: 1-dimensional binary classification problem. The shaded circles are from the class. The added information that the circles are all of the same class helps to reduce the classification error.

Considering bags of instances rather than single instances leads to a different representation than in regular supervised learning methods. One object is represented by a bag

where is the number of instances in the bag, and which may be different for different bags. However, the disadvantage here is that supervised learning methods cannot be applied directly.

One possibility to both consider information from all instances, and to retain a single feature vector representation, is to represent a bag as a distribution in instance space [8, 14]. In [8], the so-called Minimax kernel represents each bag by the minimum and maximum values of its features, resulting in a -dimensional feature vector:

(1)

where selects the -th element of vector .

The same principle can be applied to other quantiles of the instances. When we again consider a bag with instances , we obtain the -th quantile for feature by first sorting all values of one feature, and then selecting the appropriate value from the sorted list:

(2)

In principle this quantile representation can be computed on bags of any size. When the number of instances per bag is small, many of these quantiles will coincide. Therefore the number of quantiles is typically chosen to be reasonably small, as a fraction of the average number of instances in a bag. For problems where the number of instances per bag is large, the extreme quantiles ( or , or the minimum or maximum feature value) can become noisy and in these situations these extreme quantiles should be avoided. For smaller bag sizes the minimum and maximum are often used and show good performance.

Figure 2: Scatterplot of one bag with 11 instances. From this bag the 10%-, 50% and 90%-quantiles are determined for each feature. This results in a 6D feature vector for this bag.

To illustrate the quantile representation, a scatterplot of a single bag with 11 instances is shown in Fig. 2. For each of the features the 10%-quantile the 50%-quantile (the median) and the 90%-quantile are computed. For this 2D feature space and these three different quantiles this results in a 6-dimensional feature vector for each bag.

Note that the correlation between the features is lost in this representation. The instance that has the highest feature value for the first feature is typically not the same instance with the highest feature value for the second feature. On the other hand, by considering each marginal distribution independently it is avoided that correlations and higher order statistics have to be estimated. For larger feature dimensionalities and finite datasets this is very challenging and should be avoided to not suffer from the curse of dimensionality.

3 Dataset

There are two different datasets available: a training set containing 721 cells scanned from 14 different images and a test set containing 734 cells scanned from 14 different images. All cells from one image share the same label. In six images all cells are “centromere”, in five images they are labeled “coarse speckled” and five “homogeneous”. For the classes “cytoplasmatic”, “fine speckled”, and “nucleolar” there are four images each. The images contain 52 cells on average, with a minimum of 13, and a maximum of 119 cells per image.

For each cell the green channel intensity is normalized using histogram equalization. On these normalized images Gabor features are computed. A Gabor filter is the product of a Gaussian kernel and a sinusoid:

(3)

where and . The parameter defines the scale of the filter, the parameter determines the direction in which the cosine intensities vary, and determines the wavelength of the cosine. In the experiments several values for these parameters are used: , and . In total 80 filtered images are obtained per cell.

To obtain a single feature vector for each cell, some simple statistics are derived from the image: the mean of the absolute values, the maximum, and the variance over all pixels within the mask. Next to that, the maximum, the variance and the mean of the image intensity are added. Finally, to obtain a rotational invariant representation of the cells, the outputs of the four different directions

are averaged. This results in a dimensional feature vector. The results in section 5 show that this representation is expressive enough for a reasonable cell classification performance.

4 Proposed approach

We propose to formulate the problem in the multiple instance learning setting, therefore learning on images rather than individual cells.

Our procedure is shown in Algorithm 1. In our implementation, the procedure extracts the Gabor features as described in Section 3, represents each bag by a distribution of its instances and and uses a logistic classifier. However, in principle these implementations can be replaced by other features, distributions and classifiers.

1:procedure ClassifyCells()
2:     for   do
3:         for   do
4:              
5:         end for
6:         
7:     end for
8:     
9:     
10:     return
11:end procedure
Algorithm 1 Classification procedure

We examined several quantiles for the cell data: minimum, 10%, 50%, 90% and maximum. The minimum quantile turned out too noisy, but the 10%-quantile gave very good results. On the other hand, the maximum performed better than the 90% quantile, therefore we used the 10%-, 50%- and 100%-quantiles as the basis for our approach.

Adding other quantiles to this list generally deteriorated the performance, unless a quantile was selected that was close to one of the already chosen extremes, such as 11%- or 99%. Our final implementation includes the 10%-, 11%-, 50%- and 100%-quantiles.

5 Experiments

First some results are shown to evaluate the standard feature-based approach using the Gabor features. Several classifiers have been trained on the individual cells, like the -nearest neighbor, the standard support vector classifier or the logistic classifier. The best performing classifier (with a small margin) is the support vector classifier (or the Liknon [2]). This classifier handles the multi-class problem by training six one-vs-all classifiers and combining the output using a max combiner. The performance on cell level and on image level is computed, where the label of an image is computed by majority voting over all predictions of the individual cells in that image.

We then present results of our quantile-based multiple instance learning approach, where whole bags are considered as train and test objects. Each bag is represented by the quantile representation. The logistic classifier is trained on the bags. The performance on bag level can therefore be derived directly, the performance on cell level is computed by propagating the image label to all individual cells in that image.

The evaluation is the same for both approaches. The evaluation contains two performance measures: the first is the classification error on the pre-defined test set, the second is the averaged classification error using leave-one-image-out cross-validation. Because there are 28 images in total, this means a 28-fold cross-validation. Next to the classification errors, the confusion matrices are reported as well.

All experiments are performed in MATLAB® with the PRTools toolbox [6].

5.1 Cell level

Split in predefined train and test set
cell level evaluation image level evaluation
estimated class
true Ce Co Cy Fi Ho Nu
Ce 94 0 0 1 4 50
Co 9 29 6 49 3 5
Cy 0 3 47 0 1 0
Fi 25 8 0 39 36 6
Ho 1 2 0 30 141 6
Nu 16 1 0 1 17 104
estimated class
true Ce Co Cy Fi Ho Nu
Ce 2 0 0 0 0 1
Co 0 1 0 2 0 0
Cy 0 0 2 0 0 0
Fi 0 0 0 1 1 0
Ho 0 0 0 0 2 0
Nu 0 0 0 0 0 2
61.85% correct 71.43% correct
28-fold cross-validation
cell level image level
estimated class
true Ce Co Cy Fi Ho Nu
Ce 239 36 0 37 4 41
Co 54 77 5 58 6 10
Cy 8 4 92 2 0 3
Fi 25 61 0 35 87 0
Ho 2 4 4 33 281 6
Nu 39 16 1 14 8 163
estimated class
true Ce Co Cy Fi Ho Nu
Ce 5 0 0 1 0 0
Co 1 1 0 3 0 0
Cy 0 0 4 0 0 0
Fi 0 1 0 1 2 0
Ho 0 0 0 0 5 0
Nu 0 0 0 0 0 4
60.96% correct 71.43% correct
Table 1: The confusion matrices and classification performances using the support vector classifier on the eight HEp cell classes: Ce:“centromere”, Co:“coarse speckled”, Cy:“cytoplasmatic”, Fi:“fine speckled”, Ho:“homogeneous”, Nu:“nucleolar”. Top row shows results on the predefined train and test set, the bottom row the results on 28-fold cross-validation.

The top row of Table 1

shows the confusion matrices and the classification performance on the predefined train and test set. In the top left, we evaluate the individual cells in the test set. The confusion matrix shows that there is large confusion between the classes “coarse speckled” and “fine speckled”, between “fine speckled” and “homogeneous”, and between “fine speckled” and “centromere”. The errors are not symmetric though, from the “coarse” class many cells are assigned to the “fine” class, but from the “fine” class many are also assigned to the homogeneous class. Apparently, there is a gradual change from “homegeneous”, via “fine” to “coarse”, and the decision boundaries between the classes are not exactly in between the classes. Around 62% of the cells in the test set is classified correctly. In the top right, the confusion matrix and performance is given for the image level classification. The classification accuracy improves a bit from

to , but still there is a large confusion between “coarse speckled” and “fine speckled”.

The second row of Table 1 shows the same setup with results for 28-fold cross-validation over all images. The overall performances do not differ a lot with the predefined train test split, but the confusion matrices show a similar confusion pattern between the classes.

true label corr. fraction
1 homogeneous 54 88.5%
2 fine speckled 5 10.4%
3 centromere 58 65.2%
4 nucleolar 28 42.4%
5 homogeneous 42 89.4%
6 coarse speckled 22 32.4%
7 centromere 47 83.9%
8 nucleolar 27 48.2%
9 fine speckled 1 2.2%
10 coarse speckled 9 27.3%
11 coarse speckled 23 56.1%
12 coarse speckled 20 40.8%
13 centromere 40 87.0%
14 centromere 27 42.9%
true label corr. fraction
15 fine speckled 24 38.1%
16 centromere 35 92.1%
17 coarse speckled 3 15.8%
18 homogeneous 31 73.8%
19 centromere 32 49.2%
20 nucleolar 43 93.5%
21 homogeneous 42 68.9%
22 homogeneous 112 94.1%
23 fine speckled 5 9.8%
24 nucleolar 65 89.0%
25 cytoplasmatic 13 54.2%
26 cytoplasmatic 31 91.2%
27 cytoplasmatic 36 94.7%
28 cytoplasmatic 12 92.3%
Table 2: Classification performance for each individual image. From left to right: the image number, the true class label, the number of cells that is correctly classified, and the fraction of cells that is correctly classified.

In Table 2 the output per image is shown; its true label, and the number and fraction of cells in this image that is correctly assigned to the true class. For most images the performance is good, but for some images the procedure completely fails: in particular on images 2, 9, 10, 17 and 23. These images are from the “coarse speckled” and “fine speckled” classes, so these results are in line with the confusion matrices.

We also investigated the effect of other cell label combining rules on the image level performance. The mean and product rule showed a slight improvement in performance, which is consistent with previous results that these combiners are more robust than majority voting [10]. However, the improvement was not very significant, suggesting that other measures are necessaary to further increase the performance.

5.2 Proposed: Image level

Here we show experiments for the multiple instance learning approach, where each image is represented by the set of quantile levels (10%, 11%, 50% and 100%). A logistic classifier is trained and tested on the images as a whole. Therefore, in the test phase, we first obtain image labels, and then propagate these labels to the corresponding cells.

In this task, we also evaluated several classifiers (1-norm SVM, SVM, logistic and nearest neighbor) and the logistic classifier performed the best on image level.

Split in predefined train and test set
cell level evaluation image level evaluation
estimated class
true Ce Co Cy Fi Ho Nu
Ce 84 0 0 0 0 65
Co 0 33 0 68 0 0
Cy 0 0 51 0 0 0
Fi 0 0 0 51 63 0
Ho 0 0 0 61 119 0
Nu 0 0 0 0 0 139
estimated class
true Ce Co Cy Fi Ho Nu
Ce 2 0 0 0 0 1
Co 0 1 0 2 0 0
Cy 0 0 2 0 0 0
Fi 0 0 0 1 1 0
Ho 0 0 0 1 1 0
Nu 0 0 0 0 0 2
65.0% correct 64.3% correct
28-fold cross-validation
estimated class
true Ce Co Cy Fi Ho Nu
Ce 357 0 0 0 0 0
Co 0 161 0 49 0 0
Cy 0 0 109 0 0 0
Fi 0 48 0 97 63 0
Ho 0 0 0 42 288 0
Nu 66 0 0 0 0 175
estimated class
true Ce Co Cy Fi Ho Nu
Ce 6 0 0 0 0 0
Co 0 4 0 1 0 0
Cy 0 0 4 0 0 0
Fi 0 1 0 2 1 0
Ho 0 0 0 1 4 0
Nu 1 0 0 0 0 3
81.6% correct 82.1% correct
Table 3: The confusion matrices and classification performances using the logistic classifier on image level on the six cell classes: Ce:“centromere”, Co:“coarse speckled”, Cy=“cytoplasmatic”, Fi:“fine speckled”, Ho:“homogeneous”, Nu:“nucleolar”. Top row shows results on the predefined train-test split, the bottom row the results on 28-fold cross validation.

The top row of Table 3 shows the results on the predefined train test split. Although the cell level performance improves slightly, the image level performance worsens. In other words, images with more cells are classified correctly more often, but in total less images are classified correctly. With just 14 training slides, the classifier is not able to generalize well to the predefined test set.

The second row of Table 3 shows the results of the cross-validation. Here the performances are significantly improved as opposed to the train test split. The classifier is able to form a better model of the six classes because more training slides are available. The improvement in cell level performance is especially large, over 20%. This demonstrates that cells that were previously misclassified, a classified correctly due to the other cells present in the image.

We believe that the cross-validation results are a better indicator of the performance than the results on the predefined train test split. In our investigation, we have experimented with alternative features and classifiers, and in all cases, high accuracies on the train test split corresponded with an overtrained classifier, and low performances in cross-validation.

Table 4 shows the performances of our approach on each individual image. Because the classification is done on bag level and the instance labels are propagated, it is only possible to classify 0% or 100% of the cells per image correctly. The images that are misclassified (2, 4, 12, 15 and 18) are somewhat different from the images in Table 2 that also had very poor performance: 2, 9, 10, 17 and 23. This suggests that combining the two approaches might lead to further improvement in accuracy.

true label corr. fraction
1 homogeneous 61 100.0%
2 fine speckled 0 0.0%
3 centromere 89 100.0%
4 nucleolar 0 0.0%
5 homogeneous 47 100.0%
6 coarse speckled 68 100.0%
7 centromere 56 100.0%
8 nucleolar 56 100.0%
9 fine speckled 46 100.0%
10 coarse speckled 33 100.0%
11 coarse speckled 41 100.0%
12 coarse speckled 0 0.0%
13 centromere 46 100.0%
14 centromere 63 100.0%
true label corr. fraction
15 fine speckled 0 0.0%
16 centromere 38 100.0%
17 coarse speckled 19 100.0%
18 homogeneous 0 0.0%
19 centromere 65 100.0%
20 nucleolar 46 100.0%
21 homogeneous 61 100.0%
22 homogeneous 119 100.0%
23 fine speckled 51 100.0%
24 nucleolar 73 100.0%
25 cytoplasmatic 24 100.0%
26 cytoplasmatic 34 100.0%
27 cytoplasmatic 38 100.0%
28 cytoplasmatic 13 100.0%
Table 4: Classification performance for each individual image. From left to right: the image number, the true class label, the number of cells that is correctly classified, and the fraction of cells that is correctly classified.

Even with the improved results, there is confusion between “fine speckled” and “coarse speckled”. We examined the images that are misclassified. The two “fine speckled” images 2 and 15 are of very low quality and the cells in them are nearly indistinguishable. The “coarse” speckled image 12, looks similar to “fine speckled” images in the dataset, at least to an untrained eye. Perhaps the extracted Gabor features do not sufficiently capture the differences between the two classes and features designed with the help of an expert could help to reduce the class overlap.

6 Conclusions

In this paper we propose a quantile representation to represent a set of feature vectors. This representation characterizes a collection of feature vectors, or a bag of instances in the terminology of multiple instance learning. For each bag, several quantile values for each of the features are computed. This results in a vector with fixed length for bags with a variable number of instances, a representation that can be used in any standard classifier. Although this representation ignores the correlation between the features, it characterizes the marginal feature distributions of the bags well.

This quantile representation is applicable for the classification of an image containing several HEp-2 cells. All of the cells in an image share the same label, and therefore all cells are informative for the image label. Experiments show reasonably good classification performance on image level, which also results in a better performance on cell level. Unfortunately, there is still confusion between the classes “fine speckled” and “coarse speckled”, and between “fine speckled” and “homogeneous”. For non-experts like the authors, these classes are indeed very difficult to distinguish.

A downside of the current quantile representation is that the quantile levels have to be chosen appropriately. When no prior knowledge is available on what quantile level might be informative, a quantile level selection is required. It would be interesting to investigate whether this can be done by visually inspecting the distributions of the classes for each feature.

The proposed representation is a general and intuitive representation that can be readily applied to other classification problems, and might therefore be helpful for other automatic diagnostic systems.

References

  • [1] S. Andrews, I. Tsochantaridis, and T. Hofmann. Support vector machines for multiple-instance learning. In Advances in Neural Information Processing Systems, pages 561–568. MIT Press, 2003.
  • [2] C. Bhattacharyya, L.R. Grate, A. Rizki, D. Radisky, F.J. Molina, M.I. Jordan, M.J. Bissel, and I.S. Mian. Simultaneous classification and relevant feature identification in high-dimensional spaces: Application to molecular profiling data. Signal Processing, 83:729–743, 2003.
  • [3] Y. Chen, J. Bi, and J.Z. Wang. MILES: Multiple-instance learning via embedded instance selection. IEEE Transactions on Pattern Analysis and Machine Intelligence, 28(12):1931–1947, 2006.
  • [4] V. Cheplygina, David M J Tax, and M. Loog. Does one rotten apple spoil the whole barrel? In International Conference on Pattern Recognition, pages 1156–1159, 2012.
  • [5] T.G. Dietterich, R.H. Lathrop, and T. Lozano-Perez. Solving the multiple instance problem with axis-parallel rectangles. Artificial Intelligence, 89(1-2):31–71, 1997.
  • [6] Robert P W Duin, P. Juszczak, P. Paclik, E. Pekalska, D. De Ridder, David M J Tax, and S. Verzakov. PRTools, Matlab toolbox for pattern recognition. http://www.prtools.org, 2010.
  • [7] P Foggia, G Percannella, P Soda, and M Vento. Early experiences in mitotic cells recognition on HEp-2 slides. In Computer-Based Medical Systems, pages 38–43. IEEE, 2010.
  • [8] Thomas Gärtner, Peter A Flach, Adam Kowalczyk, and Alex J Smola. Multi-instance kernels. In

    International Conference on Machine Learning

    , pages 179–186, 2002.
  • [9] Arthur Kavanaugh, Russell Tomar, John Reveille, Daniel H Solomon, and Henry A Homburger. Guidelines for clinical use of the antinuclear antibody test and tests for specific autoantibodies to nuclear antigens. Archives of pathology & laboratory medicine, 124(1):71–81, 2000.
  • [10] Josef Kittler, Mohamad Hatef, Robert P. W. Duin, and Jiri Matas. On combining classifiers. Pattern Analysis and Machine Intelligence, IEEE Transactions on, 20(3):226–239, 1998.
  • [11] O. Maron and T. Lozano-Pérez. A framework for multiple-instance learning. In Advances in Neural Information Processing Systems, pages 570–576. Morgan Kaufmann Publishers, 1998.
  • [12] Petra Perner, Horst Perner, and Bernd Müller. Mining knowledge for HEp-2 cell image classification. Artificial Intelligence in Medicine, 26(1):161–173, 2002.
  • [13] Noor A Samsudin and Andrew P Bradley. Nearest neighbour group-based classification. Pattern Recognition, 43(10):3458–3467, 2010.
  • [14] David M J Tax. Multiple instance learning using bag distribution parameters. In Benelux Conference on Artificial Intelligence, pages 226–233, Maastricht, 2012.
  • [15] David M J Tax, M. Loog, Robert P W Duin, V. Cheplygina, and Wan Jui Lee. Bag dissimilarities for multiple instance learning. Similarity-Based Pattern Recognition, pages 222–234, 2011.
  • [16] J. Wang and J.D. Zucker. Solving the multiple-instance problem: A lazy learning approach. In International Conference on Machine Learning, pages 1119–1125, 2000.
  • [17] Q. Zhang and S. Goldman. EM-DD: An improved multiple-instance learning technique. In Advances in Neural Information Processing Systems, volume 14. MIT Press, 2002.