Imaging Raman spectroscopy is a powerful tool to identify chemicals and their distribution. When monochromatic light impinges on molecules, fractions of the scattered light are wavelength shifted related to the molecular vibration states. Thus, Raman spectra are fingerprints that allow a contact- and label-free identification of chemical structures. In contrast to IR absorption spectroscopy that also measures vibration transitions, Raman spectroscopy works in an aqueous environment, which makes this method promising for biological analysis,1–3 especially in the field of medical diagnostics for identification of cancerous tissue.4–8 In surgical cancer treatment, the determination of resection margins is a much-discussed topic.9–11 Too-spacious tissue removal stresses the patient, while too-tight margins reduce the chances of recovery. Usually, a biopsy is taken and examined ex vivo by a pathologist. Based on the results, surgeons later makes use of their prior experience to determine the border between cancerous and healthy tissues. To replace this time-consuming, two-stage approach and to improve its accuracy, there is a demand for spectroscopic methods that allow spatially resolved cancer detection in situ. Meanwhile, there are medical Raman microscopes available that directly allow clinical examinations of skin areas on patients.12 However, commercial Raman microscopic systems for in vivo skin cancer detection are still single-channel systems, i.e., to receive a Raman image from an area of skin, a time-consuming, step-by-step scanning process is necessary.13 Raman intensities are very low (only a fraction in the range of of the scattered light is due to Raman scattering) and in addition often superposed by fluorescence. Thus, even for only a few hundred pixels, the measurement time typically adds up to many minutes or even hours, which is far too long to perform routine examinations on patients in vivo.
To reduce acquisition times, various methods for parallel data collection are described.14 Of particular interest are full-throughput snapshot techniques, also called “multichannel spectroscopy,” “three-dimensional spectroscopy,” or “integral field spectroscopy” (IFS). These techniques work without any serial scanning procedures and do not sacrifice light fundamentally during the recording procedure. IFS was developed in astronomy more than three decades ago15 to save scarce and expensive observation time at observatories. It is based on slicing a two-dimensional (2-D) image into single strip-like segments and stringing them together into one long row in front of a long-slit spectrograph. This can be done using, e.g., a mirror stack image slicer or with a fiber bundle converter: at the sample side, the fibers of the bundle are arranged as a 2-D matrix. In front of the spectrograph’s input, the fiber front surfaces are arranged as a straight line, lying side by side in a V-groove holder. After passing the collimator and camera optics and the dispersive element of the spectrograph, the light signals emerging from each fiber generate a family of individual spectra on a large-area detector. Data reduction software evaluates the raw signal, applies calibrations, and finally provides a data cube containing the entire spectral and spatial information. A review of IFS in astronomy is given in Ref. 16. High-end IFS spectrographs are installed in the multiunit spectroscopic explorer (MUSE) system17 at the very large telescope observatory in Paranal, Chile, since spring 2014. MUSE consists of 24 connected spectrograph modules and is capable of acquiring a total number of 90,000 spatial elements (also called spaxels) within a single exposure. From every single spaxel, the entire spectrum from 465 to 930 nm is recorded at a spectral resolution of 0.22 nm in a total of 4300 spectral bins.
To a certain degree, imaging in vivo Raman spectroscopy and astronomical observations face similar challenges, namely, to efficiently detect faint signals in the presence of bright background light. This notion was the motivation to update a spectrograph based on an MUSE design for use in medical Raman spectroscopy. A fiber bundle-based optical setup was realized to record Raman images of areas of skin, which matches common sizes of lesions suspected to be cancerous. The objective of the project was to validate the concept of a future instrument that would allow dermatologists to promptly recognize in vivo the borders of the cancerous tissue without the need of a time-consuming scanning procedure. The general capability of IFS for generating large-area Raman images was presented previously.18 Similar Raman setups are described in Refs. 1920.–21. Here, we present a discussion and a characterization of a setup using an astronomy spectrograph with regard to Raman imaging of human skin ex vivo and in vivo. Eventually, comparison measurements of healthy and cancerous human skin samples in vivo and ex vivo were performed.
Spectrograph and Optical Setup
Figure 1 shows a scheme of the experimental setup. To record Raman images, an image acquisition head based on two microlens arrays (MLA) was realized. Figure 15(a) shows the top of the image acquisition head with a skin sample and the upper MLA underneath. A detailed description of the Raman acquisition optics was given previously.18 In brief: for excitation, a 784.5-nm diode laser with tunable output power up to 500 mW was used. Its fiber optic output was connected to a square core fiber with core side length. Light emerging from square core fibers shows a top-hat intensity distribution, which is favorable with regard to a preferably homogeneous illumination.22 The excitation light passes a collimation lens, a clean-up filter to remove Raman and fluorescence background generated within the fiber and spontaneous emission of the laser, and is finally guided to the sample by a 45-deg dichroic mirror. An MLA in front of the sample generates excitation spots on the sample with a 0.5-mm pitch, i.e., a square image with 1-cm side length and a sampling of 400 pixels can be recorded. In the opposite direction, the MLA collects the Raman signal from the sample. A pair of relay lenses guides the Raman signal to a further MLA that couples the signal into the fibers of the fiber bundle. A 785-nm notch filter and a long pass (LP) filter remove the Rayleigh signal.
The fiber bundle consists of 400 fibers (114/125/155 VIS/IR, , Heracle, Germany). On the sample side, the fibers are arranged within a square plate containing microholes at center distances of 0.5 mm (Fig. 2). On the spectrograph input, the fibers are arranged side by side within a V-groove holder forming a pseudoslit at a length of 118 mm. There are 421 V-grooves with a pitch of . The fibers of the bundle are arranged in groups of 20. Additional V-grooves between the groups are occupied with fibers that form a fan-out cable that is intended for other applications (calibration, tests). Due to the gaps, the groups of 20 are easily distinguishable in the raw data, thus facilitating a quick inspection by eye before starting the data reduction process.
A detailed description of the spectrograph was given previously.23 The array of optical fibers is attached to a first plane-concave silica lens. The plane side of the lens was initially intended to use index matching gel to minimize coupling losses. However, it was far more advantageous to allow an easy back and forward sliding of the fiber holder. Pushed back, the fibers could be illuminated with a high-power white-light LED array. On the basis of the resulting light spots, the optical components of the head, especially the positions of the MLAs, could be adjusted. During the alignment procedure, the 785-nm long pass filter was flipped. To allow a routine check of the alignment, the index matching gel was waived. Instead, two small pieces of adhesive tape were put on the borders of the fiber holder. The tape pieces served as distance holders and prevented scratching when sliding the fiber holder toward the lens. As a dispersive element, a volume phase holographic grating (VPHG) was used. The shape of the diffraction grating made of gelatin is 118-mm diameter circular, while the entire VPHG component appears with a square shape of 122 mm side. The grating is optimized to cover a wavelength range from 350 to 900 nm. Finally, the diffracted transmitted light is focused on the image plane of a custom-made CCD (charge coupled device) camera. The detector is a large area back-illuminated chip (CCD213, e2v, Chelmsford, United Kingdom) with and pixel size. In the range from 400 to 800 nm, the detector shows a quantum efficiency of , at 900 nm it is . The spectrograph is capable of covering a wavelength range from 350 to 900 nm with a linear dispersion of . However, for sources with a spectral energy distribution larger than one octave, the overlap of second-order signals must be suppressed by use of order separating filters. A description of the camera detector and its readout configuration can be found in Ref. 24. In the standard configuration, is needed to read out the complete CCD chip. The control software saves the raw signal as an FITS file25 and visualizes it using the free available SAOImage DS9 software (Smithsonian Astrophysical Observatory, available in Ref. 26). This viewer is useful to perform quick-look checks, e.g., to verify the presence of a suitable Raman signal by plotting a profile along the wavelength axis. For full data reduction and visualization, the raw data file is processed with the open source software P3D27,28 (available in Ref. 29). The main features of the software are: (I) the exact trace of the spectrum corresponding to each fiber is determined. Due to optical aberrations, the signal traces of the fibers are not parallel to pixels of the detector. The related deviations are determined by continuous signal traces on the detector generated when illuminating all fibers of the bundle with white light. (II) The dispersion across the detector has an arc-shaped characteristic. This effect is commonly known as spectral smile or keystone.30 By use of a Ne lamp or another light source with discrete emission lines, a wavelength calibration is applied, i.e., the wavelength corresponding to each pixel is determined as a polynomial solution. (III) Every signal trace is assigned to the spatial position of the corresponding fiber in the matrix. Finally, a data cube with two spatial axes and one spectral axis is provided as an FITS file. (IV) A viewer module of the software allows to inspect the cube in different ways (display of spatial maps at a chosen wavelength, plot of spectra of selected regions) and to perform simple data analysis tasks such as flux measurements and spectral line fits.
Epoxy Skin Phantoms
To characterize the setup, phantoms based on epoxy resin (epoxy casting resin “waterclear,” R&G Faserverbundwerkstoffe, Germany) were prepared. Since the setup is intended to measure skin samples, the absorption and scattering properties of the phantoms were aligned to human dermis31 by adding proper amounts of and black epoxy paste as scatter and absorber, respectively. After curing, the surface was polished, finally using a sand paper with grain size. The scattering and absorption coefficients of different phantoms were measured (Lambda 900 UV/VIS/NIR Spectrometer with integrating sphere; Perkin Elmer) to and , respectively, at 785-nm wavelength.
Characterization of the Experimental Setup
Image Acquisition Head
Raman imaging with IFS ideally requires a homogeneous laser illumination of the entire image field during the recording process. To verify the conditions of our setup, the illumination of a sample image was simulated by software (OpticStudio, Zemax LLC Delaware). The pseudocolor image of the simulation [Fig. 3(a)] shows that most of the excitation energy is contained in diameter spots at the focus positions of every microlens. Because of the filling-factor of 65%, spaces between the microlenses are also weakly illuminated. At least for nonscattering samples, the remaining 35% do not contribute to excitation, especially since Raman signals arising from the gaps between the lenses are not coupled into the fibers at the other end of the image acquisition head. Due to optical aberrations, the illumination geometrical spot sizes are not exactly constant across the focal plane of the MLA but increase slightly near the edge of the field of view. The simulation reveals that excitation energy within the spots drops from the MLA center to the penultimate rows by . From the penultimate row to the last row, there is an additional drop of 40%.
For comparison, Fig. 3(c) shows a real camera image of the excitation spot received by switching on the laser and putting a piece of scale paper on the top of the image acquisition head. Due to scattering of the paper, the spots appear larger than predicted by the simulation. To a certain extent, the intensity impression is an artifact of the camera optics. The actual intensities were measured using a laser power meter in combination with a shadow mask. The total power of the excitation laser at the exit of the square fiber was 400 mW. Except for the four border rows, intensities of were measured. Obviously no significant loss of laser light occurs within the excitation pathway. However, at the border rows, the intensities drop to , which is in accordance with the simulation. Obviously, the shaft of the MLA holder causes vignetting. In conclusion, an image field of is approximately homogeneously illuminated to within . It should be noted that the edges of the MLAs do not run exactly between mircrolenses. Thus, the central microlenses are not precisely in the center of the holder and hence, the intensity gradients are not perfectly symmetric.
Spatial Resolution for Skin Samples
If a scattering sample such as skin is examined, homogeneous illumination of the surface causes a decrease in spatial resolution because multiple scattering in the sample leads to crosstalk between the pixels: photons originating from illumination of a specific pixel also contribute to the Raman signal at neighboring pixels. To estimate this effect, a Monte Carlo simulation using scattering and absorption coefficients of human skin and 31 has been performed; see Fig. 4. Each fiber of the matrix captures light arising from certain sections of the sample surface. From the magnification factor of the image acquisition head, the aberrations, and the fiber core diameters, it was estimated that these sections are spots with diameters.
Figure 4(b) shows the contributions of adjacent pixels in a 0.5-mm raster. A signal strength of 1.00 is assumed if only one pixel at detection position (pixel at center of the insert) would be illuminated. If also illuminating the neighboring pixels, additional Raman-scattered photons will reach the pixel in the center. Four neighboring pixels in 0.5-mm distance add to 0.23 signal strength. Four neighboring pixels in 0.71-mm distance add to 0.18 signal strength. At 1.1-mm distance, there are 8 pixels contributing in total to 0.25 signal strength. Adding up all neighboring contributions, it turns out that approximately half of the signal detected on a certain position comes from scattered light originating from adjacent pixels.
The resulting effect on the lateral resolution is shown in Fig. 5 for the same human-skin model, where a homogeneously distributed Raman active component is restricted laterally to a half space [ for ]. This example is intended to correspond to an experiment where the boundary of a large-area tumor region is examined. It is assumed that the cancerous region has a specific Raman signal that is not present in the healthy part. This assumption is a simplification. In reality, cancer is not indicated by a clear presence or absence of a definite Raman peak, but rather by slight changes of the spectrum. Nevertheless, the assumption is meaningful to depict the consequences of channel crosstalk. Figure 5 shows two situations: (I) all spots are illuminated simultaneously, as it is the case for our setup (black circles) and (II) only one spot located at the measurement position is illuminated, as it is the case for a pointwise scanning Raman microscope. For our setup, the transition of the Raman signal at the boundary () is remarkably broadened to more than 1 mm, whereas a single spot setup would result in a resolution far better than 0.5 mm.
The Monte Carlo simulation also provides information about the signal distribution from different depths. For measurements on skin tissue, it is important to know the spatial distribution of Raman-scattering processes detected by the instrument. The Raman signal can be described with an integral over the probe volume which contains a weighting function32 The function is a measure of the sensitivity of the setup for Raman scattering at location .
Figure 6 shows the depth dependence of the weight function for Raman signal detected at one spot. The thin curve shows the signal contributions from Raman scattering at different depths as they would be received if only the detection pixel (spot) would be illuminated. In this case, most of the signal would arise in a range from the surface to depth. The bold curve shows the corresponding contributions if all pixels (spots) are illuminated. In comparison to single pixel illumination, the simultaneous illumination of all pixels results in the detection of Raman-scattering processes in larger depths due to crosstalk as discussed for Fig. 4. In the case of real skin, the most intense Raman signals arise from the epidermis and the top of the papillary dermis () but can be enhanced down to using the optical clearing technique.33
Homogeneity of the Raman Images
The epoxy phantoms were used to examine the spatial signal intensity distribution of the setup. The homogeneity of the samples, more precisely of the signals arising from the phantoms, was verified by using a scanning single channel Raman spectrometer (Laser- und Medizin-Technologie GmbH, Germany; excitation wavelength 785 nm).5 At various positions of a phantom, the spectra were measured within a square of 2-cm border length. The variation of all unprocessed spectra was less than , proving the homogeneity of the phantoms (Fig. 7).
Although homogeneous phantoms are used as samples, the received signal distribution was not uniform. Figure 8 shows the average intensity distributions of the raw spectra of a phantom, which appears as a dome-like shape with center intensities that are approximately four times higher than at the corners.
The intensity differences can be associated with (I) the spectrograph’s properties, (II) the characteristics of the image acquisition measurement head, and (III) the scattering properties of the sample.
The spectrograph’s sensitivity (I) depends on the positions of the fibers located in the V-grooves. Due to the design of the spectrograph, light emerging from fibers at the border positions of the V-groove holder is vignetted to a certain amount. The flat-field correction uncovers this characteristic. To perform the flat-field measurement, the fiber bundle matrix (Fig. 2) is removed from the image acquisition head and homogeneously illuminated using an integrating sphere. The fibers in the matrix are sorted with regards to their positions in the V-groove. Fibers at the positions and (20,20) are located at the borders of the V-groove holder, whereas the fibers at (10,20) and (11,1) are in the center. As expected, the spatial sensitivity results in a half-pipe shape showing intensity. Minima at the corners are probably caused by some vignetting due to the limited size of the available integrating sphere (data not shown). A further potential source of inhomogeneities is the distribution of the excitation intensity (II). As shown above, the excitation intensity declines from the center to the border rows. A further reason for the dome-shape is the scattering of the samples (III). As described above, scattering leads to significant signal contributions from neighboring spots. Indeed, fibers in the center of the matrix receive the largest contribution from neighbors. When moving from the center to the borders, the numbers of neighbors and, accordingly, the signal, decrease. Assuming homogeneous illumination and absorption and scattering coefficients of skin, a calculation shows that the relative Raman intensity drops from 1 to 0.6 when moving 10 fibers diagonally from the corner to a center. Finally, a combination of (I), (II), and (III) leads to the dome-like shape shown in Fig. 7. Since the actual absorption and scattering coefficients, and hence the intensity distribution, differ from sample to sample, a universal correction function cannot be applied. Therefore, it was decided to include normalization into the preprocessing of spectra. Fluorescence background was removed using a sixth-order polynomial fitted to each spectrum.34 Subsequently, the standard normal variate normalization was applied.35 It should be mentioned that the laser unit used includes two excitation lasers with 1-nm gap. Thus, the setup is capable of removing background using shifted excitation differential Raman spectroscopy (SERDS).18,36 However, the SERDS option was not used in this work, since Raman spectra reconstructions based on SERDS curves considerably differ from common Raman spectra, making a comparison to previously published results difficult. Further drawbacks of SERDS are the negative impact of photobleaching and the extension of the measurement time. It was therefore decided to apply a “classic” polynomial-based background removal algorithm. As shown in Fig. 9, highly uniform spectra result from the preprocessing. From neighboring pixels, it was estimated that the variance between the spectra is approximately only twice the variance due to noise.
To visualize the variations between the spectra of the Raman image, a principal components analysis (PCA) of the set of 400 spectra of an image was performed. Each spectrum of the image is represented as a linear combination of the orthonormal set of principal components :Figure 10 shows the first three principal components (, 2, 3) for an epoxy phantom. There are variations between the spectra in the wavelength region of 810 to 840 nm, which are probably due to slight variations of filter-transmission spectra. Near the filter edge, the transmission curve shows fringes. The angle of incidence and there with the effective filter response varies slightly with the position of each pixel.
Except for two small peaks at 825 and 900 nm, the principal components do not contain Raman bands of the epoxy sample, therefore it seems reasonable that they only describe the instrument’s variation of spectral sensitivity.
Comparison of Multichannel Setup with Single Channel Spectrograph
The measurements of epoxy phantoms and the PCA described above were also used to evaluate the influence of the Raman-image inhomogeneities in a medical application. A single-channel spectrograph (constructed at the Laser- und Medizin-Technologie GmbH, Germany) has been successfully applied to examine biopsy tissue samples.5 For this setup, a Raman image was realized using a motorized XY stage. The spectral resolution was 0.25 nm, i.e., similar as for our instrument. Partial least squares discriminant analysis (PLS-DA) has been used to discriminate normal from precancerous tissue. PLS-DA results in a discriminant function that assigns each tissue spectrum a scalar value by means of a scalar product that contains a weight function and an average spectrum (Fig. 11), the tissue type is inferred from the sign of37 is reduced.
In principle, each tissue spectrum can occur on any position of the Raman-image. Therefore, the coefficient can attain any value found in the PCA of the Raman-image inhomogeneity. It follows that for each a distribution of values with a standard deviation given asFig. 12 as a function of . The standard deviation follows from the PCA of the epoxy phantom. With Eq. (5), a value of follows. One can conclude from Fig. 12 that the Raman-image inhomogeneity will have no significant effect on the tissue discrimination in this application.
Results and Discussions
Nevus at Forearm In Vivo
The setup was already tested on porcine skin samples ex vivo.18 To verify the capability for human skin samples in vivo, a forearm of a volunteer was placed on the top of the measurement head. Figure 13(a) delineates the test area as a skin imprint of the measurement head. A nevus with diameter can be seen in the upper left. The recording time was 2 min. Due to the fluorescence of the melanin, the position of the nevus can easily be seen on the basis of the signal strength of the raw data (not shown). However, the aim of the experiment was to identify the nevus only with the aid of Raman spectra without any fluorescence background. Thus, all fluorescence background was removed by a polynomial fit followed by a normalization of the spectra. To the resulting background-free Raman spectra PCA was applied. Figure 13(b) shows the Raman spectra inside and beyond the position of the nevus in vivo. Figure 13(c) shows the difference of the principal components 2 and 3 (PC2 and PC3) received from the PCA.
The nevus can be clearly seen as peak. Its diameter of the peak’s footprint is diameter, which corresponds to 3.5 mm; i.e., the Raman image reflects the actual size of the nevus at the camera image. Thus, the achieved resolution of 0.5 mm matches the distance of the fibers in the matrix.
Discrimination of Skin Regions In Vivo
As a test case, Raman-spectral maps of six skin regions (forearm, volar forearm, palm, thumb, leg, and foot) have been measured on four volunteers. The aim of this experiment was to verify whether the system is able to distinguish between different body parts. Following preprocessing as described above, PCA was applied to the total set of spectra.38 The scores of the first two principal components are shown in Fig. 14.
As can be seen from Fig. 14, the skin regions differ in their spectral features. For example, leg (cyan points) and foot (light brown points) are well separated and can be easily discriminated.
To investigate whether one can discriminate the skin regions using more spectral information, an unsupervised clustering algorithm was applied.39 Each spectrum was reduced to the scores of four principal components. Inclusion of more than four components did not improve the result. The distribution of spectra on six clusters show that the spectra of forearm, palm, and thumb and spectra of forearm and leg are not clearly separated by the cluster algorithm (Table 1). Only the spectra of foot skin are clearly distinct, probably due to the very large thickness of the stratum corneum in this case. The main reason for the missing discriminability of the skin of different body sites is the interindividual variance of the skin. The thickness of skin layers as well as microcirculation of blood in the papillary dermis and the melanin concentration varies for different body sites. The concentration of carotenoids, which are most concentrated in the stratum corneum, is an individual parameter, depending on the skin area40 and volunteers’ lifestyle.41 This is apparent when the cluster algorithm is applied to all spectra of a single person in vivo, as shown in Table 2 for an example. Here most of the skin spectra of a body site are assigned to a single cluster.
Partition of skin spectra in vivo from four persons into six clusters obtained by hierarchical clustering. Unsupervised hierarchical clustering using a minimum variance method.39
|Dorsal forearm||Volar forearm||Palm||Thumb||Leg||Foot|
Partition of skin spectra in vivo into six clusters obtained by hierarchical clustering. Spectra from one person.
|Dorsal forearm||Volar forearm||Palm||Thumb||Leg||Foot|
For discrimination of skin of different body sites without reference to individual properties, a more refined discrimination analysis including, for example, supervised learning and feature extraction is necessary.42 Eventually, interarea differences (Fig. 14) exist and should be considered by comparison between cancerous and healthy skin samples. However, this is beyond the scope of this pilot study and would require the investigation of many more samples: clearly the objective of a detailed follow-up study.
Validation of Multiplex Raman for Cancer Diagnostics
Raman spectra of human skin cancerous and healthy tissue were recorded ex vivo and investigated as a possible future diagnostic tool for skin cancer in vivo. Pairs of healthy and cancerous skin biopsy samples were supplied from the Department of Dermatology, Charité. The samples were taken from various patients after surgery and are listed in Table 3. The experiments were approved by the Ethics committee of the Charité-Universitätsmedizin, Berlin (EA1/340/16) and conducted according to the declaration of Helsinki. All volunteers gave their written informed consent. Most of the samples were cylindrical since taken by core biopsy. After removal from the patient, samples 4a and 4b were stored at 6°C and measured one day later. All other samples were at first frozen at and stored. For the transportation from the hospital to the location of the Raman spectrograph, an ice-filled styrofoam box was used. Before starting the Raman measurements, the samples were slowly defrozen in a refrigerator at 6°C.
Biopsy skin samples. BCC, basal cell carcinoma; SCC, squamous cell carcinoma; AK, actinic keratosis. Two consecutive samples originate from one volunteer, where the first sample is normal and the second sample is cancerous.
|Sample #||Origin||Suspicion||Finding||Size||Sex of patient||Age||Frozen|
|1a||Nasal wing||Normal||—||Ø 10 mm||f||75||Yes|
|1b||″||BCC||Not confirmed||Ø 4 mm||″||″||Yes|
|2a||Forehead||Normal||—||Ø 10 mm||m||75||Yes|
|2b||″||SCC||AK grade 2-3||Ø 4 mm||″||″||Yes|
|3a||Cheek||Normal||—||Ø 4 mm||m||74||Yes|
|3b||″||BCC||SCC||Ø 10 mm||″||″||Yes|
|4a||Head||Normal||—||Ø 4 mm||m||76||No|
|4b||″||SCC||SCC||Ø 4 mm||″||″||No|
|5a||Upper arm/shoulder||Normal||—||Ø 4 mm||m||79||Yes|
|5b||″||BCC||BCC||Ø 4 mm||″||″||Yes|
|6a||Shoulder||Normal||—||Ø 10 mm||m||73||Yes|
|6b||External ear||BCC||BCC||Ø 4 mm||″||″||Yes|
|7a||Lower leg||Normal||—||Ø 10 mm||m||77||Yes|
|7b||″||SCC||Basosquamous carcinoma||Ø 6 mm||″||″||Yes|
|8a||Temple||Normal||—||Ø 2 mm||w||91||Yes|
|8b||″||BCC||BCC||Ø 2 mm||″||″||Yes|
Figure 15(a) shows the biopsy sample 6b (Table 3) placed at the top of the image acquisition head. Two parallel needles were used to hold the sample by clamping or skewering. Before and during the measurement, the sample was covered by a lid containing a wet paper tissue to avoid drying. For some of the samples, it was necessary to increase the excitation intensity. The maximum output power of the used laser source is 500 mW. Taking into account the MLA fill factor and losses in the excitation pathway limits the excitation intensity to when illuminating all 400 lenses, which is safe for in vivo application. For the samples replacing the square core fiber by a square core fiber, the excitation power is concentrated to 100 lenses, i.e., the excitation power increases fourfold.
Since the number of available biopsy samples was too low to perform any reliable refined discrimination analysis, it was decided to assess the results on the basis of the average spectra. Biopsy samples of normal skin from four patients and samples of basal-cell carcinoma (BCC) from three patients were measured, summarized, and normalized. The available samples were either entirely affected with cancer or entirely healthy. Thus, within a sample, the spectral patterns were quite similar at all positions. Figure 16 shows the resulting ranges of Raman-spectral values for both types of skin. The comparison of the Raman spectra reveals the following differences between normal skin and BCC:
• A narrower line shape for BCC at the amide I – region (1640 to )
• The intensity ratio of region (1220 to ) (amide III – region) to region (1290 to ) decreases for BCC.
• Decreasing bands in the region (830 to ) for BCC.
These differences are in accordance with published results received from a single channel Raman setup.43 However, with our setup, far more pixels could be measured in a shorter time. In Ref. 43, a 10-min measurement time was applied to record the Raman spectrum of a single diameter spot. In this work, the Raman spectra of 100 spots were recorded within 2 min. Even though the experimental conditions differ to a certain extent (up-to-dateness of the hardware, excitation wavelength, and excitation power), the achieved measurement speed increase from to is striking and documents the capability of IFS even in the field of optical cancer diagnosis.
Conclusion and Outlook
By transferring IFS from astronomy to imaging Raman spectroscopy, a setup was realized that is capable to measure the Raman spectra of 400 pixels of a sample simultaneously without any scanning procedure. Within this work, the applicability of this setup for examination of human skin patches was investigated. Monte Carlo simulation yielded that a spatial resolution of can be achieved when illuminating the whole sample area with the excitation light. This result was experimentally confirmed by recording the Raman image of a nevus in vivo. Finally, biopsy samples of cancerous (BCC, SCC, and AK) and healthy parts of skin were examined. The changes in Raman spectra match literature values received with a classic single channel spectrometer. Although the few small biopsy samples available for this experiment did not allow us to directly localize distinct malign and benign tissue on a single sample, as we would require to actually prove the capability of this method to determine resection margins as the ultimate goal for surgery, we were able to demonstrate the plausibility by means of in vivo measurements of a nevus. Clearly, the next step must be a detailed clinical study to validate this finding on the basis of statistically meaningful samples and larger pieces of tissue. However, the results of our pilot study already show great promise for the development a minimal-invasive optical medical device for detecting cancerous tissue in vivo in the future. Currently, a follow-up project is being undertaken to examine the implementation of IFS in a medical endoscope.
The authors acknowledge support from the German Federal Ministry of Education and Research (BMBF) VIP program, “Multiplex-Raman-Spektroskopie” (MRS), Grant No. 03V0843.
Elmar Schmälzlin is a senior scientist at Leibniz Institute for Astrophysics. He received his diploma in physical chemistry and his PhD in the field of laser spectroscopy from the Ludwig-Maximilians-Universität, Munich, in 1996 and 2000, respectively. He is the author of more than 30 journal papers. His current research interests include fiber optics, imaging Raman spectroscopy, and optical sensing of oxygen.
Ingo Gersonde is a senior scientist at the innoFSPEC Potsdam, University of Potsdam, Germany. He received his PhD in physics from Freie Universität Berlin, Germany. His current work is concerned with photon density wave spectroscopy and fiber spectroscopy. His work also includes chemometrics and classification methods for analysis of spectra and sensor data.
Martin M. Roth is a professor at the University of Potsdam, Institute of Physics and Astronomy. He is a speaker for the innoFSPEC Potsdam Innovation Center and head of the innoFSPEC branch at the Leibniz Institute for Astrophysics Potsdam (AIP). He received his diploma degree in physics and his PhD in astrophysics from the Ludwig-Maximilians-Universität Munich in 1986 and 1993, respectively. He is the author of more than 80 journal papers and has written two book chapters. His current research interests include resolved stellar populations in nearby galaxies, astrophotonics, and biomedical imaging.