Pressure-induced near infrared spectra response as a valuable source of information for soft tissue classification

Blaž Cugmas; Miran Bürmen; Maksimilijan Bregar; Franjo Pernuš; Boštjan Likar

doi:10.1117/1.JBO.18.4.047002

3 April 2013 Pressure-induced near infrared spectra response as a valuable source of information for soft tissue classification

Blaž Cugmas, Miran Bürmen, Maksimilijan Bregar, Franjo Pernuš, Boštjan Likar

Author Affiliations +

Journal of Biomedical Optics, Vol. 18, Issue 4, 047002 (April 2013). https://doi.org/10.1117/1.JBO.18.4.047002

Abstract

Acquiring near infrared spectra in vivo usually requires a fiber-optic probe to be pressed against the tissue. The applied pressure can significantly affect the optical properties of the underlying tissue, and thereby the acquired spectra. The existing studies consider these effects to be distortions. In contrast, we hypothesize that the pressure-induced spectral response is site- and tissue-specific, providing additional information for the tissue classification. For the purpose of this study, a custom system was designed for dynamic pressure control and rapid acquisition of spectra. The pressure-induced spectral response was studied at three proximate skin sites of the human hand. The diffuse reflectance and scattering were found to decrease with the applied contact pressure. In contrast, the concentrations of chromophores, and consequently the absorption, increased with the applied contact pressure. The pressure-induced changes in the tissue optical properties were found to be site-specific and were modeled as a polynomial function of the applied contact pressure. A quadratic discriminant analysis classification of the tissue spectra acquired at the three proximate skin sites, based on the proposed pressure-induced spectral response model, resulted in a high (90%) average classification sensitivity and specificity, clearly supporting the working hypothesis.

1. Introduction

Near infrared spectroscopy (NIRS) is a rapid, noninvasive spectroscopic technique that has considerable medical diagnostic potential.¹ The numerous biomedical applications of NIRS include performing optical biopsies, detecting breast² and skin³ cancer, making cancer treatment decisions, and defining therapeutic drug levels.¹ Most of these applications utilize the unique ability of NIRS to quantify the tissue chromophores, in particular hemoglobin, lipids, and water, which can be used to differentiate between normal and diseased tissue.¹ However, acquiring NIR spectra in vivo usually requires a fiber-optic probe to be pressed against the tissue. In general, even light pressure can significantly affect the optical properties of the underlying tissue, and thereby the acquired NIR spectra.

One of the first studies on the optical properties of the soft tissue under pressure focused on the reflectance, reduced scattering, and absorption coefficient in the spectral range from 400 to 1800 nm.⁴ The soft tissue samples, including human skin, bovine aorta, and bovine and porcine sclera, were exposed to constant contact pressures of 9.8, 98, and 196 kPa. The measured reflectance changes in the NIR region were not always consistent; however, under 1400 nm, reflectance decreased with the applied pressure. In contrast, the reduced scattering and absorption coefficient were reported to increase with the applied pressure. The increase in absorption was explained by the local tissue compression, which increased the chromophore concentration; leakage of some extracellular tissue fluids; and the volumetric water content. Likewise, the increase in scattering was explained by the pressure-induced increase in the scatterer concentration.

The study of Chen et al.⁵ confirmed that diffuse reflectance of the palm skin in the spectral range from 1100 to 1700 nm decreases with the applied pressure. The spectra were acquired continuously for 60 s after applying a contact pressure from 32 to 75 kPa. The diffuse reflectance spectra exhibited significant changes during the first 30 s of measurements and became stable afterward. Reflectance changes were attributed to the distortion of the internal tissue structure and chromophore distribution. Human skin of the inner forearm was investigated by Atencio et al.⁶ The spectra were acquired in the spectral range from 530 to 680 nm by applying contact pressure from 0 to 10 kPa. The results were inconsistent under 600 nm, while above 600 nm, the diffuse reflectance decreased with the applied pressure. However, the results were not in agreement with the research of Randeberg,⁷ where the diffuse reflectance increased with the applied pressure.

Another study reported that the changes in the soft tissue optical properties under the applied pressure are site-specific.⁸ Human skin was investigated by reflectance spectroscopy in the visible spectral range. The reflectance of the neck skin decreased with the applied pressure, indicating decrease in the scattering. In contrast, the applied pressure increased the scattering of the forehead skin. The different scattering changes were explained by the fact that skull lies close to the skin surface, preventing the dermis from collapsing into the hypodermis, which results in more light reaching the dermis layers, rich in highly scattering collagen. The absorption coefficient decreased with increasing pressure for all measurement sites. The decrease was attributed to the blood being pushed out by the applied pressure, lower perfusion, and modified tissue morphology. It was concluded that short-term ( $< 2 s$ ), low-contact pressure ( $< 9 kPa$ ) effects do not significantly affect the acquired spectra.

The existing studies used simple measurement setups during the acquisition of spectra, which exposed the sample to a constant pressure for an extended period of time. Moreover, the effects of contact pressure on the soft tissue spectra were primarily considered to be a distortion that needs to be analyzed and removed to obtain repeatable and accurate measurements. In contrast, the objective of this study was to analyze the contact pressure-induced spectral response and use the information for tissue classification. Namely, pressure-induced spectral response is tissue- and site-specific, presenting an additional useful dimension in the acquired data (Fig. 1). For the purpose of this study, a custom real-time system for rapid in vivo acquisition of NIR diffuse reflectance spectra of tissue under controlled contact pressure was developed. Such a measurement system enabled the study of the contact pressure-induced dynamic effects in the soft tissue. As a result, a novel model of the pressure-induced soft tissue spectral response was developed and successfully employed for the classification of three proximate skin sites on the human hand.

Fig. 1

Contact pressure as an additional data dimension.

2. Materials and Methods

2.1.

Samples

Measurements were conducted on the skin surface of four human volunteers (Caucasian males). Three measurement sites were chosen: namely, the palm skin above the abductor pollicis brevis muscle, skin in the middle of the wrist crease above the veins, and skin at the lateral part of the wrist above the ulnar styloid process. While the soft muscle tissue and the tissue above veins provided relatively smooth and gradually increasing resistance to the applied pressure, the tissue above the bone exhibited only limited compression. In addition, the soft tissues at the three selected sites are differently perfused and were expected to provide relevant insight into the unique tissue properties observed in the contact pressure-induced spectral response.

2.2.

Instrumentation

The measurement system used in this study is shown in Fig. 2. Spectra were acquired by a commercial NIR spectrometer (NIR-512L-1.7T1, 901 to 1685 nm, Control Development Inc., South Bend, Indiana) employing a broadband halogen light source (AvaLight-Hal LS, Avantes, Apeldoornseweg, the Netherlands) and a stainless steel, fiber-optic diffuse reflectance probe (Avantes, FCR-7IR400-2-ME) consisting of one detection and six illumination fibers. The diameter of the probe was 6.35 mm, so the contact pressure area was $31.7 {mm}^{2}$ . The probe was fixed to a metal tube, which smoothly slid through the metal rings. A force sensor (Honeywell S&C, FS01, Honeywell, Inc., Morristown, New Jersey) was attached between the end of the metal tube and the linear position sensor (9615R5.1KL2.0, 52.3 MM, Bei Sensors, Goleta, California). The travel length of the shaft spring return mechanism was 4 cm. A motorized linear stage provided precise and accurate control of the applied contact pressure.

Fig. 2

Illustration of the employed measurement system.

2.3.

Measurements

The hands of the four volunteers were fixed during the measurement process, which started by activating the motorized linear stage with a predefined motion speed. The motion was stopped when the linear position sensor exhibited full contraction, resulting in around 100 kPa of contact pressure. The spectral and pressure data were acquired synchronously at a rate of 35 Hz. Nine measurements were made at each of the three selected sites (muscle: the palm skin above the abductor pollicis brevis muscle; veins: skin in the middle of the wrist crease above the veins; bone: skin at the lateral part of the wrist above the ulnar styloid process), allowing about 5 min for the tissue to recover from the previous measurement. The 108 measurement sets obtained for the four volunteers were divided into two equal independent set (namely, the training and the test set), each containing the data collected on two volunteers.

2.4.

Spectra Processing

All the acquired spectra were calibrated using a standard diffuse reflectance tile (Spectralon, Labsphere, North Sutton, New Hampshire). The reflectance spectra $R$ were calculated according to the following equation:

Eq. (1)

R = \frac{I - D}{I_{0} - D},

where

I

is the acquired spectrum,

D

is the dark response of the sensor array, and

I_{0}

is the spectrum of the standard diffuse reflectance tile. Due to the poor sensitivity of the employed sensor array at the lower and the upper end of the spectral range, only the spectral range from 950 to 1600 nm was used. Two different modeling approaches were employed in order to study the contact pressure-induced spectral response (Fig. 3).

Fig. 3

Outline of the employed pressure-induced soft tissue spectral response modeling and classification.

The scattering and absorption properties of the tissue were estimated from the steady-state diffuse reflectance measurements in a spectrally constrained manner.⁹^,¹⁰ Briefly, a forward model of the diffuse reflectance $R$ utilizing prior knowledge on the scattering and absorption properties of chromophores was used:

Eq. (2)

R (λ) = \frac{1}{4 π} [z_{0} (μ_{eff} + \frac{1}{ρ_{1}}) \frac{e^{- μ_{eff} ρ_{1}}}{ρ_{1}^{2}} + (z_{0} + 2 z_{b}) (μ_{eff} + \frac{1}{ρ_{2}}) \frac{e^{- μ_{eff} ρ_{2}}}{ρ_{2}^{2}}],

where

z_{0}

,

μ_{eff}

,

ρ_{1}

,

ρ_{2}

, and

z_{b}

are the wavelength-dependent internal variables calculated from the absorption

μ_{a} (λ)

and reduced scattering

μ_{s}^{'} (λ)

coefficients. The absorption coefficient was modeled as a linear combination of the chromophore concentrations:

Eq. (3)

μ_{a} (λ) = c_{water} μ_{a}^{water} (λ) + c_{Hb} [α μ_{a}^{oxy Hb} (λ) + (1 - α) μ_{a}^{deoxy Hb} (λ)] + c_{lip} μ_{a}^{lip},

where

μ_{a}^{water}

,

μ_{a}^{oxy Hb}

,

μ_{a}^{deoxy Hb}

, and

μ_{a}^{lip}

are the wavelength-dependent absorption coefficients of water, oxyhemoglobin, deoxyhemoglobin, and lipids, respectively,¹¹^–¹⁴

c_{water}

,

c_{Hb}

, and

c_{lip}

are the concentrations of water, hemoglobin, and lipids, respectively, and

α

is the oxygenation fraction (i.e., saturation). The reduced scattering coefficient was modeled according to the wavelength-dependent power law:

Eq. (4)

μ_{s}^{'} = u λ^{- v},

where

u

is the scattering magnitude and

v

is the scattering power.⁹^,¹⁰^,¹⁵

To gain insight on the most prominent NIR absorption bands governing the pressure-induced spectral response, an additional analysis of the acquired spectra was performed. The acquired spectra were preprocessed by the extended form of multiplicative scatter correction (EMSC)¹⁶ to roughly separate the tissue scattering and absorption properties. In contrast to the basic form of the multiplicative scatter correction (MSC),¹⁷ EMSC considers light scattering effects as wavelength-dependent and models each acquired spectrum ( $r_{i}$ ) according to

Eq. (5)

r_{i} = a_{i} + b_{i} r_{ref} + k_{i} λ + l_{i} λ^{2},

where

a_{i}

is the baseline,

b_{i}

is the path length,

r_{ref}

is the reference spectrum, and

k_{i}

and

l_{i}

model the wavelength-dependent spectral variations. The EMSC coefficients (

a

,

b

,

k

,

l

) estimated by a least-squares solution were used to calculate the preprocessed spectrum

r_{EMSC}

as

Eq. (6)

r_{EMSC} = (r_{i} - a_{i} - k_{i} λ - l_{i} λ^{2}) / b_{i} .

In order to extract the spectral response induced by the applied contact pressure, the preprocessed reflectance spectra were subjected to the principal component analysis (PCA) followed by a detailed analysis of the spectral bands highlighted by the calculated PCA loading.

The estimated concentrations of chromophores ( $c_{water}$ , $c_{Hb}$ , and $c_{lip}$ ), the scattering magnitude $u$ and the scattering power $v$ , the oxygenation fraction $α$ , the estimated EMSC coefficients ( $a$ , $b$ , $k$ , $l$ ), and the calculated PCA scores $s$ of the EMSC preprocessed spectra were modeled by a polynomial function of the applied contact pressure $p$ :

Eq. (7)

m = \sum_{j = 0}^{N} h_{m, j} p^{j} + e, m \in {c_{water}, c_{Hb}, c_{lip}, u, v, α, a, b, k, l, s_{1}, s_{2}, \dots, s_{N_{S}}},

where

e

is the fit error. The calculated polynomial coefficients

h_{m}

(without the

h_{m, 0}

) were used to form the feature vectors, employed by the subsequent quadratic discriminant analysis (QDA)¹⁸—based classification into three classes corresponding to the three selected measurement sites. The optimal feature subsets for the two modeling approaches were selected by the sequential forward-floating selection algorithm (SFFS),¹⁹ minimizing the average classification error rate of the training set. The classification sensitivity and specificity were estimated by the independent test set. It should be stressed that only the relative spectral changes were used in this study. Therefore, the polynomial coefficients

h_{m, 0}

, which were not affected by the applied contact pressure, were not used for the classification.

In order to compare the contact pressure-induced spectral response with the findings of other studies, the mean reflectance, the reduced scattering, and absorption coefficients (1450 nm) were modeled by a first-order polynomial function of the applied contact pressure [Eq. (7)].

3. Results and Discussion

3.1.

General Findings

The chromophore concentrations, oxygenation fraction, and the scattering magnitude and power were estimated by fitting each acquired reflectance spectrum to the model [Eq. (2)] in a least-squares sense [Fig. 4(a)], while the PCA scores $s$ were calculated after applying EMSC preprocessing. The spectra exhibited significant and frequently inconsistent changes [Fig. 4(b)] for the contact pressure under 16 kPa. In general, the mean reflectance dropped significantly after the probe touched the skin. However, some of the acquired spectra did not exhibit any obvious change. The diversity of observed spectral changes may be attributed to two main mechanisms. First, the initial light coupling improves with the increasing pressure, gradually leading to a firm contact between the probe and the skin. Second, the contact pressure under the normal blood pressure of 120 mm Hg (16 kPa) substantially affects the blood flow in the upper tissue layers. At higher contact pressure, the light coupling does not improve anymore. However, the pressure is gradually transferred to the deeper tissue layers. Because of that, more and more blood is pushed out of the tissue. Consequently, the spectra exhibit substantial linear and higher-order dependence on the applied contact pressure.

Fig. 4

(a) Example of a measured and modeled [Eq. (2.)] reflectance spectrum. (b) Example of a mean reflectance and PCA score as a function of the applied contact pressure (solid line) and the corresponding polynomial fit (dashed line) calculated for the contact pressure from 20 to 100 kPa.

For the above-explained reasons, only the contact pressure range from 20 to 100 kPa, where the spectral changes were found consistent [Fig. 4(b)], was used by the subsequent analysis. The chromophore concentrations, oxygenation fraction, and the reduced scattering magnitude $u$ and power $v$ were modeled by a first-order polynomial, while the PCA scores exhibited slight nonlinearity, requiring a second-order polynomial model.

3.2.

Classification

Five of the six features describing the tissue scattering and absorption properties (Fig. 3) were selected by the SFFS method for the classification. Only the polynomial coefficient corresponding to the hemoglobin concentration was not. This was somewhat expected, as the hemoglobin absorption in the NIR spectral range is relatively low compared to the absorption of water and lipids. Moreover, the applied contact pressure in excess of 16 kPa removes most of the blood from the underlying tissue. In contrast, only three features, the two polynomial coefficients estimating the first PCA score and the quadratic polynomial coefficient estimating the EMSC parameter $k$ , were selected for the classification based on the second modeling approach. The classification results are summarized in Table 1. The average classification sensitivity and specificity and posterior probability for the first model were 81% and 73%, respectively. The second model exhibited similar classification posterior probability (74%) and slightly better average sensitivity and specificity (90%). The classification specificity was the highest for the skin above the bone (100%), while the highest classification sensitivity was observed for the wrist (100%).

Table 1

Classification sensitivity and specificity and the corresponding posterior probability obtained by polynomial modeling of the tissue absorption and scattering properties (A) and PCA scores (B).

	Average		Skin above muscle		Skin above veins		Skin above bone
	Sensitivity and specificity	Posterior probability	Sensitivity	Specificity	Sensitivity	Specificity	Sensitivity	Specificity
A	0.81	0.73	0.61	0.92	0.89	0.81	0.72	0.89
B	0.90	0.74	0.78	0.94	1.00	0.86	0.83	1.00

3.3.

Tissue Scattering and Absorption Properties as a Linear Function of the Applied Contact Pressure

For all 108 acquired measurement sets, the diffuse reflectance and reduced scattering coefficient at 1450 nm decreased with the applied contact pressure. In contrast, the absorption coefficient at 1450 nm consistently increased with the applied contact pressure. The average reflectance and absorbance changes were in accordance with the majority of the existing studies.⁴^–⁶ However, Lim et al.⁸ reported that the reduced scattering coefficient of the forehead skin increased with the applied pressure. The explanation was that the skull prevents the dermis from collapsing into the hypodermis, and therefore light can penetrate deeper into the tissue rich with the highly scattering collagen. Skin on the forehead can be considered similar to the skin above the ulnar styloid process, which was used in this study. However, no increase in the scattering was observed.

The slopes of the linear polynomial model corresponding to the selected features are listed in Table 2. The variability of the calculated slopes was the smallest for the skin above the muscle tissue and the highest for the skin above the veins. At the wrist, the tendons and bones are close to veins, so the tissue structure can vary significantly; hence, the measurement repeatability is lower. On the other hand, the muscle tissue structure is more homogeneous, resulting in higher measurement repeatability.

Table 2

Median value and the corresponding range between the 25th and 75th percentile (brackets) of the polynomial slope.

Feature	Skin above muscle	Skin above veins	Skin above bone
Mean reflectance ( $10^{- 4} {kPa}^{- 1}$ )	$- 1.20$ $[\begin{matrix} - 1.44 & - 0.94 \end{matrix}]$	$- 2.49$ $[\begin{matrix} - 2.69 & - 2.09 \end{matrix}]$	$- 1.80$ $[\begin{matrix} - 2.18 & - 1.20 \end{matrix}]$
Abs. coef. ( ${μ_{a}}^{1450 nm}$ ) ( $10^{- 3} {mm}^{- 1} {kPa}^{- 1}$ )	1.45 $[\begin{matrix} 0.90 & 1.77 \end{matrix}]$	1.17 $[\begin{matrix} 0.41 & 2.34 \end{matrix}]$	1.22 $[\begin{matrix} 0.94 & 1.96 \end{matrix}]$
Red. scat. coef. ( ${μ_{s}}^{' 1450 nm}$ ) ( $10^{- 2} {mm}^{- 1} {kPa}^{- 1}$ )	$- 0.51$ $[\begin{matrix} - 0.77 & - 0.33 \end{matrix}]$	$- 1.56$ $[\begin{matrix} - 5.80 & - 0.73 \end{matrix}]$	$- 0.94$ $[\begin{matrix} - 1.67 & - 0.67 \end{matrix}]$
Water conc. ( $c_{water}$ ) ( $10^{- 3} {kPa}^{- 1}$ )	3.51 $[\begin{matrix} 2.11 & 4.41 \end{matrix}]$	2.98 $[\begin{matrix} 1.03 & 5.87 \end{matrix}]$	3.05 $[\begin{matrix} 2.38 & 4.84 \end{matrix}]$
Saturation ( $α$ ) ( $10^{- 3} {kPa}^{- 1}$ )	$- 0.59$ $[\begin{matrix} - 0.85 & - 0.39 \end{matrix}]$	$- 2.78$ $[\begin{matrix} - 5.07 & - 1.66 \end{matrix}]$	$- 2.18$ $[\begin{matrix} - 3.05 & - 1.20 \end{matrix}]$
Lipid conc. ( $c_{lip}$ ) ( $10^{- 3} {kPa}^{- 1}$ )	7.00 $[\begin{matrix} 3.78 & 11.03 \end{matrix}]$	5.45 $[\begin{matrix} 1.16 & 12.09 \end{matrix}]$	8.38 $[\begin{matrix} 4.95 & 11.94 \end{matrix}]$
Sc. magnitude ( $u$ ) ( $10^{- 3} {kPa}^{- 1}$ )	1.31 $[\begin{matrix} 0.72 & 2.42 \end{matrix}]$	3.13 $[\begin{matrix} 2.29 & 3.99 \end{matrix}]$	3.85 $[\begin{matrix} 2.62 & 5.63 \end{matrix}]$
Sc. power ( $v$ ) ( $10^{- 4} {kPa}^{- 1}$ )	3.09 $[\begin{matrix} 1.92 & 3.93 \end{matrix}]$	5.46 $[\begin{matrix} 2.72 & 8.59 \end{matrix}]$	4.74 $[\begin{matrix} 3.58 & 6.84 \end{matrix}]$
Score 1 ( $s_{1}$ ) ( $10^{- 5} {kPa}^{- 1}$ )	$- 3.67$ $[\begin{matrix} - 5.52 & - 1.07 \end{matrix}]$	$- 7.23$ $[\begin{matrix} - 8.62 & - 5.35 \end{matrix}]$	1.54 $[\begin{matrix} 0.42 & 2.57 \end{matrix}]$
EMSC coef. ( $k$ ) ( ${kPa}^{- 1}$ )	1.87 $[\begin{matrix} 0.39 & 2.74 \end{matrix}]$	1.99 $[\begin{matrix} 1.25 & 2.80 \end{matrix}]$	0.45 $[\begin{matrix} - 0.05 & 0.87 \end{matrix}]$

The highest reflectance decrease was observed above the veins, while the reflectance decrease observed above the muscle tissue was more than two times smaller (Fig. 5). The skin above the muscle tissue also exhibited the highest increase in absorption at 1450 nm, indicating the dominant impact of scattering on the reflectance. The absorption increase observed above the muscle tissue was attributed to the increase in the concentration of water and lipids. During the measurements, the muscle provided less resistance to the applied pressure; therefore, the skin under the probe sagged deeper. As a result, the relative change in the concentration of water (Fig. 5) and lipids was higher compared to the tissue above the veins and bone, where the contact pressure was distributed to a smaller tissue volume. The differences in the tissue structure are also responsible for the higher oxygen saturation drop (Fig. 5) observed above the veins and bone. The reduced scattering coefficient consistently decreased with the contact pressure, most prominently above the veins (Fig. 5). The scattering decrease could be explained by the decrease in the refractive index mismatch due to the water displacement in the dermis and epidermis.⁸

Fig. 5

Slopes ( $h_{1}$ ) of the selected features (Table 2) modeled as a polynomial function of the applied contact pressure for the three measurement sites: skin above the muscle (muscle), skin in the middle of the wrist crease above the veins (veins) and skin above the ulnar process (bone).

3.4.

Analysis of the PCA Loadings

A detailed analysis of the PCA loadings confirmed the findings based on the results obtained by modeling the tissue scattering and absorption properties and gave additional insight on the most prominent NIR absorption bands governing the pressure-induced spectral response. The first PCA loading of the EMSC preprocessed skin pressure-induced spectral response (shown in Fig. 6) exhibits three distinct absorption bands. Two absorption bands, centered at 1210 and 1170 nm, can be attributed to the absorption of lipid mixture,¹³ which was shown to have a significant absorption band centered at 1210 nm, and a minor absorption band centered at 1170 nm, both in excellent agreement with the absorption bands highlighted by the first PCA loading. With respect to the wavelength dependence of the reduced scattering coefficient and the EMSC preprocessing of the acquired spectra, it is reasonable to assume that the calculated loadings predominantly reflect the pressure-induced changes in the absorption properties of the underlying tissue, and thereby concentration of the chromophores. Under the assumption that the lipids cannot be pushed out of the tissue by the applied pressure, the concentration of lipids can change only as a result of the increase or decrease in the relative volume fractions of the chromophores. This supposition is clearly supported by the third, water-related¹¹ absorption band of the first PCA loading observed at 1450 nm. Additionally, a decrease in the value of the first PCA loading at the locations (depicted by the two inclined arrows in Fig. 6) was observed. According to Kuenstner and Norris,¹² different hemoglobin species have absorption bands located in the spectral range from 800 to 1100 nm and in the spectral range from 1500 to 1600 nm (oxyhemoglobin at 924 nm, deoxyhemoglobin at 903 nm, and methemoglobin at 1017 nm), all in excellent agreement with the two depicted locations.

Fig. 6

First PCA loading of the skin spectra (full line). Absorption spectrum of water (dash-dotted line) and absorption spectrum of lipids (dashed line from 900 to 1350 nm). The absorption bands of lipids (1170 nm and 1210 nm) and water (1450 nm) are depicted by the vertical arrows. The absorption bands of hemoglobin are depicted by the inclined arrows.

4. Conclusions

This study showed for the first time that pressure-induced spectral response of the soft tissue in the NIR spectral range can be used as a valuable source of additional information for tissue classification. Even though the measurements were acquired at three proximate skin locations, the subcutaneous tissue provided unique spectral response to the applied pressure. In general, the applied contact pressure increased the concentrations of chromophores (water, hemoglobin, and lipids), resulting in higher absorbance and lower diffuse reflectance and scattering. The proposed classification scheme, utilizing solely a pressure-induced spectral response by polynomial modeling of the tissue scattering and absorption properties, resulted in a high (90%) average classification sensitivity and specificity. Based on the presented results, the NIR, pressure-induced spectral response of the soft tissue indicates considerable potential for improving the sensitivity and specificity of soft tissue classification.

Acknowledgments

This work was supported by the Ministry of Higher Education, Science and Technology, Republic of Slovenia under Grants L2-4072, J7-2246, and L2-2023.

References

1.

V. R. KondepatiH. M. HeiseJ. Backhaus, “Recent applications of near-infrared spectroscopy in cancer diagnosis and therapy,” Anal. Bioanal. Chem., 390 (1), 125 –139 (2008). http://dx.doi.org/10.1007/s00216-007-1651-y ABCNBP 1618-2642 Google Scholar

2.

B. J. Tromberget al., “Imaging in breast cancer—Diffuse optics in breast cancer: detecting tumors in pre-menopausal women and monitoring neoadjuvant chemotherapy,” Breast Cancer Res., 7 (6), 279 –285 (2005). http://dx.doi.org/10.1186/bcr1358 BCTRD6 0167-6806 Google Scholar

3.

L. M. McIntoshet al., “Towards non-invasive screening of skin lesions by near infrared spectroscopy,” J. Invest. Dermatol., 116 (1), 175 –181 (2001). http://dx.doi.org/10.1046/j.1523-1747.2001.00212.x JIDEAE 0022-202X Google Scholar

4.

E. K. Chanet al., “Effects of compression on soft tissue optical properties,” IEEE J. Sel. Top. Quant., 2 (4), 943 –950 (1996). http://dx.doi.org/10.1109/2944.577320 IJSQEN 1077-260X Google Scholar

5.

W. L. Chenet al., “Influence of contact state on NIR diffuse reflectance spectroscopy in vivo,” J. Phys. D-Appl. Phys., 38 (15), 2691 –2695 (2005). http://dx.doi.org/10.1088/0022-3727/38/15/022 JPAPBE 0022-3727 Google Scholar

6.

J. A. D. Atencioet al., “Influence of probe pressure on human skin diffuse reflectance spectroscopy measurements,” Opt. Mem. Neural Networks, 18 (1), 6 –14 (2009). http://dx.doi.org/10.3103/S1060992X09010020 OMNNE8 1060-992X Google Scholar

7.

L. L. Randeberg, “Diagnostic applications of diffuse reflectance spectroscopy,” Norwegian University of Science and Technology, Faculty of Information Technology, Mathematics and Electrical Engineering, (2005). Google Scholar

8.

L. Limet al., “Probe pressure effects on human skin diffuse reflectance and fluorescence spectroscopy measurements,” J. Biomed. Opt., 16 (1), 011012 (2011). http://dx.doi.org/10.1117/1.3525288 JBOPFO 1083-3668 Google Scholar

9.

J. Sunet al., “Influence of fiber optic probe geometry on the applicability of inverse models of tissue reflectance spectroscopy: computational models and experimental measurements,” Appl. Opt., 45 (31), 8152 –8162 (2006). http://dx.doi.org/10.1364/AO.45.008152 APOPAI 0003-6935 Google Scholar

10.

A. KimB. C. Wilson, “Measurement of ex vivo and in vivo tissue optical properties: methods and theories,” Optical-Thermal Response of Laser-Irradiated Tissue, Springer, Netherlands (2011). Google Scholar

11.

D. J. Segelstein, “The complex refractive index of water,” University of Missouri, (1981). Google Scholar

12.

J. KuenstnerK. Norris, “Spectrophotometry of human hemoglobin in the near infrared region from 1000 to 2500 nm,” J. Near. Infrared Spec., 2 (1), 59 –65 (1994). http://dx.doi.org/10.1255/jnirs.32 JNISEI 0967-0335 Google Scholar

13.

C. L. TsaiJ. C. ChenW. J. Wang, “Near-infrared absorption property of biological soft tissue constituents,” J. Med. Biol. Eng., 21 (1), 7 –14 (2001). IYSEAK 0021-3292 Google Scholar

14.

AltshulerG. B.AndersonR. R.MansteinD., “Method and apparatus for the selective targeting of lipid-rich tissue,” U. S. Patent No. 6,605,080 B1.

15.

A. N. Bashkatovet al., “Optical properties of human skin, subcutaneous and mucous tissues in the wavelength range from 400 to 2000 nm,” J. Phys. D-Appl. Phys., 38 (15), 2543 (2005). http://dx.doi.org/10.1088/0022-3727/38/15/004 JPAPBE 0022-3727 Google Scholar

16.

H. MartensJ. P. NielsenS. B. Engelsen, “Light scattering and light absorbance separated by extended multiplicative signal correction. Application to near-infrared transmission analysis of powder mixtures,” Anal. Chem., 75 (3), 394 –404 (2003). http://dx.doi.org/10.1021/ac020194w ANCHAM 0003-2700 Google Scholar

17.

P. GeladiD. MacDougallH. Martens, “Linearization and scatter-correction for near-infrared reflectance spectra of meat,” Appl. Spectrosc., 39 (3), 491 –500 (1985). http://dx.doi.org/10.1366/0003702854248656 APSPA4 0003-7028 Google Scholar

18.

T. Næset al., User Friendly Guide to Multivariate Calibration and Classification, NIR Publications, Chichester (2002). Google Scholar

19.

P. PudilJ. NovovičováJ. Kittler, “Floating search methods in feature selection,” Pattern Recogn. Lett., 15 (11), 1119 –1125 (1994). http://dx.doi.org/10.1016/0167-8655(94)90127-9 PRLEDG 0167-8655 Google Scholar

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation Download Citation

Blaž Cugmas, Miran Bürmen, Maksimilijan Bregar, Franjo Pernuš, and Boštjan Likar "Pressure-induced near infrared spectra response as a valuable source of information for soft tissue classification," Journal of Biomedical Optics 18(4), 047002 (3 April 2013). https://doi.org/10.1117/1.JBO.18.4.047002

Published: 3 April 2013

Access the abstract

JOURNAL ARTICLE
8 PAGES

DOWNLOAD PAPER SAVE TO MY LIBRARY

GET CITATION

CITATIONS

Cited by 14 scholarly publications and 5 patents.

Explore citations on Lens.org

KEYWORDS

Tissues

Skin

Absorption

Scattering

Reflectivity

Near infrared spectroscopy

Principal component analysis

1.

Introduction

Fig. 1

2.

Materials and Methods

2.1.

Samples

2.2.

Instrumentation

Fig. 2

2.3.

Measurements

2.4.

Spectra Processing

Eq. (1)

Fig. 3

Eq. (2)

Eq. (3)

Eq. (4)

Eq. (5)

Eq. (6)

Eq. (7)

3.

Results and Discussion

3.1.

General Findings

Fig. 4

3.2.

Classification

Table 1

3.3.

Tissue Scattering and Absorption Properties as a Linear Function of the Applied Contact Pressure

Table 2

Fig. 5

3.4.

Analysis of the PCA Loadings

Fig. 6

4.

Conclusions

Acknowledgments

References

Show All Keywords

Keywords/Phrases

Search In:

Publication Years