Open Access
5 December 2018 Portable, low-cost multispectral imaging system: design, development, validation, and utilization
Author Affiliations +
Optical spectral images can be used to estimate the amount of bulk absorbers in tissues, specifically oxy- and deoxyhemoglobin, as well as scattering parameters. Most systems that capture spectral image data are large, heavy, and expensive. This paper presents a full end-to-end analysis of a low-cost reflectance-mode multispectral imaging system operating in the visible and near-infrared spectra. The system consists of 13 LEDs mounted on a printed circuit board, a monochrome machine vision camera, and a tablet computer to control the hardware. The bill of materials for the system is less than $1000. Hardware design and implementation are detailed. Calibration, image capture, and preprocessing are also discussed. In validation experiments, excellent agreement is observed in diffuse reflectance measurements between the spectral camera setup and a spectrometer. To demonstrate that such spectral image data can yield meaningful optical measurements in vivo, the forearms of eight volunteers are imaged in the system. Their data are then analyzed to estimate the tissue optical properties of different skin layers using a Monte Carlo lookup table. In three volunteers, spectral images are captured before and after inducing erythema using a warm wet towel. Across the three subjects, a clear increase in the blood content of the superficial plexus layer was observed as a result of the erythema. Collectively, these findings suggest that a low-cost system can capture accurate spectral data and that clinically meaningful information can be derived from it.



Diffuse reflectance spectroscopy, where a broadband source illuminates a sample and the reflected spectrum is measured, is a well-trusted method for noninvasively characterizing biological tissues.1,2 The diffuse reflectance signal is a function of the wavelength-dependent bulk optical properties of the tissue. With proper modeling, it is possible to estimate the concentration of specific chromophores in the tissue. The scattering properties, which also affect the diffuse reflectance signal, can also be estimated.

In the visible and near-infrared (NIR) wavelength ranges, the main absorbers in tissue are water, oxy- and deoxyhemoglobin, and melanin. Assessing the concentrations of these chromophores can be used to identify pathological changes in tissues. Such measurements have been shown to identify lesions in several tissues, including breast,3,4 skin,5 oral cavity,6 bladder,7 esophagus,8 and the uterine cervix.9,10 Reflectance mode spectroscopy has also be used to characterize other chromophores in tissue, such as bilirubin when trying to assess jaundice in neonates.11 The main disadvantage of spectroscopy is that it is based on point measurements that average over a large area. Thus, a visual follow-up is often needed, particularly when clinical workflow requires selecting sites to biopsy.

The challenges posed by point measurements can be overcome by spectral imaging, in which images of the tissue of interest are captured at different wavelengths, resulting in a two-dimensional (2-D) array of spectral measurements. Each pixel has an associated spectrum that can be analyzed in search of pathological abnormalities, which allows clinicians to identify suspicious regions in the tissue for further examination.1214 Multispectral imaging has been primarily used in dermatological applications, to monitor portwine stain treatment,13 burns,15 and diabetic foot ulcers.16

Despite the advantages of multi- and hyperspectral imaging in obtaining clinically relevant information, few commercial systems that can assess oxygen saturation in images exist, and those that are FDA-cleared cost upward of $40,000 and are often difficult to use.1719 In contrast to these expensive, cart-filled systems from 15 years ago,20,21 the revolution in consumer electronics has made low-cost cameras and LEDs ubiquitous, bringing down the costs significantly. Several groups and companies have attempted to build a low-cost multispectral imaging system that can quantify scattering and absorption in the tissue,2224 including some based on LED illumination at various wavelengths,2528 in which spectral data were validated against spectral measurements from existing systems. However, a full end-to-end analysis that shows both the ability to reproduce a broad spectrum in an accurate way and yields reasonable in vivo scattering and absorption measurements with low-cost hardware has been elusive, and there have been questions about the overall feasibility of using low-cost hardware to make clinically relevant optical measurements.

In this work, we present the end-to-end functionality of a low-cost multispectral imaging system that captures 13 high-resolution images at different illumination wavelengths. The system is based on an area-scan camera and low-cost LEDs, such that the bill of materials is less than $1000. Fitting algorithms that analyzed the spectral data estimated the amount of each absorber (oxy- and deoxyhemoglobin, water, and melanin), as well as the overall scattering properties of the tissue. To validate our measurements, we performed a series of standard experiments used to characterize spectral systems. First, we compared spectral data measured with our system to white-light spectrometer readings of the 24 different color tiles from a color calibration card. We then captured in vivo images of skin from healthy volunteers, and reconstructed whole spectra using a lookup table based on Monte Carlo simulations. Finally, we captured images of a normal skin before and after erythema, demonstrating the ability of the measurements to noninvasively assess a change in tissue properties. Our results show that building a low-cost imaging system that can measure scattering and absorption is indeed feasible and that with proper calibration and modeling, such a system can be adapted to a specific clinical application.



A diffuse reflectance measurement on skin can be described by Ref. 29:

Eq. (1)

where S is the wavelength-dependent intensity of the source, Rskin is the reflectance of the skin, Cskin is the fraction of reflected light collected by the camera, and D is the wavelength dependent response of the detector (the camera). Similarly, the measurement on a reflectance standard (e.g., where Rstd=0.99) is given by Ref. 29:

Eq. (2)

where Cskin equals Cstd in these experiments as both skin and the reflectance standard were placed at the same aperture during measurements. The ratio M=Mskin/Mstd cancels the effects of S,D, and C,29 such that:

Eq. (3)

The parameter Rskin depends on the absorption and scattering properties of the different tissue layers in the skin (epidermis, papillary dermis, reticular dermis, etc.). The absorption properties of any tissue type are specified as

  • B = the blood volume fraction (B=1 for whole blood, 150 g hemoglobin/liter),

  • S = the oxygen saturation of hemoglobin in blood,

  • W = the water volume fraction (W=1 for pure water),

  • M = the volume fraction of typical cutaneous melanosomes in the epidermis,30

  • fi = the volume fraction of any additional (arbitrary) absorber in the tissue,

which are used to calculate the absorption coefficient:13

Eq. (4)


The absorption coefficients of oxy- and deoxyhemoglobin in whole blood [μa.oxy(λ) and μa.deoxy(λ)] are from Prahl.31 The absorption coefficient of water is from Hale and Querry,32 converted from the reported imaginary refractive index n into μa.water=4πn/λ[cm1]. The reported absorption coefficient of the interior of a cutaneous melanosome30 is described as μa.melanosome(λ)=(679.16  cm1)(λ/500  nm)3.33.

The scattering properties of any tissue type are given as:

  • μs500  nm = the reduced scattering at 500 nm,

  • f = the fraction of scattering at 500 nm due to structures much less than 500  nm (Rayleigh scattering),

  • (1f) = the fraction of scattering at 500 nm due to structures comparable to or >500  nm (Mie scattering),

  • bMie = the scattering power for such Mie scattering (typically <1),

  • g = the anisotropy of scattering,

which are used to calculate the reduced scattering coefficient:13

Eq. (5)

μs(λ)=μs500  nm[f(λ500  nm)4+(1f)(λ500  nm)bMie].

This μs is related to the scattering coefficient μs=μs1g. The values μa,μs, and g are used in a light transport model (Sec. 3.7).


Materials and Methods

Altogether, three sets of experiments were performed. To characterize the performance of the system, we captured spectral image data of different tiles from a color calibration card and compared them with point spectroscopic measurements. We then captured spectral images of skin (forearm) from eight healthy volunteers and reconstructed spectra that best fit the data. Finally, we captured spectral images of skin from three healthy volunteers before and after inducing erythema by placing a warm wet towel on the skin. In our analyses, we used a GPU version of the multilayer Monte Carlo model (CUDAMCML33) to estimate tissue composition, showing how local changes in optical properties of tissue layers can be determined from the system.


System Design

A schematic of our imaging system is shown in Fig. 1(a). Pictures of the entire system and optical head are shown in Figs. 1(b) and 1(c), respectively. Briefly, the imaging system can be subdivided into an illumination module, a detection module, and a tablet computer that controlled illumination and imaging by the system.

Fig. 1

(a) Imaging system block diagram, (b) full system setup with tablet, and (c) optical head.


In the illumination module, spectral illumination was provided by a set of LEDs across the visible and NIR regions mounted on a printed circuit board (PCB). Visible LEDs were LUXEON Z from Lumileds, NIR LEDs were LX1 from LED ENGIN. A list of the LED central wavelength and bandwidth [full width at half maximum (FWHM)] is shown in Table 1. Electrical current to drive the LEDs was provided by a simple low-side field effect transistor (FET) switch with current limited by a series resistor. As each LED has a unique voltage drop, resistor values were selected to match each LED. A microcontroller developer board (Arduino Nano) was used to interface the LED drive board to the computer system. Beam uniformity across LEDs was achieved using LEDs with a wide illumination angle without any focusing optics in the illumination path.34 In addition to spectral illumination, there was also a white LED (Cree XLamp) used to illuminate the tissue when the spectral LEDs were turned off (i.e., in preview mode), providing a view of the tissue to position the system properly; a white LED was also integrated into the multispectral PCB but not used in the analyses. This preview mode white LED was mounted separately, with its own PCB. A linear polarizer was installed in front of all the LEDs (Edmund Optics).

Table 1

Comparison of central wavelength (and FWHM) of LEDs reported by a vendor to spectrometer measurements. Note the white LED (a) was not used in the analyses described in Sec. 3.6.

LED name (internal)Reported by the vendor [nm]Measured by spectrometer [nm]
UV1395 (11)403 (12)
UV2425 (14)426 (14)
Royal blue447 (20)453 (20)
Blue470 (20)477 (27)
Cyan505 (30)505 (31)
Green530 (30)526 (38)
Lime568 (100)542 (89)
Amber590 (80)600 (18)
Red1627 (20)637 (16)
Red2655 (20)663 (19)
IR1740 (20)732 (38)
IR2850 (30)855 (34)
IR3940 (40)948 (66)

The detection module consisted of a low cost monochrome camera (Basler puA2500–14  μm) coupled mechanically to a 25-mm machine vision lens (Basler) and a linear polarizer (Midwest Optics). The camera polarizer was oriented orthogonally to the illumination module polarizer to eliminate glare. The camera and lens were selected such that 4.7  cm was imaged along the camera’s short axis at an appropriate working distance of 30 cm.

The entire system was controlled by a tablet computer (Acer) that runs custom application software, which allows the operator to preview the image scene, perform image capturing sequences, and store the data both locally and at a remote location for subsequent analyses. The bill of materials of the entire system, including the computer, was under $1000.


Spectral Image Acquisition

The spectral camera (Fig. 1) was set up at normal incidence to the target, at a distance of 30 cm. The LEDs were placed next to the camera port, where the center of the LEDs was 4 cm from the center of the LEDs, 7.5  deg. To ensure no specular glare is collected, the illumination and imaging ports were polarized orthgonally to each other. In total, five individual spectra were captured per sample and analyzed by the model for each sample/imaging condition. All the experiments were conducted in a dark room.

During an image capture procedure initiated by the user, the white LED (for preview mode) turned off, and the PCB LEDs were turned on and off sequentially, with image data captured at each illumination wavelength. Each LED had a different exposure time, determined by imaging the calibration standard and scaling appropriately. Dark frames were also captured. Altogether, four images were captured at each wavelength, with the exposure time varied linearly across the sequence to assess LED output and verify the measurements. A total sweep of all the LEDs (four sweeps of 13 LEDs) took 30  s.

For each of the experiments described below, spectral images and spectrometer readings were captured from the sample (Isample), as well as a dark frame (Idark). In each imaging session, we also captured a calibration image (Icalib), off a 99% spectralon reflectance standard (Labsphere). The reflectance was calculated using Eq. (6) for each sample:

Eq. (6)



Reflectance Spectroscopy Measurements

To validate our multispectral imaging setup, we mounted reference samples against a 3-cm diameter aperture and took measurements both using the spectral camera and using a white-light spectrometer (Flame-S 2000 spectrometer, Ocean Optics). Spectrometer-captured spectra (Fig. 2) were used for comparison. Complete spectra were captured with a 4.75 W Tungsten Halogen bulb white-light illumination light source (Ocean Optics). Light diffusely reflected from the sample was collected by a 400  μm multimode fiber, placed at normal incidence at a distance of 2 cm from the target. The light source was placed at a distance of 8 cm and at an angle of 30 deg from the normal such that the sample was illuminated obliquely, with the specular reflection pointing away from the probe. The probe only acquired the central part of the aperture to minimize effects of lateral light leakage, with an acquisition time of 300 mS.

Fig. 2

Spectrometer system for generating ground truth spectral data.



Color Card Imaging

Spectral data were captured from different tiles of a color calibration card (MacBeth chart, X-Rite ColorChecker Classic). Each of the 24 colored tiles were mounted against a 3-cm diameter aperture and imaged with our multispectral setup (Fig. 1) and with the spectrometer (Fig. 2). Pixels corresponding to a region of 2-×2-cm were averaged and calibrated according to Eq. (6); a similar process was carried out with the spectrometer measurements.


In Vivo Skin Imaging

Two sets of skin imaging experiments on healthy volunteers are presented here. First, to characterize the variation of tissue properties in normal skin, we captured multispectral images (Fig. 1) of forearm skin from eight healthy volunteers. All volunteers were Caucasian from European and Middle Eastern ethnicities and were 37±8.7-years old. The forearm was mounted against the aperture on the target. Calibration images of a white standard were acquired once at the beginning of the procedure.

In the second experiment, to verify that the system can assess changes to the skin’s tissue properties, we captured multispectral data before and after inducing erythema, which temporarily increases the blood content of the superficial skin. First, a “before” spectral image set was captured. To induce erythema, we placed a warm wet towel on the forearm for 2  min, until the skin turned visibly red. We then captured a second spectral image set immediately after the cloth was removed (“erythema”). Several pen marks were made on the participants’ forearms that matched up to an aperture that fixed the imaging distance, thus ensuring that the position and orientation were repeatable for before and after erythema images. A corresponding calibration image set was captured once for all three volunteers.


Spectral Data Analysis

Spectral data captured using the multispectral imaging system and the spectrometer were analyzed using scripts written in Python on a remote workstation. In each spectral dataset, the 2-×2-cm region of interest (ROI) at the center of the aperture was delineated for further analysis. Lateral (x,y) pixels in the ROI were averaged, and then fit to a six-layer Monte Carlo model described in Sec. 3.7, yielding a matrix of specified absorption parameters (B,S,W,M,μs500,f,bMie, and g) at each layer.


Monte Carlo Simulations

To relate the spectral reflectance measurements to tissue properties, we fit the spectral camera results using a Monte Carlo lookup table. The Monte Carlo simulations were based on a refined version of the CUDAMCML33 running on a GPU (GTX1060, Nvidia). To create the table, we modeled the tissue as a six-layer medium, representing the epidermis, basal layer, papillary dermis, dermal venous plexus, reticular dermis, and subcutaneous tissue. Within this six-layer tissue model, we considered the following parameters in each layer: B*,S*,W,M*, μs500  nm, f, bMie, and g,2 those marked with an asterisk were varied, details are shown in Table 3. Layer thickness was fixed for all analyses; epidermis (50 μm), basal layer (10 μm), papillary dermis (250 μm), dermal venous plexus (50 μm), reticular dermis (1 cm), and subcutaneous dermis (10 cm). A Python script was used to convert the table of tissue parameters to MCML input files files for 63 wavelengths, from 380- to 1000-nm in steps of 10 nm. Altogether, we simulated 1400 permutations at each of the 63 wavelengths.

To compare the measured spectra with the simulations, we converted the 63-point spectra into a 13-point spectra by multiplying (dot product) the simulations with each of the LED measured spectra, where each LED spectrum was normalized such that the sum of values was equal to unity. We then calculated the sum of squared error for each of the permutations and found the one with the lowest value. Note that because of the complexity of parameter space, in the erythema experiment, we fixed most of the parameters in the model, restricting the parameters that could change to blood content and oxygen saturation in the superficial plexus layer, as well as epidermal/basal melanin (keeping Mbas=3*Mepi35).




Multispectral Imaging System Validation

To validate the spectral camera system, spectral data were collected from the color calibration card using both the multispectral imaging system and the white-light spectrometer. A picture of the color calibration card is shown in Fig. 3 (top), and the comparison between the two spectral methods for each of the 24-color squares is shown in Fig. 3 (bottom), with the red curve representing the spectrometer, and the blue dots representing measurements from individual LEDs. As can be seen, there is a good agreement between the two systems on each of the 24 samples. We took this agreement as an indication that our LED spectral camera yields reliable spectral data.

Fig. 3

Photo of Colorchecker classic card (top), with a comparison of spectral camera to a spectrometer (bottom).



Multispectral Imaging of Healthy Skin

Using our multispectral imaging system, we imaged skin on the forearms of eight healthy subjects and fit the data to a theoretical model. Figure 4 compares LED data (green dots), the best fit spectrum (blue line), and the dot product of the best fit spectrum and the LEDs’ spectral shape (blue exes), to account for the finite LED bandwidth following the method described in Sec. 3.7. Good overall agreement can be seen point-wise between the spectral data and the best fit curve for the eight subjects (Fig. 4). Moreover, the estimated parameter concentrations are within normal physiological levels (Table 2). In assessing the performance of individual LEDs, the Lime LED (568 nm) and to a lesser degree the IR3 LED (940 nm) seem to be consistently higher than the fitted curves while maintaining a good match with the pointwise fit. We believe that this is due to the 89- and 66-nm FWHM of the LEDs (Table 1) and the spectral shape of the tissue, which begins to dramatically increase partway through the LED’s spectrum.

Fig. 4

Spectral camera measurements (green dots) versus estimated fit spectra (blue line) and the product of the LED spectrum and the estimated fit spectra (blue exes) for eight volunteers participants.


Table 2

Estimated average tissue parameters for eight volunteers (with standard deviation). For completeness, values of μa and μs at 500 nm are also shown, determined from Eqs. (4) and (5).

Epidermis0 (—)0.75 (—)0.75 (—)0.03 (0.003)40.00 (—)0 (—)1.00 (—)0.90 (—)20.38400.0
Basal layer0 (—)0.75 (—)0.75 (—)0.122 (0.01)40.00 (—)0 (—)1.00 (—)0.90 (—)82.86400.0
Papillary dermis0.0005 (—)0.55 (—)0.65 (—)0 (—)30.00 (—)0.62 (—)0.91 (—)0.90 (—)0.11300.0
Dermal venous plexus0.048 (0.028)0.819 (0.066)0.75 (—)0 (—)10.00 (—)0.62 (—)0.91 (—)0.90 (—)5.38100.0
Reticular dermis0.0005 (—)0.55 (—)0.65 (—)0 (—)30.00 (—)0.62 (—)0.91 (—)0.90 (—)0.11300.0
Subcutaneous tissue0.05 (—)0.55 (—)0.75 (—)0 (—)10.00 (—)0.60 (—)0.90 (—)0.90 (—)5.60100.0


Multispectral Imaging with Erythema

To determine if our spectral camera setup can identify changes in optical properties, we also captured additional image sets from three healthy volunteers before and after inducing erythema. The tissue properties corresponding to the spectral data were then determined from a Monte Carlo lookup table, assuming the tissue as a six-layer model (Sec. 3.7). The spectral data (dots), as well as reconstructed curves for the best fit lines (solid line) as captured before (blue) and after (red) erythema, are shown in Fig. 5.

Fig. 5

Spectral camera measurements (points) versus estimated spectral fit spectrum (line) for normal skin (blue, dots and solid line) versus induced erythema (red, exes and dashed line).


In analyzing spectral data from this experiment, we kept most of the parameters in the various layers fixed, and only varied the blood content and oxygen saturation in the superficial plexus layer, as well as the melanin content in the epidermal and basal layers (the ratio of which was fixed at 1:3).35 Numerical estimates for the fixed parameters in each layer are shown in Table 3. Plots of the sum of errors as a function of blood content at various oxygen saturation and melanin values are shown in Fig. 6. It can be seen that the blood content in the plexus layer is more than doubled as a result of inducing erythema. Additionally, the similarities between the curves representing different oxygen saturation levels show that our model and system are less sensitive to changes in oxygen saturation than to changes in blood content. Considering we have a few LEDs in the 600- to 700-nm range, which were identified to be problematic in Sec. 4.1, this result makes intuitive sense. Additionally, we also see that there is a coupling between the optimal blood content as a function of melanin concentration. Because both melanin and oxy- and deoxyhemoglobin absorb more at lower wavelengths,2 these observations also make intuitive sense.

Table 3

Base tissue parameters used to construct the tissue lookup table used in both experiments. Superscript v represents parameters that were varied. For completeness, both μa and μs at 500 nm are also given, calculated from Eqs. (4) and (5).

Basal layer00.750.750.12v40010.981.50400.0
Papillary dermis0.00050.550.650300.620.910.90.56300.0
Dermal venous plexus0.15v0.9v0.750100.620.910.916.81100.0
Reticular dermis0.00050.550.650300.620.910.90.06300.0
Subcutaneous tissue0.050.550.750100.

Fig. 6

Measurement of fitting error as a function of blood content, for different oxygen saturation levels (curves) and melanin values (a)–(l). Normal skin (blue) versus induced erythema (red). Separate plots are shown for each of the three volunteers. The vertical lines represent in each curve show where the error function was at a minimum.




The results presented here show the end-to-end characterization of a low-cost multispectral imaging system built with a $1000 bill of materials. The system was based on a PCB with 13 LEDs and an area scan camera that are controlled by a tablet. Altogether, we ran three sets of experiments to provide evidence that this system captures useful spectral data and allows for estimation of tissue chromophore content. The first experiment validated that the spectral data generated by our system yield comparable measurements to spectra measured with a spectrometer by comparing data from a color calibration card. Next, we used the system to capture image sets of the forearm of eight volunteers and fitted a tissue model (via a Monte Carlo lookup table) to the readings, demonstrating that the fitted values were within physiological ranges. Finally, we showed that our setup is capable of assessing a clear change in tissue composition (blood content of the superficial plexus layer) in spectral images captured before and after inducing erythema in three volunteers. Taken together, these results suggest that our portable, low-cost system is capable of reliably measuring spectral differences in tissue layers.

For several years, it has been known that spectral imaging systems do need to not be large or expensive, and indeed, multiple groups have built such systems.22,24,27,36,37 However, to our knowledge, these are the first end-to-end analyses that show it is possible to estimate tissue scattering and absorption properties with a low-cost system and that expensive hardware is not necessary for making clinically relevant measurements. Until now, validation of low-cost spectral analysis systems has focused primarily on the hardware, with less emphasis on the analysis and measuring optical properties. The tissue parameters measured by our system—blood content, oxygen saturation, and scattering from different tissue layers—are fundamental descriptors of tissue that provide a snapshot of its physiology on the micron scale. They provide information that clinicians would love to access but have not been able to, outside of a research setting at large university hospitals. Indeed, most systems that make such measurements that are in clinical use today are built on a cart.20,21,38 The small-form factor, light weight, and low power nature of our system allow it to be portable, carried by physicians who move between patients beds or between clinics. If successfully translated and productized, it can have a profound effect on clinical care.

The big unknown question with our setup is how does it compare with more expensive spectral imaging systems, specifically in terms of the ability to make scattering and absorption measurements. Here, Fig. 3 shows that reflectance from the LEDs matched very well with the spectrometer. Errors of 10% reflectance were only found on one LED (590 nm) on three of the 24 squares, representing <1% of the total measurements. Figures 46 show that even with inexpensive components, it is still possible to build a system capable of generating data that enables quantitative measurements of tissue scattering and absorption. The main compromise in our system was in the time required to prepare the data analysis. The time it took to iteratively run simulations to analyze a single spectral measurement was over two days using a modern GPU. This is because the process involved running tens of thousands of CUDAMCML simulations that took time and drive space. Comparing these data with spectra we measured took only a few seconds. The Monte Carlo lookup table method, while computationally expensive in our experiments, can be optimized to run much more efficiently using methods derived from artificial intelligence (AI),13,39 which will ultimately allow deployment on computers with less computational power. AI is particularly promising in this respect, as it utilizes a large Monte Carlo lookup table for training, yet the resulting classifier is not very large.

The main challenge with getting meaningful information out of the spectral measurements made by our setup is fitting it to the right model of light transport in tissue. Earlier attempts at fitting experimental data to a model were often based on a diffusion theory approximation.12 Although diffusion theory is straightforward for a semi-infinite medium, diffusion math gets complicated when there are multiple layers to the tissue. Because skin is a multilayer tissue, fitting algorithms based on diffusion theory can only be trusted so-far, and Monte Carlo approaches are preferred. We saw this with our spectral data (data not shown). Specifically, within the context of experiments like the one shown in Fig. 5, the analysis becomes challenging: inducing erythema on a semi-infinite tissue, the average increase in blood content to the entire tissue is not very large. However, within the context of a layered tissue model, it becomes very apparent that increasing (doubling) the blood content in the superficial plexus layer causes the change in spectral diffuse reflectance measured by the spectral camera. Thus, Monte Carlo-based approaches are required to make meaningful measurements using our system.

In assessing the accuracy of the in vivo analysis, Fig. 6 shows that the model and system are much more sensitive to changes in blood content and melanin content than to oxygen saturation. The reason for this is because the number and quality (stability, FWHM) of the specific LEDs needed to estimate oxygen saturation was limited. There were larger variations in output from the relevant LEDs relative to the other LEDs in the setup, specifically the Amber (590 nm). This can be seen in the color calibration data where this LED has a lower reflectance than measured by the spectrometer in the white, neutral 8, and neutral 6.5 squares but not the others (Fig. 3). The other important LED for this measurement, lime (568 nm), has a FWHM of 89 nm, and this large bandwidth exerted a non-negligible influence from nearby wavelengths on the measured reflectance, which reduced our sensitivity to measuring changes affecting that part of the spectrum.

In addition to challenges caused by the hardware, there are also some inherent challenges in working with a multilayer Monte Carlo model, as the fitting analysis becomes an underdetermined problem, with more unknowns than measurements. Such a problem requires making assumptions on tissue layers when oftentimes there is insufficient information in the literature. For example, differences in scattering properties between the papillary and reticular dermis significantly lowered our fitting error for one of the volunteers (data not shown), but making such educated guesses is not trivial, as there are not many reliable measurements of scattering by different dermal layers. Similarly, variations in the water content across the various layers have not been properly assessed. Rigorous studies of the actual values of many of these tissue parameters are necessary to improve measurement accuracy.

An alternative approach for assessing blood content and oxygen saturation in tissues is spatial frequency-domain imaging (SFDI)40 and its derivatives. SFDI holds two key advantages to unmodulated imaging (e.g., continuous wave) approaches: the ability to measure the absorption and scattering separately from each other and being able to control for the depth in which the measurement is made using the modulation frequency. However, the SFDI systems built to date are based primarily on red and NIR light and are thus less suitable for layered tissues, such as skin or the cervix, in which the layer of interest is superficial and spans only a couple of hundred microns in depth. Indeed, higher spatial frequencies and shorter wavelengths that have been demonstrated to date in vivo are required for making such measurements accurately. Our system uses shorter wavelength light (purple, blue, and green) that is more suitable for assessing changes in layered tissues. Moreover, while there have been recent advances in making low-cost SFDI devices (and variants),41,42 to our knowledge there has not yet been a noncontact, low-cost, small-form factor device validated in a clinical settings. Thus, while both methods make similar measurements and rely on similar light transport models for data analysis, they are suited for different clinical applications.

It is important to note that Monte Carlo models are specific to a particular tissue (i.e., skin), even though they can be easily adapted for alternative geometries of additional tissues. This is another issue to consider when designing such a system. Our system was designed for cervical imaging using a colposcopy imaging geometry. A system like the one described here was deployed in a cervical cancer clinic to look for cervical dysplasia.43 An exploratory trial has already began.

Additional future work will critically assess the combination of LEDs used for the spectral illumination to maximize performance at reduced cost. The current implementation is limited by assumptions made in the fitting method, which employs a limited set of tissue parameters. More extensive Monte Carlo methods will be implemented to better optimize the method, as a larger lookup table can provide finer estimates of tissue parameters.



We presented an end-to-end analysis of a low-cost multispectral imaging system, from design to development and validation, and measuring optical properties of tissue in vivo. Our results show that images captured with a low-cost spectral camera setup yield comparable fitting results as measurements captured by a spectrometer. Such measurements were made on 24 calibration samples, showing overlap in measured diffuse reflectance. We then demonstrated that it is possible to measure optical properties of the imaged tissue using a multilayer Monte Carlo model. Finally, using our system and analysis, we were able to show physical changes between normal forearm skin and erythema. We consider the experiments presented here as a proof of principle that it is possible to capture meaningful spectral measurements using a low-cost spectral camera set up. We believe our method holds considerable promise for clinical utility if implemented correctly.


This research was funded by MobileODT. FB, AB, and DL all own stock in MobileODT.


The authors wish to acknowledge the staff of MobileODT for volunteering to have their forearms imaged as part of the research and development process.



S. Svanberg, Atomic and Molecular Spectroscopy: Basic Aspects and Practical Applications, 6 Springer Science & Business Media, Berlin, Heidelberg (2012). Google Scholar


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


I. J. Bigio et al., “Diagnosis of breast cancer using elastic-scattering spectroscopy: preliminary clinical results,” J. Biomed. Opt., 5 (2), 221 –229 (2000). JBOPFO 1083-3668 Google Scholar


A. E. Cerussi et al., “In vivo absorption, scattering, and physiologic properties of 58 malignant breast tumors determined by broadband diffuse optical spectroscopy,” J. Biomed. Opt., 11 (4), 044005 (2006). JBOPFO 1083-3668 Google Scholar


A. A. M. Morales and S. V. Y. Montiel, “Retrieving the optical parameters of biological tissues using diffuse reflectance spectroscopy and Fourier series expansions. I. Theory and application,” Biomed. Opt. Express, 3 (10), 2395 –2404 (2012). BOEICL 2156-7085 Google Scholar


C.-W. Kan et al., “Comparison of linear and non-linear classifiers for oral cancer screening by optical spectroscopy,” in Proc. of AMIA Annual Symp., 1003 –1003 (2007). Google Scholar


J. R. Mourant et al., “Spectroscopic diagnosis of bladder cancer with elastic light scattering,” Lasers Surg. Med., 17 (4), 350 –357 (1995). LSMEDI 0196-8092 Google Scholar


D. S. Gareau et al., “Optical fiber probe spectroscopy for laparoscopic monitoring of tissue oxygenation during esophagectomies,” J. Biomed. Opt., 15 (6), 061712 (2010). JBOPFO 1083-3668 Google Scholar


Y. N. Mirabal et al., “Reflectance spectroscopy for in vivo detection of cervical precancer,” J. Biomed. Opt., 7 (4), 587 –595 (2002). JBOPFO 1083-3668 Google Scholar


N. M. Marín et al., “Diffuse reflectance patterns in cervical spectroscopy,” Gynecol. Oncol., 99 (3), S116 –S120 (2005). GYNOA3 Google Scholar


S. L. Jacques, D. G. Oelberg and I. Saidi, “Method and apparatus for optical measurement of bilirubin in tissue,” US Patent 5, 353, 790 (1994).


T. J. Farrell, M. S. Patterson and B. Wilson, “A diffusion theory model of spatially resolved, steady-state diffuse reflectance for the noninvasive determination of tissue optical properties in vivo,” Med. Phys., 19 (4), 879 –888 (1992). MPHYA6 0094-2405 Google Scholar


S. L. Jacques, R. Samatham and N. Choudhury, “Rapid spectral analysis for spectral imaging,” Biomed. Opt. Express, 1 (1), 157 –164 (2010). BOEICL 2156-7085 Google Scholar


J. Angelo et al., “Ultrafast optical property map generation using lookup tables,” J. Biomed. Opt., 21 (11), 110501 (2016). JBOPFO 1083-3668 Google Scholar


A. Mazhar et al., “Noncontact imaging of burn depth and extent in a porcine model using spatial frequency domain imaging,” J. Biomed. Opt., 19 (8), 086019 (2014). JBOPFO 1083-3668 Google Scholar


S. Anand et al., “Diffuse reflectance spectroscopy for monitoring diabetic foot ulcer-a pilot study,” Opt. Lasers Eng., 53 1 –5 (2014). Google Scholar


, “HyperView hypermed,” (2018) September 2018). Google Scholar


, “Near-infrared, non-invasive quantitative tissue oximeter–OxiplexTS200–ISS,” (2018) September 2018). Google Scholar


M. Imaging, “Ox-Imager CTM–modulated imaging, Inc.,” (2018) September 2018). Google Scholar


M. J. DeWeert et al., “Analysis of spatial variability in hyperspectral imagery of the uterine cervix in vivo,” Proc. SPIE, 4959 67 –77 (2003). PSISDG 0277-786X Google Scholar


H. Lange et al., “Reflectance and fluorescence hyperspectral elastic image registration,” Proc. SPIE, 5370 335 –346 (2004). PSISDG 0277-786X Google Scholar


S. Kim et al., “Smartphone-based multispectral imaging: system development and potential for mobile skin diagnosis,” Biomed. Opt. Express, 7 (12), 5294 –5307 (2016). BOEICL 2156-7085 Google Scholar


F. Vasefi et al., “Quantifying the optical properties of turbid media using polarization sensitive hyperspectral imaging (skinspect): two-layer optical phantom studies,” Proc. SPIE, 9328 93280A (2015). PSISDG 0277-786X Google Scholar


J. Spigulis, “Multispectral, fluorescent and photoplethysmographic imaging for remote skin assessment,” Sensors, 17 (5), 1165 (2017). SNSRES 0746-9462 Google Scholar


M. Goel et al., “Hypercam: hyperspectral imaging for ubiquitous computing applications,” in Proc. of the ACM Int. Joint Conf. on Pervasive and Ubiquitous Computing, 145 –156 (2015). Google Scholar


G. Belušič et al., “A fast multispectral light synthesiser based on LEDs and a diffraction grating,” Sci. Rep., 6 32012 (2016). SRCEC3 2045-2322 Google Scholar


R. Shrestha and J. Y. Hardeberg, “How are LED illumination based multispectral imaging systems influenced by different factors?,” Lect. Notes Comput. Sci., 8509 61 –71 (2014). LNCSD9 0302-9743 Google Scholar


A. S. Bernat et al., “Cloud-based processing of multi-spectral imaging data,” Proc. SPIE, 10055 1005505 (2017). PSISDG 0277-786X Google Scholar


S. L. Jacques, “Spectral imaging and analysis to yield tissue optical properties,” J. Innovative Opt. Health Sci., 2 (02), 123 –129 (2009). Google Scholar


S. L. Jacques and D. J. McAuliffe, “The melanosome: threshold temperature for explosive vaporization and internal absorption coefficient during pulsed laser irradiation,” Photochem. Photobiol., 53 (6), 769 –775 (1991). PHCBAP 0031-8655 Google Scholar


S. A. Prahl, “Optical absorption of hemoglobin,” (1996) September 1999). Google Scholar


G. M. Hale and M. R. Querry, “Optical constants of water in the 200 nm to 200 um wavelength region,” Appl. Opt., 12 (3), 555 –563 (1973). APOPAI 0003-6935 Google Scholar


E. Alerstam, T. Svensson and S. Andersson-Engels, CUDAMCML-User Manual and Implementation Notes, 1 –33 Department of Physics, Lund University, Sweden (2009). Google Scholar


F. J. Bolton et al., “Development and bench testing of a multi-spectral imaging technology built on a smartphone platform,” Proc. SPIE, 9699 969907 (2016). PSISDG 0277-786X Google Scholar


T. Fitzpatrick, “Conversation about melanin concentration,” Private Communication, (1984). Google Scholar


N. MacKinnon et al., “Melanoma detection using smartphone and multimode hyperspectral imaging,” Proc. SPIE, 9711 971117 (2016). PSISDG 0277-786X Google Scholar


J. Spigulis et al., “Smartphone snapshot mapping of skin chromophores under triple-wavelength laser illumination,” J. Biomed. Opt., 22 (9), 091508 (2017). JBOPFO 1083-3668 Google Scholar


D. Medical, “0354 01214 US DySIS product brochure FINAL,” (2015) November 2017). Google Scholar


A. Doronin et al., “Determination of human skin optical properties from hyper spectral data with deep-learning neural networks (conference presentation),” Proc. SPIE, 10489 104890U (2018). PSISDG 0277-786X Google Scholar


D. J. Cuccia et al., “Quantitation and mapping of tissue optical properties using modulated imaging,” J. Biomed. Opt., 14 (2), 024012 (2009). JBOPFO 1083-3668 Google Scholar


S. Nandy et al., “Quantitative multispectral ex vivo optical evaluation of human ovarian tissue using spatial frequency domain imaging,” Biomed. Opt. Express, 9 (5), 2451 –2456 (2018). BOEICL 2156-7085 Google Scholar


R. B. Saager et al., “Portable (handheld) clinical device for quantitative spectroscopy of skin, utilizing spatial frequency domain reflectance techniques,” Rev. Sci. Instrum., 88 (9), 094302 (2017). RSINAK 0034-6748 Google Scholar


A. S. Bernat et al., “Scattering and absorption measurements of cervical tissues measures using low cost multi-spectral imaging,” Proc. SPIE, 10485 104850D (2018). PSISDG 0277-786X Google Scholar


Frank J. Bolton is a biomedical engineer at Mobile ODT and is also currently a PhD student at Tel Aviv University. He received his bachelor’s degree in electrical engineering from the University of Johannesburg and his MS degree in biomedical engineering from the University of Cape Town in 2007 and 2011, respectively. Currently, his research interests include spectral imaging, Monte Carlo modeling, and cognitive aspects of computing. He is a member of SPIE.

Amir S. Bernat is an algorithm engineer at MobileODT. He received his MSc in physics from Ben-Gurion University and has experience performing hardware development at leading tech companies. He performs R&D on MobileODT’s multispectral LED-based camera and the optical characterization of MobileODT’s EVA-system, a smartphone based colposcope. His main professional interest is combining biophotonics with off-the-shelf imaging devices and integrating analysis algorithms focused on spectral imaging. He is a member of SPIE.

Biographies for the other authors are not available.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Frank J. Bolton, Amir S. Bernat, Kfir Bar-Am, David Levitz, and Steven Jacques "Portable, low-cost multispectral imaging system: design, development, validation, and utilization," Journal of Biomedical Optics 23(12), 121612 (5 December 2018).
Received: 8 July 2018; Accepted: 9 October 2018; Published: 5 December 2018 Logo
Cited by 21 scholarly publications and 1 patent.

Light emitting diodes

Imaging systems


Multispectral imaging


Monte Carlo methods

Back to Top