Observation of Atherosclerotic Plaque Phantoms through Saline or Blood Layers by Near-Infrared Hyperspectral Imaging ()
1. Introduction
Atherosclerosis is a disease in which plaque builds up inside the arterial wall. Many acute ischemic events result from the rupture of asymptomatic vulnerable plaques, which can cause thrombotic occlusion and myocardial infarction. Identification of potentially lethal plaques before such disruption will facilitate the development of therapeutic strategies for preventing these events [1] . The most common modalities are currently applied for clinical intravascular diagnostic imaging of atherosclerotic plaque [2] [3] including intravascular ultrasound (IVUS) [4] , optical coherence tomography (OCT) [5] , photoacoustic tomography [6] , angioscopy [7] , and combination of these methods [8] [9] . Spectroscopic techniques are also being applied, based on advances in spectroscopic photoacoustic imaging [10] , fluorescence spectroscopy [11] , Raman spectroscopy [12] , and near-infrared (NIR) spectroscopy [13] [14] . Such methods can provide chemical information about the artery and plaque.
Hyperspectral imaging (HSI) [15] is a spectroscopic imaging method that combines chemical specificity with spatial resolution. This technique can detect the presence of lipid-rich plaques. Previous work demonstrated the applicability of HSI to detection of atherosclerotic plaques in the visible range [16] . Hemoglobin in blood strongly absorbs visible light. This necessitates the use of a saline flush to deliver sufficient light to the vessel wall if HSI is implemented in vivo using an intravascular probe in a manner analogous to IVUS or OCT.
NIR light efficiently penetrates biological tissues, and the NIR region contains the characteristic absorption range of lipid-rich plaques [13] [17] . Thus, NIR spectroscopy in the wavelength range of 1000 - 2500 nm is being evaluated as a potential tool for identifying lipid-rich plaques [10] -[14] [17] [18] . In addition, the absorption of NIR by blood is low, and is dominated by water rather than hemoglobin, potentially obviating the need for a saline flush in an intravascular implementation. Consequently, NIR spectroscopy may enable lipid-rich plaques to be imaged in the presence of blood [18] .
The objective of this study was to evaluate the potential of NIR-HSI for imaging of lipid-rich plaques through saline or blood in an in vivo environment. We demonstrated that NIR-HSI could detect and distinguish lipids in a biological tissue model of atherosclerotic phantoms. We also determined the optimal effective NIR range for imaging of lipids through blood or saline.
2. Materials and Methods
2.1. Exposed Plaque Phantom
An atherosclerotic phantom was prepared using a biological tissue model, bovine fat, and saline or heparinized porcine blood. The biological tissue model (tissue-mimic model) consisted of 0.5 mg/mL hemoglobin (H7379- 10G, Sigma-Aldrich, USA), 0.1 g/mL gelatin (G2500-500G, Sigma-Aldrich, USA), and 2% Intralipid (FB- 01IL20, Terumo, Japan) [19] . Bovine fat was used as a model for lipid-rich plaques in wavelength bands around 1200, 1700, and 2300 nm. Figure 1(a) shows a schematic illustration of an exposed plaque phantom. A 1-mm layer of bovine fat was set at the center of the biological tissue model. Before measurement, a 0 - 10-mm thick layer of saline or a 0 - 2-mm thick layer of blood was overlaid on the sample.
2.2. Vulnerable Plaque Phantom
Vulnerable plaques are usually characterized by a thin fibrous cap overlaying a soft lipid pool. The thickness of the thin fibrous caps is under 65 μm, whereas that of the soft lipid pool is nearly 800 μm [20] . As shown in Figure 1(b), a vulnerable plaque phantom was prepared with a 0.1-mm thick biological tissue model overlaying the exposed plaque phantom. Then, a 0 - 4-mm thick layer of saline or blood was overlaid on the sample.
2.3. NIR-HSI Set Up
Figure 2 shows a schematic diagram of the NIR-HSI system. Supercontinuum (SC) light (prototype, Sumitomo Electric Industries, Japan) was used as a light source. The wavelength range and maximum average power were 1100 - 2400 nm and ~40 mW, respectively. Samples were irradiated with light that was wavelength-selected light using a monochromator (SPG-120IR, Shimadzu, Japan). The diffuse reflected light on the sample was detected by a mercury-cadmium-telluride (MCT) camera that collected images of 320 × 256 pixels (XEVA-2.5, Xenics, Belgium).
2.4. Data Acquisition
Spectral absorption values were obtained at wavelengths of 1150 - 1240, 1700 - 1900, and 2300 - 2400 nm at 10-nm intervals. The groups of hyperspectral data (called hyperspectral cubes) collected at those wavelengths consisted of 10, 21, and 11 images, respectively. Hyperspectral data were calibrated using black-and-white images. The white image cube was obtained from a diffuse reflection standard (WS-1, Ocean Optics, USA). The corresponding image cube was acquired under conditions identical to those for the measurement sample. The
Figure 1. Schematic view of measurement samples. Image (a) is a cross-sec- tion view of an exposed plaque phantom; Image (b) is a cross-section view of a vulnerable plaque phantom.
Figure 2. Schematic diagram illustrating the NIR-HSI system.
black images were obtained with the light source off to prevent any light from penetrating into the camera and were used to estimate the dark current noise. The reflectance factor Rij(λ) at wavelength λ and the spatial coordinates (i, j) were calculated using the following equation:
(1)
where Iij(λ) is the spectral reflectance of the sample captured by the MCT camera, Dij is the intensity of the black images, and Wij(λ) is the spectral intensity of the white images.
Finally, the so-called apparent absorbance at each pixel was calculated using the modified Beer-Lambert law [21] .
(2)
where Aij(λ) is the apparent absorbance of the measurement sample.
2.5. Image Processing
The spectral angle mapper algorithm (SAM) was used for imaging analysis [22] . SAM is a classification method that permits rapid mapping by calculating the spectral similarity between the imaging spectra and a reference spectrum. This method treats both spectra as vectors in n-dimensional space and calculates the spectral angle between them. Small angles means high similarity. This spectral angle is relatively insensitive to changes in pixel illumination because increasing or decreasing illumination does not change the direction of the vector.
In a hyperspectral cube, each pixel (i, j) has each its own vector. The angle between each vectors can be calculated as follows:
(3)
where A is a reference vector, B is the vector of each pixel, and θ is the angle between A and B. In order to obtain the hyperspectral image, the value of cosθ is divided from 0 (black) to 255 (white), expressed as image values. The spectrum of bovine fat, which was measured by Fourier transform NIR spectroscopy (FT-NIR), was adopted as the reference vector [23] .
3. Results
3.1. Optimal Wavelength Band for Lipid Observation
Figure 3 shows the hyperspectral images of the exposed plaque phantom through a 0.5-mm thick layer of saline or blood. Bovine fat was observed through a saline layer at wavelengths of 1150 - 1240 (around 1200), 1700 - 1900 (around 1700), and 2300 - 2400 (around 2300) nm. The images around 1200 and 1700 nm were especially clear. Bovine fat was observed through the blood layer only in the wavelength band around 1200 nm, but the image was blurry.
3.2. Effect of the Thickness of the Saline or Blood Layer
Figure 4 shows hyperspectral images of the exposed plaque phantoms through a 0 - 10-mm thick layer of saline
Figure 3. Images of exposed plaque phantoms taken through 0.5 mm saline (left) or blood (right). (a) and (b) are visible images; (c) and (d) are hyperspectral images at 1150 - 1240 nm; (e) and (f) are hyperspectral images at 1700 - 1900 nm; (g) and (h) are hyperspectral images at 2300 - 2400 nm.
Figure 4. Images of exposed plaque phantoms taken through saline. The three upper panels show visible images taken through saline layers with thicknesses of (a) 0 mm; (b) 2 mm; and (c) 10 mm. The five lower panels show hyperspectral images taken at wavelengths around 1200 nm through saline layers with thicknesses of (d) 0 mm; (e) 2 mm; (f) 4 mm; (g) 8 mm; and (h) 10 mm.
at wavelengths around 1200 nm. Bovine fat was observed under all conditions, but images taken through saline layers over 8 mm thick exhibited weak blurring. Figure 5 shows hyperspectral images of the exposed plaque phantoms through a 0 - 2-mm thick layer of blood at wavelengths around 1200 nm. Lipid could be observed clearly under a layer of blood with a thickness of 1.5 mm or less. Images of bovine fat through blood layers over 1.5 mm thick exhibited strong blurring, and images could not be collected at all through a blood layer more than 2 mm thick. The error rate of an image of bovine fat covered by a 0.5-mm thick layer of blood was 6.2%, whereas the error rate of an image covered by a 1.0-mm thick layer of blood was 22.3%.
Figure 6 shows the hyperspectral images of the vulnerable plaque phantoms through a 0 - 4-mm thick layer of saline with a wavelength around 1200 nm. Bovine fat was observed under all conditions, but images collected through saline layers over 4 mm thick exhibited weak blurring. Figure 7 shows the hyperspectral images of the vulnerable plaque phantoms taken through a 0 - 4-mm thick layer of blood in the wavelength band around 1200 nm. The lipid could be clearly observed when the blood layer was less than 1.0 mm thick, but images could not be collected when the blood layer was 4 mm thick. The error rate of an image of bovine fat covered by a 0.5-mm thick layer of blood was 0.7%, whereas the error rate of an image covered by a 1.0-mm thick layer of blood was 7.7%.
4. Discussion
Figure 8 shows absorption spectra of an atherosclerotic plaque and a normal artery, measured by FT-NIR [23] . The spectra were normalized at a wavelength of 1950 nm. Atherosclerotic plaques have absorption bands around 1200, 1450, 1700, 2000, and 2300 nm. The bands around 1200 and 1700 nm correspond to the second and first overtones of the CH stretching vibration mode, respectively. The band around 2300 nm corresponds to the second overtone of the CH deformation vibration mode. The band around 1450 nm corresponds to the first overtone of the OH stretching vibration mode. The band around 2000 nm corresponds to the combination tone of OH stretching and OH deformation vibration modes [24] .
As shown in Figure 3, in exposed plaque phantoms covered by saline or blood, bovine fat could not be observed at wavelengths around 2300 nm, because water has a strong absorption band at wavelengths over 1400 nm. Through a blood layer, bovine fat could not be observed at wavelengths around 1700 nm, because water in blood absorbs strongly in this range, and red blood cells strongly scatter light [25] . Figure 9 shows the optical penetration depth of blood. The spectrum was obtained by integrating sphere measurements in combination with an inverse Monte Carlo simulation [26] . As shown in Figure 9, the penetration depth was higher in the wave-
Figure 5. Images of exposed plaque phantoms taken through blood at wavelengths around 1200 nm. The upper three panels show visible images taken through blood layers with thicknesses of (a) 0 mm; (b) 0.5 mm; and (c) 2.0 mm. The five lower panels show hyperspectral images taken at wavelengths around 1200 nm through blood layers with thicknesses of (d) 0 mm; (e) 0.5 mm; (f) 1.0 mm; (g) 1.5 mm; and (h) 2.0 mm blood layer.
Figure 6. Images of vulnerable plaque phantoms taken through saline. The upper two panels show visible images taken through saline layers with thicknesses of (a) 0 mm and (b) 0.5 mm. The lower four panels show hyperspectral images taken at wavelengths around 1200 nm through saline layers with thicknesses of (c) 0 mm; (d) 0.5 mm; (e) 1.0 mm; and (f) 4.0 mm.
Figure 7. Images of vulnerable plaque phantoms taken through blood. The upper two panels show visible images taken through blood layers with thicknesses of (a) 0 mm and (b) 0.5 mm. The lower four panels show hyperspectral images taken at wavelengths around 1200 nm through blood layers with thicknesses (c) 0 mm; (d) 0.5 mm; (e) 1.0 mm; and (f) 4.0 mm.
Figure 8. Absorption spectra of arterial tissues. Solid line, atherosclerotic plaque; dashed line, normal artery.
length band around 1200 nm than in the band around 1700 nm. Thus, we considered the wavelength band around 1200 nm to be optimal for the observation of atherosclerotic plaques through saline and blood. However, as shown in Figure 5 and Figure 7, images collected through blood exhibited prominent and blurring at wavelength band around 1200 nm. The blurred images were indistinct due to light scattering through the blood layer,
Figure 9. Optical penetration depth of blood.
making accurate detection of the atherosclerotic plaque difficult. To solve this problem, a scatter correction method is required.
NIR-HSI has the potential to be applied as an angioscopic technique for identifying lipid-rich plaques without clamping and saline injection. The minimum focus of an angioscope is about 1 mm. These results support the idea that NIR-HSI at a wavelength band around 1200 nm could be used to detect lipid-rich plaques through saline or blood layers with an angioscope. However, because our measurement sample was an optical phantom, it will be necessary to perform future studies on actual human atherosclerotic plaques.
5. Conclusion
Atherosclerotic plaque phantoms were observed through saline or blood layers with the NIR-HSI system. In saline and blood environments, wavelengths around 1200 nm were optimal for lipid detection. Hyperspectral images collected through saline layers were clearer than those collected through blood layers in the phantoms of both types. NIR-HSI has the potential to be used in angioscopic diagnostic techniques aimed at identifying lipid- rich plaques without clamping and saline injection.
NOTES
*Corresponding author.