Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Single snapshot spatial frequency domain imaging for risk stratification of diabetes and diabetic foot

Open Access Open Access

Abstract

Diabetic foot is one of the major complications of diabetes. In this work, a real-time Single Snapshot Multiple-frequency Demodulation (SSMD) - Spatial Frequency Domain Imaging (SFDI) system was used to image the forefoot of healthy volunteers, diabetes, and diabetic foot patients. A layered skin model was used to obtain the 2D maps of optical and physiological parameters, including cutaneous hemoglobin concentration, oxygen saturation, scattering properties, melanin content, and epidermal thickness, from every single snapshot. We observed a strong correlation between the measured optical and physiological parameters and the degree of diabetes. The cutaneous hemoglobin concentration, oxygen saturation, and epidermal thickness decrease, whereas the melanin content increases with the progress of diabetes. The melanin content further increases, and the reduced scattering coefficient and scattering power are lower for diabetic foot patients than those of both healthy and diabetic subjects. High accuracies (AUC) of 97.2% (distinguishing the diabetic foot patients among all subjects), 95.2% (separating healthy subjects from the diabetes patients), and 87.8% (classifying mild vs severe diabetes), respectively, are achieved in binary classifications in sequence using the SSMD-SFDI system, demonstrating its applicability to risk stratification of diabetes and diabetic foot. The prognostic value of the SSMD-SFDI system in the prediction of the occurrence of the diabetic foot and other applications in monitoring tissue microcirculation and peripheral vascular disease are also addressed.

© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Diabetes is a common disease that affected 425 million people worldwide in 2017, and the number is estimated to reach 700 million in 2045. According to the International Diabetes Federation (IDF), one in eleven adults are expected to have diabetes among whom 2/3 are at the prime time. The low awareness and poor control of the disease for diabetes patients can easily lead to various complications, including cardiovascular, retinal-related complications, and diabetic foot [1]. A quarter of diabetes patients will suffer from diabetic foot eventually. The Global Lower Limb Amputation Study Group estimated that between 25 and 90 percent of amputations are related to diabetes, and 46% of diabetes deaths occur at the prime time [1,2].

Diabetic foot is generally classified into neuropathy, peripheral vascular disease, and mixed type. Diabetic foot diagnosis is foremost based on neurologic examination of skin temperature and perfusion pressure, transcutaneous oxygen pressure, the dorsal foot pulses, and the measurement of ankle or toe blood pressure [3]. Vascular ultrasound and angiography have been widely used in the clinic to assist the diagnosis. Non-invasive detection methods such as X-ray, CT imaging, MRI, and radionuclide scanning are also used for diabetic foot, in particular, with infection and severe ischemic ulcers [3,4]. However, current non-invasive methods have several deficiencies. First, the measurement may not necessarily reflect ischemia at the site of the ulcer. With different blood vessels supplying blood to the feet, low blood flow in one vessel may not accurately reflect blood supply at the location of the ulcer [5,6]. Second, international guidelines indicate that there is a grey area where healing probabilities are unclear [7]. For example, toe pressure values between 30 and 50mmHg indicate a 20-80% chance of healing without vascular intervention [6,8,9]. This condition is common in diabetic patients whose microcirculation is affected but has not altered the ankle pressures nor toe pressures used to assess ischemia. Such uncertainty may cause delayed vascular intervention and prolonged healing times. Overcoming these shortcomings of the existing non-invasive diagnostics remains one major challenge in the field of diabetic foot ulcers [7]. There is an unmet need for a non-invasive real-time continuous imaging technique for mapping local micro-circulation over a large field of view and timely risk stratification of diabetes and diabetic foot ulcers [10,11].

In this article, we present a real-time non-contact Single Snapshot Multiple Frequency Demodulation - Spatial Frequency Domain Imaging (SSMD-SFDI) platform implemented with a digital mirror device to map the optical properties of the subsurface of the foot skin. The SSMD-SFDI platform is suitable for imaging and monitoring dynamic turbid medium and processes over a large field of view. In the context of foot imaging, this method not only enables rapid mapping of both structural and physiological parameters of the layered foot skin but also identifies the stage of diabetes and assesses the risk of diabetes and diabetic foot. The potential applications of the technique in early diagnosis, disease tracking, and efficacy evaluation of tissue microcirculation and cardiovascular disease complications, in general, are discussed at the end.

2. Materials and methods

2.1 SSMD-SFDI

SSMD-SFDI has been presented elsewhere [12,13]. Here the principle of SSMD-SFDI is briefly outlined. SFDI is a non-contact optical imaging technology with the unique capability of allowing wide-field quantitative mapping of tissue optical properties (absorption and scattering coefficients). In SFDI, the incident spatially modulated light comprises a DC component, $I_{DC}^{(0)}$, and an AC component, $I_{AC}^{(0)}\cos (2\pi fx + \phi )$, modulated at a specific frequency f. The reflectance of the spatially-modulated light is demodulated to obtain the tissue modulation transfer function (MTF). Conventional SFDI utilizes a demodulation technique referred to as the three-phase method and requires three images ${I_{{0^o}}}$, ${I_{{{120}^o}}}$, ${I_{{{240}^o}}}$ acquired at different phase delays (0, 2π/3, and 4π/3) to compute MTF (the ratio of ${I_{AC}}$ and $I_{AC}^{(0)}$) at one spatial frequency [12,23]:

$${I_{AC}} = \frac{{\sqrt 2 }}{3}\sqrt {{{({I_{{0^o}}} - {I_{{{120}^o}}})}^2} + {{({I_{{{120}^o}}} - {I_{{{240}^o}}})}^2} + {{({I_{{{240}^o}}} - {I_{{0^o}}})}^2}}$$

In contrast, the recently proposed single snapshot multiple frequency demodulation (SSMD) method [13,14] can extract multiple modulation transfer functions from a single structured light image containing multiple components of different spatial modulation frequencies. The amplitude of the $i$-th AC component is given by

$${I_{AC,i}} = \frac{{\sqrt {{{[\int\!\!\!\int\limits_\sigma {I(x,y)\cos (2\pi {f_{x,i}}x + 2\pi {f_{y,i}}y)dxdy} ]}^2} + {{[\int\!\!\!\int\limits_\sigma {I(x,y)\sin (2\pi {f_{x,i}}x + 2\pi {f_{y,i}}y)dxdy} ]}^2}} }}{{\int\!\!\!\int\limits_\sigma {{{\cos }^2}(2\pi {f_{x,i}}x + 2\pi {f_{y,i}}y)dxdy} }}$$
where (${f_{x,i}},{f_{y,i}}$) is the spatial frequency of the $i$-th spatial modulation pattern and σ represents the kernel window of an appropriate size over which the integration is being performed [16].

SSMD can extract both AC and DC components from a single structured light image, whereas conventional three phase demodulation requires three images to obtain AC and DC components. The main advantage of SSMD over other single snapshot SFDI approaches includes the superiority in SNR, the robustness against measurement imperfections, and the capability of extracting multiple AC components of different frequencies from one single snapshot [1315].

2.2 Mapping a two-layer structure to an equivalent homogeneous medium in SFDI

Biological tissue such as the skin is highly structured. The foot skin consists of three main layers from the surface: the epidermis (∼800µm thick, the blood-free layer), dermis (1–4mm thick, vascularized layer), and subcutaneous fat (from 1 to 6mm thick, dependent on the body site). Light penetration is mostly confined within the top two (epidermis and dermis) layers in reflection measurements [17,18], such as SFDI. Accounting for the two-layer structure has been found to be critical to the accurate recovery of cutaneous hemoglobin concentration and oxygen saturation [16,19]. To this end, we have developed a novel approach that maps the two-layer structure to a homogeneous medium for spatial frequency domain imaging (Fig. 1).

 figure: Fig. 1.

Fig. 1. A layered structure is mapped to the equivalent homogeneous medium of the same reduced scattering coefficient and an effective absorption coefficient.

Download Full Size | PDF

The principle of mapping a two-layer structure to an equivalent homogeneous medium is briefly outlined here [12,16,19]. Ignoring the difference in the scattering property between the epidermis and dermis and noting that light absorption in the epidermis is dominated by melanin and that in the dermis by oxy- and deoxyhemoglobins, the mapped homogeneous medium has the same reduced scattering coefficient $\mu _s^{\prime}(\lambda )$ as the original medium and an equivalent absorption coefficient $\mu _a^{}(\lambda )$ determined by requiring total light absorption is equivalent within the two systems, i.e.,

$${\mu _a}L = {\mu _{a,epidermis}}(\lambda )h + {\mu _{a,dermis}}(\lambda )(L - h)$$
where h is the epidermal thickness and L is the mean penetration depth for the modulated light at the spatial frequency f and the wavelength λ, given by
$$L(q,\lambda ) = \frac{{{{(1 + Ql)}^2}{{(2\mu _t^{\prime})}^{ - 2}} + {{(1 + \mu _t^{\prime}l)}^2}{{(2Q)}^{ - 2}} - 2(1 + Ql)(1 + \mu _t^{\prime}l){{(Q + \mu _t^{\prime})}^{ - 2}}}}{{{{(1 + Ql)}^2}{{(2\mu _t^{\prime})}^{ - 1}} + {{(1 + \mu _t^{\prime}l)}^2}{{(2Q)}^{ - 1}} - 2(1 + Ql)(1 + \mu _t^{\prime}l){{(Q + \mu _t^{\prime})}^{ - 1}}}}$$
where $\mu _t^{\prime} \equiv {\mu _a} + \mu _s^{\prime}$, $Q \equiv \sqrt {{q^2} + 3{\mu _a}({\mu _a} + \mu _s^{\prime})}$, $q \equiv 2\pi f$, and l is the extrapolation length for diffusing light.

The fractal refractive index fluctuation dominates light scattering by the epidermis and dermis. Consequently, the reduced scattering coefficient is described by a power law:

$$\mu {^{\prime}_s}(\lambda ) = \mu {^{\prime}_s}({\lambda _0}){({\lambda /{\lambda_0}} )^{ - b}}$$
where b is the scattering power and λ and λ0 are the wavelengths in nanometers [12].

Ignoring the layered structure will lead to significant errors in the recovered physiological parameters and underestimation of the penetration depth of spatially modulated light. The above mapping of a layered structure to an equivalent homogeneous medium not only resolves this deficiency but also obtains both optical and physiological parameters at the same time from one single snapshot [16,20].

2.3 Experiment procedure and data analysis

The experimental setup is shown in Fig. 2. The sinusoidal fringe pattern from the digital micromirror device (DMD, DLP Light Crafter 4500, Texas Instruments) is projected onto the specimen. Backscattered light from the sample is then imaged onto the CCD camera (Point Gray Grasshopper3 GS3-U3-51S5C-C). The illumination uses white light, which comprises three monochromatic components of wavelengths: 623 nm, 540 nm, and 460 nm, respectively. One millimeter on the specimen surface corresponds to 23 pixels on the CCD. The full usable image window is 33mmx33mm.

 figure: Fig. 2.

Fig. 2. Schematic diagram of the SSMD-SFDI imaging system.

Download Full Size | PDF

The SSMD-SFDI system was used to image the bottom of the foot of healthy subjects, diabetes, and diabetic foot subjects. The spatial modulation frequency of the white illumination light is 0.2 mm−1. The CCD exposure time was fixed at 55000µs, and the sampling rate was 6fps. Forefoot (the front part of the foot) was imaged by the SSMD-SFDI system after resting in a supine position for 5-10 minutes [21].

A region of interest (ROI) of size 8.5mmx8.5mm at the center of the image window was selected for data analysis. Single Snapshot Multiple Frequency Demodulation (SSMD) was used to compute the modulation transfer functions, $\textrm{MT}{\textrm{F}_{\textrm{DC}}} = {I_{DC}}/I_{DC}^{(0)}$ and $\textrm{MT}{\textrm{F}_{\textrm{AC}}} = {I_{AC}}/I_{AC}^{(0)}$, at the spatial frequencies f = 0 and 0.2 mm−1 for red, green, and blue light from the recorded color image over the ROI after color correction [12]. The properties of the two-layer skin including oxy- and deoxyhemoglobin concentrations, the reduced scattering coefficient at 540 nm, the scattering power, the melanin concentration, alpha (skin roughness), and the epidermal thickness are then fitted using the function fmincon in MATLAB by minimizing the least squared error:

$$\textrm{error} = \sum\limits_{i = 1}^3 {[{{(\textrm{MT}{\textrm{F}_{\textrm{AC}}}({\lambda _i}) - \textrm{mt}{\textrm{f}_{\textrm{AC}}}({\lambda _i}))}^2} + {{(\textrm{MTF}{}_{\textrm{DC}}({\lambda _i}) - \textrm{mt}{\textrm{f}_{\textrm{DC}}}({\lambda _i}))}^2}]}$$
where i = 1,2,3 represent the three colors, and the theoretical values of the modulation transfer functions mtfDC and mtfAC are computed with the enhanced diffusion model for the equivalent homogeneous medium [20,22].

Monte Carlo simulations have been performed to determine the dependence of the penetration depth on the wavelength and the spatial modulation frequency of spatially modulated light into layered skin [12]. We selected the three visible wavelengths of 623 nm, 540 nm, and 460 nm and the spatial modulation frequency of 0.2 mm−1 in this study as (1) it probes the epidermis and dermis predominantly and can be imaged with ordinary RGB cameras; and (2) it offers high sensitivity for separating the three major chromophores (melanin, deoxy-hemoglobin, and oxyhemoglobin) in the skin. Furthermore, the accuracy of the diffusion model deteriorates with the increase of the spatial modulation frequency higher than 0.2 mm−1.

In this study, a total of 219 ROIs were obtained for 26 healthy volunteers, 163 diabetes patients (87 mild diabetes and 76 severe diabetes subjects), and 30 diabetic foot patients. The diabetic foot patients are mainly of mixed type. Some patients with diabetic foot have severe ulceration at the base of the foot. The imaging window excludes the scabbing area of the foot.

3. Results

3.1 2D maps for optical and physiological properties of the foot skin

The photo of the imaged area and the 2D ROI images of a typical diabetes subject are shown in Fig. 3, including the photograph and the maps for total hemoglobin concentration (THb), tissue oxygen saturation (SO2), melanin concentration (Melanin), reduced scattering coefficient ($\mu _s^{\prime}$), scattering power, alpha, and epidermal thickness. These images show the spatial variations of the skin property. The features in Fig. 3 reflect both the surface texture and the underlying spatial distribution of structures and chromophores of the skin. The ring shape appearing at the lower right part of most images is attributed to a capillary loop. Hereafter, the properties averaged over the ROI will be used as the basis for diagnosis and risk stratification.

 figure: Fig. 3.

Fig. 3. The photo of the imaged area and 2D images for optical and physiological properties of typical foot skin. The red rectangle outlines the ROI within the image window.

Download Full Size | PDF

3.2 Characteristics and differentiation of diabetes and diabetic foot

Figure 4 compares the optical and physiological parameters between healthy, diabetes, and diabetic foot subjects. The total hemoglobin concentration and oxygen saturation are lower, and melanin concentration is higher for diabetes patients than the healthy control, consistent with the fact that the diabetes patients have symptoms of local microcirculation ischemia and hypoxia at the extremity. The increase of melanin concentration reflects that the body endocrine is abnormal for diabetes patients. In contrast to diabetes subjects, the melanin content continues to increase in diabetic foot patients. The reduced scattering coefficient and scattering power are significantly lower, and the epidermal thickness is higher in diabetic foot patients, indicating a skin structural alteration with the further deterioration of the disease. The anti-trend increase of the total hemoglobin in diabetic foot patients is likely due to drug intervention (Metformin and Gliclazide), which enhances microcirculation.

 figure: Fig. 4.

Fig. 4. Comparison of healthy, diabetes, and diabetic foot subjects. The significance is presented by “*” (p < 0.05), “**” (p < 0.01), and “***” (p < 0.001). The values of THb have been multiplied by 10 for clarity.

Download Full Size | PDF

Glycated hemoglobin (HbA1c) reflects the average level of blood glucose control in the past 2-3 months and is not affected by occasional changes in blood glucose [24]. According to IDF, the range of control is when HbA1c is less than 6.5%, and HbA1c is a risk factor for the occurrence of chronic complications leading to diabetic nephropathy, atherosclerosis, cataract, etc. when HbA1c is higher than 9%. We hence further divided the experimental data into four groups according to HbA1c: healthy subjects, mild diabetes patients (HbA1c: 5.5%-9%), severe diabetes patients (HbA1c: >9%), and diabetic foot patients. It is worth mentioning that HbA1c of diabetic foot subjects ranged from 5.7% to 9.2% in our data set, which was lower than that in severe diabetes subjects (HbA1c: 9.5%-16.5%). The same trends observed in Fig. 4 are also evident across the four groups (see Fig. 5).

 figure: Fig. 5.

Fig. 5. Structural properties and physiological parameters of healthy, mild diabetes, severe diabetes, and diabetic foot subjects. The significance is presented by “*” (p < 0.05), “**” (p < 0.01), and “***” (p < 0.001). The values of THb have been multiplied by 10 for clarity.

Download Full Size | PDF

Referring to both Fig. 4 and Fig. 5, several important observations are in order. First, the reduction in blood supply to the patient’s extremities (feet) leads to a decrease in intracellular thiols, which ultimately contributes to the increase in melanin production of the body, as shown in the steady rise of the melanin content from healthy, diabetes, to diabetic foot subjects. Second, microvascular ischemia and neuropathy lead to tissue nutritional dystrophia, which is manifested in the reduced epidermal thickness with the progress of diabetes. The thinned skin barrier is more prone to damage by fungal infection and trauma, eventually resulting in foot ulcers. The patients presented with foot ulcers, after healing, tend to have a thicker epidermis due to scabbing.

The characteristics of the optical and physiological parameters for healthy, mild diabetes, severe diabetes, and diabetic foot subjects are summarized in Table 1. The number in “()” shows the standard deviation computed across the subjects of a common type.

Tables Icon

Table 1. Mean and Standard Deviation of Optical and Physiological Parameters for Healthy, Mild Diabetes, Severe Diabetes, and Diabetic Foot Subjects.

Based on the Wilcoxon rank sum test, the oxygen saturation, epidermal thickness, and alpha are significantly different between mild and severe diabetes patients, attributed to the changes in microcirculation such as ischemia, hypoxia, and nutritional dystrophia (see Table 2).

Tables Icon

Table 2. Key Properties Differentiating Mild and Severe Diabetes.

The total hemoglobin concentration, melanin concentration, epidermal thickness, and scattering power are significantly altered between severe diabetes and diabetic foot patients, as the outcome of the deterioration of microvascular ischemia and neuropathy leading to altered skin structure (see Table 3).

Tables Icon

Table 3. Key Properties Differentiating Severe Diabetes and Diabetic Foot.

3.3 Risk stratification of diabetes and diabetic foot

Three binary classifications using support vector machines (SVM) [25,26] were then conducted in sequence to classify the cohort of subjects into diabetic foot patients, healthy subjects, severe diabetes, or mild diabetes patients, as outlined in Fig. 6. A weighted SVM scheme with nested 5-fold cross-validation was used to suppress any potential biases and account for the imbalanced data of a limited sample size explicitly [25,26]. The ROC (Receiver Operating Characteristic) curves for the three successive binary classifications were shown in Fig. 7, Fig. 8, and Fig. 9, respectively. The AUC (the area under the ROC curve) for the three binary classifications is 0.972, 0.952, and 0.878, respectively.

 figure: Fig. 6.

Fig. 6. Diagram for the diagnosis of disease progression.

Download Full Size | PDF

 figure: Fig. 7.

Fig. 7. ROC curve for distinguishing the diabetic foot patients from the rest.

Download Full Size | PDF

 figure: Fig. 8.

Fig. 8. ROC curve for distinguishing the healthy subjects from the rest (mild and severe diabetes patients).

Download Full Size | PDF

 figure: Fig. 9.

Fig. 9. ROC curve for distinguishing mild and severe diabetes subjects.

Download Full Size | PDF

At present, there are no clinical methods that can be used to predict the occurrence of diabetic foot [9]. The optical and physiological parameters measured by the SSMD-SFDI system may potentially be used to assess the risk of developing into the diabetic foot from diabetes. Among the measured parameters differentiating diabetic foot patients from other subjects, the melanin content reflects the internal endocrine condition of diabetes patients and is least likely to be affected by drug intervention. Other parameters, including the epidermal thickness, in particular, will also be of prognostic values but cannot be included as they are either altered by drug intervention or scabbing in the diabetic foot in the current data set. The melanin content has been found to increase from healthy, diabetes, and diabetic foot subjects steadily (see Table 1). We thus used a single-factor logistic regression with 5-fold cross-validation on the melanin content to predict the risk of the diabetic foot [27]. Figure 10 shows the ROC curve for the logistic regression model for predicting the diabetic foot in severe diabetes subjects (AUC = 93.0%). The optimal cutoff point for the melanin content using Youden’s index is 0.79 mM, with the sensitivity and specificity at 90.0% and 84.2%, respectively [28].

 figure: Fig. 10.

Fig. 10. ROC curve for the logistic risk prediction model for the diabetic foot.

Download Full Size | PDF

4. Discussion

We have presented SSMD-SFDI, a non-contact imaging technique, for risk stratification of diabetes and diabetic foot. SSMD-SFDI is capable of mapping optical properties (reduced scattering coefficient and scattering power) and physiological parameters including hemoglobin concentration, oxygen saturation, melanin content, and epidermal thickness over a large field of view for foot skin in real-time. The melanin concentration increases, whereas the hemoglobin concentration decreases with the occurrence of mild diabetes. The hemoglobin concentration, oxygen saturation, alpha, and epidermal thickness further decrease when the patient’s condition deteriorates into severe diabetes. The reduced scattering coefficient and scattering power alter for diabetes foot patients, resulting from permanent tissue structural damage. The measured optical and physiological parameters were shown to be capable of risk stratification of diabetes and diabetic foot.

The hyperglycemia toxicity of diabetes is related to local microcirculation complications. Microvascular dysfunction can be attributed to endothelial dysfunction, increased brittleness of vessels, novel branches of arteries, and platelet aggregation, causing prolonged healing of ulcers in diabetic patients [29]. Microvascular ischemia and neuropathy lead to tissue dystrophia, fungal infections, and trauma, ultimately resulting in the occurrence of diabetic foot. These physiological alterations manifest themselves in the reduction of cutaneous hemoglobin concentration and oxygen saturation, and the increase of melanin concentration. The structural parameters such as the epidermal thickness, the reduced scattering coefficient, and the scattering power mark the alterations in local tissue structure [3033]. The oxygen saturation drops significantly from mild diabetes to severe diabetes (see Fig. 5). No significant difference in SO2 between diabetes and the diabetic foot is observed, likely due to drug interventions on diabetic foot patients to enhance microcirculation. The reduction in blood supply to the patient’s extremities (feet) results in a decrease in intracellular thiols, which ultimately leads to an increase in melanin production of the body.

The optical and physiological parameters measured by the SSMD-SFDI system can accurately differentiate between healthy subjects, mild diabetes patients, severe diabetes patients, and diabetic foot patients. In our study of 26 healthy volunteers, 163 diabetic patients (87 mild diabetes and 76 severe diabetes subjects), and 30 diabetic foot patients, the accuracies (AUC) of the successive binary classifications were 97.2% (distinguishing the diabetic foot patients among all subjects), 95.2% (separating the healthy subjects from diabetes patients), and 87.8% (classifying mild vs severe diabetes), respectively. We have further shown the potential of the melanin content to assess the risk of diabetes developing into the diabetic foot with an accuracy of 93.0%. We note that a weighted SVM scheme with nested 5-fold cross-validation has been adopted to suppress potential biases and address the limited sample size in our analysis. A study of a larger sample size and, in particular, monitoring the same subjects over the full span of disease progression will be warranted to corroborate the above results.

As it directly interrogates local microcirculation, metabolic, and structural alterations, the SSMD-SFDI approach can complement or even replace the existing diagnostic assessment of microcirculation for diabetes patients [5,6,34]. User-friendliness, low cost, and noninvasiveness are the advantages of the SSMD-SFDI system as well as one portable device convenient for bedside or point-of-care measurements. The SSMD-SFDI system is hence suitable for the early diagnosis and treatment tracking of the diabetic foot and peripheral vascular diseases [34,35].

Finally, the study of tissue changes in microstructure and local blood circulation will be helpful to understand pathogenesis, the progress of the disease, and guide clinical treatment of the diabetic foot. The SSMD-SFDI system can be useful in monitoring and efficacy tracking of patients for the long-term in the ward. Future work is planned to study the application of this system in areas such as monitoring anesthesia, neurological function detection, and image-guided surgery.

5. Conclusion

In summary, we have demonstrated a method of using the SSMD-SFDI system and a layered skin model to assess the local microcirculation in diabetes and diabetic foot. The optical and physiological parameters were obtained, including cutaneous hemoglobin concentration, oxygen saturation, scattering properties, melanin content, alpha, and epidermal thickness. The parameters such as cutaneous hemoglobin concentration, oxygen saturation, epidermal thickness, and reduced scattering coefficient show a strong correlation with the progress of diabetes and have been demonstrated to be applicable to risk stratification and prognosis of diabetes and diabetic foot. The melanin content has been shown as one potential risk predictor for the occurrence of diabetic foot. As a non-contact technique, SSMD-SFDI is suitable for the efficacy tracking and long-term monitoring of diabetes. This system also has potential applications in monitoring microcirculatory and peripheral vascular diseases, in general.

Funding

Key Research Program of Zhejiang Natural Science Foundation of China (LZ16H180002); US National Science Foundation (NSF) (1607664); Wenzhou Science and Technology Major Project (ZS2017022); National Natural Science Foundation of China (81670777); Natural Science Foundation of Zhejiang Province (LZ19H020001).

Acknowledgments

We thank Chao Z for clinical support.

Disclosures

We report a patent application related to the technology and analysis methods described in this study. The authors declare that there are no other conflicts of interest related to this article.

References

1. N. H. Cho, J. E. Shaw, S. Karuranga, Y. Huang, J. D. da Rocha Fernandes, A. W. Ohlrogge, and B. Malanda, “IDF Diabetes Atlas: Global estimates of diabetes prevalence for 2017 and projections for 2045,” Diabetes Res. Clin. Pract. 138, 271–281 (2018). [CrossRef]  

2. N. Maria, K. Tatjana, C. Heiner, D. Sigrid, S. Bjorn, M. Stephan, R. Gerhard, V. A. Kristien, I. Andrea, and G. Alena, “Incidence of lower extremity amputations in the diabetic compared with the non-diabetic population: A systematic review,” PLoS One 12(8), e0182081 (2017). [CrossRef]  

3. B. A. Lipsky, A. R. Berendt, P. B. Cornia, J. C. Pile, and E. J. G. Peters, “2012 infectious diseases society of America clinical practice guideline for the diagnosis and treatment of diabetic foot infections,” J. Am. Podiatr. Med. Assoc. 103(1), 2–7 (2013). [CrossRef]  

4. C. C. L. M. Naves, “The diabetic foot: a historical overview and gaps in current treatment,” Adv Wound Care 5(5), 191–197 (2016). [CrossRef]  

5. O. A. Mennes, J. J. van Netten, R. Slart, and W. Steenbergen, “Novel optical techniques for imaging microcirculation in the diabetic foot,” Curr. Pharm. Des. 24(12), 1304–1316 (2018). [CrossRef]  

6. B. E. Sumpio, R. O. Forsythe, K. R. Ziegler, J. G. van Baal, M. J. Lepantalo, and R. J. Hinchliffe, “Clinical implications of the angiosome model in peripheral vascular disease,” J. Vasc. Surg. 58(3), 814–826 (2013). [CrossRef]  

7. R. J. Hinchliffe, J. R. Brownrigg, J. Apelqvist, E. J. Boyko, R. Fitridge, J. L. Mills, J. Reekers, C. P. Shearman, R. E. Zierler, and N. C. Schaper, “International Working Group on the Diabetic, “IWGDF guidance on the diagnosis, prognosis and management of peripheral artery disease in patients with foot ulcers in diabetes,” Diabetes/Metab. Res. Rev. 32(Suppl 1), 37–44 (2016). [CrossRef]  

8. J. R. W. Brownrigg, R. J. Hinchliffe, J. Apelqvist, E. J. Boyko, R. Fitridge, J. L. Mills, J. Reekers, C. P. Shearman, R. E. Zierler, and N. C. Schaper, “Effectiveness of bedside investigations to diagnose peripheral artery disease among people with diabetes mellitus: a systematic review,” Diabetes/Metab. Res. Rev. 32(Suppl 1), 119–127 (2016). [CrossRef]  

9. J. R. W. Brownrigg, R. J. Hinchliffe, J. Apelqvist, E. J. Boyko, R. Fitridge, J. L. Mills, J. Reekers, C. P. Shearman, R. E. Zierler, and N. C. Schaper, “Performance of prognostic markers in the prediction of wound healing or amputation among patients with foot ulcers in diabetes: a systematic review,” Diabetes/Metab. Res. Rev. 32(Suppl 1), 128–135 (2016). [CrossRef]  

10. S. Gioux, A. Mazhar, and D. J. Cuccia, “Spatial frequency domain imaging in 2019: principles, applications, and perspectives,” J. Biomed. Opt. 24(7), 1–18 (2019). [CrossRef]  

11. S. Nothelfer, F. Bergmann, A. Liemert, D. Reitzle, and A. Kienle, “Spatial frequency domain imaging using an analytical model for separation of surface and volume scattering,” J. Biomed. Opt. 24(07), 1 (2019). [CrossRef]  

12. X. Chen, W. Lin, C. Wang, S. Chen, J. Sheng, B. Zeng, and M. Xu, “In vivo real-time imaging of cutaneous hemoglobin concentration, oxygen saturation, scattering properties, melanin content, and epidermal thickness with visible spatially modulated light,” Biomed. Opt. Express 8(12), 5468 (2017). [CrossRef]  

13. M. Xu, Z. Cao, W. Lin, X. Chen, L. Zheng, and B. Zeng, “Single snapshot multiple frequency modulated imaging of subsurface optical properties of turbid media with structured light,” AIP Adv. 6(12), 125208 (2016). [CrossRef]  

14. Z. Cao, W. Lin, X. Chen, B. Zeng, and M. Xu, “Real-time spatial frequency domain imaging by single snapshot multiple frequency demodulation technique,” Proc. SPIE 10059, 100590Z (2017). [CrossRef]  

15. J. Vervandier and S. Gioux, “Single snapshot imaging of optical properties,” Biomed. Opt. Express 4(12), 2938 (2013). [CrossRef]  

16. R. B. Saager, A. Truong, D. J. Cuccia, and A. J. Durkin, “Method for depth-resolved quantitation of optical properties in layered media using spatially modulated quantitative spectroscopy,” J. Biomed. Opt. 16(7), 077002 (2011). [CrossRef]  

17. S. L. Jacques, “Optical properties of biological tissues: a review,” Phys. Med. Biol. 58(11), R37–R61 (2013). [CrossRef]  

18. S. L. Jacques, R. Samatham, and N. Choudhury, “Rapid spectral analysis for spectral imaging,” Biomed. Opt. Express 1(1), 157–164 (2010). [CrossRef]  

19. D. Yudovsky, J. Q. Nguyen, and A. J. Durkin, “In vivo spatial frequency domain spectroscopy of two layer media,” J. Biomed. Opt. 17(10), 107006 (2012). [CrossRef]  

20. M. Xu, “Low-coherence enhanced backscattering beyond diffusion,” Opt. Lett. 33(11), 1246–1248 (2008). [CrossRef]  

21. S. Kato, K. Yoshitani, Y. Kubota, Y. Inatomi, and Y. Ohnishi, “Effect of posture and extracranial contamination on results of cerebral oximetry by near-infrared spectroscopy,” J. Anesth. 31(1), 103–110 (2017). [CrossRef]  

22. M. Xu, “Diagnosis of the phase function of random media from light reflectance,” Sci. Rep. 6(1), 22535 (2016). [CrossRef]  

23. M. Reilly and M. Xu, “Analytical model for sub-diffusive light reflection and the application to spatial frequency-domain imaging,” Proc. SPIE 9319, 93191A (2015). [CrossRef]  

24. A. Sugawara, K. Kawai, S. Motohashi, K. Saito, and H. Sone, “Erratum: HbA1c variability and the development of microalbuminuria in type 2 diabetes (Tsukuba Kawai Diabetes Registry 2),” Diabetologia 55(8), 2128–2131 (2012). [CrossRef]  

25. A. Vabalas, E. Gowen, E. Poliakoff, and A. J. Casson, “Machine learning algorithm validation with a limited sample size,” PLoS One 14(11), e0224365 (2019). [CrossRef]  

26. S. Angraal, B. J. Mortazavi, A. Gupta, R. Khera, and H. M. Krumholz, “Machine learning prediction of mortality and hospitalization in heart failure with preserved ejection fraction,” JACC Heart Fail 8(1), 12–21 (2020). [CrossRef]  

27. J. C. Stoltzfus, “Logistic regression: a brief primer,” Acad. Emerg. Med. 18(10), 1099–1104 (2011). [CrossRef]  

28. E. F. Schisterman, D. Faraggi, B. Reiser, and J. Hu, “Youden Index and the optimal threshold for markers with mass at zero,” Stat Med 27(2), 297–315 (2008). [CrossRef]  

29. R. O. Forsythe, J. Brownrigg, and R. J. Hinchliffe, “Peripheral arterial disease and revascularization of the diabetic foot,” Diabetes, Obes. Metab. 17(5), 435–444 (2015). [CrossRef]  

30. M. Xu, “Plum pudding random medium model of biological tissue toward remote microscopy from spectroscopic light scattering,” Biomed. Opt. Express 8(6), 2879–2895 (2017). [CrossRef]  

31. Z. Xu, M. Reilley, R. Li, and M. Xu, “Mapping absolute tissue endogenous fluorophore concentrations with chemometric wide-field fluorescence microscopy,” J. Biomed. Opt. 22(6), 066009 (2017). [CrossRef]  

32. T. T. Wu, J. Y. Qu, and M. Xu, “Unified Mie and fractal scattering by biological cells and subcellular structures,” Opt. Lett. 32(16), 2324–2326 (2007). [CrossRef]  

33. M. Xu, T. T. Wu, and J. Y. Qu, “Unified Mie and fractal scattering by cells and experimental study on application in optical characterization of cellular and subcellular structures,” J. Biomed. Opt. 13(2), 024015 (2008). [CrossRef]  

34. R. O. Forsythe and R. J. Hinchliffe, “Assessment of foot perfusion in patients with a diabetic foot ulcer,” Diabetes/Metab. Res. Rev. 32(Suppl 1), 232–238 (2016). [CrossRef]  

35. J. Allen and K. Howell, “Microvascular imaging: techniques and opportunities for clinical physiological measurements,” Physiol. Meas. 35(7), R91–R141 (2014). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (10)

Fig. 1.
Fig. 1. A layered structure is mapped to the equivalent homogeneous medium of the same reduced scattering coefficient and an effective absorption coefficient.
Fig. 2.
Fig. 2. Schematic diagram of the SSMD-SFDI imaging system.
Fig. 3.
Fig. 3. The photo of the imaged area and 2D images for optical and physiological properties of typical foot skin. The red rectangle outlines the ROI within the image window.
Fig. 4.
Fig. 4. Comparison of healthy, diabetes, and diabetic foot subjects. The significance is presented by “*” (p < 0.05), “**” (p < 0.01), and “***” (p < 0.001). The values of THb have been multiplied by 10 for clarity.
Fig. 5.
Fig. 5. Structural properties and physiological parameters of healthy, mild diabetes, severe diabetes, and diabetic foot subjects. The significance is presented by “*” (p < 0.05), “**” (p < 0.01), and “***” (p < 0.001). The values of THb have been multiplied by 10 for clarity.
Fig. 6.
Fig. 6. Diagram for the diagnosis of disease progression.
Fig. 7.
Fig. 7. ROC curve for distinguishing the diabetic foot patients from the rest.
Fig. 8.
Fig. 8. ROC curve for distinguishing the healthy subjects from the rest (mild and severe diabetes patients).
Fig. 9.
Fig. 9. ROC curve for distinguishing mild and severe diabetes subjects.
Fig. 10.
Fig. 10. ROC curve for the logistic risk prediction model for the diabetic foot.

Tables (3)

Tables Icon

Table 1. Mean and Standard Deviation of Optical and Physiological Parameters for Healthy, Mild Diabetes, Severe Diabetes, and Diabetic Foot Subjects.

Tables Icon

Table 2. Key Properties Differentiating Mild and Severe Diabetes.

Tables Icon

Table 3. Key Properties Differentiating Severe Diabetes and Diabetic Foot.

Equations (6)

Equations on this page are rendered with MathJax. Learn more.

I A C = 2 3 ( I 0 o I 120 o ) 2 + ( I 120 o I 240 o ) 2 + ( I 240 o I 0 o ) 2
I A C , i = [ σ I ( x , y ) cos ( 2 π f x , i x + 2 π f y , i y ) d x d y ] 2 + [ σ I ( x , y ) sin ( 2 π f x , i x + 2 π f y , i y ) d x d y ] 2 σ cos 2 ( 2 π f x , i x + 2 π f y , i y ) d x d y
μ a L = μ a , e p i d e r m i s ( λ ) h + μ a , d e r m i s ( λ ) ( L h )
L ( q , λ ) = ( 1 + Q l ) 2 ( 2 μ t ) 2 + ( 1 + μ t l ) 2 ( 2 Q ) 2 2 ( 1 + Q l ) ( 1 + μ t l ) ( Q + μ t ) 2 ( 1 + Q l ) 2 ( 2 μ t ) 1 + ( 1 + μ t l ) 2 ( 2 Q ) 1 2 ( 1 + Q l ) ( 1 + μ t l ) ( Q + μ t ) 1
μ s ( λ ) = μ s ( λ 0 ) ( λ / λ 0 ) b
error = i = 1 3 [ ( MT F AC ( λ i ) mt f AC ( λ i ) ) 2 + ( MTF DC ( λ i ) mt f DC ( λ i ) ) 2 ]
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.