Next Article in Journal
An Efficient Parallel Algorithm for Detecting Packet Filter Conflicts
Next Article in Special Issue
Short Text Classification with Tolerance-Based Soft Computing Method
Previous Article in Journal
Vector Fitting–Cauchy Method for the Extraction of Complex Natural Resonances in Ground Penetrating Radar Operations
Previous Article in Special Issue
Exploring the Suitability of Rule-Based Classification to Provide Interpretability in Outcome-Based Process Predictive Monitoring
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Automatic Classification of Foot Thermograms Using Machine Learning Techniques

1
School of Science and Technology, University of Trás-os-Montes e Alto Douro, Quinta de Prados, 5001-801 Vila Real, Portugal
2
INESC TEC—INESC Technology and Science, 4200-465 Porto, Portugal
3
Mathematics Centre CMAT, Pole CMAT—UTAD, Quinta de Prados, 5001-801 Vila Real, Portugal
*
Author to whom correspondence should be addressed.
Algorithms 2022, 15(7), 236; https://doi.org/10.3390/a15070236
Submission received: 31 May 2022 / Revised: 30 June 2022 / Accepted: 4 July 2022 / Published: 6 July 2022
(This article belongs to the Special Issue Algorithms for Machine Learning and Pattern Recognition Tasks)

Abstract

:
Diabetic foot is one of the main complications observed in diabetic patients; it is associated with the development of foot ulcers and can lead to amputation. In order to diagnose these complications, specialists have to analyze several factors. To aid their decisions and help prevent mistakes, the resort to computer-assisted diagnostic systems using artificial intelligence techniques is gradually increasing. In this paper, two different models for the classification of thermograms of the feet of diabetic and healthy individuals are proposed and compared. In order to detect and classify abnormal changes in the plantar temperature, machine learning algorithms are used in both models. In the first model, the foot thermograms are classified into four classes: healthy and three categories for diabetics. The second model has two stages: in the first stage, the foot is classified as belonging to a diabetic or healthy individual, while, in the second stage, a classification refinement is conducted, classifying diabetic foot into three classes of progressive severity. The results show that both proposed models proved to be efficient, allowing us to classify a foot thermogram as belonging to a healthy or diabetic individual, with the diabetic ones divided into three classes; however, when compared, Model 2 outperforms Model 1 and allows for a better performance classification concerning the healthy category and the first class of diabetic individuals. These results demonstrate that the proposed methodology can be a tool to aid medical diagnosis.

1. Introduction

Technological advances in infrared cameras improved the speed of response and resolution, which allowed for the application of this technology to different areas, namely, to medical use [1]. With this technique, it is possible to visualize the temperature distribution of the human body and detect thermal variations associated with some diseases, such as breast cancer, vascular disorders, muscle pain, complications related to diabetic foot, among others [1,2].
The application of infrared thermography has increased significantly over the years, especially with regard to the study of problems related to diabetic foot [3,4], that is, one of the main complications observed in Diabetes Mellitus (DM) patients, who have a 12% to 25% risk of developing foot ulcers during life [5,6].
Currently, medical information systems that support detection and decision making are changing, evolving from a traditional manual data analysis to more modern and specialized computer methods, known as Computer-Assisted Diagnosis (CAD), thus providing tools that allow health professionals to achieve more accurate diagnosis [7].
The evaluation of medical images by humans is susceptible to errors due to several factors, namely: negligence, fatigue, and even sensory overload, as consequence of the large amount of information that needs to be analyzed. Other factors to be considered are that human visual capacities contain limitations, resulting from the visual perception of optical illusions, which affects the accuracy of the diagnosis. It also happens that health institutions may not have enough specialists to perform the diagnostic task [8]. Therefore, it is important to bridge the gap between the user’s perception and the accuracy of the diagnosis; as the literature review shows, the development of CAD systems helps to overcome these problems, with image processing algorithms with greater detection performance based on automated segmentation, image enhancement, and restoration and feature extraction and classification approach [3,4,9].
There are several studies that focused on the diagnosis of diabetic foot diseases, using thermal images, based on the measurement of skin temperature variation. Many of them seek to classify foot thermograms through asymmetric analysis, which consists of comparing the temperature of the foot with the contralateral one, in order to define a limit that allows for categorizing between healthy and diabetic [10,11,12]. This type of approach has shown good results in several studies; however, it has some limitations, for example, when the patient has similar complications in both feet, as it is not possible to detect the risk areas and incorrectly classify the thermogram. Nevertheless, as techniques to detect areas of risk without contralateral analysis were still needed, as well as to have a consistent accurate algorithm that can effectively detect, analyze and classify the minimum thermal variations in the diabetic foot, machine learning techniques have gained an increasing impact in medical areas. The support vector machine algorithm was used by Saminathan et al. [13] to apply asymmetric analysis to foot thermograms and classify them as normal or prone to ulcer.
To combat the limitations previously mentioned, in [14] the authors propose a thermal change index, which represents the degree of variation of the plantar temperature, to group the thermograms of diabetic individuals in five classes; for that, a public plantar thermogram database [15] was used. This public dataset has been used by several researchers with the aim of analyzing the distribution of the temperature of the feet of (non)diabetic individuals. One of them was Cruz Vega et al. [9], who proposed a deep learning technique to classify the thermograms of the feet of diabetic individuals into five classes, by comparing two classes at a time. Another group was Khandakar et al. [16], who proposed a machine learning-based scoring technique with feature selection, optimization techniques and learning classifiers to perform a binary classification of the thermograms of the feet into diabetic and nondiabetic classes. The same investigators also use the public dataset [15] to propose a machine-learning framework, based on a thermal change index, to classify the thermograms of diabetic individuals into five classes [17]. Isaza and Diaz [18] also used the dataset in [15] to propose the implementation of deep learning architectures under data augmentation methods, using the Fourier transform, to classify each foot thermogram separately as belonging to a diabetic or nondiabetic individual. Despite the great progress, there is still a need to develop techniques to classify thermal images that not only present a robust, accurate and detailed classification of the thermograms, but that are also practical enough to be used by health professionals in their clinical practices.
With this in mind, with this work we aim to develop a functional methodology to analyze and classify a diversity of thermal changes in the plantar region of diabetic and healthy individuals; the public dataset [15] was used. Additionally, based on the experience acquired when working with the data from the aforementioned dataset and on the information collected during the literature review, a three-level scale was adopted to classify the feet of a DM individual into three categories. Thus, in this paper, two models, using different machine learning algorithms, for classifying thermographic images of the feet are proposed. In the first one, the foot thermograms are classified into four classes: healthy and three categories for diabetics. In the second, the classification process has two stages: Initially, each foot is classified as belonging to a DM or healthy individual. Then, if the foot is considered an individual DM foot, a classification refinement is made in a second stage.
In order to extract the image descriptors used for foot classification, a cluster algorithm divides the thermogram into different regions. For each region, a set of regional descriptors are calculated, being used as features to train and test both models. Five machine learning classification algorithms were evaluated in both models: logistic regression (LR), support vector machine (SVM) with two variants, and K-nearest neighbor (k-NN), with two variants.
This paper is organized as follows: In Section 2 the used dataset is introduced and the proposed methodology is described. In Section 3 the obtained results are presented. Finally, Section 4 discusses the results and presents conclusions.

2. Materials and Methods

2.1. Thermographic Image Dataset

In our study, we used a public dataset of thermal images (thermograms) of the plantar region, available online.
In their work, Contreras et al. [15] described the acquisition protocol to collect the thermograms of both feet of 167 individuals (122 diabetics and 45 healthy) recruited from four medical institutions located in the city of Puebla, Mexico, over a period of 3 years. Before the image acquisition procedure, each participant had a resting period of 15 min, in a room with controlled temperature, in order to reach a state of thermodynamic equilibrium and to improve the accuracy of temperature variation detection. Two different infrared cameras were used for thermogram acquisition (FLIR E60 and FLIR E6). The emissivity parameter was set to 0.98, which is the recommended emissivity for measurement of human skin with an IR camera.
The camera was fixed by a vertical tripod placed one meter away from the feet. The raw thermogram dimensions are 320 × 240 pixels, and each one has a corresponding temperature file of equal dimensions, in which the temperature value of each point (pixel in the thermogram) is stored.
After acquisition, the raw images were processed to separate the left and right foot in two different sub-images. These sub-images are segmented and spatial transformed to position the foot in a vertical alignment. The publicly available dataset is composed of individual thermograms of each foot, stored in a csv format.

2.2. Methodology

With the aim of detecting and classifying abnormal changes in the plantar temperature, a methodology having two processing stages (Figure 1), feature extraction and classification, is proposed.
The first stage consists of the extraction of relevant features from the thermograms. After applying machine learning classification algorithms, the subject’s foot is classified as healthy or in one of the three classes of pathology severity.

2.3. Feature Extraction

The extraction of features is one of the most important steps in computer-aided diagnostic systems. It intends to extract the most significant and non-redundant information from the images to facilitate the learning process.
Previous studies showed that an increase in temperature in the plantar region of diabetic patients is related with foot ulcer appearance [19]. Based on this fact, we calculate the global temperature features average, maximum and minimum from the entire plantar region to retain the general foot condition. However, these parameters are not enough to characterize the status of the foot, because it does not have a uniform temperature. Regional features should also be extracted to describe the pattern of temperatures in the different sections of plantar region. Following a previous work of the same authors, the K-means algorithm is applied to divide the foot into five different clusters with different values of temperature [20]. The clustering algorithm allows the distribution of image pixels between the different thermal clusters, as shown in Figure 2.
For each cluster, the average, maximum and minimum temperatures are calculated, obtaining this way a total of 15 regional features of plantar temperatures. The clustering algorithm allows for the detection of small hot spots revealing plantar regions that are prone to the appearance of ulcers. In addition to the features already presented, the Cluster Thermal Index (CTI) is also included in the list of features to provide a quantitative estimation of the thermal changes of the foot caused by DM. The CTI provides a measure of temperature deviation between a subject and the control group, taking into account not only the temperature difference between the clusters but also the range of temperature in the control group. The index is based on the average temperature of each cluster in relation to reference temperatures obtained from the healthy individuals (control group). To establish the reference temperatures in the control group, firstly the clusters are sorted in ascending order according with the centroid value. After, the average temperature of each cluster ( C G i ¯ ) is calculated. Table 1 shows the reference temperatures obtained in the control group, with 45 subjects.
The CTI is calculated as the average of the positive differences between the average temperature in each cluster of a subject and the correspondent reference temperatures, obtained in the control group, as is (1).
CTI = i = 1 k | C G I i I N D i | k ,
where k represents the number of clusters, (CGi) is the reference temperature of cluster i and I N D i is the average temperature in cluster i for the subject.
Table 2 illustrates the result of feature extraction in thermograms of a non-diabetic and a diabetic subject. A total of 19 features are extracted from each thermogram, which are used to feed the classifiers algorithms of the next stage. It is observed that the global temperature values in the diabetic subject and their respective cluster temperatures are further away from the reference values, as shown in Table 1. This results in a CTI value of 3.23. In the case of the non-diabetic subject, we observe smaller deviations between the features values and the references temperatures, resulting in an CTI value close to 0 (0.24). The extraction of features in the different clusters allows for the detection of temperature deviations in foot sections that may not be evident when analyzing the global features.

2.4. Classification Algorithms

Image classification involves the extraction of information from an image and then associating it with a class label. Classifying images in the machine learning domain can be approached as a supervised learning task. In this type of learning, the system must “learn” inductively an expression of a model that describes the data, called the target function. This function is used to predict the value of a variable, called a dependent variable or output variable, from a set of variables, called independent variables or input variables, also known with features.
There are several supervised machine learning algorithms that deal with classification. In this study, Logistic Regression, Support Vector Machines and k-Nearest Neighbors are used.
In general, the logistic regression classifier uses a linear combination with more than one input variable as an argument for the sigmoid function, obtaining a value between 0 and 1 as the output of the function. In particular, an input that produces a result greater than 0.5 is considered to belong to class 1. On the other hand, if the output is less than 0.5, the corresponding input will be classified as belonging to class 0.
In the SVM algorithm, the objective is to find a hyperplane in an N-dimensional space (N represents the number of input variables), which classifies the data. In order to separate the classes, there may be several possible hyperplanes; however, the ideal is to find a plane that maximizes the distances between the class data, called the maximum margin. Maximizing the distance from the margin provides some reinforcement so that the points of future data can be classified with more confidence [21]. The classifier has different kernel functions, allowing it to adjust to different types of data. In this study the classifier will be tested with the linear and the quadratic kernel.
The k-NN algorithm consists of assigning a new example for the most common class among its closest k-neighbors. To classify a point x as belonging to one of the classes, the k observed data points closest to x are found. That is, the classification rule is to assign x to the population that has the most observed data points out of the k- nearest neighbors [22]. For example, if k is 1, then x is simply assigned to your nearest neighbor’s class. This classifier allows for changing the number of neighbors (k), in order to improve the classification. In this work, 3-NN and Weighted k-NN will be applied.

2.5. Proposed Models

In this work, two models for classifying the thermographic images are proposed. The LR algorithm and the variants of SVMs, k-NN were applied to both models.
In the first model, the foot thermograms are classified into four classes: healthy and three categories for diabetics, which represent progressive severity of the disease.
In the second model, the classification consists of two stages: initially a binary classification is performed and the foot thermogram is cataloged as belonging to a healthy or diabetic individual; then, and only for cases classified as diabetic, a second classifier is used to obtain a refinement of the classification for diabetic cases.

2.6. Performance Evaluation

To evaluate and compare the performance of classification algorithms, several experiments were carried out using the metrics: Accuracy, Sensitivity, Specificity, Precision, F-score, and Area Under the Curve (AUC). The first five metrics are defined by (2) to (6), where TP, FP, TN and FN represent the number of cases of True Positive, False Positive, True Negative and False Negative, respectively.
Accuracy = T P + F N T P + T N + F P + F N
Sensitivity = T P T P + F N
Specificity = T N F P + T N
Precision = T P T P + F P
F - score = 2 T P 2 × T P + F P + F N
A Receiver Operating Characteristic (ROC) curve was used to capture the trade-off between the true positive rate (sensitivity) and false positive rate (1-specificity) at thresholds over a continuous range.
The AUC summarizes the information obtained using the ROC curve and, instead of depending on a specific operational point, it is an effective and combined measure of sensitivity and specificity that describes the inherent validity of diagnostic tests [23], with an AUC of 1.0 representing a perfect test, while an AUC of 0.5 defines a random test. The higher this metric (closer to 1), the better the discriminatory power of the model.
In binary classification, accuracy, sensitivity, specificity, precision, and F-score are calculated using the expressions (2) to (6). In the case of multi-class classification, those metrics can be calculated for each class in a one-vs.-all method and determine how often each class makes accurate predictions.
As in the used dataset, the classes are not balanced [15], that is, the percentage of samples in each one of the classes is not similar; it is important to pay attention to the impact of class imbalance on the performance metrics of the classification. With this in mind, the metrics that should be taken into greater consideration are F-score and AUC.

3. Results and Discussion

This section presents the results of experiments with the two models described in Section 2. In Section 3.2, the results of the multi-classification concerning to Model 1 are described, while in Section 3.3 the results of Model 2 are analyzed.

3.1. Experiments

In order to find the features that allow one to obtain better performance, the different classification algorithms were tested and evaluated in different subsets of features. The subset that obtains the best results contains the features referring to the average, maximum, and minimum of the general thermogram of the foot and the respective five clusters. The CTI index is also a relevant resource, so 19 features will be used as inputs to train the different classification algorithms.
To assess the generalizability and validation of the both models, the k-fold cross-validation technique was adopted. In this technique, the dataset D is randomly divided into k subsets, D1; D2;…; Dk, of approximately equal size. The model is trained and tested k times and, for each time I = 1,2,⋯, k, it is trained in all subsets, except Di, which is used for testing. In this work, the 10-fold cross-validation technique is used to divide the 80% of the whole sample into 10 equal-sized subsamples. While preparing the folds, the percentage of samples for each class was maintained in every fold, named the stratified-fold technique. This way, the model gets equally distributed data for training and test folds.
After finding the best classifier algorithm for each one of the two models presented in Section 2.4 and in order to be able to compare them, the remaining 20% of the dataset not used to train/validate the models was used for testing.

3.2. Classification with Model 1

In this subsection, the results of Model 1 are presented. In Model 1, the thermogram is classified into one of four classes, healthy and three different categories for diabetics, according to the degree of variation of the plantar regions, with Class 1 being an earlier stage of the disease revealing small thermal variations and Class 3 a more advanced stage of the disease with large thermal variations. These categories are based on the ones proposed in [14], where the authors use the Thermal Change Index (TCI), which represents the degree of change in plantar temperature, to group the thermograms into five classes. Class 1 corresponds to a TCI ≤ 2, Class 2 to a 2 < TCI ≤ 3, Class 3 to a 3 < TCI ≤ 4, Class 4 to a 4 < TCI ≤ 5 and Class 5 to a TCI > 5. An example of this classification can be seen in Figure 3.
Given the low number of samples for each class, when compared to the number of healthy examples, in the present work the five categories were regrouped into three. Additionally, as the classification algorithms used are more robust and accurate when applied to balanced datasets, when performing this reduction in the number of classes, a more balanced distribution of the thermograms per class was obtained.
The performed class reduction consists of joining the categories with the closest TCI indexes, thus grouping the feet that have more similar temperatures. Applying this methodology, three classes were defined: Class 1 for TCI ≤ 3, Class 2 for 3 <TCI ≤ 5 and Class 3 for TCI > 5. Figure 4 presents the reclassification of the thermograms of Figure 3. Thus, considering the new categorization, the dataset used will have 90 healthy thermograms (27% of the total sample), and 244 diabetic ones, divided as follows: Class 1 has 92 thermograms (28%), Class 2 has 64 (19%) and Class 3 has 88 (27%).
The evaluation of each classification algorithm is measured, separately, for each one of the four classes, using a one-vs.-all method. The results of the performance obtained for each class with the different algorithms are presented in Table 3. The comparison of the classification algorithms is based on the average values, of the four classes, for each one of the metrics used; these values are presented in the Average row in Table 3.
Although the classes considered are more balanced than the ones in [14], the thermograms are not equally divided between the classes; thus, special attention should be paid to the F-score and AUC metrics.
Analyzing all the results obtained with the different algorithms, presented in Table 3, and focusing, in particular, on the Average row, both SVM Linear and LR algorithms present similar metric results, obtaining a difference lower than 1% between each metric, and the same happens with the 3-NN and Weighted k-NN algorithms; however, the first group performs slightly better than the second group, mainly in the AUC metric, which is almost 2% higher. The SVM Quadratic is the algorithm that presents the best results, obtaining in some metrics a difference of 4% when compared to the algorithm that presents worse metric results.
In Figure 5, the ROC curves obtained with the Quadratic algorithm, for each class, are illustrated.
From Figure 5, it is perceptible that, although all the AUC results are higher than 91%, the classes with the best classification are Classes 2 and 3, with the two highest AUC values.

3.3. Classification with Model 2

In this subsection, the results of Model 2 are presented. This model is proposed in order to improve the performance of the classification of the healthy and diabetic classes; a very important step as it contributes to the division of the classification between healthy and diabetic thermograms.
In Model 2, the classification is divided in two stages: initially a binary classification is performed, in order to catalog the foot thermogram as belonging to a healthy or a diabetic individual; for the feet classified as diabetic, a second classifier is used to obtain a refinement of the classification.
The evaluation of the classification algorithms is measured considering as a true positive when the model correctly predicts the diabetic individual and as a true negative when the model correctly predicts the healthy individual. The performances of the different algorithms are shown in Table 4.
It can be seen that LR is the algorithm with the highest values in the different metrics, obtaining differences greater than 10% compared to the algorithm 3-NN, which has the worst results, while the SVM algorithms are the second-best ones, with similar results, followed by Weighted k-NN; however, this was higher than the 3-NN algorithm. In Figure 6, the ROC curves for all of the five algorithms are illustrated.
In the second stage of Model 2, a multi-classification, similar to the one in Model 1, is applied. As the aim of this stage is to refine the diabetic classification, only thermograms of diabetic individuals are used and the same three classes used in Model 1 are considered.
The evaluation of each classification algorithm is measured, separately, for each one of the three classes, using a one-vs.-all method. The results of performance obtained for each class with the different algorithms are shown in Table 5. The comparison of the algorithms is based on the average values, of the three classes, for each of the metrics; these values are shown in the Average row of Table 5.
Analyzing all the results obtained with the different algorithms, presented in Table 5, and focusing, in particular, on the Average row, both SVM Linear and Weighted k-NN algorithms present similar performance, obtaining a difference lower than 1% between each metric, and the same happens with the LR and 3-NN algorithms; however, the first group performs slightly better than the second group. The Quadratic SVM is the algorithm that presents the best results, obtaining in some metrics a difference of 3% when compared to the algorithm that presents worse metric results. In Figure 7, the ROC curves obtained with the SVM Quadratic algorithm, for each one of the three class, are illustrated.
From Figure 7 it is perceptible that, although all the AUC results are higher than 91%, Class 1 is the one whose performance excels, with the highest AUC value. Classes 2 and 3 present very similar AUC values, but considering the values of other metrics, the results in Class 3 are slightly better than in Class 2.

3.4. Comparison of Model 1 with Model 2

In order to compare the performance of both models, the remaining 20% of the dataset not used to train and validate the models was used.
For this comparison, the classification algorithms that showed best performance for each one of the two models were used: as shown in Section 3.2 and Section 3.3, the SVM Quadratic algorithm for Model 1, while for Model 2 is used the Logistic Regression algorithm, for stage 1 (binary classification), and the SVM Quadratic algorithm for stage 2 (multi classification).
The results of the performance obtained with both models are presented in Table 6. To compare them, as before, the average of the values of each metric, presented in the Average row of Table 6.
As only 20% of the data were used to perform the comparison between the two models, which results in a low amount of thermograms per class, it appears that the results of Table 6 are lower when compared to those presented in Table 5 (results of stage 2 in Model 2). However, it must be noted that Model 2 is composed of two stages, and the errors in the binary classification of stage 1 harm the global performance of the model.
Nevertheless, Model 2 performs better than Model 1, for all of the metrics evaluated. It can be verified that the values of the metrics for Classes 2 and 3 are the same for both models, while for the Classes 1 and Healthy the performance is better with Model 2. This result suggests that Model 2 is the one that is able to best discriminate between Healthy individuals and those with slight signs of pathology (Class 1). Additionally, on average, F-score is almost 10% higher for Model 2 and Accuracy 5%; thus, this approach outperforms Model 1 and allows for a more accurate performance classification concerning the healthy category and the first class of diabetic individuals.

4. Discussion and Conclusions

As previously mentioned, there are research works, reported in the literature, that apply machine learning or deep learning techniques with good performance results (taking into account the presented metrics) to analyze the temperature distribution on the soles of the feet of diabetic and healthy individuals, using thermograms. These works either aim to classify the thermograms into only two classes—(1) normal/healthy and (2) diabetics—or proceed with a more differentiated classification of the diabetic class. The works of Saminathan et al. [13], Khandakar et al. [16] and Isaza and Diaz [18] use techniques that only allow for a binary classification of thermograms (fitting the first case), while the works of Cruz Vega et al. [9] and of Khandakar et al. [16] (that only use part of the database [15], as they only involve the classification of thermograms of diabetic individuals into five classes) are in the second group, although in [9] only two by two classes comparisons are performed.
The work developed in this paper aimed to obtain a thermogram classification tool that would allow a response with clinical interest, in the sense of being suitable for being used by health professionals in a clinical environment. Thus, it was intended to develop a tool that is prepared to receive a thermogram at the entrance and that gives as output its classification as healthy, mild diabetic, moderate diabetic or severe diabetic. To this end, the 334 thermograms in the public dataset [15] were used to develop two approaches that allow for the classification of foot thermograms into four classes, either directly, using Model 1, or in two phases, with Model 2. Both models use a one-vs.-all method for class comparison and present good performance (given the obtained values of the metrics); therefore, with this work, we think that we have taken a step towards obtaining a tool that can be used in clinical practice to efficiently classify a diabetic foot.
To synthetize what has been conducted, it is observed that, in this article, a classification methodology is proposed, comparing the performance of two different models, and different machine learning algorithms to classify foot thermograms of healthy and diabetic individuals are presented. For the extraction of the features, the concept of cluster was used both to divide the thermograms into different regions and to extract different regional parameters that are used as features to train and test the models used.
The two proposed models proved to be efficient, allowing us to classify a foot thermogram as belonging to a healthy or diabetic individual, with the diabetic ones divided into three classes. It can be concluded that Model 2, which firstly classifies the feet as healthy or diabetic and then obtains a refinement of the diabetic ones in three different classes, outperforms Model 1, which classifies in just a step a thermogram into four classes (healthy and three other classes for diabetics), allowing for a more precise classification regarding the healthy category and the first class of diabetics, corresponding to an early stage of the disease.
The illustrated results show that the proposed methodology can be a useful tool in helping health professionals in the classification of thermograms and in the early detection of the risk of injuries and in the prevention of ulcerations in the feet.
A limitation of this work is the fact that the used public dataset is unbalanced, thus resulting in the low representation of some classes. As future work, we intend to extend the dataset, in order to balance the number of thermograms of healthy individuals with DM and also the number of diabetics in each class. By using balanced datasets with a significant number of samples, classification models can achieve a more accurate prediction. The extension of the dataset will not only imply the inclusion of more detailed clinical information on the subjects in study (which will allow for a better characterization of the obtained thermograms) but will also involve the application of the concept of active thermography. We observe that active thermography consists of temperature measurement in dynamic states as a result of a certain thermal provocation and allows one to obtain more accurate quantitative information regarding the variation of the temperature [24,25,26]; in our case, external stimulation will be applied by using a cold bag of gel to slightly cool down the soles of the feet, and then a thermographic camera will be used to record the temperature recovery to the initial conditions.

Author Contributions

Conceptualization, V.F. and A.T.; methodology, V.F. and P.T.; software, P.T.; validation, V.F., A.T. and P.T.; formal analysis, V.F. and A.T.; investigation, V.F. and A.T.; resources, V.F., A.T. and P.T.; data curation, V.F., A.T. and P.T.; writing—original draft preparation, A.T. and P.T.; writing—review and editing, V.F. and P.T.; visualization, P.T.; supervision, V.F. and A.T.; project administration, not applicable; funding acquisition, not applicable. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The used dataset is public [15] and is available at https://ieee-dataport.org/open-access/plantar-thermogram-database-study-diabetic-foot-complications (accessed on 2 January 2022).

Acknowledgments

The research of the Ana Teixeira was partially financed by Portuguese Funds through FCT (Fundação para a Ciência e a Tecnologia) within the Projects UIDB/00013/2020 and UIDP/00013/2020.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jones, B.F. A Reappraisal of the Use of Infrared Thermal Image Analysis in Medicine. IEEE Trans. Med. Imaging 1998, 17, 1019–1027. [Google Scholar] [CrossRef]
  2. Lahiri, B.B.; Bagavathiappan, S.; Jayakumar, T.; Philip, J. Medical Applications of Infrared Thermography: A Review. Infrared Phys. Technol. 2012, 55, 221–235. [Google Scholar] [CrossRef]
  3. Adam, M.; Ng, E.Y.K.; Tan, J.H.; Heng, M.L.; Tong, J.W.K.; Acharya, U.R. Computer Aided Diagnosis of Diabetic Foot Using Infrared Thermography: A Review. Comput. Biol. Med. 2017, 91, 326–336. [Google Scholar] [CrossRef]
  4. Madhava Prabhu, S.; Verma, S. A Systematic Literature Review for Early Detection of Type II Diabetes. In Proceedings of the 2019 5th International Conference on Advanced Computing and Communication Systems (ICACCS 2019), Coimbatore, India, 15–16 March 2019; pp. 220–224. [Google Scholar] [CrossRef]
  5. Apelqvist, J.; Larsson, J. What Is the Most Effective Way to Reduce Incidence of Amputation in the Diabetic Foot? Diabetes/Metab. Res. Rev. 2000, 16, S75–S83. [Google Scholar] [CrossRef]
  6. Leung, P. Diabetic Foot Ulcers—A Comprehensive Review. Surgeon 2007, 5, 219–231. [Google Scholar] [CrossRef]
  7. Faust, O.; Rajendra Acharya, U.; Ng, E.Y.K.; Hong, T.J.; Yu, W. Application of Infrared Thermography in Computer Aided Diagnosis. Infrared Phys. Technol. 2014, 66, 160–175. [Google Scholar] [CrossRef]
  8. Suissa, S.; Ernst, P. Optical Illusions from Visual Data Analysis: Example of the New Zealand Asthma Mortality Epidemic. J. Clin. Epidemiol. 1997, 50, 1079–1088. [Google Scholar] [CrossRef]
  9. Cruz-Vega, I.; Hernandez-Contreras, D.; Peregrina-Barreto, H.; Rangel-Magdaleno, J.d.J.; Ramirez-Cortes, J.M. Deep Learning Classification for Diabetic Foot Thermograms. Sensors 2020, 20, 1762. [Google Scholar] [CrossRef] [Green Version]
  10. Van Netten, J.J.; Prijs, M.; Van Baal, J.G.; Liu, C.; Van Der Heijden, F.; Bus, S.A. Diagnostic Values for Skin Temperature Assessment to Detect Diabetes-Related Foot Complications. Diabetes Technol. Ther. 2014, 16, 714–721. [Google Scholar] [CrossRef] [Green Version]
  11. Liu, C.; van Netten, J.J.; van Baal, J.G.; Bus, S.A.; van der Heijden, F. Automatic Detection of Diabetic Foot Complications with Infrared Thermography by Asymmetric Analysis. J. Biomed. Opt. 2015, 20, 026003. [Google Scholar] [CrossRef] [Green Version]
  12. Eid, M.M.; Yousef, R.N.; Mohamed, M.A. A Proposed Automated System to Classify Diabetic Foot from Thermography. Int. J. Sci. Eng. Res. 2018, 9, 371–381. [Google Scholar]
  13. Saminathan, J.; Sasikala, M.; Narayanamurthy, V.B.; Rajesh, K.; Arvind, R. Computer Aided Detection of Diabetic Foot Ulcer Using Asymmetry Analysis of Texture and Temperature Features. Infrared Phys. Technol. 2020, 105, 103219. [Google Scholar] [CrossRef]
  14. Hernandez-Contreras, D.; Peregrina-Barreto, H.; Rangel-Magdaleno, J.; Gonzalez-Bernal, J.A.; Altamirano-Robles, L. A Quantitative Index for Classification of Plantar Thermal Changes in the Diabetic Foot. Infrared Phys. Technol. 2017, 81, 242–249. [Google Scholar] [CrossRef]
  15. Hernandez-Contreras, D.A.; Peregrina-Barreto, H.; Rangel-Magdaleno, J.D.J.; Renero-Carrillo, F.J. Plantar Thermogram Database for the Study of Diabetic Foot Complications. IEEE Access 2019, 7, 161296–161307. [Google Scholar] [CrossRef]
  16. Khandakar, A.; Chowdhury, M.E.H.; Ibne Reaz, M.B.; Md Ali, S.H.; Hasan, M.A.; Kiranyaz, S.; Rahman, T.; Alfkey, R.; Bakar, A.A.A.; Malik, R.A. A Machine Learning Model for Early Detection of Diabetic Foot Using Thermogram Images. Comput. Biol. Med. 2021, 137, 104838. [Google Scholar] [CrossRef]
  17. Khandakar, A.; Chowdhury, M.E.H.; Reaz, M.B.I.; Ali, S.H.M.; Abbas, T.O.; Alam, T.; Ayari, M.A.; Mahbub, Z.B.; Habib, R.; Rahman, T.; et al. Thermal Change Index-Based Diabetic Foot Thermogram Image Classification Using Machine Learning Techniques. Sensors 2022, 22, 1793. [Google Scholar] [CrossRef]
  18. Anaya-Isaza, A.; Zequera-Diaz, M. Fourier Transform-Based Data Augmentation in Deep Learning for Diabetic Foot Thermograph Classification. Biocybern. Biomed. Eng. 2022, 42, 437–452. [Google Scholar] [CrossRef]
  19. Fraiwan, L.; AlKhodari, M.; Ninan, J.; Mustafa, B.; Saleh, A.; Ghazal, M. Diabetic Foot Ulcer Mobile Detection System Using Smart Phone Thermal Camera: A Feasibility Study. Biomed. Eng. Online 2017, 16, 117. [Google Scholar] [CrossRef] [Green Version]
  20. Filipe, V.; Teixeira, P.; Teixeira, A. A Clustering Approach for Prediction of Diabetic Foot Using Thermal Images. In Lecture Notes in Computer Science (Including Subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics); Springer: Berlin/Heidelberg, Germany, 2020; Volume 12251, pp. 620–631. [Google Scholar] [CrossRef]
  21. Kotsiantis, S.B. Supervised Machine Learning: A Review of Classification Techniques. Informatica 2007, 31, 249–268. [Google Scholar]
  22. Wettschereck, D.; Aha, D.W.; Mohri, T. A Review and Empirical Evaluation of Feature Weighting Methods for a Class of Lazy Learning Algorithms. Artif. Intell. Rev. 1997, 11, 273–314. [Google Scholar] [CrossRef]
  23. Hajian-Tilaki, K. Receiver Operating Characteristic (ROC) Curve Analysis for Medical Diagnostic Test Evaluation. Casp. J. Intern. Med. 2013, 4, 627–635. [Google Scholar]
  24. Nowakowski, A.; Kaczmarek, M.; Ruminski, J.; Hryciuk, M.; Renkielska, A.; Grudzinski, J.; Siebert, J.; Jagielak, D.; Rogowski, J.; Roszak, K.; et al. Medical Applications of Model-Based Dynamic Thermography. Thermosense XXIII 2001, 4360, 492–503. [Google Scholar] [CrossRef]
  25. Strakowska, M.; de Mey, G.; Wiecek, B.; Strzelecki, M. A Three Layer Model for the Thermal Impedance of the Human Skin: Modeling and Experimental Measurements. J. Mech. Med. Biol. 2015, 15, 1550044. [Google Scholar] [CrossRef] [Green Version]
  26. Strąkowska, M.; Strąkowski, R.; Strzelecki, M.; de Mey, G.; Więcek, B. Thermal Modelling and Screening Method for Skin Pathologies Using Active Thermography. Biocybern. Biomed. Eng. 2018, 38, 602–610. [Google Scholar] [CrossRef]
Figure 1. Methodology overview.
Figure 1. Methodology overview.
Algorithms 15 00236 g001
Figure 2. Thermogram division in clusters. First row: RGB representation of thermal image and the clusters image representation, in non-diabetic (left) and diabetic (right); Second row: the 5 individual clusters.
Figure 2. Thermogram division in clusters. First row: RGB representation of thermal image and the clusters image representation, in non-diabetic (left) and diabetic (right); Second row: the 5 individual clusters.
Algorithms 15 00236 g002
Figure 3. Classification of 5 thermograms of diabetic individuals [14].
Figure 3. Classification of 5 thermograms of diabetic individuals [14].
Algorithms 15 00236 g003
Figure 4. Reclassification of thermograms of diabetics in 3 classes.
Figure 4. Reclassification of thermograms of diabetics in 3 classes.
Algorithms 15 00236 g004
Figure 5. ROC Curves for Model 1, with SVM Quadratic algorithm.
Figure 5. ROC Curves for Model 1, with SVM Quadratic algorithm.
Algorithms 15 00236 g005
Figure 6. ROC Curves for stage 1 in Model 2 for the different algorithms: binary classification.
Figure 6. ROC Curves for stage 1 in Model 2 for the different algorithms: binary classification.
Algorithms 15 00236 g006
Figure 7. ROC curves for stage 2 in Model 2, with SVM Quadratic algorithm.
Figure 7. ROC curves for stage 2 in Model 2, with SVM Quadratic algorithm.
Algorithms 15 00236 g007
Table 1. Average and standard deviation of the temperature (°C) for the Control Group, per cluster.
Table 1. Average and standard deviation of the temperature (°C) for the Control Group, per cluster.
Cluster12345
T c ¯ 23.9325.5226.5927.6128.88
S t d .   d e v . 1.761.541.571.510.92
Table 2. Extraction of features in the thermograms of two subjects: Non-diabetic and Diabetic.
Table 2. Extraction of features in the thermograms of two subjects: Non-diabetic and Diabetic.
Features Non-DiabeticDiabetic
Global Avg.26.7529.75
Max.29.9033.10
Min.22.9426.35
RegionalCluster 1Avg.23.8327.85
Max.24.8028.23
Min.22.9426.35
Cluster 2Avg.25.8228.67
Max.26.3228.24
Min.24.8129.14
Cluster 3Avg.26.8229.66
Max.27.3729.16
Min.26.3330.18
Cluster 4Avg.27.9430.73
Max.28.5530.20
Min.27.3831.21
Cluster 5Avg.29.1331.76
Max.29.9031.23
Min.28.5533.10
Deviation from
reference
Cluster Thermal Index0.243.23
Table 3. Metrics scores for Model 1, with 4 classes.
Table 3. Metrics scores for Model 1, with 4 classes.
AlgorithmClassAccuracySensitivitySpecificityPrecisionAUCF-Score
SVM LinearHealthy0.8510.6530.9230.7580.9270.701
10.8430.7970.8610.6860.9230.738
20.9370.8330.9590.8160.9390.825
30.9810.9460.9950.9860.9840.966
Average0.9030.8070.9350.8120.9440.807
SVM QuadraticHealthy0.8540.7780.8830.7090.9150.742
10.8770.7300.9330.8060.9300.766
20.9370.8330.9590.8160.9730.825
30.9810.9590.9900.9730.9830.966
Average0.9120.8250.9410.8260.9500.825
3-NNHealthy0.8170.7080.8570.6460.8910.675
10.8320.6490.9020.7160.8660.681
20.9370.8130.9640.8300.9520.821
30.9740.9590.9790.9470.9890.953
Average0.8900.7820.9260.7850.9240.783
Weighted k-NNHealthy0.8170.7080.8570.6460.9070.675
10.8320.6620.8970.7100.9060.685
20.9370.7920.9680.8440.9680.817
30.9810.9730.9850.9600.9890.966
Average0.8920.7840.9270.7900.9430.786
LRHealthy0.8660.7640.9030.7430.9240.753
10.8770.7570.9230.7890.9390.772
20.9180.8130.9410.7500.9640.780
30.9590.9050.9790.9440.9940.924
Average0.9050.8100.9370.8060.9550.807
Table 4. Metrics scores for stage 1 in Model 2—binary classification.
Table 4. Metrics scores for stage 1 in Model 2—binary classification.
AlgorithmAccuracySensitivitySpecificityPrecisionAUCF-Score
SVM Linear0.8580.8980.7500.9070.9020.903
SVM Quadratic0.8580.8880.7780.9160.9080.902
3-NN0.7870.8320.6670.8720.8720.851
Weighted k-NN0.8250.8670.7080.8900.9200.879
LR0.8920.9230.8060.9280.9540.926
Table 5. Metrics scores for stage 2 in Model 2—multi-classification.
Table 5. Metrics scores for stage 2 in Model 2—multi-classification.
AlgorithmClassAccuracySensitivitySpecificityPrecisionAUCF-Score
SVM Linear10.9640.9460.9750.9590.9960.952
20.9390.9170.9460.8460.9790.880
30.9750.9460.9920.9860.9830.966
Average0.9590.9360.9710.9300.9860.933
SVM Quadratic10.9800.9860.9750.9610.9980.973
20.9490.8750.9730.9130.9830.894
30.9690.9590.9750.9590.9810.959
Average0.9660.9400.9750.9440.9870.942
3-NN10.9590.9460.9670.9460.9810.946
20.9180.8330.9460.8330.9400.833
30.9590.9460.9670.9460.9840.946
Average0.9460.9080.9600.9080.9680.908
Weighted k-NN10.9640.9460.9750.9590.9960.952
20.9390.8750.9590.8750.9810.875
30.9740.9730.9750.9600.9890.966
Average0.9590.9310.9700.9310.9890.931
LR10.9540.9190.9750.9580.9910.938
20.9180.8960.9260.7960.9530.843
30.9640.9320.9840.9720.9820.952
Average0.9460.9160.9620.9090.9750.911
Table 6. Comparison between Model 1 and Model 2.
Table 6. Comparison between Model 1 and Model 2.
ModelClassAccuracySensitivitySpecificityPrecisionF-Score
Model 1Healthy0.8330.6670.8960.7060.686
10.8030.6110.8750.6470.629
20.9390.9380.9400.8330.882
30.9700.9290.9810.9290.929
Average0.8860.7860.9230.7790.781
Model 2Healthy0.9240.8330.9580.8820.857
10.8940.7780.9380.8240.800
20.9390.9380.9400.8330.882
30.9700.9290.9810.9290.929
Average0.9320.8690.9540.8670.867
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Filipe, V.; Teixeira, P.; Teixeira, A. Automatic Classification of Foot Thermograms Using Machine Learning Techniques. Algorithms 2022, 15, 236. https://doi.org/10.3390/a15070236

AMA Style

Filipe V, Teixeira P, Teixeira A. Automatic Classification of Foot Thermograms Using Machine Learning Techniques. Algorithms. 2022; 15(7):236. https://doi.org/10.3390/a15070236

Chicago/Turabian Style

Filipe, Vítor, Pedro Teixeira, and Ana Teixeira. 2022. "Automatic Classification of Foot Thermograms Using Machine Learning Techniques" Algorithms 15, no. 7: 236. https://doi.org/10.3390/a15070236

APA Style

Filipe, V., Teixeira, P., & Teixeira, A. (2022). Automatic Classification of Foot Thermograms Using Machine Learning Techniques. Algorithms, 15(7), 236. https://doi.org/10.3390/a15070236

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