5 November 2013 Assessment of human burn scars with optical coherence tomography by imaging the attenuation coefficient of tissue after vascular masking
Author Affiliations +
Abstract
The formation of burn-scar tissue in human skin profoundly alters, among other things, the structure of the dermis. We present a method to characterize dermal scar tissue by the measurement of the near-infrared attenuation coefficient using optical coherence tomography (OCT). To generate accurate en face parametric images of attenuation, we found it critical to first identify (using speckle decorrelation) and mask the tissue vasculature from the three-dimensional OCT data. The resulting attenuation coefficients in the vasculature-masked regions of the dermis of human burn-scar patients are lower in hypertrophic (3.8±0.4  mm −1 ) and normotrophic (4.2±0.9  mm −1 ) scars than in contralateral or adjacent normal skin (6.3±0.5  mm −1 ). Our results suggest that the attenuation coefficient of vasculature-masked tissue could be used as an objective means to assess human burn scars.

1.

Introduction

Wound healing after burn injury involves the sequential phases of inflammation, tissue formation, and tissue remodeling.1 Scar formation during this process is a natural facet of healing, and hypertrophic scars are one pathology arising from an excessive response to injury during healing.2 Such scars are characterized by prolific neovascularization; and excessive content34.5 and abnormal arrangement of collagen.6 In normal skin, collagen is arranged in a “basket-weave”-like structure, whereas collagen in hypertrophic scars is more unidirectionally aligned.6,7 Previous studies have shown that, compared with normal skin, hypertrophic scars contain collagen fibrils with reduced diameter6,8 and have higher water content.9,10

Various techniques are available or under investigation for the treatment of pathological scars including surgery, pressure therapy, corticosteroid injections, silicone gel sheeting, radiotherapy, and laser therapy.11 Evaluating the efficacy of these clinical treatments typically requires the longitudinal assessment of a scar’s severity. Current scar assessment methods are primarily based on the direct observation and palpation of the features of scars. For example, the Vancouver Scar Scale uses observer assessment of pigmentation, vascularity, pliability, and height to characterize the scar.12 Assessment with such methods is inevitably somewhat subjective, leading to high inter-observer variability. There is reported research on alternative, less subjective techniques. Electron microscopy has been used to study the characteristic morphological profiles of pathological scars but requires histological slices of excised scar tissue.6,13 Multiphoton14,15 and confocal microscopies16 have been applied to scar tissue or wounds, ex vivo and in vivo, including for high-resolution imaging of collagen and elastic fiber components. Ultrasound imaging has been investigated as a means of assessing the thickness of burn scars in pediatric patients noninvasively.17 Other techniques, such as laser Doppler perfusion imaging and laser speckle perfusion imaging, have used the degree of vascularization as a surrogate indicator of scar status.18

Changes in tissue morphology in pathological scarring can result in changes in the optical properties of the tissue in the near-infrared wavelength range. This has driven a growing body of work exploring the use of optical coherence tomography (OCT) in scar assessment.1920.21.22.23.24.25.26 OCT is a noninvasive imaging modality capable of acquiring three-dimensional (3-D) scans of tissue at resolutions from 1 to 20 μm from superficial tissues at depths of up to 1 to 2 mm.19 The morphological features of normal and diseased skin have been studied in vivo with OCT imaging.20,21 Wound healing has been evaluated noninvasively by high-resolution OCT,22 and polarization-sensitive (PS) OCT has been used to measure the birefringence of both normal and burn-damaged human skin.23,24 OCT has also been used to noninvasively assess vascularity as a measure of normal skin and scar status.25,26

For objective scar assessment, absolute quantities that characterize the tissue are required to monitor the progression of healing. One such quantity that may be extracted from OCT scans is the attenuation rate (parameterized by the attenuation coefficient), which is the rate at which the OCT signal decreases with depth in the tissue. This rate can be extracted by fitting the measured OCT signal versus depth (an A-scan) to a specific model.2728.29 Previous studies have shown that the attenuation coefficient measured with OCT can be used to distinguish between different types of tissue.3031.32.33 Previous research has also explored the attenuation coefficient of normal human skin tissue.27,34,35 A very effective means of visualizing attenuation coefficients is through the presentation of en face maps, i.e., parametric images of the attenuation coefficient.30,31 We hypothesize that differences in tissue microstructure between scarred and normal skin give rise to differences in the attenuation coefficient of near-infrared light. In this article, we extend the measurement of OCT attenuation coefficients and their presentation in parametric images to the case of pathological scarred tissue.

A difficulty arises when applying the attenuation coefficient method to in vivo OCT scans of skin. The vasculature of scarred and normal skin generates strong forward scattering at near-infrared wavelengths due, in part, to the high anisotropy and high scattering cross-section of red blood cells,36 causing artificially high attenuation of the OCT signal that is not representative of the surrounding nonvascular component of the tissue. Here, we present an algorithm that automatically removes areas of vasculature from the calculations, allowing more accurate quantification of the attenuation coefficient of scarred and normal skin tissue. We present results from six clinically assessed scars that demonstrate a systematically lower attenuation coefficient of scarred tissue as compared with normal skin tissue and propose a mechanism explaining this observation.

2.

Materials and Methods

2.1.

OCT Imaging

A polarization-sensitive swept-source OCT system (PSOCT-1300, Thorlabs, Newton, New Jersy) with a central wavelength of 1325 nm was utilized to perform the scanning on patients. The measured axial and transverse resolutions (full-width-at-half-maximum) of the system are, respectively, 17 μm (in air) and 16 μm. Light illuminated the skin through an objective lens with a working distance of 25 mm to deliver a scanning beam with a numerical aperture of 0.03 and a measured optical power of 3.2 mW to the skin.

The tissue-imaging setup has been described in detail previously.26 Briefly, a 10×10mm outline was drawn on the skin to guide the attachment of a metal fiducial marker, which had a central hole (3-mm diameter) for imaging. The marker was used as a guide to remove motion artifacts, as previously described.37 An adjustable spacer between the OCT scanner and the skin was used to maintain a suitable distance from the objective lens of the OCT scanner to the skin surface. A layer of ultrasound gel was applied to reduce the refractive index mismatch at the skin surface, thereby reducing imaging artifacts and enhancing the imaging depth.38 The scanner was positioned on the skin and handheld during data acquisition.

Scans were acquired on six patients (age: 19 to 58 years, three females, three males) undergoing follow-up examination 12 or 24 months after initial treatment for burns. The scanning protocol was approved by the Human Research Ethics Committee of Royal Perth Hospital and The University of Western Australia. Written consent was obtained from all patients prior to scanning. For each patient, locations on the scar and on the contralateral or adjacent normal skin region were selected for 3-D OCT imaging. Each OCT data volume covered a region of 4×1.5×3mm (1088×1088×512pixels) in the x, y (lateral), and z (axial into the skin) directions, respectively, acquired within the center of the outlined region, including a portion of the fiducial marker to enable motion correction. With these scanning parameters, over-sampling was used in order to record the speckle decorrelation for vasculature segmentation. The system was operated at an A-scan rate of 4.9kHz, leading to a B-scan rate of 4.5 frames per second and 3-D data volume acquisition time of 4min.

2.2.

Vasculature Masking

The flow chart for vasculature masking is shown in the left column of Fig. 1. We denote by i2 the sum of the squared OCT signals from the two detection channels of the PS-OCT system. This is equivalent to the square of a conventional OCT signal and is sometimes also referred to as the OCT signal power. In the remainder of this article, the term OCT data refers to i2. Blood vessels were automatically segmented from this OCT data using speckle decorrelation.39 In brief, the correlation between OCT data within a small moving window was calculated between pairs of adjacent B-scans, A and B, using the standard normalized cross-correlation equation:

(1)

corrMap(x,z)=p=1Mq=1N{[iA2(x+p,z+q)iA2(x,z)¯][iB2(x+p,z+q)iB2(x,z)¯]}p=1Mq=1N[iA2(x+p,z+q)iA2(x,z)¯]2p=1Mq=1N[iB2(x+p,z+q)iB2(x,z)¯]2,
where the window size is M×N, and iA2 and iB2 are OCT data from a pair of neighboring subimages in adjacent B-scans in an OCT data volume that are defined by the window. iA2¯ and iB2¯ are, respectively, the mean values of iA2 and iB2.

Fig. 1

Data processing flow for the calculation of the attenuation coefficient of vasculature-masked tissue in scarred and normal skin (see text for details). MIP: maximum intensity projection.

JBO_19_2_021111_f001.png

Regions of low cross-correlation in the resulting correlation volume image were identified as blood vessels. Motion artifacts caused by patient movement during scanning were then removed from this correlation volume by the use of a registration algorithm utilizing the metal fiducial marker. Details of the motion correction are described elsewhere.37 By identifying the voxel of the lowest correlation (i.e., the highest decorrelation) along each A-scan (z-direction) at depths in the skin with adequately high OCT signal-to-noise ratio, a two-dimensional (2-D) en face maximum intensity projection (MIP) image of the vasculature was produced, showing the (x,y) locations of blood vessels. Further details of the algorithm were given previously.26 This MIP image was then thresholded to form a binary mask of blood vessels, indicating (x,y) locations to be removed from subsequent calculation of the attenuation coefficient.

2.3.

Calculation of the Attenuation Coefficient

Assuming that single scattering is the dominant scattering process, the rate of attenuation of the OCT signal with depth in a homogeneous sample follows a negative exponential function.27 In addition, this signal is modulated by the confocal function of the focusing optics, and the sensitivity fall-off of our swept source OCT system. The skin OCT data, i2, as a function of depth z can be expressed as follows:

(2)

i2(z)F(z)S(z)ρe2μtz,
where F(z) and S(z) are, respectively, the confocal and sensitivity fall-off functions, and ρe2μtz represents the single-scattering component, where ρ is the initial value of the reflectance and μt is the attenuation coefficient.31

To extract the attenuation coefficient, we first corrected our in vivo scans for the confocal function and sensitivity fall-off of the system using a calibration scan of a low-scattering suspension of polystyrene microspheres (Polybead®, Polysciences, Inc., Warrington, Pennsylvania), as described by Scolaro et al.31 The calibration OCT data can be written as

(3)

i02(z)F(z)S(z)ρ0e2μt0z,
where ρ0 is the initial value of the reflectance, and μt0 is the attenuation coefficient of the calibration suspension. μt0 was calculated to be 0.1mm1 at our OCT wavelength using Mie theory40 and assuming negligible absorption. Each voxel value of the skin scans was corrected by division by the corresponding voxel value in the calibration scan volume. After taking the logarithm of the corrected data, the resulting logarithmic OCT data at each (x,y) location as a function of depth (z) is given by

(4)

ln[i2(z)i02(z)]=2(μtμt0)z+a,
where a is a constant for each A-scan dependent on the ratio ρ/ρ0.

The aforementioned fiducial marker-based registration algorithm was applied to the corrected OCT data to reduce the motion artifacts and to maintain the spatial consistency between the vasculature mask and the corrected OCT data. A Canny edge detector41 was used to identify the surface of the skin. For each dataset, the average depth to the dermis from the tissue surface within the scan field-of-view was empirically estimated. Speckle was reduced using a moving window averaging function with a uniform kernel of 40×40μm within the xy plane over regions not masked out. The attenuation coefficient μt was calculated from a linear least-squares fit to the averaged logarithmic OCT data over a depth of 200 μm from the average depth into the dermis [taking the refractive index of skin to be 1.43 (Ref. 42) to correct for the difference between physical and optical path length]. The goodness-of-fit of the least-squares regression was calculated to assess the fitting quality.

After performing the calculation on each A-scan at vasculature-masked locations covered by the central 1×1mm (x×y) region, a 2-D image in the xy plane was generated with the attenuation coefficient represented graphically as the pixel intensity for each location. In addition, a one-dimensional normalized histogram was computed, comprising all attenuation values in the image (not including locations masked by the vasculature). Each histogram contained typically 100,000 to 170,000 values, corresponding to the number of vasculature-free A-scans in each dataset. The histogram shows the distribution of attenuation coefficients for scar tissue and the contralateral or adjacent normal skin tissue. All data processing was performed using MATLAB® (vR2012a, The MathWorks, Inc., Natick, Massachusetts).

3.

Results

3.1.

Vasculature-Induced Artifacts

Figure 2 shows the calculated attenuation coefficient and related data for normal skin on the lower right leg of a 19-year-old Caucasian male patient. The OCT scanning area (4×1.5mm) was in the center of the region marked with a square outline (10×10mm) on the skin, as shown in Fig. 2(a). Figure 2(b) shows the 2-D en face vasculature MIP of the scanning region. The blue square (1×1mm) delineates the region over which the attenuation coefficient map in Fig. 2(c) was extracted using the corrected, registered, and averaged structural OCT data. The attenuation coefficient determined for all (vascular and nonvascular) tissue regions is shown. The corresponding goodness-of-fit (0poor fitting, 1good fitting) is presented in Fig. 2(d). By comparing Figs. 2(b) and 2(c), it is apparent that the attenuation coefficient values in the circled vascular regions are much lower than in the surrounding tissues, even giving rise to nonphysical negative values. Figure 2(e) compares the normalized attenuation coefficient distribution between the vascular (i.e., blood vessels) and nonvascular (i.e., all tissue excepting blood vessels) regions and highlights this problem. The attenuation coefficient distribution in the vascular region is artificially broadened and extended to low values.

Fig. 2

Attenuation coefficient of normal skin showing the corrupting effect of the vasculature: (a) Photograph and (b) vasculature MIP of normal skin. (c, d) En face attenuation coefficient map and corresponding goodness-of-fit of the vascular and nonvascular tissues in the region outlined in blue in (b) (each image 1×1mm). Examples of vascular regions are circled by the dotted lines in (b), (c), and (d). (e) Comparison of normalized attenuation coefficient distribution between vascular and nonvascular tissue regions. (f, g) Fitting examples, respectively, from nonvascular and vascular regions marked by the small purple squares and indicated by arrows in (c).

JBO_19_2_021111_f002.png

Figures 2(f) and 2(g) show two representative fitting examples, respectively, for nonvascular and vascular tissues, marked by the two small purple squares in Fig. 2(c). The region of OCT signal used for the fitting (200 μm in length) is shown in blue, while the remainder of the A-scan is rendered partially transparent. The attenuation coefficient is extracted from the line of best fit (red) for each A-scan. Figure 2(f) demonstrates the applicability of the single-scattering model to the nonvascular dermal tissue with a high goodness-of-fit value of 0.91. In Fig. 2(g), the presence of two blood vessels at different depths gives rise to the sharp drops in the observed signal. The subsequent increase in the measured OCT signal below the first blood vessel gives rise to an apparent negative attenuation coefficient value with a poor goodness-of-fit value of 0.06. It is also evident from Fig. 2(d) that the low attenuation in the outlined vascular regions is accompanied by poor goodness-of-fit to the single-scattering, log-linear model. Since scarred skin has been shown to have larger blood vessel diameters and higher vasculature density than normal skin,26 this vascular artifact is expected to be still more prominent in scar tissue. Consequently, we incorporated a vasculature-masking step in the determination of the attenuation coefficient in order to minimize this corrupting effect.

3.2.

Case Study of a Hypertrophic Scar

Figure 3 shows a 12-month-old hypertrophic scar originating from a hot water burn on the medial left inner forearm of a 30-year-old Caucasian female patient. Figures 3(a) and 3(b) are photographs of the contralateral normal skin and scar tissue, respectively. OCT scans were performed in the central 4×1.5mm of the outlined regions. Figures 3(c) and 3(d) show the vasculature MIPs for the normal skin and scar, respectively, with the central blue squares (1×1mm) identifying the regions for attenuation coefficient calculation. As noted by Liew et al.,26 the vasculature is observed to be much more prolific in the scar tissue. Figures 3(e) and 3(g) depict the masks generated (black: vascular region; white: nonvascular region) after thresholding the vasculature MIPs. Figures 3(f) and 3(h) present the 2-D en face attenuation coefficient maps, respectively, for normal skin and scar. The attenuation coefficient is visibly lower in the hypertrophic scar than in the contralateral normal skin. From the distribution of attenuation coefficient values for each scan plotted in Fig. 3(i), we observe that the attenuation coefficients in scar tissue have a distribution of similar width to normal skin tissue but are on average much lower, with the mean attenuation coefficients of scar and normal skin being 3.8 and 6.4mm1, respectively.

Fig. 3

Case study of the attenuation coefficient of a hypertrophic scar: (a, b) Photographs of the contralateral normal skin and scar. (c, d) Vasculature MIPs of the normal skin and scar in (a) and (b), respectively. (e, g) Vasculature masks for normal skin and scar (1×1mm), respectively. (f, h) En face attenuation coefficient (mm1) maps of, respectively, normal skin and scar (1×1mm). (i) Histogram of attenuation coefficients for normal skin and scar.

JBO_19_2_021111_f003.png

3.3.

Case Study of a Normotrophic Scar

The results of a case study on a 12-month-old normotrophic scar are shown in Fig. 4. The scar, outlined in red in Fig. 4(a), was caused by a wood fire burn on the left lateral forearm of a 58-year-old Caucasian female. Figures 4(b) and 4(c) show MIPs of the vasculature of the contralateral normal skin (photograph not shown) and the scar, respectively. After masking of the vasculature using Figs. 4(d) and 4(f), the calculated attenuation coefficient maps are shown in Figs. 4(e) and 4(g). Figure 4(h) plots the normalized attenuation coefficient distribution of the scar and normal skin, which have mean values of 5.4 and 6.7  mm1, respectively. Note the increased degree of overlap of the distributions between this (nonpathological) normotrophic scar (red) and the contralateral normal skin (green), in comparison to the (pathological) hypertrophic scar in Fig. 3(i).

Fig. 4

Case study of the attenuation coefficient of a normotrophic scar. (a) Photograph of the scar (which is outlined in dotted red). (b, c) Vasculature MIPs of the normal skin and scar, respectively. (d, f) Vasculature masks for normal skin and scar (1×1mm), respectively. (e, g) En face attenuation coefficient (mm1) maps of, respectively, normal skin and scar (1×1mm). (h) Histogram of attenuation coefficients for normal skin and scar.

JBO_19_2_021111_f004.png

3.4.

Results from All Subjects

The mean (±standard error) attenuation coefficient of the vasculature-masked dermal tissue in hypertrophic scars (n=3), normotrophic scars (n=3), and normal skin (n=6), calculated from data on all subjects, is summarized in Fig. 5. The hypertrophic scars give rise to the lowest attenuation coefficient (3.8±0.4mm1); the normotrophic scars values are greater (4.2±0.9mm1); and the contralateral or adjacent normal skin has the largest values (6.3±0.5mm1). Student’s t-test shows that the difference between the attenuation coefficients of the scar and normal skin is statistically significant (p<0.001).

Fig. 5

Mean (±standard error) attenuation coefficient of hypertrophic scars (n=3), normotrophic scars (n=3), and normal skin (n=6) of all subjects.

JBO_19_2_021111_f005.png

4.

Discussion

In this study, we specifically focused on characterization of the dermis, since it undergoes very significant changes in content and structure during the wound healing process2 and occupies most of the OCT imaging depth range. As the epidermis is only about 100-μm thick and is largely restored with new epidermal tissue during the wound healing process, we chose to exclude it. Our results show on average 36% lower dermal attenuation in scarred skin compared with normal skin. This strong effect can potentially be used as the basis for an objective method of scar characterization. As pathological scarring commonly leads to excessive production of both collagen fibers and blood vessels, the present method complements the scar assessment method solely based on vasculature that we have recently reported.26 The integration of both methods could potentially be used to track the efficacy of various treatment strategies targeting both collagen fibers and blood vessels.

To explain the contrast that we observe in the vasculature-masked tissue, we have employed a model proposed by Jacques which has been shown to accurately describe the attenuation of light in the dermis.43 This model was verified against measurements of the reduced scattering coefficient, (1g)μs, where μs is the scattering coefficient and g is the anisotropy factor, of ex vivo samples of dermis.44 Jacques’ model treats scattering in the dermis as being equivalent to single scattering by a collection of infinite cylinders, each of which represents a collagen fiber. Jacques also considered the contribution of Rayleigh scattering by subwavelength scatterers, but it has been shown that Rayleigh scattering is insignificant compared with scattering by the collagen fibers at the wavelength of interest.44 The key parameters in the model are, thus, the refractive index of the collagen fibers and the background material, the distribution of fiber diameters, and the volume fraction of collagen fibers. The refractive index of the collagen fibers was calculated by assuming that the fibers are composed of a combination of water and collagen. Mie theory, describing the scattering of light by infinite cylinders, was used to calculate μs and g.

In relating Jacques’ model to the mathematical form of the OCT signal in Eq. (2), we note that μt is proportional to μsa(g), where a(g) represents the directional nature of scattering.45 We, however, found that g varies insignificantly compared with μs when contrasting normal skin and scarred tissue, enabling a(g) to be neglected when assessing the relative magnitudes of μt for the normal and scar cases. We did not have access to data on the distribution of fiber diameters for the particular cases presented in this article. We, thus, approximate the two ensembles of fibers, i.e., scar and normal, by their mean diameter and spacing from the literature.7,46 In particular, the mean spacing of collagen fibers has been found to be 10.37 μm and 8.53 μm for normal and scar tissues, respectively, while the mean diameter was found to be 5.39 and 5.77 μm, respectively. Thus, if all other parameters remained equal, the higher density of fibers would result in μs for scar tissue being larger than that for normal tissue. It has, however, been noted that the percentage of water content is greater in hypertrophic scar tissue than in normal tissue (approximately 64% in normal skin).9,10 The water content of hypertrophic scar tissue has been found to vary with scar age and, in one study, varied monotonically between 85.1% for a 0.5-year-old scar to 63.8% for a 2-year-old scar.9 Another study noted that the water content of hypertrophic scar was 82%, although the age of the scar was not noted.10 We cannot suppose that this water content is identical to that likely to be found in the collagen fibers. We note, however, that Jacques43 assumed a water content of 65% in his model, which is in agreement with the measured figure for normal skin tissue.

Our approach was, thus, to calculate the ratio μsscar/μsnormal, where μsscar and μsnormal are the scattering coefficients for scar and normal skin tissues, respectively, for different percentages of water content in scar collagen fibers, while the water content percentage of collagen fibers in normal skin was fixed at 64% consistent with the literature. We used refractive indices for water and collagen of 1.3241 and 1.4367 at a wavelength of 1325 nm, respectively.47 The scattering coefficients were calculated using Mie theory for scattering by cylinders, at normal incidence, using the mathematical formulation presented by Bohren and Huffman.40 Figure 6 shows the simulated ratio of scattering coefficient between scar and normal skin, μsscar/μsnormal, as a function of water content in the scar collagen fibers. The green triangle and red circle highlight simulated scattering coefficient ratios between scar and normal skin of particular interest: respectively, for the water content percentage of 64% (2-year-old hypertrophic scars) and 85% (0.5-year-old hypertrophic scars).9 The purple square corresponds to the ratio derived from our calculated attenuation coefficients of the scars and normal skin, and the resulting inferred water content (76%) lies between that of the 0.5- and 2-year-old hypertrophic scars. This level is consistent with the 1.2-year-old mean age of the scars in this study but not with the lower reported water content of 1-year-old hypertrophic scars (66.5%).9 Although requiring further study, our model and previous reports suggest that the contrast in attenuation coefficient may come from higher water content in scars. Thus, results of the type plotted in Fig. 6 could be used to study the water content in scar tissue noninvasively from measurement of the OCT attenuation coefficient ratio.

Fig. 6

Simulated μsscar/μsnormal ratios versus water content in scar collagen fibers. The triangle and circle represent water content values reported in the literature,9 and the square represents the value calculated from the data in this study.

JBO_19_2_021111_f006.png

The masking of blood vessels removed a significant artifact from the calculation of the attenuation coefficient caused by the strong scattering by blood vessels. To the best of our knowledge, this is the first reported incorporation of vasculature masking into the calculation of attenuation coefficients from OCT data. Blood comprises plasma and other constituents, among which erythrocytes dominate the optical properties. Erythrocytes have a biconcave, disk-like shape with a size considerably larger than our OCT wavelength, which leads to strong forward scattering. This, together with the high-scattering cross-section and strong absorption of light,36 generates a strong drop in the OCT signal, causing blood to exhibit a much higher attenuation coefficient than the nonvascular tissue in the skin. A previous study on fresh porcine blood, using 1300 nm OCT light, reported an attenuation coefficient of 12.15mm1,48 which is almost twice the mean value of the vasculature-masked skin tissue measured in this study. This is particularly important in the evaluation of hypertrophic scars, which are known to have enhanced vasculature.26 Analysis of the attenuation in these masked regions of blood flow may potentially provide additional information on metabolism such as the noninvasive monitoring of blood glucose concentration.49

Skin comprises a collection of microstructures with varying scattering and absorption properties. Quantification of scar and normal skin attenuation coefficients is complicated by tissue heterogeneity. Even after masking out the vasculature, heterogeneity in depth can invalidate the single-scattering exponential decay model and results in modified attenuation coefficient values, even negative values when the fitting window initiates in a low-signal region and extends into a high-signal region. We empirically found a fitting range of 200 μm in the dermis to be an appropriate compromise between including sufficient data in the estimate of the attenuation coefficient and minimizing the effects of tissue heterogeneity. The average attenuation coefficient of normal skin we determined using this fitting length is 6.3±0.5mm1, which lies between the results reported by Schmitt et al.27 (4.6 to 4.7mm1) and Kholodnykh et al.34 (10 to 13mm1) for forearm at 1300-nm wavelength. As discussed by Kholodnykh et al.,34 the variation in attenuation coefficient could be caused by factors such as differences in experimental protocol, the use of a clearing agent, and variations in the skin tissue among individuals. One different feature of our experimental protocol from that of others is that we removed the contribution of blood vessels to target the nonvascular tissue.

In vivo OCT skin scanning can be affected by the motion between the scanner head and the skin. Such motion introduces image artifacts that distort the structural OCT images.37 We used a fiducial marker fixed to the skin to guide a feature-based registration algorithm to remove such artifacts. However, motion artifacts also artificially increase speckle decorrelation, resulting in errors in blood vessel segmentation. An example of such an artifact appears as a pattern of horizontal lines in the vasculature MIP of Fig. 2(b). Alternative solutions have been proposed, such as rigidly affixing the scan head to the patient’s skin,50,51 although this can present additional difficulties when scanning acute burn wounds. Another feasible method is to reduce the 3-D data acquisition time, although this may reduce the sensitivity of vessel detection in areas of slow blood flow due to reduced speckle decorrelation between rapidly acquired B-scans.

5.

Conclusions

In this article, we have presented a method for in vivo assessment of near-infrared attenuation coefficients of the dermis using OCT. We corrected for artifact caused by the presence of blood vessels using an automated method to segment and mask vasculature. We demonstrated the method of assessment on several in vivo human burn scars. Our results indicate that scars of mean age 1.2 years have on average 36% lower attenuation coefficients than contralateral or adjacent normal skin, and our modeling suggests that the lower attenuation coefficients arise from the previously reported increased water content in scars. This method could potentially be integrated with an automated vascularity quantification method in an objective scar assessment protocol to characterize both the vasculature and the connective tissue in scars. Future work will involve the investigation of this objective assessment protocol for the longitudinal study of human burn scars.

Acknowledgments

The authors thank Loretta Scolaro, Lixin Chin, Blake R. Klyen, Rodney W. Kirk, Dirk Lorenser, and Xiaojie Yang for discussions regarding the attenuation coefficient calculation and OCT data processing. The authors also thank Sharon Rowe at the Telstra Burn Reconstruction and Rehabilitation Unit, Royal Perth Hospital, for help in patient management. P. Gong acknowledges a joint scholarship from The University of Western Australia and the China Scholarship Council. R. A. McLaughlin is supported by a fellowship from Cancer Council Western Australia. P. R. T. Munro is supported by a Discovery Early Career Research Award from the Australian Research Council (DE120101331). This study is partially supported by the National Breast Cancer Foundation, Australia. Y. M. Liew was supported by a SIRF and UIS scholarship from The University of Western Australia.

References

1. A. J. SingerR. A. F. Clark, “Cutaneous wound healing,” N. Engl. J. Med. 341(10), 738–746 (1999).NEJMAG0028-4793 http://dx.doi.org/10.1056/NEJM199909023411006 Google Scholar

2. D. T. NguyenD. P. OrgillG. F. Murphy, “The pathophysiologic basis for wound healing and cutaneous regeneration,” Chapter 4 in Biomaterials for Treating Skin Loss, D. OrgillC. Blanco, Eds., pp. 25–57, Woodhead Publishing Limited, Cambridge, England (2009). Google Scholar

3. M. B. WitteA. Barbul, “General principles of wound healing,” Surg. Clin. North Am. 77(3), 509–528 (1997).SCNAA70039-6109 http://dx.doi.org/10.1016/S0039-6109(05)70566-1 Google Scholar

4. F. B. Niessenet al., “On the nature of hypertrophic scars and keloids: a review,” Plast. Reconstr. Surg. 104(5), 1435–1458 (1999).PRSUAS0032-1052 http://dx.doi.org/10.1097/00006534-199910000-00031 Google Scholar

5. T. L. TuanL. S. Nichter, “The molecular basis of keloid and hypertrophic scar formation,” Mol. Med. Today 4(1), 19–24 (1998).MMTOFK1357-4310 http://dx.doi.org/10.1016/S1357-4310(97)80541-2 Google Scholar

6. H. A. Linareset al., “The histiotypic organization of the hypertrophic scar in humans,” J. Invest. Dermatol. 59(4), 323–331 (1972).JIDEAE0022-202X http://dx.doi.org/10.1111/jid.1972.59.issue-4 Google Scholar

7. P. D. H. M. Verhaegenet al., “Differences in collagen architecture between keloid, hypertrophic scar, normotrophic scar, and normal skin: an objective histopathological analysis,” Wound Repair Regen. 17(5), 649–656 (2009).WREREU1067-1927 http://dx.doi.org/10.1111/wrr.2009.17.issue-5 Google Scholar

8. K. J. Stewart, “A quantitative ultrastructural study of collagen fibrils in human skin, normal scars, and hypertrophic scars,” Clin. Anat. 8(5), 334–338 (1995).CLANE81098-2353 http://dx.doi.org/10.1002/(ISSN)1098-2353 Google Scholar

9. M. Babuet al., “Differentiation of keloid and hypertrophic scar; correlation of the water proton relaxation times with the duration of the scar,” Physiol. Chem. Phys. Med. NMR 25(2), 113–120 (1993).PCPNER0748-6642 Google Scholar

10. C. W. KischerM. R. Shetlar, “Collagen and mucopolysaccharides in the hypertrophic scar,” Connect. Tissue Res. 2(3), 205–213 (1974).CVTRBC0300-8207 http://dx.doi.org/10.3109/03008207409152245 Google Scholar

11. T. A. Mustoeet al., “International clinical recommendations on scar management,” Plast. Reconstr. Surg. 110(2), 560–571 (2002).PRSUAS0032-1052 http://dx.doi.org/10.1097/00006534-200208000-00031 Google Scholar

12. Z. Tyacket al., “A systematic review of the quality of burn scar rating scales for clinical and research use,” Burns 38(1), 6–18 (2012).BURND80305-4179 http://dx.doi.org/10.1016/j.burns.2011.09.021 Google Scholar

13. M. Santucciet al., “Keloids and hypertrophic scars of Caucasians show distinctive morphologic and immunophenotypic profiles,” Virchows Archiv 438(5), 457–463 (2001).VARCEM0945-6317 http://dx.doi.org/10.1007/s004280000335 Google Scholar

14. J. Chenet al., “Multiphoton microscopy study of the morphological and quantity changes of collagen and elastic fiber components in keloid disease,” J. Biomed. Opt. 16(5), 051305 (2011).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.3569617 Google Scholar

15. V. Da Costaet al., “Nondestructive imaging of live human keloid and facial tissue using multiphoton microscopy,” Arch. Facial Plast. Surg. 10(1), 38–43 (2008).1521-2491 http://dx.doi.org/10.1001/archfacial.2007.18 Google Scholar

16. S. Lange-Asschenfeldtet al., “Applicability of confocal laser scanning microscopy for evaluation and monitoring of cutaneous wound healing,” J. Biomed. Opt. 17(7), 076016 (2012).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.17.7.076016 Google Scholar

17. X. Q. Wanget al., “Ultrasound assessed thickness of burn scars in association with laser Doppler imaging determined depth of burns in paediatric patients,” Burns 36(8), 1254–1262 (2010).BURND80305-4179 http://dx.doi.org/10.1016/j.burns.2010.05.018 Google Scholar

18. C. J. Stewartet al., “A comparison of two laser-based methods for determination of burn scar perfusion: Laser Doppler versus laser speckle imaging,” Burns 31(6), 744–752 (2005).BURND80305-4179 http://dx.doi.org/10.1016/j.burns.2005.04.004 Google Scholar

19. W. DrexlerJ. G. Fujimoto, Optical Coherence Tomography: Technology and Applications, Springer, Berlin (2008). Google Scholar

20. J. Welzelet al., “Optical coherence tomography of the human skin,” J. Am. Acad. Dermatol. 37(6), 958–963 (1997).JAADDB0190-9622 http://dx.doi.org/10.1016/S0190-9622(97)70072-0 Google Scholar

21. J. Welzel, “Optical coherence tomography in dermatology: a review,” Skin Res. Technol. 7(1), 1–9 (2001).0909-752X http://dx.doi.org/10.1034/j.1600-0846.2001.007001001.x Google Scholar

22. M. J. Cobbet al., “Noninvasive assessment of cutaneous wound healing using ultrahigh-resolution optical coherence tomography,” J. Biomed. Opt. 11(6), 064002 (2006).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.2388152 Google Scholar

23. M. C. Pierceet al., “Birefringence measurements in human skin using polarization-sensitive optical coherence tomography,” J. Biomed. Opt. 9(2), 287–291 (2004).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.1645797 Google Scholar

24. K. H. Kimet al., “In vivo imaging of human burn injuries with polarization-sensitive optical coherence tomography,” J. Biomed. Opt. 17(6), 066012 (2012).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.17.6.066012 Google Scholar

25. F. Bazant-Hegemarket al., “Multi-beam resolution video-rate swept-source optical coherence tomography (OCT) provides endogenous contrast for in vivo blood flow independent of flow direction,” Proc. SPIE 7554, 75542Z (2010).PSISDG0277-786X http://dx.doi.org/10.1117/12.841906 Google Scholar

26. Y. M. Liewet al., “In vivo assessment of human burn scars through automated quantification of vascularity using optical coherence tomography,” J. Biomed. Opt. 18(6), 061213 (2013).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.18.6.061213 Google Scholar

27. J. M. SchmittA. KnüttelR. F. Bonner, “Measurement of optical properties of biological tissues by low-coherence reflectometry,” Appl. Opt. 32(30), 6032–6042 (1993).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.32.006032 Google Scholar

28. J. M. SchmittA. Knüttel, “Model of optical coherence tomography of heterogeneous tissue,” J. Opt. Soc. Am. A 14(6), 1231–1242 (1997).JOAOD60740-3232 http://dx.doi.org/10.1364/JOSAA.14.001231 Google Scholar

29. L. ThraneH. T. YuraP. E. Andersen, “Analysis of optical coherence tomography systems based on the extended Huygens-Fresnel principle,” J. Opt. Soc. Am. A 17(3), 484–490 (2000).JOAOD60740-3232 http://dx.doi.org/10.1364/JOSAA.17.000484 Google Scholar

30. R. A. McLaughlinet al., “Parametric imaging of cancer with optical coherence tomography,” J. Biomed. Opt. 15(4), 046029 (2010).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.3479931 Google Scholar

31. L. Scolaroet al., “Parametric imaging of the local attenuation coefficient in human axillary lymph nodes assessed using optical coherence tomography,” Biomed. Opt. Express 3(2), 366–379 (2012).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.3.000366 Google Scholar

32. C. Xuet al., “Characterization of atherosclerosis plaques by measuring both backscattering and attenuation coefficients in optical coherence tomography,” J. Biomed. Opt. 13(3), 034003 (2008).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.2927464 Google Scholar

33. F. J. van der Meeret al., “Localized measurement of optical attenuation coefficients of atherosclerotic plaque constituents by quantitative optical coherence tomography,” IEEE Trans. Med. Imaging 24(10), 1369–1376 (2005).ITMID40278-0062 http://dx.doi.org/10.1109/TMI.2005.854297 Google Scholar

34. A. I. Kholodnykhet al., “Accurate measurement of total attenuation coefficient of thin tissue with optical coherence tomography,” IEEE J. Sel. Top. Quant. Electron. 9(2), 210–221 (2003).IJSQEN1077-260X http://dx.doi.org/10.1109/JSTQE.2003.814194 Google Scholar

35. P. LeeW. GaoX. Zhang, “Performance of single-scattering model versus multiple-scattering model in the determination of optical properties of biological tissue with optical coherence tomography,” Appl. Opt. 49(18), 3538–3544 (2010).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.49.003538 Google Scholar

36. A. Rogganet al., “Optical properties of circulating human blood in the wavelength range 400–2500 nm,” J. Biomed. Opt. 4(1), 36–46 (1999).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.429919 Google Scholar

37. Y. M. Liewet al., “Motion correction of in vivo three-dimensional optical coherence tomography of human skin using a fiducial marker,” Biomed. Opt. Express 3(8), 1774–1786 (2012).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.3.001774 Google Scholar

38. Y. M. Liewet al., “Reduction of image artifacts in three-dimensional optical coherence tomography of skin in vivo,” J. Biomed. Opt. 16(11), 116018 (2011).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.3652710 Google Scholar

39. J. EnfieldE. JonathanM. Leahy, “In vivo imaging of the microcirculation of the volar forearm using correlation mapping optical coherence tomography (cmOCT),” Biomed. Opt. Express 2(5), 1184–1193 (2011).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.2.001184 Google Scholar

40. C. F. BohrenD. R. Huffman, Absorption and Scattering of Light by Small Particles, Wiley, New York (1983). Google Scholar

41. J. Canny, “A computational approach to edge detection,” IEEE Trans. Pattern Anal. Mach. Intell. PAMI-8(6), 679–698 (1986).ITPIDJ0162-8828 http://dx.doi.org/10.1109/TPAMI.1986.4767851 Google Scholar

42. T. Gambichleret al., “In vivo data of epidermal thickness evaluated by optical coherence tomography: effects of age, gender, skin type, and anatomic site,” J. Dermatol. Sci. 44(3), 145–152 (2006).JDSCEI0923-1811 http://dx.doi.org/10.1016/j.jdermsci.2006.09.008 Google Scholar

43. S. L. Jacques, “Origins of tissue optical properties in the UVA, visible and NIR regions,” in OSA TOPS on Advances in Optical Imaging and Photon Migration, R. R. AlfanoJ. G. Fujimoto, Eds., Vol. 2, pp. 364–371, Optical Society of America, Washington, DC (1996). Google Scholar

44. I. S. SaidiS. L. JacquesF. K. Tittel, “Mie and Rayleigh modeling of visible-light scattering in neonatal skin,” Appl. Opt. 34(31), 7410–7418 (1995).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.34.007410 Google Scholar

45. D. Levitzet al., “Quantitative characterization of developing collagen gels using optical coherence tomography,” J. Biomed. Opt. 15(2), 026019 (2010).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.3377961 Google Scholar

46. P. D. Verhaegenet al., “Adaptation of the dermal collagen structure of human skin and scar tissue in response to stretch: an experimental study,” Wound Repair Regen. 20(5), 658–666 (2012).WREREU1067-1927 http://dx.doi.org/10.1111/wrr.2012.20.issue-5 Google Scholar

47. A. N. Bashkatovet al., “Estimation of wavelength dependence of refractive index of collagen fibers of scleral tissue,” Proc. SPIE 4162, 265–268 (2000).PSISDG0277-786X http://dx.doi.org/10.1117/12.405952 Google Scholar

48. D. P. PopescuM. G. Sowa, “In vitro assessment of optical properties of blood by applying the extended Huygens-Fresnel principle to time-domain optical coherence tomography signal at 1300 nm,” Int. J. Biomed. Imaging 2008, 591618 (2008).IJBIBD1687-4196 http://dx.doi.org/10.1155/2008/591618 Google Scholar

49. R. O. Esenalievet al., “Noninvasive monitoring of glucose concentration with optical coherence tomography,” Opt. Lett. 26(13), 992–994 (2001).OPLEDP0146-9592 http://dx.doi.org/10.1364/OL.26.000992 Google Scholar

50. C. Blatteret al., “In situ structural and microangiographic assessment of human skin lesions with high-speed OCT,” Biomed. Opt. Express 3(10), 2636–2646 (2012).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.3.002636 Google Scholar

51. K. Königet al., “Clinical optical coherence tomography combined with multiphoton tomography of patients with skin diseases,” J. Biophotonics 2(6–7), 389–397 (2009).JBOIBX1864-063X http://dx.doi.org/10.1002/jbio.v2:6/7 Google Scholar

© 2014 Society of Photo-Optical Instrumentation Engineers (SPIE) 0091-3286/2014/$25.00 © 2014 SPIE
Peijun Gong, Peijun Gong, Robert A. McLaughlin, Robert A. McLaughlin, Yih Miin Liew, Yih Miin Liew, Peter R. T. Munro, Peter R. T. Munro, Fiona M. Wood, Fiona M. Wood, David D. Sampson, David D. Sampson, "Assessment of human burn scars with optical coherence tomography by imaging the attenuation coefficient of tissue after vascular masking," Journal of Biomedical Optics 19(2), 021111 (5 November 2013). https://doi.org/10.1117/1.JBO.19.2.021111
JOURNAL ARTICLE
11 PAGES


SHARE
Back to Top