Translator Disclaimer
14 January 2016 Spatially resolved, diffuse reflectance imaging for subsurface pattern visualization toward development of a lensless imaging platform: phantom experiments
Author Affiliations +
A portable, spatially resolved, diffuse reflectance lensless imaging technique based on the charge-coupled device or complementary metal-oxide semiconductor sensor directly coupled to the fiber optic bundle is proposed for visualization of subsurface structures such as superficial microvasculature in the epithelium. We discuss an experimental method for emulating a lensless imaging setup via raster scanning a single fiber-optic cable over a microfluidic phantom containing periodic hemoglobin absorption contrast. To evaluate the ability of the technique to recover information about the subsurface linear structures, scattering layers formed of the Sylgard® 184 Silicone Elastomer and titanium dioxide were placed atop the microfluidic phantom. Thickness of the layers ranged from 0.2 to 0.7 mm, and the values of the reduced scattering coefficient (μs) were between 0.85 and 4.25  mm1. The results demonstrate that fiber-optic, lensless platform can be used for two-dimensional imaging of absorbing inclusions in diffuse reflectance mode. In these experiments, it was shown that diffuse reflectance imaging can provide sufficient spatial sampling of the phantom for differentiation of 30  μm structural features of the embedded absorbing pattern inside the scattering media.



Structural and functional analysis of superficial vasculature provides important diagnostic information about microcirculation pathogenesis.1 Dense plexus vessels participate in oxygen delivery to the dermal layer via a network of microvessels, which rise close to the surface of the epithelium and form intrapapillary capillary loops (IPCLs).2,3 In mucosal layers, the capillary network is located 250 to 350  μm below the surface.4,5 Imaging of such loops requires high spatial sampling at the depth range of 0.2 to 0.6 mm. Individual IPCLs are periodically distributed in the arches of the dermal papilla space. The distance between oxygenated and deoxygenated vessels of a capillary is minimal, and sometimes there is no distance at all. The mean range of distances between two successive papilla is 120  μm (diffusion length of oxygen from a vessel into the interstitial tissue is about 100  μm).2,6 Structural alterations of capillaries’ periodicity indicate some degree of a medical condition, such as rheumatic diseases or cancer.6 Architectural disorganization, giant capillaries, hemorrhages, angiogenesis, and avascular areas characterize >95% of patients with overt scleroderma (autoimmune rheumatic diseases).7 Dilation, elongation, and tortuousness are the features of developing carcinomas.6,8 As invasion of cancer progresses, organization of the IPCLs is destroyed and the network is characterized by large vessels with no hairpin structure at the terminal branches.

There are several approaches to imaging of the superficial vascular network. Diffuse reflectance tomography (DOT) is used to evaluate tissue within a volume sampled by the light traveling from a source to a detector.5,911 Source-to-detector distance is one of the aspects determining the resolution of the resultant image in DOT. The technique recovers local optical properties and quantifies changes in relative concentrations of the chromophores of interest. The most traditional method of microvascular imaging is capillary magnification through a microscope or by means of cameras. This technique is termed capillaroscopy. It is mainly employed to analyze microcirculation of the vasculature from the nailfold of a finger or a toe because the scattering layer of the epidermis is very thin at those locations.7,12 When the IPCLs are located beneath thick and highly scattering epithelium layers (most of mucosa and skin), spatial architecture of capillaries is generally invisible for reflectance magnification imaging. Furthermore, сapillaroscopy is highly dependent on conditions of data collection and quality of the imaging optics, restricting imaging to the depth of field of a thin area within a narrow field of view (FOV).

Other imaging techniques applicable for reflectance imaging of microcapillaries, such as optical coherence tomography (OCT), broadband and narrowband endoscopic imaging/spectroscopy, and intravital microscopy also suffer from various limitations.13 OCT can visualize the IPCL, but the equipment is complex, requires mechanical scanning, and is relatively expensive.14 Intravital microscopy has a limited scope for clinical applications toward human subjects.15 Overall, most of these imaging methods rely on conventional lens-based imaging approaches, which are mainly characterized by their FOVs, magnification factors, spatial resolution, and focusing. The spatial resolution and FOV are inversely proportional for the same sensor pixel density, and achieving a higher spatial resolution (through optical magnification) results in diminished FOV. Also, due to the arbitrary three-dimensional (3-D) geometries of vascular networks and living tissue motion artifacts, focusing becomes challenging under mobile conditions in-vivo. In case of microvascular pathogenesis, acquisition of large FOVs with high spatial resolutions as well as increased depth of field to study the progression of abnormal conditions is critical. Therefore, to visualize subsurface patterns effectively, a technique has to possess sufficient spatial sampling while preserving instantaneous, large FOV necessary for wide-area analysis in real time of deeply lying microvascular structures.16 Portability, reproducibility, and the ability to provide optical sectioning of depths under highly scattering conditions can be fulfilled by lensless, fiber-optic imaging platforms capable of performing high sampling with multiple source–detector geometries.17,18 Lensless imaging encompasses a broad group of applications based on a complementary metal-oxide semiconductor (CMOS)/charge-coupled device chip used for imaging of a sample with no lenses attached in front of the sensor plane to achieve a large FOV with high spatial sampling.19 This includes x-ray imaging, lensless compressive imaging, diffractive (phase) imaging, high-resolution lensless endoscopy, and super-resolution imaging techniques.2024

The design of the fiber coupled, lens-free imaging platform proposed by our group is based on a fiber-optic plate (FOP) directly mounted onto a CMOS sensor. The FOP is placed onto, or close to, the object’s surface, and reflected light is transferred to the imaging array without the use of conventional lenses. Such imaging architecture can provide a sufficiently high sampling of diffusely reflected signal and incorporate multiple source–detector configurations. Commercially available FOPs have unit fiber-core sizes of 3 to 5  μm and can be tapered to achieve magnification/demagnification effects.25 At the current stage of its development, the fiber-optic, lensless imaging platform cannot demonstrate a clear advantage over lens-based systems in terms of spatial resolution and quality of the image. However, it offers a number of advantages for specific applications where diffuse reflectance imaging is valuable for its ability to recover information about attenuation of light by the tissue chromophores and where magnification is not potent.

The design described in the previous paragraph was tested and verified via performing spatially resolved imaging through raster scanning fiber-optic detectors across the absorption perturbations, under different scattering conditions and depths (thicknesses). We aimed to establish the applicability of the lensless, FOP coupled imaging sensor for diffuse reflectance imaging. To emulate a 2-D fiber-optic array arranged in a grid, a single fiber (detector) was scanned across an FOV for continuous data acquisition. Broadband illumination allowed for selection of the optimal spectral bands of wavelengths at each scan point. Through collection of sufficient statistics on parameters such as scanning step size, thickness of the scattering layers, scattering optical properties, and choice of the detector size, the method is evaluated for further development.




Spatially Resolved, Diffuse Reflectance Imaging System

The spatially resolved, diffuse reflectance imaging setup is described in Fig. 1. Widefield, oblique collimated illumination of the entire sample was achieved with a 0.40-mm core optical fiber (Silica fiber, NA 0.22, Thorlabs Inc.) coupled to a broadband white light source, HPLS-30-02 (Thorlabs Inc.), and a fiber optic collimator (OZ Optics, Canada). For image acquisition, regions of interest of 1  mm×0.060  mm and 0.50  mm×0.050  mm were chosen to minimize the time of mechanical scanning. Two individual optical fibers were selected to collect the reflected signal: a 10-μm multimode silica detector fiber (NA 0.1, Thorlabs Inc.), and a 50  μm (NA 0.22, Thorlabs Inc.). One fiber at a time was connected to a spectrometer (Ocean Optics USB4000), which had a spectral range of 200 to 1100 nm with a resolution of 1.5 to 2.3 nm full width at half maximum. Integration time range was between 5 and 200 ms. The spectrometer data were processed using a custom script in MATLAB (MathWorks) written in our laboratory.

Fig. 1

(a) The experimental setup illustrates the raster scan setup for collection of the reflectance signal from the phantom with microchannels (area C) filled with hemoglobin through reservoirs A and B. Optical fiber was raster scanned in a pattern depicted above. The y step size was kept constant at 5  μm, and the x step size was varied while keeping a constant field of view (FOV). The illuminating fiber provides oblique, collimated illumination (angle was kept close to 45 deg with the y axis pointing downward) of the microfluidic phantom covered by scattering layers (D) of different thicknesses and optical properties. The detector fiber then transfers the reflected photons to the spectrometer. (b) Image of microfluidic platform taken with 10× objective at ≈ unit magnification. The FOV for the scanning experiments with 10-μm core fiber is outlined by the rectangle (1×0.060  mm2). The area of microfluidic is 5.75×5  mm2.


For comparison with the traditional lens-based setup, several images were acquired with an RGB CMOS chip (Motic, China) using a 10× objective (0.25 NA, Motic, China). The sensor was placed close to the back focal plane of the objective, permitting a magnification ratio of 1. The CMOS matrix had physical dimensions of 6.8  mm×5.3  mm and a maximum resolution of 3 megapixels (MP; 2048×1536), which corresponds to a pixel pitch 3.2  μm, neglecting packaging area losses at the chip boundary. Spectral resolution of the CMOS was estimated as the spectral width of the red–green–blue channels (100  nm). The Software Motic Images Plus 2.0 for PC software package was used for image acquisition, and the images had a digital resolution of 8 bits. The visualization of images and the image processing were performed by an algorithm created in MATLAB (MathWorks).

The microfluidic platform is composed of narrow channels covered with various scattering layers (prepared separately) placed on the top (D in Fig. 1). The microchannel plate was manufactured from a clear polydimethylsiloxane (PDMS, silicone) material; n=1.56 in the visible spectrum without scattering particles. The parallel grooves were designed to approximate a nailfold microvascular spatial pattern with mean physical dimensions: 34-μm wide groove, separated by 82-μm spacing, and a groove height of 30  μm.2,26 The solution of the 90% to 95% oxygenated, human hemoglobin (Agat-Med, Russia, 160  g/l) was injected in reservoir A in Fig. 1 (2  mm×2  mm) by means of 18-gauge needle syringes. Hemoglobin was then redistributed from reservoir A to reservoir B (2  mm×2  mm) to obtain a close to uniform distribution of hemoglobin in the microfluidic channels (area C in Fig. 1).


Fabrication of Top Scattering Layers

To approximate the optical properties of biotissue, the microfluidic channels were covered by layers of variable thickness and scattering properties, produced with Sylgard® 184 Silicone Elastomer DOW/Corning, a two-component silicone product with a refractive index of 1.41. This silicone allows construction of complex geometries and easy incorporation of scatterers. The scatterer concentrations were based on weight percentage of TiO2 (Sigma Aldrich) reported in the literature. About 0.005%, 0.010%, 0.015%, 0.020%, and 0.025% concentrations were selected to match the scattering properties (μs, mm1) of the epithelium.2729 The geometric thickness of the layers was measured with a high-precision gauge. The results were found to be consistent with similar silicone phantom thickness measurements reported in Ref. 30. The averaged thickness of the silicone layers with standard deviation and their scattering properties are given in Table 1.

Table 1

Physical values of the μs′ (mm−1) at λ=545  nm for a range of TiO2 concentrations (columns), and thicknesses (rows) of the upper diffusive layers. Thus, B2 will be the layer with μs′=1.564  mm−1 and thickness of 0.32 mm.

Thickness (mm)Scattering coefficient [μs′ (mm−1)]

Individual scattering blocks were made by mixing the desired concentration of scatterer with the curing agent component of the silicone elastomer. To obtain a homogeneous mixture, it was placed in an ultrasonic bath (Branson, 1510R-DTH) for 10 min at 30 kHz to break residual clusters. The curing agent was then mixed with the silicone (19 weight ratio) under careful stirring. Remaining air bubbles were removed by using a vacuum pump (Barnant) that kept the mixture under low pressure for 30 min. A small portion (0.5 ml) of the final mixture was placed between two standard microscopic glass slides, separated by two spacers within the range of 0.200 to 0.800 mm. Finally, a thin single-layer scattering phantom block was produced by curing at 30°C for 6 h or at room temperature for 24 h. Then this block was gently placed over the microfluidic channel area before starting measurements. Electrostatic forces kept the phantom scattering layer and the microfluidic together. To verify the scattering properties of the phantom blocks in the visible range (400 to 700 nm), the integrating sphere technique (UV-VIS-NIR Spectrophotometer Shimadzu UV-3600, Japan) was used. An inverse adding doubling software program was used to compute the optical properties of the samples.31


Image Processing

Spectral data acquired at each point resulted in a dataset of the dimensions {x,y,wavelengths}={500(1000)  μm/xstep size,50(60)  μm/ystep size,3648wavelengths}. The raster scan was performed for x and y step sizes of 2, 5, 10, 20, and 50  μm, for two FOVs 0.50×0.050  mm2 and 1×0.60  mm2, as shown in Fig. 1. For choosing a spectral band with the highest image contrast, spatial variance analysis was performed on a dataset obtained without any scattering layers. The spectra were divided into 10-nm bands across 450 to 800 nm. The mean of each band was used to calculate the normalized variance in the spatial dimensions. Before calculating variance, the data were mean centered and normalized to the mean.

Due to the oblique setting of light, background intensity varied across the FOV. In order to account for variations, we assumed that the fluctuations in the background illumination were of low spatial frequency compared to the fluctuations exhibited by the embedded hemoglobin contrast obtained from the microfluidic platform. As shown in Fig. 2(b) (smoothed version denoted by dashed line), a lowpass filter was employed with an adaptive moving average kernel having a size (N) [Eq. (2)] of half the scan line length in the x direction of the image. The filtered image (Iavg(x,y)) was then subtracted from the original image I(x,y) [Eq. (1)]. This resulted in an image with background correction [Eq. (3)], showing pure contrast of the microfluidic channels filled by hemoglobin. Because the scan was performed line by line and an image was formed by a collection of such lines, the moving average filter was estimated to correct for the background for each line.

Eq. (1)


Eq. (2)

Iavg(i,y)=1N(IiN2,y+IiN2+1,y+IiN2+N,y)where  i=1tox,foralinearscan,y=1

Eq. (3)


Fig. 2

Image processing steps: (a) a raw reflectance image formed through the raster scan is shown. (b) A line plot of the intensity variation across the white line at 25  μm in (a). The dashed line in plot (b) represents a moving average filter with a kernel size of half the scan line width emulating the background variation in light intensity due to the illumination variance. (d) A background corrected version of the scan line and (c) corresponding background corrected image.




The alternating dark and light bands observed on all images represent an absorption contrast within narrow FOVs after the light has traveled through the hemoglobin grooves of the phantom under a scattering layer [Fig. 1(b)]. Properties of all diffusive layers for a selected wavelength, λ=545  nm, are listed in Table 1. For labeling purposes, scattering properties of the layers were coded with letters A to E and thickness was coded from 1 to 5. All together, five thickness and five scattering coefficients were tested experimentally, resulting in 25 datasets for image analysis. Images were presented with a modified aspect ratio in order to better visualize the patterns.

Dependence of image resolution on the scanning step size, thickness of the upper layers, scattering property, and detector size can be observed in Figs. 35. Quality of an image was significantly affected by the scanning step for the 50-μm detector (Fig. 3). Images formed from 2-μm step size scans demonstrated superior contrast despite the interrogated volume sampled by the fairly large detector diameter. As anticipated, increasing the reduced scattering coefficient and scanning step size decreased the quality of the image due to blurring from “A1-2” (least scattering and scanning step size) to “E1-50” (largest scattering and scanning step size).

Fig. 3

Grayscale images of the microfluidic platform for a fixed thickness (d=0.24  mm), different scattering properties of the upper layer at λ=545  nm. The scanning step (from left to right) tested for 50  μm detector core (a), FOV=0.05×0.5  mm2 is 2, 5, 20, 50  μm and (b) 10× magnification images of the same scattering layers acquired with microscope camera. The horizontal bars on the bottom right corner of the images represent 100  μm.


Fig. 4

(a) Images for five thicknesses (x-axis) and five scattering properties (y-axis) scanned with 10  μm core fiber (NA=0.1) are presented. The red line demarcates scattering properties, the yellow line divides one image into two halves: lower part (NL: no layer) is an image scanned without any scattering layer from the same location as the image collected via a top layer (image above the yellow line). (b) The same experiments were repeated for 50-μm core fiber detector. Scanning step is 2  μm. (c) The layout of the microfluidic phantom and the corresponding scanned FOV. The white channels in the microfluidic are the dark channels in the scanned images (a and b) as these channels were filled with hemoglobin.


To compare deterioration of the optical signal due to the thickness (μm, from left to right) and increased scattering property of the layers [μs (mm1) from bottom to top], the matrix of the images with all 25 experimental conditions is displayed in Fig. 4 for (a) 10-μm and (b) 50-μm detectors. Analyzing the contrast for the scanning step of 2  μm, we observe that decreasing the step size (i.e., a detector size) can dramatically increase the imaging quality (at least for linear objects) due to high spatial sampling (Fig. 3). Also, to ensure that the observed pattern of dark lines is due to the absorption by the hemoglobin in the channels and not an image-processing artifact, two images of no-layer scan and scan through the scattering layer were performed on exactly the same location. The images obtained after scanning with a scattering layer were aligned with its corresponding no-layer scan (Fig. 4).

Fig. 5

(a) Scanned FOV (1×0.06  mm2) diffuse reflectance images for all thicknesses (x-axis) and (b) all scattering (y-axis) acquired with 10-μm core fiber, 2-μm step size are compared with 10× magnification images. Preprocessing algorithm is same for (a) and (b). The small bars inset in the bottom left of the figure represent 100  μm.


Diffuse reflectance images collected with both detectors were compared with 10× magnification images of a conventional microscope (Figs. 3, 5, and 6). For every matching condition, magnification images were put through the same image processing algorithm as for the scanned figures. Significant loss of visibility of the absorption pattern occurred for lens-based images for a higher scattering property of the layers as well as for the deeper-located sample.

Analysis of the amplitude of the optical signal from the scanning experiments and lens-based images are compared in Fig. 6. The 2-D image of the microfluidic platform represents a pattern of periodically absorbing bands embedded in the transparent background of the PDMS. Although it is expected that the image of parallel channels should yield the shape of the “rect” function, the manufacturing of the microfluidic and distribution of absorbers inside the grooves indicates that the shape of the amplitude is contaminated and deviates from the square amplitude of the “rect” function, even for the surface image of the 10× objective image [Fig. 6(e)]. In case of diffusive media, an ideal shape of the optical signal is traded for recovery of the attenuation amplitude from the scattering subsurface of the sample, where an advantage of the absorption-based method is evident (Figs. 5 and 6). Although the images obtained through scanning show distortions, it is estimated that an optomechanical system composed of a multifiber assembly with an instant large FOV image collection would get rid of distortions generated due to the uncertainty of the scanning stage.

Fig. 6

Fiber scanned images of large FOVs obtained with two detectors: (a) 10-μm fiber with scanning step=10  μm×50  μm and (b) 50-μm fiber with scanning step=10  μm×20  μm are compared with lens-based images: (c) 10× magnification objective of the microscope, and (d) a 50-mm planoconvex lens. Amplitude of the signal for both (a and c) conditions are plotted in (e). Dashed lines represent a line extracted from (a), and the solid line indicates a line extracted from (c). No scattering layers were present.


To evaluate the reproducibility of images and stability of the experimental setup, we performed 10 repeated trials of raster scans on layer A3 (Table 1). The resultant mean and standard deviation are shown in Fig. 7. The interexperiment variability appears to be quite low and is presumably caused mainly by uncertainty of placement of the fiber by the xy stage (0.5  μm) and temporal fluctuations of light intensity. Thus, major variations in the diffuse reflectance signal can be mostly attributed to the absorption contrast by the hemoglobin channels of the microfluidic platform.

Fig. 7

The plot shows standard deviation for a sample size of 10 trials for different scan step sizes. The FOV and the scattering layer were kept constant for all these trials and the mean standard deviation (MSTD) shows low signal variability. The insets represent the scanned image corresponding to the structure of the periodic absorbing channels.



Discussion and Conclusions

Imaging of microvasculature in vivo is one of the core research directions in the diffuse reflectance imaging field. The main limiting factor of diffuse reflectance imaging is the scattering in biological tissues, which decreases image contrast and restricts the sampling depth. Diffuse optical tomography is a process of mapping local optical properties of tissues with multiple source–detector configurations and reconstruction of the attenuation coefficients by means of mathematical algorithms. The resultant spatial sampling is quite poor, and the need to improve resolution of microvascular images persists.11,18 We propose a method to overcome some of the limitations of diffuse reflectance imaging by employing the latest advances in sensor and fiber-optic technologies. In this study, a fiber-optic scanning system was employed to collect spatially resolved diffuse reflectance images from a phantom for testing technical parameters for a design of an imaging architecture. Overall, the results have shown that sampling microscale subsurface patterns with fiber optics can provide a low-cost solution for clinical imaging through scattering media.

As a proof of principle of the device, it was demonstrated that images formed via mechanical scanning of the pattern maintain intensity variations of the absorbing structure through scattering layers of up to μs=4  mm1 and 0.72 mm thick (Figs. 35). Selection of the diameters of the detector-fibers was based on the physical dimensions of the microfluidic phantom: an individual channel of 30-μm width has sufficient sampling points to be detected with a 10-μm detector, whereas the entire period between two neighboring grooves is probed by a 50-μm detector. According to the Nyquist theorem, dimensions of the detectors and step sizes are within the limits of necessary sampling and should not show any image aliasing. Even so, the mechanical scanning system was set for overlapping sampling of the region of interest for better visualization of absorption contrast. While keeping the area of interrogation constant, we collected light in virtual subpixel arrangements to compensate potentially limited sampling of the detectors and to confirm that the experiments were not affected by spatial frequency artifacts. This is especially evident from a full FOV figure of the entire microfluidic [Figs. 6(a) and 6(b)], where one can see all the features of the phantom as compared with the lens-based images.

Traditional diffuse reflectance imaging systems suffer from low resolution (mm)/large sampling volumes due to the source–detector configurations of the probes.11 Resolution of our fiber-optic–based system is the function of an individual core dimension of the FOP and the interrogation volume (λ, μs); i.e., the smaller the core of a fiber, the closer the source–detector distance and the higher the sampling. Although further minimization of fiber core (detector) size results in diminished light collection efficiency, leading to longer integration times, it can be alleviated by exploiting high sensitivity sensors (e.g., electron multiplying charge coupled devices) and by binning the pixels (similar to mosaic pattern binning of RGB images). As an example, modern endoscopes rely on fiber-optic image conduits with small cores (2 to 5  μm) for deep cavity imaging. We built our prototype to operate under similar light intensity levels.

In order to verify the performance of our lensless, fiber-optic approach against standard lens-based imaging, we acquired reflectance images of the microfluidic phantom with scattering layers using a 10× (0.25NA) microscopic objective and a 50-mm planoconvex (0.26NA) lens (Figs. 3, 5, and 6). Figure 5(b) depicts limits of contrast visibility that demonstrates the diminishing capability of a lens-based system to image contrast under higher scattering and larger depths, whereas the lens-free imaging technique performs better at recovering absorption contrasts under similar conditions. Although focusing errors of the lens-based setup can be overcome by physical scan of the focal spot and stitching together the obtained images to gather a large depth of field (and FOV), it still requires mechanical motion, and for a portable device, it limits practicality and raises the cost. To emphasize, the proposed method of fiber-optic imaging is not in direct competition with the lens-based technologies in terms of spatial sampling; rather, it can offer a focus-independent interrogation depth for a specific clinical demand. In particular, it can be used to assess optical density of the epidermal layer (to estimate its thickness) and calculate density of IPCLs and their pattern periodicities, especially in case of horizontal capillary loops, which are close to surface (e.g., nailfold, lips, etc.). Furthermore, it is advantageous under mobile conditions as it is not affected by loss of focusing on uneven surfaces, collects spectral information, and has a sampling depth that can be adjusted by implementing narrow-band imaging. By contrast, any lens-based approach requires a working distance between the image and object planes governed by the optics used. Typically, this is an unwanted trait when designing portable point-of-care devices, which need to have a small form factor (e.g., a device the size of a pulse oximeter). Finally, it has been found that lensless imaging is able to detect absorption inhomogeneities on a depth of more than one transport mean free path (1/μs), which outperforms lens-based systems in our previous estimates.32

In summary, imaging of subsurface linear patterns through a diffusive environment based on the absorption contrast was demonstrated using a fiber-optic scanning system. Spatially resolved diffuse reflectance images collected from the microfluidic phantom demonstrate that the fiber-optic imaging method can provide sufficient sampling of microstructures to recover information about subsurface lines as small as 30  μm. Such an imaging architecture was tested as a proof of principle experiment toward the development of a lensless fiber-optic imaging platform. The compact design of the FOP-based camera can be further improved with implementation of smaller diameter core fibers and smaller pixels of CMOS sensors.


Alexandre Douplik gratefully acknowledges support of the Erlangen Graduate School in Advanced Optical Technologies (SAOT) by the German National Science Foundation (DFG) in the framework of the excellence initiative. This work was also supported in part by the Ryerson University Health Research Grant and the NSERC Personal Discovery Grant, Canada. The research was also partially supported by the Ontario Graduate Scholarship, Ryerson University, Toronto, Canada.



J. H. Takano et al., “Detecting early oral cancer: narrowband imaging system observation of the oral mucosa microvasculature,” Int. J. Oral Maxillofac. Surg., 39 208 –213 (2010). Google Scholar


I. M. Braverman and A. Yen, “Ultrastructure of the capillary loops in the dermal papillae of psoriasis,” J. Invest. Dermatol., 68 (1), 53 –60 (1977). Google Scholar


A. Yen and I. M. Braverman, “Ultrastructure of the human dermal microcirculation: the horizontal plexus of the papillary dermis,” J. Invest. Dermatol., 66 131 –142 (1976). Google Scholar


I. M. Braverman and A. Yen, “Ultrastructure of the human dermal microcirculation: II. The capillary loops of the dermal papillae,” J. Invest. Dermatol., 68 (1), 44 –52 (1977). JIDEAE 0022-202X Google Scholar


R. A. Bennett, R. N. Pittman and S. M. Sullivan, “Capillary spatial pattern and muscle fiber geometry in three hamster striated muscles,” Am. J. Physiol., 260 (2), H579 –H585 (1991). Google Scholar


Y. Kumagai, M. Toi and H. Inoue, “Dynamism of tumour vasculature in the early phase of cancer progression: outcomes from oesophageal cancer research,” Lancet Oncol., 3 (10), 604 –610 (2002). LOANBN 1470-2045 Google Scholar


M. Cutolo et al., “Nailfold videocapillaroscopy assessment of microvascular damage in systemic sclerosis,” J. Rheumatol., 27 155 –160 (2000). Google Scholar


H.-Y. Li et al., “Current clinical applications of magnifying endoscopy with narrow band imaging in the stomach,” Diagn. Ther. Endosc., 2012 1 –9 (2012). Google Scholar


A. J. Welch and M. J. C. Van Gemert, Optical-Thermal Response of Laser-Irradiated Tissue, Springer, London (2011). Google Scholar


A. Kienle et al., “Spatially resolved absolute diffuse reflectance measurements for noninvasive determination of the optical scattering and absorption coefficients of biological tissue,” Appl. Opt., 35 (13), 2304 –2314 (1996). APOPAI 0003-6935 Google Scholar


D. A. Boas et al., “Imaging the body with diffuse optical tomography,” IEEE Signal Process. Mag., 18 (6), 57 –75 (2001). ISPRE6 1053-5888 Google Scholar


J.-M. Sainthillier et al., “Geometrical capillary network analysis,” Skin Res. Technol., 9 312 –320 (2003). Google Scholar


R. Kiesslich et al., “New imaging techniques and opportunities in endoscopy,” Nature Rev. Gastroenterol. Hepatol., 8 547 –553 (2011). Google Scholar


S. Marschall et al., “Optical coherence tomography-current technology and applications in clinical and biomedical research,” Anal. Bioanal. Chem., 400 (9), 2699 –2720 (2011). ABCNBP 1618-2642 Google Scholar


M. J. Pittet and R. Weissleder, “Intravital imaging,” Cell, 147 983 –991 (2011). Google Scholar


J. Allen and K. Howell, “Microvascular imaging: techniques and opportunities for clinical physiological measurements,” Physiol. Meas., 35 (7), R91 –R141 (2014). PMEAE3 0967-3334 Google Scholar


G. Saiko and A. Douplik, “Real-time optical monitoring of capillary grid spatial pattern in epithelium by spatially resolved diffuse reflectance probe,” J. Innov. Opt. Health Sci., 05 (02), 1250005 (2012). Google Scholar


E. M. C. Hillman and S. A. Burgess, “Sub-millimeter resolution 3D optical imaging of living tissue using laminar optical tomography,” Laser Photonics Rev., 3 (1–2), 159 –179 (2009). Google Scholar


O. E. D. Online, “Oxford English Dictionary Online,” Oxford English Dictionary, (2010) Google Scholar


H. N. Chapman and K. A. Nugent, “Coherent lensless x-ray imaging,” Nature Photonics, 4 833 –839 (2010). Google Scholar


H. Jiang, G. Huang and P. Wilford, “Multi-view in lensless compressive imaging,” in Proc. 2013 Picture Coding Symposium (PCS 2013), 41 –44 (2013). Google Scholar


J. Spence, “Diffractive (lensless) imaging,” Science of Microscopy, 1196 –1227 Springer, New York (2007). Google Scholar


O. Mudanyali, W. Bishara and A. Ozcan, “Lensfree super-resolution holographic microscopy using wetting films on a chip,” Opt. Express, 19 17378 (2011). OPEXFF 1094-4087 Google Scholar


I. N. Papadopoulos et al., “High-resolution, lensless endoscope based on digital scanning through a multimode optical fiber,” Biomed. Opt. Express, 4 (2), 260 –270 (2013). BOEICL 2156-7085 Google Scholar


K. K. Hamamatsu Photonics, “FOP (fiber optic plates),” (2015) September ). 2015). Google Scholar


K. Kretsos and G. B. Kasting, “A geometrical model of dermal capillary clearance,” Math. Biosci., 208 (2), 430 –453 (2007). MABIAR 0025-5564 Google Scholar


F. Ayers et al., “Fabrication and characterization of silicone-based tissue phantoms with tunable optical properties in the visible and near infrared domain,” Proc. SPIE, 6870 687007 (2008). Google Scholar


J. Pekar and M. S. Patterson, “Fabrication and characterization of phantoms with tissue-like optical properties from 500 to 700 nm,” Med. Laser Appl., 25 147 –153 (2010). Google Scholar


R. C. Chang et al., “Fabrication and characterization of a multilayered optical tissue model with embedded scattering microspheres in polymeric materials,” Biomed. Opt. Express, 3 1326 (2012). BOEICL 2156-7085 Google Scholar


D. M. de Bruin et al., “Optical phantoms of varying geometry based on thin building blocks with controlled optical properties,” J. Biomed. Opt., 15 (2), 025001 (2010). JBOPFO 1083-3668 Google Scholar


S. A. Prahl, M. J. C. van Gemert and A. J. Welch, “Determining the optical properties of turbid mediaby using the adding–doubling method,” Appl. Opt., 32 (4), 559 –568 (1993). APOPAI 0003-6935 Google Scholar


G. Saiko et al., “Optical detection of a capillary grid spatial pattern in epithelium by spatially resolved diffuse reflectance probe: Monte Carlo verification,” IEEE J. Sel. Top. Quantum Electron., 20 (2), 1 –9 (2014). Google Scholar


Irina Schelkanova is an MSc and PhD student in the Physics Department of Ryerson University (Toronto, Canada). His area of research and interests covers a broad range of topics in photonics, biomedical optics, microcirculation, and cancer detection. Currently, he is working on the thesis: “Development of a Novel Miniature Fiber-Optics Technology for Quantification of Local Microvasculature Density Based on Hemoglobin Spatial Distribution.”

Aditya Pandya is a PhD student at the Physics Department in Ryerson University (Toronto, Ontario, Canada). His areas of interest include novel variants of modeling light propagation in tissues using Monte Carlo (with GPUs, parallelized, and efficient approaches), surface enhanced Raman spectroscopy (SERS), nano- and biophotonics, cancer detection, data mining, machine learning algorithms, astrophotography, and gesture and voice recognition.

Guennadi Saiko obtained his MSc and PhD degrees and is currently associated with the Physics Department at Ryerson University. Prior to that, he worked in the Erlangen Graduate School in Advanced Optical Technologies at Friedrich-Alexander University Erlangen-Nuremberg, Erlangen, Germany, and before that at Ontario Cancer Institute, Toronto, Canada. His interests are related to biotissue optics and bioimaging with a primary focus on cancer detection.

Duoaud Shah is currently a project manager with the Techna Institute at the University Health Network in Toronto, Canada. He earned his MSc in medical biophysics from the University of Toronto, where his work was focused on designing, fabricating, and using microfluidic devices for single cell analysis.

Lothar Lilge obtained his education in physics from the WJ Goethe University in Germany, and from the West, Wilhelms University in Germany. His training included 3 years at Mass General Hospital, Boston, and 2 years at McMaster University Canada. He is a faculty in medical biophysics at the UofT and senior scientist at the PMH (Canada). His research interests are photodynamic therapy, optical spectroscopy for cancer risk assessment and development of tools for microfluidics.

Alexandre (Sasha) Douplik obtained his MSc, PhD, Habilitationsfach (Doctor of Science) and is an associate professor of the Physics Department at Ryerson Technical University and an affiliate scientist of Keenan Research Centre of the LKS Knowledge Institute at St. Michael Hospital, both in Toronto, Canada. He is a guest professor of the Erlangen Graduate School in Advanced Optical Technologies at Friedrich-Alexander University Erlangen-Nuremberg, Erlangen, Germany. His interests are related to the field of laser therapy/surgery and spectral/fluorescence bioimaging.

Biographies for the other authors are not available.

© 2016 Society of Photo-Optical Instrumentation Engineers (SPIE) 1083-3668/2016/$25.00 © 2016 SPIE
Irina Schelkanova, Aditya Pandya, Guennadi Saiko, Lidia Nacy, Hannan Babar, Duoaud F. Shah, Lothar Lilge, and Alexandre Douplik "Spatially resolved, diffuse reflectance imaging for subsurface pattern visualization toward development of a lensless imaging platform: phantom experiments," Journal of Biomedical Optics 21(1), 015004 (14 January 2016).
Published: 14 January 2016

Back to Top