Novel chromaticity similarity based color texture descriptor for digital pathology image analysis

Pathology images are color in nature due to the use of chemical staining in biopsy examination. Aware of the high color diagnosticity in pathology images, this work introduces a compact rotation-invariant texture descriptor, named quantized diagnostic counter-color pattern (QDCP), for digital pathology image understanding. On the basis of color similarity quantified by the inner product of unit-length color vectors, local counter-color textons are indexed first. Then the underlined distribution of QDCP indexes is estimated by an image-wise histogram. Since QDCP is computed based on color difference directly, it is robust to small color variation usually observed in pathology images. This study also discusses QDCP’s extraction, parameter settings, and feature fusion techniques in a generic pathology image analysis pipeline, and introduces two more descriptors QDCP-LBP and QDCP/LBP. Experimentation on public pathology image sets suggests that the introduced color texture descriptors, especially QDCP-LBP, outperform prior color texture features in terms of strong descriptive power, low computational complexity, and high adaptability to different image sets.


Introduction
Pathology is a medical sub-specialty that studies and practices the diagnosis of disease through examining biopsy samples or surgical specimens under microscopes by pathologists. It serves as the golden truth of cancer diagnosis. To address subjectivity in pathology examination [1,2], digital pathology exploits image analysis techniques and pattern recognition algorithms for histological information understanding in tissue images, and merges as a promising approach owing to its time-efficiency, consistency, and objectivity. Essentially, a digital pathology diagnosis system is a pattern recognition system. Given a query pathology image, a machine understands it by comparing a set of quantitative features from the image against the stored feature sets in the database. Hence, extraction of discriminative features from color pathology images is important.
construction seems similar to certain texture features, its motivation and focus on countercolor texture representation are novel and distinct.
In addition, we elaborate the computation of the QDCP descriptor in a generic pathology image analysis framework. Parameter settings and texture feature fusion techniques are discussed. Since QDCP contains complementary information to classic grayscale texture descriptors, two numerical features, QDCP-LBP and QDCP/LBP, are introduced for complete color texture representation. Experimentation demonstrates that the proposed descriptors greatly boost the performance of texture analysis in pathology images and always achieve top performance in different image sets. It is noteworthy that the high adaptability of the proposed QDCP-LBP feature to different datasets is very attractive in image description. It is true that one can use a specially-designed analysis tool for good performance if images' property and content are known beforehand. However generally, such information is not well known. Therefore, an analysis tool having such high adaptability is highly demanding in image analysis.
To summarize, the main contribution of this work are: 1. Motivated by high color diagnosticity in pathology images, we propose the use of counter color's spatial arrangement to characterize histological content composed of different stained tissue substances. To the best of our knowledge, this innovation is first explicitly explained and presented in color texture analysis literature, and provides a new vision to counter-color analysis in pathology images.

2.
A rotation-invariant color texture descriptor, QDCP, is introduced. It is advanced in counter-color content description for pathology image analysis. In a very compact form, QDCP is able to achieving better, or at least comparable, classification performance to the state-of-the-arts in different pathology image sets.
3. Two numerical color texture descriptors, QDCP-LBP and QDCP/LBP, are introduced based on different feature fusion techniques. With smaller computation complexity, experimentation indicates that the proposed texture features, especially QDCP-LBP, are capable of always achieving top performance in different pathology image analysis tasks. Such high adaptability enables QDCP-LBP to be a competitive descriptor in various image analysis tasks.
The rest of this paper is organized as follows. In Section of Prior Arts, the state-of-the-art color texture analysis methods are briefly reviewed. Section Methods specifies the construction of the proposed QDCP descriptor, and elaborates the computation of QDCP from pathology images in a general digital pathology pipeline. Experimental Design with applications of pathology images and the results are presented in Section Experimental Design and Results and Discussion, respectively. Finally, conclusions are drawn in the last section. correlation that exists between color channels. Hence the obtained texture descriptors may have abundant information redundancy, resulting in less compact descriptors. On the other hand, this paradigm may lose information about interactions between channels, leading to ineffective features. To obtain a compact color texture descriptor, the LBP operator is applied to the hue channel directly [24]. The obtained feature achieves good performance on the PAS-CAL visual object classes challenge 2007 image benchmark. To further address the circular nature of hue in color texture representation, a circular Hue-LBP (in short, CHLBP) [16] is proposed based on the angles between hue values. As the hue component is ill-defined for achromatic colors, the proposed Hue-LBP is unreliable in image areas containing a large amount of achromatic pixels.
The second category of algorithms considers color inter-channel dependency in texture feature extraction. In the opponent color (OC) texture description method [25], authors propose to apply scalar texture analysis methods to each pair of color channels for color texture quantification. For example, to compute OC-LBP from an RGB-format image, 6 LBP histograms are extracted from the red-red, green-green, blue-blue, red-green, red-blue, and green-blue channel pairs. But this method is computationally intensive. Instead of examining texture patterns within channel pairs, an image indexing study introduces a 3-D co-occurrence matrix to summarize the joint distribution of LBP in the red, green, and blue channels for an RGB image [26]. Although the so-called joint-LBP is demonstrated outperforming OC-LBP in endoscopic image classification [27], the disadvantage of joint-LBP is the large size of the color co-occurrence matrix whose elements may become sparse and unstable.
The last category of methods treats a color signal as a vector and quantifies color textures using vector processing. Norm-LBP is extracted from color vector's magnitude in the RGB or CIELAB color space [27]. Shortly, the local color vector binary patterns (LCVBP), composed of 3 angular texture patterns and the Norm-LBP, is proposed for face recognition [14]. In this work, the color angular patterns are quantified by applying the LBP operator to the relative phases within channel pairs in the YIQ color space. Due to high correlations between color channels, the three angular pattern descriptors may have information redundancy. To obtain a compact descriptor, a face expression recognition study introduces J-LCVBP, which adopts the sin distance [28] to measure similarity between two color vectors and uses of the J-th largest sin distance as a threshold to construct a LBP descriptor within a neighborhood [29]. Though J-LCVBP is good at addressing micro skin color difference associated with different expressions on human faces, the over-emphasis on small color variation in image smooth areas makes it less effective to describe counter-color textures in pathology images. Local color contrast (LCC) makes the use of angular difference to quantize color contrast statistics between a center pixel and the local mean color derived from its neighborhood [30]. LCC is an analogy to the local intensity standard deviation in a grayscale image, and irrelevant to local color texton description. To resist the color changes of an illuminant, local angular patterns (LAP) are extracted from the red-green, red-blue, and green-blue channel pairs in the RGB color space [31]. For each channel pair, angular differences between center pixels and their local means are quantized and fed to the LBP framework to generate a LAP. Due to high correlations between color channels, the 3 LAP descriptors are not compact and may have information redundancy.
The state of the art of color texture analysis algorithms with the proposed texture descriptor QDCP are summarized in Table 1. Note that though the vector-based color texture analysis methods are demonstrated more effective for color image analysis, they have different limitations as discussed above. Besides, as existing methods are based on color signal order statistics, which is sensitive to small color variation, they are not stable to represent medical information in terms of counter-color histological substances in pathology images.

Methods
In this section, we first introduce the novel color similarity based texture descriptor, called quantized diagnositc counter-color patters (QDCP). It is proposed for effectively coding medical information conveyed by counter colors in pathology images. Then related issues about applying QDCP to pathology image analysis are discussed. To clarify, vector variables in this work are denoted in boldface.

Quantized diagnostic counter-color patterns
In brief, given a color image I, its orientation image I O is generated using color's magnitudedirection (MD) representation [19] (or the so-called brightness-chromaticity (CB) model in literature). Based on predetermined neighborhood O N and threshold T O , various local counter-color texton in I O are indexed, forming a scalar image I QDCP . Finally, the statistics of local color texton are summarized by an image-wise histogram H QDCP , which is then used as a color texture feature in analysis. It is noteworthy that the introduced descriptor QDCP combines both structural and statistical texture analysis approaches. On one hand, a texton which represents spatial placement of counter-color elements in a small neighborhood is analyzed and indexed. On the other hand, the obtained histogram approximates the distribution of texton patterns in an image. Details about the introduced color texture descriptor is represented as follows.
MD representation of color images. A color image is represented by a function I: Z 2 ! Z 3 that maps a pixel p = (x, y) in the 2-dimensional image plane to a 3-dimensional vector I(p) = [i 1 , i 2 , i 3 ] T in a color space. For instance, in the RGB color model, i 1 , i 2 , and i 3 correspond to the red, green, and blue components, respectively. Since color vector's elements are highly correlated, chroma information is more closely related to the relative values among color channels, and color vector's orientation/direction correlates well with color description [18,[32][33][34]. Hence, to describe spatial arrangement of image elements in counter colors, we propose to extract the descriptor from color vectors' orientations. Specifically, any color vector I(p) = [i 1 , i 2 , i 3 ] T can be uniquely described by its magnitude I M (p) and orientation I O (p) in color magnitude-directional representation [19]: where I M ðpÞ ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi In (1), different from image magnitude I M (p) which is a scalar variable, the direction of a color vector is defined by a unit-length vector is a self-contained directional variable on the chromatic unit sphere, analysis of color texture based on it is also simplified and restricted on the unit sphere.
where To avoid the inverse trigonometric operation which is computational intensive, the inner product I O (p)I O (i) T is used to measure color similarity in this work. For ease of reference, in this work A small m p,i indicates two different color orientations, and m p,i = 1 implies that two color vectors are co-linear. Local texton indexing based on color similarity. We first consider color spatial distribution in a small image patch. For clarification, given a color orientation image I O , a small patch centered at a pixel p is composed by its N neighborhood pixels within a predefined area O N . In this work, to obtain rotation-invariant descriptor, the classical circular neighborhood is taken to define O N . For instance, given a center pixel, O 8 consists of 8 equally spaced pixels on a circle of unit radius around it, under which there are 36 different rotation-invariant patterns. It should be noted that in the circle neighborhood, there must be some neighbor points not on integer coordinates. Color vectors for these point are computed using the bilinear interpolation based on the nearest 4 pixels on image integer coordinates. The assumption behind this operation is that image color is smooth and without abrupt change.
For each pixel p in I O whose neighbors are all within an image, a vector M(p) = {m p,i , i 2 O N } is generated to quantify local color differences between the center pixel I O (p) and its neighbor pixels I O (i) for i 2 O N . That is, M(p) provides information on local color similarity with respect to the center pixel.
Note that the real-value quantity m p,i is sensitive to small color variation which may be introduced by inconsistent staining depths, imaging noise, or uneven illumination. To obtain a robust descriptor against small color variation, a binarization b p,i 2 {0, 1} with respect to a threshold T O is introduced to determine whether color difference m p,i is smaller enough to distinguish image colors. To be more precise, b p,i = 1 for m p,i < T O , which suggests a major color change due to different stained substances in images. With the obtained similarity vector M(p) and threshold T O , a value I QDCP (p) is generated to represent the rotation-invariant local counter-color pattern around pixel p as follows: Note the min(�) together with the mod (�) operations in (5) ensure QDCP's rotation-invariant property. Since a tissue slide may be placed along any orientation during imaging, the property of rotation invariance is of paramount significance in pathology image analysis. The binarization threshold T O can be either predetermined or adaptively set based on applications. It is noteworthy that for the same image, different T O leads to distinct color texture descriptors. Briefly, fine textures associated with small color change can be extracted with a large T O , and a smaller threshold on the color similarity vector results in a coarse color texture descriptor. We discuss the parameter setting of T O particularly in pathology image analysis in the next section.
Summary of texton statistics. From the viewpoint of texture statistics analysis approaches, texture is considered as a probabilistic generator of texton and the underlying probability distribution of texton can be used for image abstract representation. Following the histogram of equivalent patterns (HEP) paradigm [35], an image-wise histogram H QDCP , the approximation of color textons's distribution, is produced from the scalar image I QDCP for image characterization. Mathematically, the value in the j th bin in H QDCP is where δ(�) is the Dirac delta function. It is noteworthy that since QDCP is computed based on color difference which is less sensitive to image luma components, it is robust to vignette and inconsistent illumination in pathology images.

Pathology image analysis using QDCP
Overview of color texture based analysis pipeline. Fig 2 depicts pathology image analysis workflow in this work. Given a color pathology image I, it is first converted to the YCbCr color space, a linear transformation of the RGB domain in (7).
We select the YCbCr model to extract the QDCP descriptor for reasons as follows. First, in the YCbCr domain, the color similarity metric in (4) is more effective to distinguish stain colors in pathology images. For instance, in Masson's trichrome stained pathology images, collagen and mucus appear in green-blue colors, while muscle and cytoplasm are in red. In the RGB color space, all colors in the green-blue plane are orthogonal to red, which implies that obtained color similarity metric is 0 between green-blue colors and red. In contrast, after above nonorthonormal transformation in (7), this issue is solved and the similarity metric has more different values, which facilitates distinguishing histological objects in pathology images using QDCP. Second, the performance of QDCP depends on the selection of color difference threshold T O . In the RGB color space, a color vectors is valid with positive vector components in a 3-D Cartesian coordinate system, as depicted in Fig 1(a) Then the YCbCr-format image is decomposed into orientation image I O and magnitude image I M following (1), from which QDCP and LBP are extracted to form a long color texture feature f ct using different feature fusion techniques. Note, as QDCP describes color texture patterns in the orientation image only, texture patterns in the magnitude image should also be included in analysis. Given I M (p) being a linear variable, we select LBP to characterize content in the magnitude image due to the descriptive power of LBP. We want to point out that in the YCbCr domain, I M /I O decomposition is superior to the Luminance/Chromaticity separation (corresponding to Y-axis/CbCr-plane) in pathology image analysis for following reason. In the YCbCr color space, CbCr components are ill-defined when the vector is close to an achromatic color. Consequently, color analysis on the CbCr plane may be problematic, especially for pathology images which usually contains achromatic background areas. The decomposition of I M /I O can avoid this problem as an achromatic color in the YCbCr color space has a welldefined orientation, which can be denoted by a unit vector [1, 0, 0].
Finally, a classifier is applied to the long feature vector f ct to obtain image analysis result. Note that color variation is often observed in pathology images [36] and color texture features are reliable when examined pathology images have limited color variation [11,12]. Hence, in this paper, we assume that such appearance variation, if it exists, has been removed beforehand by color normalization approaches [37][38][39][40].
Selection of threshold T O for QDCP. To quantitatively represent diagnostic countercolor patterns by QDCP, an appropriate threshold T O is vital. Generally, a small threshold leads to a less-descriptive feature, while a very large T O makes QDCP sensitive to small color change introduced by imaging noise or non-uniform imaging illumination. It should be noted that selection of the threshold for a QDCP descriptor is application dependent. It is non-trivial to set the optimal threshold by a generic algorithm. So in this study, we provide several strategies to select T O for QDCP's construction in pathology image analysis. First, if prior knowledge about chemical dyes is known, a value which is slightly larger than the inner product of corresponding stains' color vectors can be used as T O . By using this threshold, QDCP is capable of extracting information about diagnostic counter-color patterns composed of stained histological components. Otherwise, we believe T O 2 (0.75 − 0.85), which corresponds to colors separating around 30-40 degrees, is a good start to compute QDCP from the YCbCr domain. Experimentation later demonstrates that T O 2 (0.75 − 0.85) is capable of generating considerably good classification in different datasets. In addition, the various QDCP descriptors generated under different thresholds can be combined together, so that both coarse and fine color textures are taken into consideration for image understanding.
QDCP and LBP feature fusion. QDCP and LBP represent image textures in the orientation image and the magnitude image, respectively. To form a complete color texture feature, we can combine QDCP and LBP by concatenation or joint manner.
• In the concatenation manner, we compute individual histograms of QDCP and LBP indexes.
Then the two histograms are concatenated together. In this work, the concatenation fusion is represented by QDCP-LBP.
• In the joint fusion method, a joint 2-dimensional histogram is generated, where co-occurrence of QDCP and LBP indexes is recorded. Then the 2-dimensional histogram is vectorized to form a long feature vector. The joint fusion manner is represented by QDCP/LBP in this paper.
Compared to the joint fusion method which generates a long feature vector, the concatenation manner is computationally efficient with the penalty of discarding the co-occurrence information between QDCP and LBP.

Experimental design
In this section, experimental design with applications in pathology image analysis, which includes information about testing pathology image sets and evaluation methodology, is specified.

Pathology image sets
This study takes two public pathology image sets, the GlomDB glomeruli dataset [12] and the ALL-IDB2 dataset [42], as experimental data to evaluate the proposed texture descriptor.
The GlomDB glomeruli dataset [12] is published for color and texture descriptor evaluation. The image set contains 1976 16-by-16 non-overlapping square patches of textures selected by manual segmentation from 15 Masson's trichrome stained renal biopsy samples. Among the 1976 sub-images, half textures correspond to glomeruli while the other half are non-glomeruli patches. Fig 4 shows a portion of one renal biopsy image in the GlomDB set, where two glomeruli are observed. Compared to other tissue substances in the renal image that have relatively small color variations, color texture patterns in glomerulus areas are more complicated.
The ALL-IDB2 image set [42] is published for testing the performances of classification systems on blood pathology images. Acute Lymphoblastic Leukemia (ALL) is a serious blood cancer that can be fatal for children. In pathology, identification of blast cells in microscopic images of blood samples is essential in ALL diagnosis. Though morphology characteristics of white blood cells are considered the classical features for distinguishing normal lymphocytes cell and lymphoblasts in ALL images, texture feature LBP is demonstrated having good performance in this classification scenario [9,10]. Hence, we include this image set in this study to examine our new color texture descriptor. Specifically, the ALL-IDB2 dataset contains 260 cropped area of interest of normal and blast cells from 108 blood sample images generated by Canon PowerShot G5 camera and stored with 24 bit color depth in the RGB format. Among the 260 cell images, 130 images contain lymphoblasts from ALL patients and the rest images have normal white blood cells. Fig 5 shows examples of lymphoblasts and normal lymphocytes cell images in the ALL-IDB2 image set.

Evaluation protocols
Evaluation of the proposed descriptors is based on pathology images classification in either the GlomDB image set or the ALL-IDB2 dataset. Following the procedure described in previous section, texture features (QDCP, QDCP-LBP, and QDCP/LBP) are extracted and passed to the classifier. To obtain a solid conclusion, fisher's linear discrinminant (FLD) classifier, the support vector machine (SVM) with a radial basis function (RBF) kernel, and K-nearest neighbors (KNN) are exploited to evaluate features' discriminative power. Specifically, we randomly partition the image set into training and testing set following the 10-fold nested cross-validation methodology. In the inner 5-fold validation loop of the training stage, parameters of the proposed descriptor T O and hyper-parameters of classifiers (σ in the RBF kernal and the box constraint for declassification penalty of SVM, and the number of neighbors k and the distance metric of KNN) are automatically selected to optimize classification accuracy. In the outer testing evaluation loop, the obtained threshold is applied to extract QDCP features which are then passed to the trained classifier.
The agreement between groundtruth and classification results is estimated using two metrics, which are classification accuracy (ACC) and the receiver operating characteristic (ROC) curve analysis. ACC is intuition and easy to compute, representing the probability of a correct classification for a query image. ROC analysis is a more comprehensive measurement than ACC, as it is more statistically consistent in classification evaluation [43]. To quantitatively compare ROC curves of various color texture descriptors, the area under the ROC curve (AUC) is calculated. Both ACC and AUC are in the range of [0, 1], and large values indicate better classification. The 10-fold cross-validation is repeated 10 times, and the statistics of ACC/AUC are summarized based on the 100 performance indexes.
Based on above evaluation procedure, two comparison experiments are performed. In the first experiment, the QDCP descriptor is extracted from the RGB domain, the YCbCr domain, Color texture descriptor in pathology image analysis the CIELab domain, and the I1H2H3 domain. The first three color spaces are the most common color spaces used in texture analysis literature, and the I1H2H3 color space is recommended for the Masson's trichrome stained images in the GLombDB set [12].
In the second experiment, QDCP-based descriptors are compared to prior arts of color texture features. In specific, the grayscale LBP is applied to image luma component only and used as comparison baseline. The non-vector processing methods, ind-LBP and OC-LBP [25], are included in our comparison experiment. Ind-LBP is obtained by concatenating LBP histograms independently extracted from the three channels in the YCbCr color space, and the OC-LBP computes 6 LBP histograms from channel pairs in the YCbCr domain. Finally, 4 state-of-the-art vector-processing based color texture features are examined. Norm-LBP applies the LBP operator to color vector's magnitude [27]. LCVBP extracts 3 color angular patterns from inter-channel pairs and concatenates them with norm-LBP [14]. LCC-LBP combines the local contrast histogram to LBP descriptor extracted from the image luma component [30]. And LBP-LAP combines 3 angular feature patters computed from inter-channel pairs and 3 LBP histograms from the P1P2P3 space [31]. Note in this work the circular neighborhood O 8 is applied to all descriptors.

Results and discussion
Quantitative results for each image set are presented and compared to prior works. The sensitivity analysis of QDCP to T O is examined. All the algorithms and evaluations in this study are implemented in Matlab and accessible to audience. Table 2 summarizes glomerulus image classification in the fist comparison experiment that evaluates discriminative power of various color spaces. In this pathology texture analysis image set, QDCP is most descriptive in the YCbCr domain. Table 3 summarizes glomerulus image classification using different color texture features over the GlomDB image set, where the top four classification results are highlighted. To visualize the glomerulus classification, ROC curves associated with SVM are depicted in Fig 6(a). Because we observe similar trends in the results associated with FLD and KNN, their ROC curves are omitted here.

Glomerulus texture classification on GlomDB images
The proposed QDCP-based descriptors (QDCP, QDCP-LBP, and QDCP/LBP) achieve the top performance. This is because QDCP is very descriptive for glomerulus textures that are composed of different stained histological substances as shown in Fig 6(a). In contrast, as color texture descriptors in prior arts are usually based on signal order statistics which are sensitive to small color variation, the resulting features are not reliable in such a scenario. In addition, as expected, QDCP-LBP and QDCP/LBP improve the classification performance compared to QDCP. Four extra observations are obtained in the GlomDB image set. First, ind-LBP and OC-LBP take color information into account, and thus improve classification performance compared to LBP. Second, norm-LBP and LAP achieves no better result than the comparison baseline, grayscale LBP. For one hand, both norm-LBP and LBP focus on numerical representation of textures associated with color brightness; On the other hand, since LAP is obtained based on ordering of color vectors' similarity, counter-color texture is prone to being overwhelmed by textures composed by small color variations. Third, LCC has the worst performance, as LCC is an analogy to the local intensity standard deviation in a grayscale image, and less descriptive to color textures. Fourth, both LCVBP and LAP-LBP include LBP and the angular texture patterns extracted from inter-channel pairs, they boost the classification performance compared to LBP. Table 4 reports lymphoblast classification in the color space comparison experiment. In the ALL-IDB2 image set, the YCbCr color space and the I1H2H3 color space achieve similar performance. The CIELAB color domain is in the middle, and the RGB color space lags behind.

Lymphoblast classification in ALL-IDB2 images
Lymphoblast classification in the ALL-IDB2 image set is reported in Table 5 and Fig 6(b). Again, the top four classification results are highlighted in Table 5. In the ALL-IDB2 dataset, QDCP, QDCP-LBP and LCVBP obtain the top performance for all three classifiers, but the proposed descriptors are more compact than LCVBP, with less than 1/2 size of LCVBP. We also notice that performance of the joint QDCP/LBP is about 8% worse than QDCP-LBP. Theoretically, QDCP/LBP should outperform QDCP-LBP, since QDCP/LBP also maintains the co-occurrence information between QDCP and LBP indexes. However, it should be noted that the information gain of QDCP/LBP is with a penalty of relatively large feature dimension. If the experimental image is not big enough, the resulted QDCP/LBP histogram may be sparse and sensitive to noise and image artifact, which degrades the performance of QDCP/LBP. We  believe this is the major reason for the contradictory observation on performance of QDCP/ LBP and QDCP-LBP in the glomerulus texture classification and the lymphoblast classification. It should be noted that it is not also the case in practice that the properties of targeting images are known. Hence, we suggest to try both QDCP-LBP and QDCP/LBP for better results. Different from observations in the GlomDB image set, LBP obtains comparable classification results to other textures. This may be because lymphoblast image color in ALL-IDB2 dataset is much smoother and has less color texture patterns. Another interesting observation is about OC-LBP's classification results. There is a big performance gap between FLD and the other two classifiers. We believe the reason is that OC-LBP features associated with normal and lymphoblast images are not linearly separable and both SVM and KNN are capable of better handling this nonlinear situation.

Sensitivity analysis of QDCP against parameter setting
The effectiveness of the proposed QDCP descriptor depends on the selection of threshold T O in (5). To examine the sensitivity of QDCP against T O , after representing a query image in the YCbCr domain, we manually vary the value of the parameter when calculating QDCP and feed the obtained QDCP descriptor to a classifier. The average performances in the two testing data sets are report in   Color texture descriptor in pathology image analysis Specifically, in Fig 7(a) which corresponds to GlomDB color texture classification, QDCP is most descriptive when threshold T O is set around 0.85 (which is about 30 degree of color difference) . Fig 7(b) depicts the effects of T O on lymphoblast classification. QDCP obtains strong discriminative power with T O = 1 and T O = 0.8, which corresponds to 0 degree and 36 degree of color difference in the orientation domain. An interesting observation is that in lymph cell image classification, T O = 1 leads to the best performance. This may be because image color in the ALL-IDB2 set is very smooth, and all color change in images are very informative for down-streaming classification.
It is noteworthy that the golmerulus images contains many counter-color content, whereas lymph cell images, in contrast, are relatively dull. On one hand, though the two image sets have distinct image properties, generally, QDCP is descriptive when T O 2 [0.8, 0.9]. We believe this observation is applicable to other images. On the other hand, through all experiments in this study, we notice that QDCP-based descriptor, especially QDCP-LBP, have high adaptability to different image sets. We believe that the flexibility of T O contributes to such adaptability. In other words, different T O makes the proposed features capable of describing different textures in images. This adaptability is very attractive in image analysis, because usually image properties are not well known beforehand.

Conclusion
Color texture patterns in pathology images provide insightful information for disease diagnosis. Motivated by the close relation between color diagnosticity and counter-color histological structures in pathology images, this paper introduced a novel compact numerical descriptor, QDCP, for image counter-color texture representation. After projecting image colors to a unit-radius sphere, angular difference between center pixels and their neighbor pixels were summarized in an image-wise histogram. On the basis of QDCP, two new color descriptors, QDCP-LBP and QDCP/LBP, were introduced. Experimentation on publicly available pathology image sets suggested that the proposed descriptors, especially QDCP-LBP, were very descriptive and outperformed state-of-the-art color texture descriptors in terms of discriminative power, computational efficiency and adaptability to different image sets.
Due to the promising performance of the proposed descriptors, our work can be extended in two directions in future. First, observing the high adaptability to different pathology image sets, we want to apply these descriptors to nature texture images, examining their performance in terms of description and adaptability. Second, many algorithms proposed to improve LBP can be adopted by the QDCP-based descriptors. For instance, we may increase the neighborhood dimensionality, introduce multi-thresholding and fuzzy-thresholding descriptors, and examine the uniform patterns on the basis of QDCP.