Email updates

Keep up to date with the latest news and content from Diagnostic Pathology and BioMed Central.

This article is part of the supplement: New trends in digital pathology: Proceedings of the 9th European Congress on Telepathology and 3rd International Congress on Virtual Microscopy

Open Access Proceedings

Colour model analysis for microscopic image processing

Gloria Bueno1*, Roberto González1, Oscar Déniz1, Jesús González2 and Marcial García-Rojo2

Author Affiliations

1 Escuela Técnica Superior de Ingenieros Industriales, Universidad de Castilla-La Mancha, Avda. Camilo José Cela, s/n, 13071 Ciudad Real, Spain

2 Hospital General de Ciudad Real. Calle Tomelloso s/n. 13004 Ciudad Real, Spain

For all author emails, please log on.

Diagnostic Pathology 2008, 3(Suppl 1):S18  doi:10.1186/1746-1596-3-S1-S18


The electronic version of this article is the complete one and can be found online at: http://www.diagnosticpathology.org/content/3/S1/S18


Published:15 July 2008

© 2008 Bueno et al; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

This article presents a comparative study between different colour models (RGB, HSI and CIEL*a*b*) applied to a very large microscopic image analysis. Such analysis of different colour models is needed in order to carry out a successful detection and therefore a classification of different regions of interest (ROIs) within the image. This, in turn, allows both distinguishing possible ROIs and retrieving their proper colour for further ROI analysis. This analysis is not commonly done in many biomedical applications that deal with colour images. Other important aspects is the computational cost of the different processing algorithms according to the colour model. This work takes these aspects into consideration to choose the best colour model tailored to the microscopic stain and tissue type under consideration and to obtain a successful processing of the histological image.

Background

A challenge still facing scientists is the efficient analysis and management of biomedical data, including images. Advances in biomedical imaging diagnosis have been possible thanks to the development of new imaging technologies. Anatomical Pathology has also benefited from these new technologies, which have provided solutions for whole slide scanning by means of motorized microscopes and scanners [1], that is, whole slide imaging (WSI). However, the image processing performed with these slides is still limited both in data processed and processing methods.

Much research has been carried out on the development of algorithms for histological image analysis. Most of them are based on the segmentation of just one region of interest (ROI), which is usually the nucleus, and its classification for diagnosis purposes. To this end, statistical information techniques, region growing algorithms, active contour models and morphological methods have been used for ROI detection and processing [2-5].

The main problem with these methods is that they are not designed to process large amounts of data, which is the case when working with WSI in pathology. Besides that, many of these methods show limited results because they are mainly focused on a single structure or a type of tissue.

There is a need to develop more general and efficient image processing methods. To this end the colour model should be analysed, as well as the distance colour model applied to the processing algorithm in order to reduce the computational cost and obtain, in an efficient way, a set of heterogeneous, complex and specific image analysis. In this work different colour models and distances have been studied and applied under a general parallel image-processing model designed and implemented with MPP (Massively Parallel Processing).

Methods

There are three main colour models:

RGB: channel Red, channel Green and channel Blue,

HSI: channel Hue, channel Saturation and channel Intensity,

L*a*b*: channel Luminance, channel a*, that is, range of channel between Red to Green and channel b* that is range of colours between Yellow to Blue.

All colour models have their advantages and drawbacks. It is necessary to identify which colour model is suitable to represent and reproduce the ROI under consideration for each tissue type and WSI modality. Analysing the distance colour formulae applied between two colours may do this, d ( x , y ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemizaq2aaeWaaeaacuWG4baEgaWcaiabcYcaSiqbdMha5zaalaaacaGLOaGaayzkaaaaaa@32A7@ , x = ( x 1 , x 2 , x 3 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmiEaGNbaSaacqGH9aqpcqGGOaakcqWG4baEdaWgaaWcbaGaeGymaeJaeiilaWcabeaakiabdIha4naaBaaaleaacqaIYaGmaeqaaOGaeiilaWIaemiEaG3aaSbaaSqaaiabiodaZaqabaGccqGGPaqkaaa@39BB@ , y = ( y 1 , y 2 , y 3 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmyEaKNbaSaacqGH9aqpcqGGOaakcqWG5bqEdaWgaaWcbaGaeGymaedabeaakiabcYcaSiabdMha5naaBaaaleaacqaIYaGmaeqaaOGaeiilaWIaemyEaK3aaSbaaSqaaiabiodaZaqabaGccqGGPaqkaaa@39C3@ .

The distances considered within this study are: the Euclidean distance for the RGB model (Equation 1), the NBS colour distance formulae for HSI model (Equation 2) and the CIEDE2000 for the CIEL*a*b*, colour model (Equation 3).

d ( x , y ) = ( x 1 y 1 ) 2 + ( x 2 y 2 ) 2 + ( x 3 y 3 ) 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemizaq2aaeWaaeaacuWG4baEgaWcaiabcYcaSiqbdMha5zaalaaacaGLOaGaayzkaaGaeyypa0ZaaOaaaeaadaqadaqaaiabdIha4naaBaaaleaacqaIXaqmaeqaaOGaeyOeI0IaemyEaK3aaSbaaSqaaiabigdaXaqabaaakiaawIcacaGLPaaadaahaaWcbeqaaiabikdaYaaakiabgUcaRmaabmaabaGaemiEaG3aaSbaaSqaaiabikdaYaqabaGccqGHsislcqWG5bqEdaWgaaWcbaGaeGOmaidabeaaaOGaayjkaiaawMcaamaaCaaaleqabaGaeGOmaidaaOGaey4kaSYaaeWaaeaacqWG4baEdaWgaaWcbaGaeG4mamdabeaakiabgkHiTiabdMha5naaBaaaleaacqaIZaWmaeqaaaGccaGLOaGaayzkaaWaaWbaaSqabeaacqaIYaGmaaaabeaaaaa@506E@ (1)

d ( x , y ) = 1.2 2 x 2 y 2 ( 1 cos ( 2 π Δ H 100 ) ) + Δ S 2 + ( 4 Δ I ) 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemizaq2aaeWaaeaacuWG4baEgaWcaiabcYcaSiqbdMha5zaalaaacaGLOaGaayzkaaGaeyypa0JaeGymaeJaeiOla4IaeGOmaiJaey4fIOYaaOaaaeaacqaIYaGmcqWG4baEdaWgaaWcbaGaeGOmaidabeaakiabdMha5naaBaaaleaacqaIYaGmaeqaaOWaaeWaaeaacqaIXaqmcqGHsislcyGGJbWycqGGVbWBcqGGZbWCdaqadaqcfayaamaalaaabaGaeGOmaiJaeqiWdaNaeuiLdqKaemisaGeabaGaeGymaeJaeGimaaJaeGimaadaaaGccaGLOaGaayzkaaaacaGLOaGaayzkaaGaey4kaSIaeuiLdqKaem4uam1aaWbaaSqabeaacqaIYaGmaaGccqGHRaWkdaqadaqaaiabisda0iabfs5aejabdMeajbGaayjkaiaawMcaamaaCaaaleqabaGaeGOmaidaaaqabaaaaa@5B4C@ (2)

Δ E 00 = ( Δ L K L S L ) 2 + ( Δ C K C S C ) 2 + ( Δ H K H S H ) 2 + R T ( Δ C K C S C ) ( Δ H K H S H ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeuiLdqKaemyrau0aaSbaaSqaaiabicdaWiabicdaWaqabaGccqGH9aqpdaGcaaqaamaabmaajuaGbaWaaSaaaeaacqqHuoarcuWGmbatgaqbaaqaaiabdUealnaaBaaabaGaemitaWeabeaacqWGtbWudaWgaaqaaiabdYeambqabaaaaaGccaGLOaGaayzkaaWaaWbaaSqabeaacqaIYaGmaaGccqGHRaWkdaqadaqcfayaamaalaaabaGaeuiLdqKafm4qamKbauaaaeaacqWGlbWsdaWgaaqaaiabdoeadbqabaGaem4uam1aaSbaaeaacqWGdbWqaeqaaaaaaOGaayjkaiaawMcaamaaCaaaleqabaGaeGOmaidaaOGaey4kaSYaaeWaaKqbagaadaWcaaqaaiabfs5aejqbdIeaizaafaaabaGaem4saS0aaSbaaeaacqWGibasaeqaaiabdofatnaaBaaabaGaemisaGeabeaaaaaakiaawIcacaGLPaaadaahaaWcbeqaaiabikdaYaaakiabgUcaRiabdkfasnaaBaaaleaacqWGubavaeqaaOWaaeWaaKqbagaadaWcaaqaaiabfs5aejqbdoeadzaafaaabaGaem4saS0aaSbaaeaacqWGdbWqaeqaaiabdofatnaaBaaabaGaem4qameabeaaaaaakiaawIcacaGLPaaadaqadaqcfayaamaalaaabaGaeuiLdqKafmisaGKbauaaaeaacqWGlbWsdaWgaaqaaiabdIeaibqabaGaem4uam1aaSbaaeaacqWGibasaeqaaaaaaOGaayjkaiaawMcaaaWcbeaaaaa@69D6@ (3)

Where KL, Kc, KH are weight factors and the rest of components, SL, SC, SH, C', H', may be calculated by means of the, {L*, a*, b*} coordinates [6].

Moreover, another aspect to be considered is how to deal with the colour coordinates, that is as a vector or in a marginal way. These aspects have been analysed within this work. To this end the 3*2 distances to the most representative colour ROIs and statistically identified on the image were calculated on different WSI, which is to prostate biopsies and lung cytology stained with hematoxiline-eosine (HEO), inmunohistochemistry and papanicolau. The images were obtained by the ALIAS II automatic microscope and the processing was done using our own libraries, implemented by the research group, running under MPI on a grid composed by 17 nodes Intel XEON (3,2 GHz) INFINIBAND net (10 GB full-duplex) architecture. The results are shown as follows.

Results

The results applied to microscopic images show that the Euclidean and NBS vector distance for the RGB and HSI model respectively distinguish between different ROIs but the vector CIEDE2000 distance for the CIEL*a*b* model reproduces in a better way, the original colour. However, the computational cost of the last one is higher than the other two colour models.

Figure 1 shows the result for a biopsy stained with Hematoxiline-Eosine and Figure 2 for a cytology stained with Papanicolau.

thumbnailFigure 1. Colour distances for ROI detection applied to biopsies. a) Original Image, b) RGB (Euclidean), c) HSI (NBS), d) CIEL*a*b (CIEDE2000). Colour distances for ROI detection applied to biopsies.

thumbnailFigure 2. Colour distances for ROI detection applied to cytology. a) Original Image, b) RGB (Euclidean), c) HSI (NBS), d) CIEL*a*b (CIEDE2000). Colour distances for ROI detection applied to cytology.

The computational cost for the three colour distance vector models against different number of ROIs is shown in Figure 3.

thumbnailFigure 3. Computational cost of the colour distances vs. number of ROIs. Computational cost for the different colour models against the number of ROIs analysed.

To quantify the goodness of the distance formulae a ROC analysis has been carried on. Figure 4 shows this analysis for two ROIs in a prostate biopsy at 10× stained with HEO. The true pixels belonging to the ROIS were indicated by experts at Hospital General Ciudad Real. Figure 4b) and 4c) show the true values for the two regions of interest, the glandular light and the nucleus, extracted from the original image (Figure 4a). Figure 4d) to i) show the different colour distance results for these two regions. Finally, Table 1 shows the ROC analysis for the Eucludian, NBS and CIEDE2000 colour distance for the RGB, HSI and CIEDEL*a*b* models. It is shown that the % of specificity is higher for the CIEDE2000 distance with lower value of FP.

Table 1. ROC analysis of the colour distance formulae for two ROIs

thumbnailFigure 4. Colour distance validation. a) Original Image, b) True section of glandular light, c) True section of nucleus, d) RGB (Euclidean), e) HSI (NBS), f) CIEL*a*b (CIEDE2000), g) RGB (Euclidean), h) HSI (NBS), i) CIEL*a*b (CIEDE2000). Colour distances validation on different ROIs.

Conclusion

This article has presented a comparative study between RGB, HSI and CIEL*a*b* colour models applied histological images. This analysis, in turn, allows both distinguishing possible regions of interest and retrieving their proper colour for further region analysis.

The results applied to prostate biopsies stained with HEO and lung cytologies stained with papanicolau show that the vector CIEDE2000 distance for the CIEL*a*b* model reproduces in a better way the original colour.

Therefore, this comparison does allow us to choose the best colour model tailored to the microscopic stain and tissue type under consideration to obtain a successful processing. Moreover, a compromise between the computational cost and the results focus always to distinguish between different colour detection and colour retrieval for further ROI analysis should be kept. The colour model should be taken into consideration when defining standards for histological images.

Acknowledgements

This work has been possible thanks to project MEC DPI2004-01346 carried out at Engineering School in collaboration with Hospital General at Ciudad Real.

This article has been published as part of Diagnostic Pathology Volume 3 Supplement 1, 2008: New trends in digital pathology: Proceedings of the 9th European Congress on Telepathology and 3rd International Congress on Virtual Microscopy. The full contents of the supplement are available online at http://www.diagnosticpathology.org/supplements/3/S1

References

  1. Rojo MG, García GB, Mateos CP, García JG, Vicente MC: Critical comparison of 31 commercially available digital slide systems in pathology.

    Int J Surg Pathol 2006, 14:285-305. PubMed Abstract | Publisher Full Text OpenURL

  2. Gao M, Bridgman P, Kumar S: Computer aided prostate cancer diagnosis using image enhancement and JPEG2000.

    Proc SPIE International Conf on Applications of Digital Image Processing 2003, 5203:323-334. OpenURL

  3. Li Y, Lu D, Lu X, Liu J: Interactive colour image segmentation by region growing combines with image enhancement based on bezier model.

    Proceedings of the Third International Conference on Image and Graphics (ICIG'04) 2004, 100:96-99. OpenURL

  4. Begelrnan G, Gur E, Rivlin E, Rudzsky M, Zalevsky Z: Cell nuclei segmentation using fuzzy logic engine.

    IEEE Intern Conference on Image Processing 2004, 2937-2940. OpenURL

  5. Hu MP, Ding XY: Automated cell nucleus segmentation using improved snake.

    Proc Intern Conference on Image Processing 2004, 4:2737-2740. OpenURL

  6. Wyszecki G, Stiles WS: Color Science: Concepts and Methods, Quantitive Data and Formulae. 2nd edition. Wiley; 2000.