Next Article in Journal
Construction of a Ginseng Root-Meristem Sensor and a Sensing Kinetics Study on the Main Nitrogen Nutrients
Next Article in Special Issue
Supervised Machine Learning Methods and Hyperspectral Imaging Techniques Jointly Applied for Brain Cancer Classification
Previous Article in Journal
Muscle Co-Activation around the Knee during Different Walking Speeds in Healthy Females
Previous Article in Special Issue
Hyperspectral Imagery for Assessing Laser-Induced Thermal State Change in Liver
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Curve-Based Classification Approach for Hyperspectral Dermatologic Data Processing

1
Department of Education and Pedagogy, UiT the Arctic University of Norway, 9019 Tromsø, Norway
2
Institute for Applied Microelectronics, Universidad de Las Palmas de Gran Canaria, 35016 Las Palmas de Gran Canaria, Spain
3
Department of Dermatology, Hospital Universitario de Gran Canaria Doctor Negrín, 35016 Las Palmas de Gran Canaria, Spain
4
Department of Dermatology, Complejo Hospitalario Universitario Insular-Materno Infantil, 35016 Las Palmas de Gran Canaria, Spain
5
Department of Electromedicine, Complejo Hospitalario Universitario Insular-Materno Infantil, 35016 Las Palmas de Gran Canaria, Spain
6
Department of Mathematics and Statistics, UiT the Arctic University of Norway, 9019 Tromsø, Norway
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(3), 680; https://doi.org/10.3390/s21030680
Submission received: 15 December 2020 / Revised: 13 January 2021 / Accepted: 15 January 2021 / Published: 20 January 2021
(This article belongs to the Special Issue Trends and Prospects in Medical Hyperspectral Imagery)

Abstract

:
This paper shows new contributions in the detection of skin cancer, where we present the use of a customized hyperspectral system that captures images in the spectral range from 450 to 950 nm. By choosing a 7 × 7 sub-image of each channel in the hyperspectral image (HSI) and then taking the mean and standard deviation of these sub-images, we were able to make fits of the resulting curves. These fitted curves had certain characteristics, which then served as a basis of classification. The most distinct fit was for the melanoma pigmented skin lesions (PSLs), which is also the most aggressive malignant cancer. Furthermore, we were able to classify the other PSLs in malignant and benign classes. This gives us a rather complete classification method for PSLs with a novel perspective of the classification procedure by exploiting the variability of each channel in the HSI.

1. Introduction

Hyperspectral (HS) imaging (HSI) combines conventional imaging and spectroscopy methods in a single imaging technique providing both spatial and spectral information of the captured display [1,2]. It is also a fitting method for medical applications due to its non-invasive, non-ionizing, and label-free nature [3]. Dermatology is one of the medical fields where HSI have shown its potential as an appropriate imaging technique [4].
Currently, the most common form of cancer, with more than 1.3 million new cases worldwide in 2018, is skin cancer [5]. There are several types of skin lesions. Pigmented skin lesions (PSLs) contain a wide variety of types, including cancerous and non-cancerous PSLs [6]. There are two types of PSLs depending on the type of growth of the tissue: benign and malignant. Nevus, which is benign, has a slow growth rate and is noncancerous, while e.g., melanomas, which is malignant, are invasive and potentially metastatic tumors [6]. Other types of skin cancer produced by different types of cells include squamous cell carcinoma and basal cell carcinoma. However, melanomas are much more dangerous than other types of skin cancer, and an early detection of this skin lesion can be extremely important in improving the patient survival [7]. The diagnosis of PSLs is performed by dermatologists employing their naked eye or dermascopic cameras, which enhance the morphological visualization of the PSL. After the inspection, an analysis of the lesion following the ABCDE (Asymmetry, Border irregularity, Color, Diameter and Evolving size, shape, or color) rule is applied in order to establish a preliminary diagnosis [7]. A suspicious lesion will accordingly get a histopathological analysis carried out to determine the definitive diagnosis. These diagnostic tools based on conventional imaging have been used to help dermatologists in preliminary diagnosis. However, conventional imaging has limitations that could be surpassed by the enriched spectral information provided by the HSIs.
HSI technology has already been explored as a target technology to aid in PSLs diagnosis. A recent study by Leon et al. showed promising results in discriminating PSLs by random forests and artificial neural networks [8]. Another recent study by Hosking et al. used digital dermascopy HSI with machine learning methods in order to detect melanomas [9]. Kazianka et al. used HSIs to differentiate between normal skin, melanomas, and moles [10]. They used unsupervised approaches to segment the images and moreover evaluated the performance of several supervised classifiers aimed to retrieve the diagnosis of each pixel. The results indicated that it was possible to differentiate melanomas from moles with high specificity and sensitivity. In another study, Nagaoka et al. [11] collected a dataset for discriminating between melanomas and other PSLs, evaluating the statistical significance of an index defined to this end. Tomatis et al. used a classifier based on a neural network model [12]. There have also been developed some commercial systems as SIAscope/SIAscopy [13] and MelaFind [14,15,16]. MelaFind is being used in several studies. In [15], Elbaum et al. conduct a leave-one-out cross-validation procedure using a database composed of 183 melanocytic nevus and 63 melanomas. Monheit et al. [16] applied it in a multicenter study. At last, Song et al. performed a paired comparison between MelaFind and a reflectance confocal microscopy system to differentiate between melanoma and non-melanoma PSLs by comparing MelaFind with a confocal microscopy system. In the study by Stamnes et al., they discriminate between malignant and benign PSLs, i.e., no specific melanoma detection [17].
The current default method of data analysis (DA) is deep learning (DL). DL emerged as a competing DA method around 2009 and solves a wide range of problems exemplified by the versatile TensorFlow-package [18]. DL is also applied in the HSI field. However, DL often needs a huge amount of data in order to train the parameters enough in the neural network [19]. Thus, in the proposed method, this presents a difficult problem due to the lack of data. We are pursuing this problem by an ongoing data collection project, aiming for future applications of DL systems.
Most of the research in the use of HSI for skin analysis is focused on the automatic diagnosis of PSLs. Clearly, a successful methodology of this type would be an extremely useful decision support tool for general practitioners (GPs) and dermatologists. In our contribution in this direction, we were able to find clear patterns for the melanoma PSLs, the other malignant PSLs and the benign PSLs. In the novel exploitation of the variability of each channel in the HSI cubes, we made curves of the mean and standard deviation of sub-images in these channels. Furthermore, using these curves and their inherent patterns, we were able to discriminate between malignant and benign PSLs, but also between the aggressive melanoma PSL and the others.
The outline of the paper is as follows. First, we describe the developed dermatological HSI acquisition system, and thereafter, we describe the datasets we are using in this curve-based classification method. Next, we point out how the data are pre-processed and labeled before they are used as input in our new classification method, which consists of a training, validation, and testing phase.

2. Materials and Methods

2.1. Hyperspectral In Vivo Dermatologic Data

The system used to acquire the HS dermatologic images is described in [20], using the same database as in Leon et al. [8]. The database consists of 76 images of PSLs from 61 patients: 36 cancerous and 40 noncancerous. The dataset is divided into a training set, a test set, and a validation set. The training set is subsequently divided into melanoma PSLs, malignant PSLs, and benign PSLs. Each mentioned set is mutually exclusive. The acquisition system is composed of an HS snapshot camera able to capture HS images in the very near infrared (VNIR) range, between 450 and 950 nm, with a spatial resolution of 50 × 50 pixels and 125 spectral bands. This acquisition system employs a customized dermascopic contact structure and a halogen source light (150 W) coupled to a fiber optic ring light guide for cold light emission. The effective capturing area of the system is 12 × 12 mm, and the acquisition time is lower than 1 s. The system was applied to create an HS database composed of 76 images of PSLs from 61 patients. The data acquisition was done at the Hospital Universitario de Gran Canaria Doctor Negrín and Complejo Hospitalario Universitario Insular-Materno Infantil (Spain). The research protocol and consent procedures were approved by the Comité Ético de Investigación Clínica-Comité de Ética en la Investigación (CEIC/CEI) at the University Hospital Doctor Negrin.

2.2. Curve-Based Classification Approach

In the pre-processing step, a calibration is done of the HSI using the white and dark reference images following Equation (1), where CI is the calibrated image, RI is the raw image obtained from the HS camera, and DI and WI are the dark and white reference images, respectively. The white reference image is obtained capturing an image of a white reference tile that reflects 99% of the incident light. The dark reference image is captured by keeping the tap in the lens of the camera.
C I = R I D I W I D I
Finally, the data is normalized in the range [0,1] to avoid differences in the spectral signature intensities caused by possible different illumination conditions.

2.3. Region of Interest Curves

The region of interest (ROI) was obtained through the extraction of 7 × 7 sub-images for all the channels. Subsequently, the mean (x-coordinate of the curves) and standard deviation (y-coordinate of the curves) of these ROIs from each HS channel were plotted against each other, thus creating ROI curves for each HSI cube. These curves were rather different, some almost linear, other more convoluted. Thus, we applied a morphological dilation filter to the more complex curves to ease the interpretation and subsequently use these filtered curves in function fits. Morphological dilation is a common way to enhance or change digital images in some way. Thus, often used in bridging gaps, our morphological filter can be seen as an inverse dilation, a pruning.
There were some issues of this method that needed to be addressed. In particular, how the ROIs are chosen. So, the first issue is where to collect this ROI in the HSI, i.e., which coordinates to choose? To shed some light on this, we established some criteria for each of the classes of PSLs and further picked the ROIs best suited for creating the curves with adequate classification properties. In each channel, we started from (x,y) = (16, 16) and went to (x,y) = (35, 35). This gave us 20 × 20 = 400 HS ROI curves for each PSL HSI cube collected from different locations. However, with some overlap and since each ROI area is 7 × 7 pixels, we cover the HSI from the start at (x,y) = (13, 13) to the end at (x,y) = (38, 38): thus, an area of 26 × 26 = 676 pixels. We tried a wider area, but the results did not improve, so we kept this size due to computational cost. Thus, we tested how the ROI extraction affected the mean square error (MSE) of the curves chosen by suitable criteria for the melanoma as well as the other malignant and benign PSLs.
The CPU used to run the code was an Intel Xeon 2.8 GHz running on a 16 GB RAM computer.

Curve-Fitting

By trying to classify the PSLs through fits of the ROI curves, all the HSI ROI curves were fitted by a quartic (fourth degree) polynomial function, where the coefficients are found through the nlinfit function of MATLAB®, which minimizes the weighted least squares equation,
i = 1 n w i ( y i f ( x i , b ) ) 2
where the weights were chosen through a weight function, w(y) = 1/(0.011 + 0.011y)2 and n = 125, the number of data points. Further is f(x, b) = b0 x4 + b1 x3 + b2 x2 + b3 x + b4, the quartic regression function, with parameters bj, j = 0,…,4. Then, the nlinfit-function uses the iteratively reweighted least-square algorithm, which is designed to deal with outliers [21]. Furthermore, the use of a quartic polynomial with a weight function proved to be the best option with regard to flexibility and robustness to outliers. We also tested higher degree polynomial functions; however, they were prone to overfitting.
Figure 1 gives an overview of the algorithm. The constants a, b, c, d, e, f, and g are found though the training phase. The features are as follows: df (first-order derivative), ddf (second-order derivative), totm = a·df + b·ddf, mabdf = mean + abs(df) and mean.

3. Results

Here, we will represent the obtained experimental results. We start the representation of the training phase.

3.1. The Training Phase of the Method

As mentioned above, we acquired an ROI from each channel of the HSIs, where the curves were formed through the calculation of the mean and standard deviation of the ROIs from each channel. First, we investigated the melanoma PSLs. The features coming most natural from earlier experiments are a high mean first-order derivative (Max df) and a high mean curvature given by the mean of the second derivative (Max ddf). However, during the training phase, we found that a high sum of mean df (multiplied by 1, a in Figure 1), mean ddf (multiplied by 0.1, b in Figure 1), and the derivative at the maximum value gave the most successful feature (Max totm > 2.86, c in Figure 1). In Table 1 and Figure 2a, we can observe that this feature produced rather good fits according to the MSEs: the MSE ddf following this measure and the MSE df lower for two melanomas. Thus, all the plots of the HSI ROI curves from the melanoma PSLs are from the fits of the Max totm curves. Moreover, the recurring theme in Figure 2c–f is that all the HSI ROIs came from the interior and/or near the edges of the melanomas. In Table 1, we can also observe the most distinct outlier, the P82C1000 PSL, which did not have the same characteristics as the other melanomas, behaving more similar to the malignant PSLs. In Figure 2a, we can observe that this PSL is the one with an MSE ddf diverging most clearly from the MSE totm.
Regarding the malignant PSLs, we inferred based on earlier experiments that the HSI ROI curves with a high mean (Max Mean > 0.05, d in Figure 1) or a high product of max x and max y (Max Prod > 0.395, e in Figure 1) from the fitted curve are plausibly originating from this type of PSL. This proved successful for some of the HSI ROI curve fits from the malignant PSLs. As seen in Table 2 and Figure 3a, this was the most difficult type PSL to fit, giving rise to rather high MSEs. The Max Mean was the best feature accordingly. Some of the curves with least MSE Mean are plotted in Figure 3b, where all have means lying above the standard deviation (y-value) of 0.05 (d in Figure 1). From these HSI ROI mean standard deviation curves, the most complex data occurred—thus, as mentioned above, giving rise to rather high MSEs. The corresponding ROIs of the plotted malignant curves are also here lying on the inside or at the edge of PSLs, as seen in Figure 3c–f.
The benign PSLs HSI ROI curves proved easier to fit than the malignant ones. From earlier experiments, two useful features came to be low mean (Min Mean < 0.05, g in Figure 1), low absolute value first order derivative (Min abs df), and the sum of these (Min mabdf < 0.109, f in Figure 1). This proved rather successful with a few exceptions, as can be observed in Table 3 and Figure 4a, with Min abs df and Min mabdf as similar rather successful features. In Figure 4b, four of the curves with lowest MSE abs df are plotted with their corresponding placements of the ROIs given in Figure 4c–f. We can here see that for these PSLs, the ROIs are somewhat different placed than for the melanomas and malignant ones—not all inside the PSLs, but at some sort of feature, nevertheless.

Classification

From the training phase, we inferred that the classification procedure can consist of two consecutive steps.
First, we identify the curves where the first derivative increases, possibly converging at maximum value of the HSI ROI curves mean. Additionally, the curves with high mean positive curvature and less double inflexion points were also identified as melanomas. In addition, the sum df + 0.1ddf proved successful (totm).
Then, we find the two clusters of the rest of the curves, such that they can be classified as benign or malignant. This is mainly done by choosing what curves have mean above y = 0.05 (malignant) and below y = 0.05 (benign). For the benign curves, an additional useful feature was the sum of the mean of the fitted curve and the mean absolute value of the first-order derivative in each point (mabdf).

3.2. Validation and Test Experiment

We also performed an experiment in order to validate the training phase with the same type of ROI curves as in the training part but different PSLs. This resulted in the usual 125 points mean and standard deviation curve; then, a quartic polynomial function was fitted, as mentioned earlier. Based on the features developed in the training phase, then, we were able to first discriminate between melanomas and the rest of the PSLs, where a typical melanoma fit is given in Figure 5a, and then discriminate between malignant and benign PSLs. Examples of their respective fits given in Figure 5b,c.
When looking at the results, we denote correctly classifying a malignant PSL as a true positive (tp) and correctly classifying a benign PSL as a true negative (tn). Furthermore, falsely classifying a malignant PSL is denoted as a false positive (fp) and at last, falsely classifying a benign PSL is denoted as a false negative (fn). By using the validated classification criteria, we were able to test the method with a test set. The results are given in Table 4, where the P23C1001 gave a false positive due to some regression difficulties regarding outliers.
The accuracy of this proposed method is 90%. The sensitivity is 100% and the specificity is 80%. Considering the melanoma discrimination, both sensitivity and specificity are 100%. In Table 5 results from other researchers are provided.
From Table 5, we can discern that most research studies have been discrimination of the aggressive melanoma PSLs, where Nagaoka et al. [11] and Leon et al. [8] had the highest sensitivity and specificity. Stamnes et al. [17] also had high such measures; however, these were only discriminating between benign and malignant PSLs. Thus, the sensitivity and specificity in the method developed here is comparable to the others.
The works listed in Table 5 mostly use an array of machine learning methods, while the proposed method only use a regression method, making it a simpler approach than the others.
The computation time was on average 0.16 s for each PSL. However, including the best ROI search, this jumps to 70 s. Thus, it is a rather slow algorithm due to rather complex calculations. Compared to the previous work in [8], in which six of the 10 patients were correctly identified with an average execution time of 0.5 s, the proposed methodology identifies nine out of 10 patients with an average of 160 s. However, this is only the core algorithm; an efficient segmentation algorithm has to be included also, raising the execution time.

4. Limitations and Future Directions

Further studies must be conducted in order to improve the limitations in this research. The first limitation is the rather low number of samples in particular for the melanoma class. In addition, in the benign and malignant class, it would have been desirable of an increase of samples. However, the number of samples was enough to find and validate the patterns for each class of PSLs considered here. Another limitation is the computing time. However, this will be decreased by the automatic choosing of ROIs by segmentation of the PSLs, which have several known potential solutions. We will also try to improve the regression of the curves by further experimentation of the weight function or other types of regression. The developed system in this paper is at an experimental level; however, it has the potential to assist in the diagnosis of skin cancer and reduce the number of biopsies relating to this.

5. Conclusions

The work presented here shows that HSIs can be used to discriminate between melanoma, malignant, and benign PSLs as a potential future diagnostic tool for dermatologists or GPs. The proposed curve-based method is able to do this discrimination by using the inherent patterns of these curves. The presented method was able to correctly classify nine out of 10 in the test set, identifying all cancerous cases accurately (100% of sensitivity). Thus, it is a rather good result. In addition, in the test phase, it had a sensitivity of 80% overall and 100% for the melanomas. Since this method is relatively experimental, we will improve it further in order to increase its accuracy and lower its computation time. The limitation of the proposed method is mainly the lack of data, which is remedied in an ongoing project of data collection. The merits of the proposed method is the utilization of the richer structure provided by the HSI technology in a novel way.

Author Contributions

Conceptualization, F.G. and G.M.C.; methodology, S.U.; software, S.U.; validation, S.U., formal analysis, S.U.; investigation, S.U.; resources, G.M.C.; data curation, I.C., G.C., A.G., P.A., J.A.H.; writing—original draft preparation, S.U.; writing—review and editing, E.Q.; visualization, S.U.; project administration, E.Q., F.G., G.M.C.; funding acquisition, F.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been supported by the Canary Islands Government through the ACIISI (Canarian Agency for Research, Innovation and the Information Society), ITHACA project “Hyperspectral Identification of Brain Tumors” under Grant Agreement ProID2017010164 and it has been partially supported also by the Spanish Government and European Union (FEDER funds) as part of support program in the context of Distributed HW/SW Platform for Intelligent Processing of Heterogeneous Sensor Data in Large Open Areas Surveillance Applications (PLATINO) project, under contract TEC2017-86722-C4-1-R.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to patient confidentiality.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hong, D.; Wu, X.; Ghamisi, P.; Chanussot, J.; Yokoya, N.; Zhu, X.X. Invariant Attribute Profiles: A Spatial-Frequency Joint Feature Extractor for Hyperspectral Image Classification. IEEE Trans. Geosci. Remote Sens. 2020, 58, 3791–3808. [Google Scholar] [CrossRef] [Green Version]
  2. Huang, H.; Luo, F.; Liu, L.; Yang, Y. Dimensionality reduction of hyperspectral images based on sparse discriminant manifold embedding. ISPRS J. Photogramm. Remote Sens. 2015, 106, 42–54. [Google Scholar] [CrossRef]
  3. Lu, G.; Fei, B. Medical hyperspectral imaging: A review. J. Biomed. Opt. 2014, 19, 010901. [Google Scholar] [CrossRef] [PubMed]
  4. Johansen, T.H.; Møllersen, K.; Ortega, S.; Fabelo, H.; Garcia, A.; Callico, G.M.; Godtliebsen, F. Recent advances in hyperspectral imaging for melanoma detection. Wiley Interdiscip. Rev. Comput. Stat. 2019, 12, e1465. [Google Scholar] [CrossRef]
  5. Bray, F.; Ferlay, J.; Soerjomataram, I.; Siegel, R.L.; Torre, L.A.; Jemal, A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [Green Version]
  6. LeBoit, P.E.; Burg, G.; Weedon, D.; Sarasin, A. Pathology & Genetics Skin Tumors; World Health Organization: Geneve, Switzerland, 2006. [Google Scholar]
  7. Tsao, H.; Olazagasti, J.M.; Cordoro, K.M.; Brewer, J.D.; Taylor, S.C.; Bordeaux, J.S.; Chren, M.-M.; Sober, A.J.; Tegeler, C.; Bhushan, R.; et al. Early detection of melanoma: Reviewing the ABCDEs American Academy of Dermatology Ad Hoc Task Force for the ABCDEs of Melanoma. J. Am. Acad. Dermatol. 2015, 72, 717–723. [Google Scholar] [CrossRef]
  8. Leon, R.; Martinez, B.; Fabelo, H.; Ortega, S.; Melian, V.; Castaño, I.; Carretero, G.; Almeida, P.; Garcia, A.; Quevedo, E.; et al. Non-Invasive Skin Cancer Diagnosis Using Hyperspectral Imaging for In-Situ Clinical Support. J. Clin. Med. 2020, 9, 1662. [Google Scholar] [CrossRef] [PubMed]
  9. Hosking, A.-M.; Coakley, B.J.; Chang, D.; Talebi-Liasi, F.; Lish, S.; Lee, S.W.; Zong, A.M.; Moore, I.; Browning, J.; Jacques, S.L.; et al. Hyperspectral imaging in automated digital dermoscopy screening for melanoma. Lasers Surg. Med. 2019, 51, 214–222. [Google Scholar] [CrossRef] [PubMed]
  10. Kazianka, H.; Leitner, R.; Pilz, J. Segmentation and classification of hyper-spectral skin data. In Studies in Classification, Data Analysis, and Knowledge Organization; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  11. Nagaoka, T.; Nakamura, A.; Kiyohara, Y.; Sota, T. Melanoma Screening System Using Hyperspectral Imager Attached to Imaging Fiberscope. In Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society, San Diego, CA, USA, 28 August–1 September 2012. [Google Scholar]
  12. Tomatis, S.; Carrara, M.; Bono, A.; Bartoli, C.; Lualdi, M.; Tragni, G.; Colombo, A.; Marchesini, R. Automated melanoma detection with a novel multispectral imaging system: Results of a prospective study. Phys. Med. Biol. 2005, 50, 1675–1687. [Google Scholar] [CrossRef] [PubMed]
  13. Moncrieff, M.D.; Cotton, S.; Claridge, E.; Hall, P. Spectrophotometric intracutaneous analysis: A new technique for imaging pigmented skin lesions. Br. J. Dermatol. 2002, 146, 448–457. [Google Scholar] [CrossRef] [PubMed]
  14. Fink, C.; Jaeger, C.; Jaeger, K.; Haenssle, H.A. Diagnostic performance of the MelaFind device in a real-life clinical setting. J. Dtsch. Dermatol. Ges. 2017, 15, 414–419. [Google Scholar] [CrossRef] [PubMed]
  15. Elbaum, M.; Kopf, A.W.; Rabinovitz, H.S.; Langley, R.G.B.; Kamino, H.; Mihm, M.C.; Sober, A.J.; Peck, G.L.; Bogdan, A.; Gutkowicz-Krusin, D.; et al. Automatic differentiation of melanoma from melanocytic nevi with multispectral digital dermoscopy: A feasibility study. J. Am. Acad. Dermatol. 2001, 44, 207–218. [Google Scholar] [CrossRef] [PubMed]
  16. Monheit, G.; Cognetta, A.B.; Ferris, L.; Rabinovitz, H.; Gross, K.; Martini, M.; Grichnik, J.M.; Mihm, M.; Prieto, V.G.; Googe, P.; et al. The performance of MelaFind: A prospective multicenter study. Arch. Dermatol. 2011, 147, 188–194. [Google Scholar] [CrossRef] [PubMed]
  17. Stamnes, J.J.; Ryzhikov, G.; Biryulina, M.; Hamre, B.; Zhao, L.; Stamnes, K. Optical detection and monitoring of pigmented skin lesions. Biomed. Opt. Express 2017, 8, 2946–2964. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Abadi, M.; Barham, P.; Chen, J.; Chen, Z.; Davis, A.; Dean, J.; Devin, M.; Ghemawat, S.; Irving, G.; Isard, M.; et al. TensorFlow: A System for Large-Scale Machine Learning. In Proceedings of the 12th USENIX Symposium on Operating Systems Design and Implementation, Savannah, GA, USA, 2–4 November 2016; pp. 265–283. [Google Scholar]
  19. Li, S.; Song, W.; Fang, L.; Chen, Y.; Ghamisi, P.; Benediktsson, J.A. Deep learning for hyperspectral image classification: An overview. IEEE Trans. Geosci. Remote Sens. 2019, 57, 6690–6709. [Google Scholar] [CrossRef] [Green Version]
  20. Fabelo, H.; Carretero, G.; Almeida, P.; Garcia, A.; Hernandez, J.A.; Godtliebsen, F.; Melian, V.; Martinez, B.; Beltran, P.; Ortega, S.; et al. Dermatologic Hyperspectral Imaging System for Skin Cancer Diagnosis Assistance. In Proceedings of the 2019 34th Conference on Design of Circuits and Integrated Systems, Bilbao, Spain, 20–22 November 2019. [Google Scholar]
  21. Holland, P.W.; Welsch, R.E. Robust regression using iteratively reweighted least-squares. Commun. Stat. Theory Methods 1977, 6, 813–827. [Google Scholar] [CrossRef]
  22. Song, E.; Grant-Kels, J.M.; Swede, H.; D’Antonio, J.L.; Lachance, A.; Dadras, S.S.; Kristjansson, A.K.; Ferenczi, K.; Makkar, H.S.; Rothe, M.J. Paired comparison of the sensitivity and specificity of multispectral digital skin lesion analysis and reflectance confocal microscopy in the detection of melanoma in vivo: A cross-sectional study. J. Am. Acad. Dermatol. 2016, 75, 1187–1192. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic overview of the algorithm.
Figure 1. Schematic overview of the algorithm.
Sensors 21 00680 g001
Figure 2. (a) Plot of mean square error (MSE) first-derivative feature (MSE df), second-derivative feature (MSE ddf) and sum feature (MSE totm) for the hyperspectral image (HSI) ROI curves from the training melanoma PSLs. (b) Plot of quartic polynomial fits of melanoma HSI ROI curves chosen according to the low MSE totm. (c) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P62C1003 PSL. (d) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P81C1005 PSL. (e) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P94C1005 PSL. (f) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P95C1000 PSL.
Figure 2. (a) Plot of mean square error (MSE) first-derivative feature (MSE df), second-derivative feature (MSE ddf) and sum feature (MSE totm) for the hyperspectral image (HSI) ROI curves from the training melanoma PSLs. (b) Plot of quartic polynomial fits of melanoma HSI ROI curves chosen according to the low MSE totm. (c) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P62C1003 PSL. (d) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P81C1005 PSL. (e) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P94C1005 PSL. (f) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P95C1000 PSL.
Sensors 21 00680 g002aSensors 21 00680 g002b
Figure 3. (a) Plot of MSE product of max x and max y (MSE Prod) and maximum mean (MSE Mean) for the fitted HSI ROI curves from the malignant PSLs. (b) Plot of quartic polynomial fits of malignant HSI ROI curves chosen according to low MSE mean. (c) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P80C1003 PSL. (d) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P88C1000 PSL. (e) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P91C1003 PSL. (f) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P112C1000 PSL.
Figure 3. (a) Plot of MSE product of max x and max y (MSE Prod) and maximum mean (MSE Mean) for the fitted HSI ROI curves from the malignant PSLs. (b) Plot of quartic polynomial fits of malignant HSI ROI curves chosen according to low MSE mean. (c) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P80C1003 PSL. (d) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P88C1000 PSL. (e) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P91C1003 PSL. (f) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P112C1000 PSL.
Sensors 21 00680 g003
Figure 4. (a) Plot of MSE mean feature (MSE Mean), absolute value of the first derivative feature (MSE abs df), and the sum of these (MSE mabdf) for the HSI ROI curves from the benign PSLs. (b) Plot of quartic polynomial fits of benign HSI ROI curves chosen according to low MSE abs df. (c) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P29C2000 PSL. (d) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P13C2000 PSL. (e) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P18C1000 PSL. (f) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P30C1000 PSL.
Figure 4. (a) Plot of MSE mean feature (MSE Mean), absolute value of the first derivative feature (MSE abs df), and the sum of these (MSE mabdf) for the HSI ROI curves from the benign PSLs. (b) Plot of quartic polynomial fits of benign HSI ROI curves chosen according to low MSE abs df. (c) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P29C2000 PSL. (d) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P13C2000 PSL. (e) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P18C1000 PSL. (f) The corresponding ROI in red drawn on the HSI for one of the channels used in the creation of the HSI ROI curve from the P30C1000 PSL.
Sensors 21 00680 g004
Figure 5. Examples of the fitted curve plots with data points from the validation phase: (a) Plot of a quartic polynomial fit of a melanoma HSI ROI curve and the data points. (b) Plot of a quartic polynomial fit of a malignant HSI ROI curve and the data points. (c) Plot of a quartic polynomial fit of a benign HSI ROI curve and the data points.
Figure 5. Examples of the fitted curve plots with data points from the validation phase: (a) Plot of a quartic polynomial fit of a melanoma HSI ROI curve and the data points. (b) Plot of a quartic polynomial fit of a malignant HSI ROI curve and the data points. (c) Plot of a quartic polynomial fit of a benign HSI ROI curve and the data points.
Sensors 21 00680 g005
Table 1. The results for the region of interest (ROI) investigation for the melanoma pigmented skin lesions (PSLs) with the corresponding curves with maximum mean first derivative, Max df, maximum curvature, Max ddf, maximum weighted sum of mean df, mean ddf and df at the maximum point, Max totm and the corresponding MSEs of the chosen curves.
Table 1. The results for the region of interest (ROI) investigation for the melanoma pigmented skin lesions (PSLs) with the corresponding curves with maximum mean first derivative, Max df, maximum curvature, Max ddf, maximum weighted sum of mean df, mean ddf and df at the maximum point, Max totm and the corresponding MSEs of the chosen curves.
PSLMax dfMax ddfMax totmMSE dfMSE ddfMSE totm
P62C10030.199595.22611.27160.007290.0223530.022353
P81C10050.9011114.4655.01230.0441380.0106910.010691
P82C10000.17383–0.1158−0.33730.013810.0217810.0026241
P94C10051.138845.8125.04850.0087420.0012140.0012136
P95C10000.6049119.0923.32120.0019830.0190060.015033
Table 2. The results for the ROI investigation for the malignant PSLs with the corresponding curves with maximum product of max x times max y, Max Prod, maximum mean, Max Mean, and the corresponding MSEs of the chosen curves.
Table 2. The results for the ROI investigation for the malignant PSLs with the corresponding curves with maximum product of max x times max y, Max Prod, maximum mean, Max Mean, and the corresponding MSEs of the chosen curves.
PSLMax ProdMax MeanMSE ProdMSE Mean
P67C10030.200580.150513.41332.0846
P32C10000.059820.0540240.022440.006436
P87C10000.334990.214311.3230.41494
P106C10000.91940.3584615.47411.997
P21C10000.137090.118650.0398190.16209
P56C10020.731140.171115.63093.7864
P66C10010.176690.115425.35950.64848
P75C10000.368190.180034.19024.234
P77C10000.315250.154021.35320.65383
P80C10030.0941590.134040.20150.033497
P88C10000.163560.142311.59060.038208
P89C10010.236280.153013.10983.8229
P90C10020.331550.196753.8720.087464
P91C10030.17630.128430.0396250.039625
P101C10000.635370.2415425.6391.0538
P104C10000.770460.203048.37248.3684
P110C10000.181870.0975680.967770.96777
P112C10000.349580.333960.192030.045852
P116C10040.263480.145932.92842.9284
P92C10041.47340.2828919.80521.99
P109C10000.381930.234670.233740.48969
Table 3. The results for the ROI investigation for the benign PSLs with the corresponding curves with minimum mean, Min Mean, minimum absolute value of the first derivative Min abs df, minimum sum of mean and absolute value of the first derivative Min mabdf, and the corresponding MSEs of the chosen curves.
Table 3. The results for the ROI investigation for the benign PSLs with the corresponding curves with minimum mean, Min Mean, minimum absolute value of the first derivative Min abs df, minimum sum of mean and absolute value of the first derivative Min mabdf, and the corresponding MSEs of the chosen curves.
PSLMin MeanMin absdfMin mabdfMSE MeanMSE absdfMSE mabdf
P60C30040.0358750.042110.0781260.432490.491820.49182
P24C10000.0477580.137510.202625.72922.97092.9709
P61C10040.0373920.0195410.0705250.115430.120840.049766
P78C30000.0540880.142070.229122.07380.988250.98825
P83C10030.0532590.0519720.155442.23840.0068060.20734
P29C10000.0179860.0377970.0671690.375930.0561220.056122
P29C20000.0164990.0286430.0509630.129630.0295770.029577
P13C20000.0187380.0250280.0686590.0044590.0185440.073912
P13C30000.0207950.0589140.0869640.0975570.0337920.020257
P16C10000.0177590.0499080.0731010.189810.0130010.013001
P17C10010.0198610.0210360.0551370.153540.0130070.028137
P17C20020.0238280.0121840.0475590.0797090.0207210.018856
P18C10000.0238870.0103990.0548890.117620.0073390.0073387
P25C20000.0533360.212510.275974.17420.912040.91204
P25C30000.0578660.139190.220132.366413.1749.6771
P26C10000.0502460.187330.26054.58071.43191.4319
P27C10000.0234080.0529520.081730.121960.0391680.047658
P27C20000.0207740.0451390.0694340.248550.160420.16042
P27C30000.020750.0580850.0860340.0290030.0295910.016113
P27C40000.0191930.0261950.0519280.0706610.0377390.037739
P29C30000.0172970.0427380.0680370.148740.010730.022491
P30C10000.020940.005740.0450260.099030.0018850.0032779
Table 4. The PSLs in the test set and their classification correctness.
Table 4. The PSLs in the test set and their classification correctness.
PSLCorrectness
P100C1000tp
P23C1001fp
P102C1000tp
P28C1000tn
P107C1003tn
P69C1003tp
P13C1000tn
P74C1002tp
P14C1000tn
P97C1004tp
Table 5. Comparison with results from other studies. (* Sensitivity for melanoma detection. ¥ Only reported sensitivity for three melanoma lesions. α Only reported sensitivity for four melanoma lesions).
Table 5. Comparison with results from other studies. (* Sensitivity for melanoma detection. ¥ Only reported sensitivity for three melanoma lesions. α Only reported sensitivity for four melanoma lesions).
ReferencesPatientsImagesBandsRange (nm)Sensitivity (%)Specificity (%)
Tomatis et al. [12]1278139115483–95080.4 *75.6
Kazianka et al. [10]-310300-95 *-
Moncrieff et al. [13]3113488400–1000100 *, ¥5.5
Fink et al. [14]11136010430–950100 *, ¥5.5
Song et al. [22]553610430–95071.4 *, α25
Monheit et al. [16]1257161210430–95098.2 *9.5
Nagaoka et al. [11]97134124380–78096.0 *87
Stamnes et al. [17]-15710365–10009797
Stamnes et al. [17]-71210365–10009993
Hosking et al. [9]1005221350–95036 *100 *
Leon et al. [8]6176116450–95087.5/100 *100
Proposed6176125450–95010080/100 *
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Uteng, S.; Quevedo, E.; M. Callico, G.; Castaño, I.; Carretero, G.; Almeida, P.; Garcia, A.; A. Hernandez, J.; Godtliebsen, F. Curve-Based Classification Approach for Hyperspectral Dermatologic Data Processing. Sensors 2021, 21, 680. https://doi.org/10.3390/s21030680

AMA Style

Uteng S, Quevedo E, M. Callico G, Castaño I, Carretero G, Almeida P, Garcia A, A. Hernandez J, Godtliebsen F. Curve-Based Classification Approach for Hyperspectral Dermatologic Data Processing. Sensors. 2021; 21(3):680. https://doi.org/10.3390/s21030680

Chicago/Turabian Style

Uteng, Stig, Eduardo Quevedo, Gustavo M. Callico, Irene Castaño, Gregorio Carretero, Pablo Almeida, Aday Garcia, Javier A. Hernandez, and Fred Godtliebsen. 2021. "Curve-Based Classification Approach for Hyperspectral Dermatologic Data Processing" Sensors 21, no. 3: 680. https://doi.org/10.3390/s21030680

APA Style

Uteng, S., Quevedo, E., M. Callico, G., Castaño, I., Carretero, G., Almeida, P., Garcia, A., A. Hernandez, J., & Godtliebsen, F. (2021). Curve-Based Classification Approach for Hyperspectral Dermatologic Data Processing. Sensors, 21(3), 680. https://doi.org/10.3390/s21030680

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop