 Research
 Open Access
 Published:
Segmentation in dermatological hyperspectral images: dedicated methods
BioMedical Engineering OnLine volume 15, Article number: 97 (2016)
Abstract
Background
Segmentation of hyperspectral medical images is one of many image segmentation methods which require profiling. This profiling involves either the adjustment of existing, known image segmentation methods or a proposal of new dedicated methods of hyperspectral image segmentation. Taking into consideration the size of analysed data, the time of analysis is of major importance. Therefore, the authors proposed three new dedicated methods of hyperspectral image segmentation with special reference to the time of analysis.
Methods
The segmentation methods presented in this paper were tested and profiled to the images acquired from different hyperspectral cameras including SOC710 Hyperspectral Imaging System, Specim sCMOS50V10E. Correct functioning of the method was tested for over 10,000 2D images constituting the sequence of over 700 registrations of the areas of the left and right hand and the forearm.
Results
As a result, three new methods of hyperspectral image segmentation have been proposed: fast analysis of emissivity curves (SKE), 3D segmentation (S3D) and hierarchical segmentation (SH). They have the following features: are fully automatic; allow for implementation of fast segmentation methods; are profiled to hyperspectral image segmentation; use emissivity curves in the model form, can be applied in any type of objects not necessarily biological ones, are faster (SKE—2.3 ms, S3D—1949 ms, SH—844 ms for the computer with Intel^{®} Core i7 4960X CPU 3.6 GHz) and more accurate (SKE—accuracy 79 %, S3D—90 %, SH—92 %) in comparison with typical methods known from the literature.
Conclusions
Profiling and/or proposing new methods of hyperspectral image segmentation is an indispensable element of developing software. This ensures speed, repeatability and low sensitivity of the algorithm to changing parameters.
Background
Image segmentation is the process of dividing an image into parts defined as areas which are homogeneous in terms of selected properties. Today, segmentation is one of the most widely used [1, 2] and overused word in the area of biomedical image analysis (PubMed database contains 23,202 results—using the search strategy “all fields”, scopus—150,147 results—using the search strategy “all fields”). Practically everywhere where a region of interest (ROI) is separated using any method, the authors call it segmentation [3, 4]. Of course, in some cases, this is justified [5]. Formally, however, segmentation of images is divided into four main groups [6–12]:

point operations,

edge methods,

region methods and

hybrid methods.
Today, the most developed are region methods which include [13–16]:

region growing,

region merging,

region splitting,

split and merge,

watershed segmentation.
Not all segmentation methods can be applied in the analysis of medical hyperspectral images [17]. The point methods based on the selection of an appropriate binarization threshold are the most common. The watershed and hybrid methods are equally popular. Their application is limited due to two main elements. The first one is segmentation of an object whose limit, contour, is a change in the brightness of pixels. It is most often calculated for one or three (RGB) brightness levels. The second element is the lack of dependence of segmentation on the reference (expected) spectral spectrum. Therefore, in comparison to the segmentation of images in grey levels or colour images, hyperspectral imaging gives much more opportunities [18]. However, this excess of data is not, in every case, used by the authors of the algorithm (the authors of the article). The main limitation is the abovementioned lack of methods dedicated for hyperspectral image segmentation. This type of dedicated methods should use, apart from the conventional 2D image analysis, spatial information for individual wavelengths. One of the significant factors facilitating this type of analysis is organization of data written to *.raw, *.dat or *.cube files by a hyperspectral camera. Typical organization of this type of data is shown in Fig. 1.
The data are saved in a file by a hyperspectral camera (Fig. 1a) sequentially, starting with the full spectral range for the first (mth) row of the matrix. In the next steps, data for subsequent rows of the matrix of the image L _{ GRAY }(m, n, i) are saved in the file. Therefore segmentation methods can be easily implemented (minimizing computational complexity) using information about the spectral amplitude values stored for each wavelength for one pixel [19, 20]. In the stored *.raw, *.dat or *.cube file this is one (nth) column (see Fig. 1a). In terms of data organization in a file, access to these data is at the beginning (Fig. 1a, b). This type of data organization has a significant influence on the access time. Moreover, an additionally favorable element is the organization of particular first rows of each image consecutively as the first read data from *.raw, *.dat or *.cube files. This specific data organization and specificity of hyperspectral imaging were further used in the dedicated methods of hyperspectral image segmentation presented in the following sections.
Materials
The presented segmentation methods are related to the images acquired from different hyperspectral cameras (SOC710 Hyperspectral Imaging System, Specim sCMOS50V10E). In total, more than 10,000 2D images were acquired, constituting a sequence of more than 700 registrations, for 20 healthy subjects (40 % of women) aged 20–55 years. The test area was the area of the left and right hand as well as the forearm. The subjects expressed their free and informed consent for the study. The images were acquired in accordance with the Helsinki Declaration. Data were obtained retrospectively and no measurements or tests were carried out on the subjects as part of this work. This work only describes new techniques for image analysis. The discussed methods were tested on hyperspectral images from public databases and images described in the authors’ earlier works—for example [5]. The processing methods presented in this article were tested on a computer with Intel^{®} Core i7 4960X CPU 3.6 GHz. The resolution of acquired images L _{ GRAY }(m, n, i), where mrow, ncolumn, iframe number, was standardized to M × N × I = 696 × 520 × 128 pixels (where Mnumber of rows, Nnumber of columns, Inumber of bands). The spectral range for 128 bands was from 0.4 to 1.0 µm, dynamic range 12 bit, line rate 33 lines/s, pixels per line 520. The test subjects were illuminated by means of halogen lamps with a power of 100 W and linear radiation in the range of the camera operation.
Methods
The new dedicated methods of hyperspectral image segmentation were divided into 3 different areas (methods):

fast analysis of emissivity curves (absorption) of an object,

3D segmentation,

hierarchical segmentation.
The methods were preceded by image preprocessing.
Image preprocessing
These methods were preceded by image preprocessing. It involved median filtering of the image L _{ GRAY }(m, n, i) with a mask sized M _{ h } × N _{ h } × I _{ h } = 3 × 3 × 3 pixels. The size of the filter mask was selected based on the maximum size of a single artefact which was not higher than 4 pixels. The resulting image L _{ MED }(m, n, i) was further calibrated on the basis of the calibration bar which, at each registration, was located in the upper part of the image [5], or on the basis of the recorded images [dark L _{ DARK }(m, n, i) and white L _{ WHITE }(m, n, i)] allowing for conventional calibration and normalization of the image. The resulting image after this preprocessing of the image L _{ C }(m, n, i) was the basis for further transformations.
It should be emphasized here that median filtering for the fixed mask size is one of the simplest filtration methods. In a more developed form, median filtering should be adaptive. The necessity of applying the adaptive method is caused by the imaging character of hyperspectral cameras. For each camera, 2D images acquired for threshold wavelength values have the biggest noise [21]. Therefore, one of the possibilities is the application of an adaptive filtration method based on enlarging the mask size of the median filter for 2D images at the band edge. As it was shown in [21], it is the mask size changed depending on the wavelength in the range from M _{ h } × N _{ h } = 7 × 7 pixels to M _{ h } × N _{ h } = 3 × 3 pixels. In extreme cases when there is no information which 2D image will be analysed, the following mask size can be adopted: M _{ h } × N _{ h } = 7 × 7 pixels, taking into account the necessary removal of some minor details for median 2D images (for median wavelengths) [22–26].
Fast analysis of emissivity curves
According to the authors, the analysis of emissivity curves should be the most commonly used technique of hyperspectral image segmentation. Due to the specific data organisation, it can be used even when reading the beginning of the data in the *.raw, *.dat and *.cube file. Therefore, it can be used for the rough, screening test of compliance of spectrum amplitudes with the model (the expected waveform). The practical realization of such segmentation relates to the analysis of the first row of the matrix L _{ GRAY }(m = 1, n, i) or, after image preprocessing, the matrix L _{ C }(m, n, i). Assuming that the reference waveform is L _{ PAT }(i), the analysis for the first column n = 1 and the subsequent ones was formulated as:
For the adopted tolerance threshold p _{ r } of the difference between the reference and measured emissivity, the number of individual wavelengths λ (number i) for which the values L _{ D }(m, n, i) exceeded the adopted threshold p _{ r } were calculated—the obtained results constitute a new matrix L _{ R }(m, n), i.e.:
where:
Depending on the needs (speed of calculations), this analysis can be narrowed down to the aforementioned wavelength values i at m = 1 and i = 1. Table 1 shows the summary of mean analysis times for different values of m and i for Intel^{®} Core i7 4960X CPU 3.6 GHz. As shown in Table 1, the mean analysis time depends largely on the value of m. It is the result of time needed to read the information from the data files, *.raw, *.dat or *.cube (see Fig. 1a). The highest values are obtained at the first reading, and they are equal to 2.3 ms. This increased time results from getting access to data on the disk (in the tested operating system Windows 7 Professional). The time of reading other data for the next values of m and i increases linearly and for i∈(1, 100) and m = 100, it is 1.6 ms. For comparison, the time of reading any image L _{ GRAY }(m, n, i) for i = 1 or i = 100 is a minimum of 10 ms.
The sum of the difference values between the reference and measured waveforms L _{ R }(m, n) is the basis for changing the colour space and, thus, segmentation. The image L _{ V }(m, n), as defined below, is reliable for the grey levels:
In this case, interpretation of the results obtained is intuitive. The value equal to “0” is the maximum value of the error L _{ V }(m, n), the values close to “1” are the minimum error values. A sample image L _{ V }(m, n) is shown in Fig. 2. Segmentation of an object (objects) present in a scene can be performed based on the matrix L _{ V }(m, n) or L _{ R }(m, n). Sample segmentation results for various binarization thresholds (values 0.2, 0.4, 0.6 and 0.8) were presented in Fig. 3. It should be underlined that the obtained result is for the segmentation based on the first row of the matrix of 2D images for specific wavelengths, while the time of segmentation for the processor Intel^{®} Core i7 4960X CPU 3.6 GHz did not exceed 2.3 ms.
3D segmentation
A different approach to segmentation is the analysis of the whole sequence of i images L _{ C }(m, n, i). For this purpose, the ROI containing a fragment of the segmented object is marked either manually or automatically in the image L _{ C }(m, n, i). On this basis, the waveform L _{ PAT }(i) is created as:
where: M _{ R } and N _{ R }—number of rows and columns of the ROI respectively. The image L _{ V }(m, n) is calculated in the next step:
The image L _{ V }(m, n) contains a lot of artefacts resulting from the noise occurring in the border images L _{ V }(m, n, i) for i ∈ {1, 2, I − 1, I}. Therefore, 3D filtration was proposed, which enabled median filtering of each pixel for the 8neighbourhood system. Schematic diagram of the filtration process for a threedimensional plane of an object (border) is shown in Fig. 4a. The process of filtration relates to the 8neighbourhood system of each pixel of the matrix L _{ v }(m, n, i) on the basis of which the median value is calculated. A median filter is used in relation to the plane formed from the object—in this case sized 30 × 30 pixels (see Fig. 4a). On this basis, the filtered image Ls(m, n), after binarization with the thresholds p _{ r1 } , p _{ r2 } and p _{ r3 } equal to 30, 20 and 10 % respectively (thresholding with the upper threshold), is shown using colours (red, green and blue respectively) in Fig. 4b. This image, L _{ S }(m, n), resulting from binarization with the properly selected (manually or automatically) threshold, is the result of segmentation.
As it can be concluded from the above description, due to the specificity of this segmentation method, the time of analysis is at the level of 1949 ms for the processor Intel^{®} Core i7 4960X CPU 3.6 GHz.
Hierarchical segmentation
The dedicated hierarchical method gives the best results in segmentation. The word hierarchical refers here to the hierarchy in the resolution analysis of the image L _{ C }(m, n, i). The analysis is associated with the analysis of the resolution of the image L _{ C }(m, n, i) reduced to M/8 × N/8. Then, the analysis of the knearest neighbours was proposed. The knearest neighbour analysis was implemented for three features, w(1), w(2) and w(3). These features are absolute values of differences between the analysed image L _{ C }(m, n, i) and the reference waveform L _{ PAT }(i) for the selected three wavelengths. The three selected wavelengths are, in this case, the local maxima of the waveform L _{ PAT }(i). In general, the number of features [local maxima of the waveform L _{ PAT }(i)] is arbitrary. It can be limited only by computational complexity (the need to determine individual local maxima). In this example, due to the facilitated visualization (3D graph), the number of sought local maxima was limited to three [three features w(1), w(2) and w(3)]. Therefore, the segmentation method starts with classification in 2, 3, 4 and 5 regions, taking into account the features w(1), w(2) and w(3)—Fig. 5.
As the number of regions, the number of classes, is not generally known (the number of the segmented objects), the value of standard deviation of the mean (std) is used as a criterion. The number of classes is increased from 2 to 5 regions. The correct number of classes is the one for which the value of the interclass mean square error (std) is the smallest. Figure 6 shows examples of segmentation results for the thumb (Fig. 2a). The results were obtained for two, three and four classes and three features [w(1), w(2) and w(3)].
For the sample results shown in Fig. 6, the smallest value of the mean square error (std = 0.21) was for 3 classes. The image L _{ HC }(m, n, i) created in this way is the basis for increasing its resolution and automatic designation of segmented areas in the source image L _{ GRAY }(m, n, i) or in the image after preprocessing of the images L _{ C }(m, n, i) (Fig. 5).
Hyperspectral images segmentation [27–29] with the hierarchical method is linked with yet another new question concerning the variability of size/shape of the object in the segmentation for next ith values taking into consideration the criterion described with the formula (4). The segmented object may therefore change its shape fulfilling the criterion (4) for the next ith images. A simple method of binarization conducted for the next ith images does not render satisfying results due to the possibility of adding next new objects (not by changing the size of the existing segmented object). Therefore, erosion L _{ E }(m, n, i) was proposed with conditional dilatation L _{ D }(m, n, i) performed for each ith binary image starting from L _{ BIN }(m, n, i = 1), i.e.:
where:
and m _{ S }, n _{ S } the row and column position of the mask SE sized M _{ S } × N _{ S } depending on the size of the segmented object. The mask changed in this case has a size of M _{ S } × N _{ S } = 17 × 17 pixels. The values of thresholds p _{ ec } and p _{ dc } are set by the user once for all images from the given type of camera. In the present case of SOC710 Hyperspectral Imaging System the value was p _{ ec } = p _{ dc } = 0.5. The results for the segmentation without erosion and conditional dilatation and with erosion and conditional dilatation are shown in Fig. 7.
Experimental and discussion
The three new proposed segmentation methods need to be tested in practice in order to confirm the accuracy and time of analysis. The basis for the performed experiments is:

The three methods described in this article:

fast analysis of emissivity curves (SKE),

3D segmentation (S3D),

hierarchical segmentation (SH),


The three known segmentation methods (for quantitative comparison of the results obtained):

method based on brightness thresholding (SPJ)—the binarization threshold is selected manually and automatically using Otsu’s formulas [18],

watershed method (SWS) preceded by filtration with an averaging filter whose mask size is in the range from 3 × 3 pixels to 9 × 9 pixels,

method based on mathematical morphology (SMM), especially erosion and conditional dilation,


Manual method of object selection considered further as the benchmark (SP).
Over 700 complete sequences of images L _{ GRAY }(m, n, i) (a total of more than 10,000 2D images), whose acquisition conditions are given in the “Materials” section, were subject to segmentation. Common measures were used to compare the quality of the detected objects—FP the number of falsepositive detected pixels, FN—false negative as well as TP, TN—true negative and true positive respectively. On their basis, sensitivity TPR = TP/(TP + FN), specificity SPC = TN/(TN + FP) and balanced accuracy ACC = (TPR + SPC)/2 were determined. The obtained results of mean values of TPR, SPC and ACC for the compared segmentation methods are presented in Table 2.
According to the presented Table 2, segmentation based on a hierarchical approach, proposed in this article, has the greatest value of ACC (ACC = 92). The other proposed segmentation methods are at a similarly high level (of ACC values). Among the known nonprofiled segmentation methods, the method based on brightness thresholding has the largest value of ACC. However, it is a semiautomatic method whose results are closely dependent on the method of selecting the binarization threshold. The other features of the discussed segmentation methods look slightly different—Table 3.
The results for a PC with Intel^{®} Core i7 4960X CPU 3.6 GHz clearly show the superiority of the method of 3D segmentation (S3D) whose analysis time does not depend on the number of detected (segmented) objects. Moreover, this method is not sensitive to image rotation. This is due to the idea of its operation shown in Fig. 4 where no direction is privileged. Of course, fast analysis of emissivity curves is the quickest—its analysis time for m = 1 and i = 1 is equal to 2.3 ms (Table 1). The other known segmentation methods consume several times more of the CPU time.
As it was mentioned in the “Background” section, numerous similar segmentation methods are known from the literature. Hence, the method utilizing emissivity curves of a given object (SKE method) is used in many applications such as specim spectral imaging [30, 31], hyperspectral imaging system [5, 18] or imageJ [32]. However, in none of such applications the time of preliminary segmentation was obtained taking into account the emissivity curve at the level of 2.3 ms. For example in [33] the authors obtain the time of analysis equal to 100 ms using Intel^{®} Core i5 CPU M460 @2.5 GHz 4 GB RAM. Usually, the time of analysis increases, e.g. in Specim systems, which results from the fact that the registered sequence of ith 2D images is preprocessed. A significant element of this processing is lowering the resolution of images. Then, this initiallyprepared 3D image is analyzed. In the case when many different objects are registered, the process becomes burdensome unnecessarily increasing the involvement of staff (operator). In the case of the second and third method in question (S3D and SH), the situation is different. These methods, especially SH based on erosion and conditional dilatation, is not used in any commercial applications concerning hyperspectral imaging and it was not presented in this use, namely for correction of hyperspectral images. The majority of authors, for instance in the articles [4, 6, 8, 11, 13], perform segmentation in the field of an image, one selected 2D image from ith images of a sequence. Another widely tested approach is the use of clustering in the space of particular wavelengths e.g. at work [6]. Both mentioned types of segmentation do not cope with two elements: variable levels of noise of hyperspectral images for band extreme values and with the correction of the shape of the segmented image for subsequent ith images. In particular, the correction of the object shape for subsequent ith 2D images of the series is especially important in terms of determining the compatibility of its emissivity curves with the calibration curve. Therefore, the methods presented in this paper are new segmentation methods which provide a new quality in a considerably shorter period of time when compared to the known methods of hyperspectral image segmentation.
Conclusion
The presented three new segmentation methods have the following advantages:

they are fully automatic,

they enable to use a fast segmentation method

they are profiled to hyperspectral image segmentation,

they use emissivity curves as the model,

they can be applied to any type of objects (not necessarily biological),

they are faster and more precise compared to conventional methods known from the literature.
Work is currently underway to implement the described segmentation methods, especially the first described fast method for digital circuits DSP. The aim is the initial, at the acquisition phase, object segmentation. This will enable to automatically limit the acquired image and, thus, reduce the hyperspectral camera operation time.
Abbreviations
 ROI:

region of interest
 SKE:

fast analysis of emissivity curves
 S3D:

segmentation 3D
 SH:

hierarchical segmentation
 SPJ:

method based on brightness thresholding
 SWS:

watershed method
 SMM:

method based on mathematical morphology
 SP:

manual method
References
Pang B, Zhang D, Wang K. The bielliptical deformable contour and its application to automated tongue segmentation in Chinese medicine. IEEE Trans Med Imaging. 2005;24(8):946–56.
Morgan P, Frankish C. Image quality, compression and segmentation in medicine. Audiov Media Med. 2002;25(4):149–54.
Suetens P, Bellon E, Vandermeulen D, Smet M, Marchal G, Nuyts J, Mortelmans L. Image segmentation: methods and applications in diagnostic radiology and nuclear medicine. Eur J Radiol. 1993;17(1):14–21.
Piqueras S, Krafft C, Beleites C, Egodage K, Eggeling F, GuntinasLichius O, Popp J, Tauler R, Juan A. Combining multiset resolution and segmentation for hyperspectral image analysis of biological tissues. Anal Chim Acta. 2015;881:24–36.
Koprowski R, Wilczyński S, Wróbel Z, BłońskaFajfrowska B. Calibration and segmentation of skin areas in hyperspectral imaging for the needs of dermatology. Biomed Eng Online. 2014;13:113.
Veganzones MA, Tochon G, DallaMura M, Plaza AJ, Chanussot J. Hyperspectral image segmentation using a new spectral unmixingbased binary partition tree representation. IEEE Trans Image Process. 2014;23(8):3574–89.
Porwik P. Efficient spectral method of identification of linear Boolean function. Control Cybern. 2004;33(4):663–78.
Fu D, Xie XS. Reliable cell segmentation based on spectral phasor analysis of hyperspectral stimulated Raman scattering imaging data. Anal Chem. 2014;86(9):4115–9.
Hennessy R, Bish S, Tunnell JW, Markey MK. Segmentation of diffuse reflectance hyperspectral datasets with noise for detection of melanoma. Conf Proc IEEE Eng Med Biol Soc. 2012;2012:1482–5.
Eches O, Benediktsson JA, Dobigeon N, Tourneret JY. Adaptive Markov random fields for joint unmixing and segmentation of hyperspectral images. IEEE Trans Image Process. 2013;22(1):5–16.
Piqueras S, Duponchel L, Tauler R, Juan A. Resolution and segmentation of hyperspectral biomedical images by multivariate curve resolutionalternating least squares. Anal Chim Acta. 2011;705(1–2):182–92.
Zhang Q, Plemmons R, Kittle D, Brady D, Prasad S. Joint segmentation and reconstruction of hyperspectral data with compressed measurements. Appl Opt. 2011;50(22):4417–35.
Tarabalka Y, Chanussot J, Benediktsson JA. Segmentation and classification of hyperspectral images using minimum spanning forest grown from automatically selected markers. IEEE Trans Syst Man Cybern B Cybern. 2010;40(5):1267–79.
Liu Z, Yan JQ, Zhang D, Li QL. Automated tongue segmentation in hyperspectral images for medicine. Appl Opt. 2007;46(34):8328–34.
Christensen MP, Euliss GW, McFadden MJ, Coyle KM, Milojkovic P, Haney MW, Gracht J, Athale RA. ACTIVEEYES: an adaptive pixelbypixel imagesegmentation sensor architecture for highdynamicrange hyperspectral imaging. Appl Opt. 2002;41(29):6093–103.
Gao L, Smith RT. Optical hyperspectral imaging in microscopy and spectroscopy—a review of data acquisition. J Biophotonics. 2015;8(6):441–56.
Foster KR, Koprowski R, Skufca JD. Machine learning, medical diagnosis, and biomedical engineering research—commentary. BioMed Eng OnLine. 2014;13:94.
Koprowski R. Hyperspectral imaging in medicine: image preprocessing problems and solutions in Matlab. J Biophotonics. 2015;8(11–12):935–43.
Mitchell T. Machine learning. New York: McGraw Hill; 1997. p. 414.
Krzanowski WJ. Principles of multivariate analysis, a user’s perspective. New York: Oxford University Press; 1988. p. 608.
Koprowski R. Processing hyperspectral medical images. Berlin: Springer; 2017. p. 140.
Lefèvre S, Aptoula E, Perret B, Weber J. Morphological template matching in color images advances in lowlevel color image processing. Berlin: Springer; 2013.
Galeano J, Jolivot R, Marzani F. Analysis of human skin hyperspectral images by nonnegative matrix factorization advances. Soft Comput. 2011;7095:431–42.
Jia H, Ding S, Meng L, Fan S. A densityadaptive affinity propagation clustering algorithm based on spectral dimension reduction Neural Computing and Applications. 2014;25:1557–67.
Carlinet E, Géraud T. MToS: a Tree of Shapes for Multivariate Images. IEEE Trans Image Process. 2015;24(12):5330–42.
Halimi A, Dobigeon N, Tourneret JY. Unsupervised unmixing of hyperspectral images accounting for endmember variability. IEEE Trans Image Process. 2015;24(12):4904–17.
Grana M, Chyzhyk D. Image understanding applications of lattice autoassociative memories. IEEE Trans Neural Netw Learn Syst. 2015.
Piqueras S, Krafft C, Beleites C, Egodage K, von Eggeling F, GuntinasLichius O, Popp J, Tauler R, de Juan A. Combining multiset resolution and segmentation for hyperspectral image analysis of biological tissues. Anal Chim Acta. 2015;881:24–36.
Banas K, Banas A, Gajda M, Pawlicki B, Kwiatek WM, Breese MB. Preprocessing of Fourier transform infrared spectra by means of multivariate analysis implemented in the R environment. Analyst. 2015;140(8):2810–4.
Lin Y, Puttonen E, Hyyppä J. Investigation of tree spectral reflectance characteristics using a mobile terrestrial line spectrometer and laser scanner. Sensors (Basel). 2013;13(7):9305–20.
Serranti S, Cesare D, Marini F, Bonifazi G. Classification of oat and groat kernels using NIR hyperspectral imaging. Talanta. 2013;103:276–84.
Pasqualin C, Gannier F, Yu A, Malécot CO, Bredeloux P, Maupoil V. SarcOptiM for ImageJ: high frequency online sarcomere length computing on stimulated cardiomyocytes. Am J Physiol Cell Physiol. 2016;311:C277–83.
Koprowski R, Wilczyński S, Wróbel Z, Kasperczyk S, BłońskaFajfrowska B. Automatic method for the dermatological diagnosis of selected hand skin features in hyperspectral imaging. Biomed Eng Online. 2014;22(13):47.
Authors’ contributions
RK suggested the algorithm for image analysis and processing, implemented it and analyzed the images. PO performed the measurements and assessment, from the medical perspective, of the usefulness of the segmentation method in question in medical practice. Both authors read and approved the final manuscript.
Acknowledgements
The Authors would like to thank Mr. Raphael Stachiewicz from ENFORMATIC Sp. z. o. o., Poland, for providing the hyperspectral camera. The Authors wish to sincerely thank Doctor Sławomir Wilczyński from Medical University of Silesia in Katowice in Poland for the possibility of using the collected hyperspectral images.
Competing interests
The authors declare that they have no competing interests.
Funding
This work was financially supported by Medical University of Silesia in Katowice, Grant number KNW1023/N/6/O.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Koprowski, R., Olczyk, P. Segmentation in dermatological hyperspectral images: dedicated methods. BioMed Eng OnLine 15, 97 (2016). https://doi.org/10.1186/s1293801602195
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1293801602195
Keywords
 Fast segmentation method
 Image processing
 Thresholding
 Conditional erosion
 Conditional dilatation
 Hyperspectral imaging