Research Papers: Imaging

Ultrafast wavelength multiplexed broad bandwidth digital diffuse optical spectroscopy for in vivo extraction of tissue optical properties

[+] Author Affiliations
Alyssa Torjesen, Raeef Istfan, Darren Roblyer

Boston University, Department of Biomedical Engineering, Boston, Massachusetts, United States

J. Biomed. Opt. 22(3), 036009 (Mar 10, 2017). doi:10.1117/1.JBO.22.3.036009
History: Received December 22, 2016; Accepted February 17, 2017
Text Size: A A A

Open Access Open Access

Abstract.  Frequency-domain diffuse optical spectroscopy (FD-DOS) utilizes intensity-modulated light to characterize optical scattering and absorption in thick tissue. Previous FD-DOS systems have been limited by large device footprints, complex electronics, high costs, and limited acquisition speeds, all of which complicate access to patients in the clinical setting. We have developed a new digital DOS (dDOS) system, which is relatively compact and inexpensive, allowing for simplified clinical use, while providing unprecedented measurement speeds. The dDOS system utilizes hardware-integrated custom board-level direct digital synthesizers and an analog-to-digital converter to generate frequency sweeps and directly measure signals utilizing undersampling at six wavelengths modulated at discrete frequencies from 50 to 400 MHz. Wavelength multiplexing is utilized to achieve broadband frequency sweep measurements acquired at over 97 Hz. When compared to a gold-standard DOS system, the accuracy of optical properties recovered with the dDOS system was within 5.3% and 5.5% for absorption and reduced scattering coefficient extractions, respectively. When tested in vivo, the dDOS system was able to detect physiological changes throughout the cardiac cycle. The new FD-dDOS system is fast, inexpensive, and compact without compromising measurement quality.

Figures in this Article

Diffuse optical spectroscopy (DOS) has emerged over the past three decades as an important noninvasive method to characterize optical scattering and absorption and, in turn, calculate chromophore concentrations in thick tissue.1,2 Diffuse optical methodologies can be broadly categorized into three distinct measurement types: time domain, frequency domain (FD), and continuous wave (CW). This work focuses on FD-DOS, which allows for the separation of optical scattering and absorption through the measurement of amplitude attenuation and phase delay of photon density waves generated by temporally intensity modulated light sources at the tissue surface. Prior FD-DOS systems have utilized either a single modulation frequency39 or a sweep of multiple frequencies57 in the radio frequency (RF) range. Most single modulation frequency systems utilize either homodyne4,7,8,10 or heterodyne3 techniques to downconvert the high frequency signal of interest by mixing it with a reference signal. Broad bandwidth DOS instruments, which utilize light sources that are swept through a range of RF frequencies, have typically utilized a benchtop network analyzer to generate and measure signals up to the gigahertz range.1,11 FD-DOS devices based around network analyzers have relatively large instrument footprints, and the cost of fabrication can range from approximately $30,000 to $60,000 (based on our internal cost estimates), limiting access to patients in the clinic. Although some analog systems have been successfully miniaturized,12,13 the difficulties involved in RF system design potentially limit the adoption of these techniques. Furthermore, FD-DOS instruments often have slow measurement acquisition speeds compared to their CW counterparts, and although recent improvements have been made in data acquisition rates for single modulation frequency measurements,14 data acquisition speeds are limited for broad modulation frequency bandwidth FD-DOS systems.1,10,15

We recently demonstrated the feasibility of FD-DOS measurements using all digital signal generation and detection using direct digital synthesizer (DDS) integrated circuits and a 3.6-gigasample/s two-channel analog-to-digital converter (ADC).16 These digital components reduced system costs, but measurement acquisition speed was limited due to data transfer bottlenecks between system components. Furthermore, the high sampling rate of the ADC resulted in large datasets and lengthy data transfer times. A second implementation of a digital FD-DOS system tested the feasibility of utilizing a slower ADC, sampling at only 25  megasamples/s, and utilizing undersampling techniques.17 This substantially reduced system complexity but measurement speed was still a major limitation due to limited data transfer rates. Here, we present a new FD-DOS system that overcomes these past limitations by fully integrating digital signal synthesis and detection to rapidly acquire multiwavelength, multiplexed, broadband frequency sweep measurements at a repetition rate of up to 97 Hz. The digital DOS (dDOS) system was validated relative to a gold-standard network analyzer based DOS system and tested in vivo. This new high-speed, broad-bandwidth, digital, miniaturized DOS system minimizes analog circuitry while reducing cost and device footprint, potentially enabling access to more patients in the clinic and expanding the capabilities of this technology, especially for monitoring fast physiological changes.

Hardware Integration

Figure 1 shows the core components and signal/data flow for the new dDOS system. The hardware components of the dDOS system are integrated into a single motherboard that both generates and measures RF signals while communicating and transferring data to the host computer. A development board (ZedBoard, MircoZed Zynq-7010 SoC) was used as the central processing and data storage core for the system. It contains an ARM Cortext-A9 processor and an Artix 7 FPGA, 1 GB of RAM to store data read from the ADC, and a Linux server to transfer data between the board and host computer. The user sets the desired frequency sweep settings in software and the Linux server configures the DDS boards accordingly. The two-channel ADC samples both channels simultaneously and transfers data to the RAM on the MicroZed module via a parallel low voltage differential signal bus. Data are then transferred from the Linux server to the host computer via an Ethernet link at a rate of up to 170 Mbps.

Graphic Jump Location
Fig. 1
F1 :

Schematic overview of dDOS system. GUI, graphical user interface; ADC, analog-to-digital converter; FPGA, field programmable gate array; RAM, rapid access memory; DDS, direct digital synthesizer; APD, avalanche photodiode.

Light Sources and Modulation

The dDOS system utilizes NIR laser diodes with wavelengths of 658, 690, 785, 808, 830, and 850 nm that are housed in thermoelectric cooling modules (Thorlabs, LDM9T). All laser diodes are coupled to 400-μm core diameter fibers bundled into a single ferrule. Each laser diode is driven with a DC current, which is mixed with an RF current using a bias-tee to intensity-modulated light. DC current is supplied by an eight-channel, high-stability laser diode controller (ILX Lightwave, LDC-390, Newport Corporation). The dDOS system replaces the RF current output functionality of the network analyzer by implementing six DDS chips (AD9910, Analog Devices), each on individual daughterboards. Each DDS can output up to 20 mA of current modulated up to 400 MHz. The RF output of each DDS board is low-pass filtered at 400 MHz and routed to a 71 directional coupler. The directional coupler sends the majority of the output power to the bias-tee, and the lesser output from each directional coupler is routed to a power combiner (Mini-Circuits, ZBSC-615+) and measured as a reference signal with the second channel of the ADC.

The dDOS system is capable of utilizing the six individual DDS chips to simultaneously modulate a bank of laser diodes at offset frequencies. The data are then demultiplexed in the FD in postprocessing. A typical dDOS measurement, for example, may sweep through frequencies between 50 and 400 MHz in steps of 1 MHz with the modulation frequency for each wavelength offset by 10 MHz. This multiplexing reduces measurement duration to less than 100  ms/frequency sweep (in the slowest case), allowing for detection of rapid physiological changes. For this work, modulation frequencies between 50 and 300 MHz or between 50 and 400 MHz were used depending on the optical attenuation of the sample and data acquisition rate requirements. The frequency step size can be set between 1 and 7 MHz. For wavelengths with modulation sweeps starting above 50 MHz, the modulation frequency was incremented up to the maximum frequency then wrapped back to 50 MHz and incremented to the starting frequency. Data were unwrapped such that modulation frequencies appear in increasing order.

Detection

Reflectance mode geometry was used for all dDOS measurements, with the source fiber bundle and active area of an avalanche photodiode (APD) placed directly on the surface of the phantom at a source detector separation of 10 to 30 mm depending on the specific measurement. An APD module (Hamamatsu C5658) with a 0.5-mm active area, 2.50×105  V/W photoelectric sensitivity, and a high band cutoff of 1 GHz was used for most phantom measurements. For in vivo experiments, data were acquired with a 3-mm active area APD.

A 14-bit, two-channel ADC (Texas Instruments, ADS62P49) with a 2 V peak-to-peak full-scale input voltage was used to sample the signal at each channel at 250 MSPS. The electrical output of the APD was high-pass filtered at 41 MHz and routed to the first channel of the ADC, while the combined reference signal from the power combiner was routed to the second channel of the ADC. The user can select the number of samples collected by each channel of the ADC at each frequency step (e.g., 4096,8192,163842n). Increasing the number of samples per step substantially increases data transfer and processing time. Both channels of the ADC run off of the same clock and are sampled simultaneously.

Signal Processing

All signal processing is performed in MATLAB (R2014b, Mathworks Natick, Massachusetts). Since the two-channel ADC has a sampling rate of 250 MSPS, all signals modulated above 125 MHz are aliased. Given that the modulation frequency is known for each wavelength at each frequency step, the signal of interest is easily located in the baseband. The processing algorithm eliminates modulation frequencies that are multiples of the Nyquist frequency and removes any measurements in which multiple wavelengths are aliased to the same frequency in the baseband (e.g., signals at 120 and 130 MHz are both read as 120 MHz and eliminated).

A fast Fourier transform (FFT) with a rectangular window was used at each frequency step to calculate the amplitude and phase for the reference and sample channels. The appropriate frequency bin was located for each wavelength at each modulation frequency step. The magnitude of the FFT of the measurement channel relative to that of the reference channel was considered the raw amplitude value, and, similarly, the phase offset between the measurement and reference channels was considered the raw phase value.

Measurement Calibration and Optical Property Extraction

The raw amplitude and phase represent the amplitude attenuation and phase delay induced both by the tissue and the instrument. In order to remove the instrument response, a measurement was taken on a silicone phantom with known optical properties, where calibration optical properties were obtained using multidistance measurements. An analytical solution to the P1 approximation of the Boltzman transport equation with boundary conditions for semi-infinite geometry was employed as a forward model for the phantom measurement.5 Calibration factors for amplitude and phase were determined by ratiometrically or differentially comparing the measured amplitude and phase, respectively, to the theoretical amplitude and phase from the forward model. These calibration factors were then applied to the subsequent raw amplitude and phase measurements in order to obtain calibrated amplitude and phase.

An iterative least-squares fitting algorithm was employed to minimize the error between the forward model calculated and measured calibrated amplitude and phase to extract absorption (μa) and scattering (μs) parameters at each wavelength.5 Chromophore extinction coefficients were obtained from Zijlstra and Buursma.18 The Beer–Lambert law was used to extract chromophore concentrations (oxy- and deoxyhemoglobin) from multiwavelength μa data.13

System Overview

The completed dDOS system is shown in Fig. 2, where the motherboard, DDS daughterboards, MicroZed module, ADC, and RF output connections are pictured. The primary hardware is housed in an enclosure measuring 10×12×6  in.

Signal to Noise

Signal to noise was measured on a breast-simulating silicone optical phantom. SNR was compared for measurements taken using either a single laser or all six lasers modulated simultaneously. A comparison was also done for measurements taken using either 4096 or 8192 samples. The noise floor was calculated as the average noise level, in dBc, relative to the frequency of interest. Figure 3 shows that when all six lasers were modulated simultaneously, the noise floor was 46.3  dBc, whereas utilizing a single laser yielded a lower noise floor of 50.1  dBc. When increasing the sample length to 8192 samples per frequency step, the noise floor was reduced to 49.2 and 54.5  dBc for multiplexed and single wavelength measurements, respectively. Although the noise floor was lower with 8192 samples per frequency step, acquisition and data transfer times of these measurements were substantially slower than of those with 4096 samples per frequency step, and there was no improvement in accuracy of optical properties (e.g., the accuracy of 10 measurements on two phantoms with 4096 samples per frequency step was 1.29% compared to 1.85% with 8192 samples per frequency step). 4096 samples per frequency step were collected for all subsequent measurements to optimize speed.

Graphic Jump Location
Fig. 3
F3 :

Noise floor when collecting 4096 samples per frequency step in (b) a single wavelength measurement or (a) a multiplexed measurement with a 10 MHz offset between wavelengths.

Accuracy of Optical Property Extractions

An example of calibrated (black) and fit (red) amplitude and phase measured from 50 to 400 MHz in steps of 1 MHz at each wavelength is shown in Fig. 4. Data were calibrated using a silicone phantom with known optical properties in order to remove the instrument response. Measurements at 658 and 850 nm have somewhat higher noise levels compared to other wavelengths, likely due to the lower modulation index achieved for these laser diodes.

Graphic Jump Location
Fig. 4
F4 :

Amplitude (in arbitrary units, A.U.) and phase measured with dDOS system at each wavelength. Red lines are calibrated data and black lines are model fits.

dDOS optical property extraction accuracy was determined by measuring a set of 10 tissue-simulating silicone phantoms with varying amounts of titanium dioxide and nigrosin, as the scattering and absorbing agents, respectively, using both the dDOS and the network analyzer-based benchtop DOS systems. For each phantom, 10 FD measurements were taken with each system with a 15-mm source–detector separation and all six wavelengths modulated at frequencies between 50 and 400 MHz with 1-MHz frequency steps. Measurements on one highly attenuating phantom at 658, 690, and 850 nm were excluded due to phase errors resulting from low detected signal levels. The mean of the absolute difference in optical properties for all measured phantoms between the two systems was 5.3% and 5.5% for μa and μs, respectively; when frequency sweeps from 50 to 300 MHz were used, the mean accuracy in optical properties measured with the two systems was 4.9% and 6.0% for μa and μs, respectively.

Figure 5 shows scatterplots (top panels) and Bland–Altman plots (bottom panels) for μa and μs measured with the dDOS and benchtop DOS systems. For the Bland–Altman plots, the differences between optical properties with the dDOS and benchtop system are plotted on the y-axis versus average optical properties for each system on the x-axis, for each wavelength and phantom. 96% of μa measurements and 98% of μs measurements fell within 1.96 standard deviations of the mean difference, indicating acceptable equivalence between the two systems.19μa and μs both tend to skew slightly lower for the dDOS system compared to the benchtop DOS system. μa values had the best agreement between systems when the mean absorption coefficient was <0.015  mm1, while μs had the best agreement between systems when the mean reduced scattering coefficient was >1  mm1. In healthy breast tissue, μa in the NIR is typically <0.01  mm1 and μs in the NIR is typically >0.8  mm1.20

Graphic Jump Location
Fig. 5
F5 :

(a and b) Scatterplots and (c and d) Bland–Altman plots of optical properties obtained with the dDOS system compared to those obtained with the benchtop system. Top and bottom dashed lines represent 1.96 standard deviations from the mean difference; center dashed lines represent the mean difference.

Precision

The stability of amplitude, phase, and optical properties was characterized during a 2-h drift test. Measurements were taken at 1-min intervals on a silicone optical phantom using 658, 690, 785, 808, 830, and 850 nm lasers modulated simultaneously in frequency sweeps from 50 to 400 MHz with a 15-mm source–detector separation. Amplitude and phase were examined for the duration of the drift test for the 850-nm laser modulated at 50 MHz and were plotted in Fig. 6. For this wavelength and modulation frequency, the standard deviation of the phase was 0.38 deg and the standard deviation of the amplitude relative to the mean amplitude was 0.81%.

Graphic Jump Location
Fig. 6
F6 :

2-h drift test showing dDOS system stability. Measurements were taken approximately every minute and show variation from value at the first time point (t=0). (a) Amplitude drift at 50 MHz and 850 nm and (b) phase drift at 50 MHz and 850 nm.

System stability was further analyzed for all wavelengths and modulation frequencies used in the drift test, and results are shown in Table 1. Amplitude and optical property stability were defined as the standard deviation divided by the mean at each modulation frequency, averaged for all wavelengths. Phase stability was defined as the standard deviation throughout the drift test averaged for all wavelengths and modulation frequencies. The stability in μa and μs was 2.0% and 0.8%, respectively.

Table Grahic Jump Location
Table 1System stability during a 2-h drift test with measurements taken approximately every minute.
Speed

Most prior FD-DOS instruments performed frequency sweeps with each wavelength modulated individually in a sequential manner.5,17 The dDOS system utilizes wavelength multiplexing, which eliminates switching time between lasers and reduces data acquisition time, resulting in the ability to measure fast physiological changes. The rate at which frequency sweep measurements can be performed is dependent on the range of modulation frequencies, the step size, and the number of samples acquired at each frequency step. When multiplexing six laser wavelengths, sweeping from 50 to 300 MHz with 7 MHz steps, and acquiring 4096 samples per frequency step, the maximum speed achieved for dDOS system was 97.2 Hz. Increasing the frequency sweep range from 50 to 400 MHz reduces the measurement rate to 68.1 Hz. When decreasing the frequency step size to 1 MHz, the 50 to 300 MHz sweep can be acquired at 13.9 Hz, and the 50 to 400 MHz sweep can be acquired at 9.8 Hz. These measurement speeds allow for FD-dDOS data acquisition that can capture fast physiological changes, including during the cardiac cycle.

To test the maximum measurement speed attainable with the dDOS system, measurements were also taken with a single modulation frequency. When modulating at only 50 MHz, a measurement rate of 2557 Hz can be achieved. Although this is not the preferred data acquisition mode, because utilizing additional modulation frequencies improves model fit accuracy, it provides a useful comparison to other FD systems, which utilize single modulation frequencies.3,4,8,10,14,21

A summary of the performance specifications is given in Table 2.

Table Grahic Jump Location
Table 2dDOS performance specifications.

All in vivo human measurements were carried out under a protocol approved by the Boston University Institutional Review Board.

Cuff Occlusion Measurement

An in vivo cuff occlusion measurement was performed using a blood pressure cuff placed on the upper arm of a male volunteer, while optical measurements were taken on the forearm. The cuff was loosely secured in order to measure baseline molar concentrations of oxyhemoglobin (ctO2Hb) and deoxyhemoglobin (ctHHb). The cuff was then was inflated to 200 mmHg for 2  min, and then released and data were acquired for 5 additional minutes. Frequency sweeps between 50 and 400 MHz, with frequency steps of 1 MHz, were used to modulate all six laser diodes simultaneously. Measurements were taken approximately every 8 s for 10 min and a source–detector separation of 10 mm was used for each measurement.

Figure 7 shows ctO2Hb and ctHHb throughout the cuff occlusion test. At baseline average, ctO2Hb and ctHHb were 41.6 and 41.2  μM, respectively. At the time of cuff release, ctO2Hb had dropped to 19.3  μM and ctHHb had increased to 87.6  μM. After cuff release, ctO2Hb rebounded to a local maximum of 67.7  μM and ctHHb dropped to 36.5  μM, before both chromophore concentrations began to return toward their baseline values. These results are consistent with reduced arterial supply and vascular drainage during cuff occlusion, followed by a rush of arterial blood into the forearm after cuff release.

Graphic Jump Location
Fig. 7
F7 :

Forearm cuff occlusion test. Cuff was occluded during time interval between dashed lines. Measurements were taken approximately every 8 s and a three-point moving average is plotted for each chromophore.

Rapid In Vivo Measurements

To test the ability of the dDOS system to measure physiological changes at the cardiac rate, frequency sweep measurements were taken on the finger of a male volunteer. Measurements were taken in reflectance mode on the index finger with a 10-mm source–detector separation. Frequency sweeps from 50 to 300 MHz with 7 MHz steps were used to modulate all six laser diodes simultaneously. About 350 frequency sweeps were taken with a measurement rate of 97.2 Hz.

An electrocardiogram (ECG) signal was measured simultaneously with the optical signal and synced using an external data acquisition board. Figure 8 shows five-point moving averages of μa, μs, ctTHb, ctO2Hb, and ctHHb from the high-speed finger measurement temporally aligned with the ECG. μa values fluctuate temporally with the ECG signal, likely due to increased blood volume in the finger after ventricular contraction. Similarly, ctTHb and ctO2Hb vary temporally with the ECG signal, while ctHHb shows minimal fluctuation throughout the cardiac cycle. Vertical dashed lines mark the peak of the ECG R-wave, the start of the upstroke of the pulsatile waveform, and the dicrotic notch. There is a 155-ms lag in the upstroke of the μa signal after the R-wave of the ECG, due to the transmission time of the pressure wave to the peripheral vasculature. The delay is often referred to as the pulse transit time and previously shown to range between 200 and 300 ms for photoplethysmographic signals measured in the fingertip.22 The dicrotic notch, which is correlated with end of the systolic ejection period, is likely correlated with the dip in μa, ctTHb and ctO2Hb marked by the rightmost vertical line. In this case, the dichrotic notch appears roughly 303 ms after the upstroke of the pulsatile waveform, in agreement with other report values in the literature.23,24 Fluctuations in μa, ctTHb, and ctO2Hb appear to occur at the same frequency as the ECG signal, suggesting that the dDOS system is capable of detecting pulsatile, oxygenated, arterial blood flow.

Graphic Jump Location
Fig. 8
F8 :

High-speed dDOS measurements were taken on the index finger. An ECG signal was obtained simultaneously with the optical measurements and temporally aligned with optical property and chromophore values. From top to bottom, plots represent the ECG, μa, μs, ctTHb, ctO2Hb, and ctHHb. The leftmost vertical dashed line represents the approximate time of the peak of the ECG R-wave, the middle dashed line represents the approximate time of the start of the upstroke of the optical signal following the arterial transit time, and the rightmost dashed line represents the approximate timing of the dicrotic notch based on visual inspection.

Chromophore Maps in Healthy Breast Tissue

Optical properties of breast tissue were measured on a healthy female volunteer. Measurements were taken with the dDOS system using 658, 690, 785, 808, 830, and 850 nm lasers with a 28-mm source–detector separation and modulation frequencies from 50 to 400 MHz. 36 point measurements were taken in which the source–detector pair was manually scanned across the right breast in a 6×6  cm grid pattern with 1 cm spacing between measurement points. A quadrant of the areolar complex was located at the top left corner of the imaging grid. Optical properties were calculated and chromophore values determined for each data point. Data from the 658- and 850-nm wavelengths were excluded due to noise in the phase measurement at these wavelengths. Figure 9 shows a chromophore map for ctTHb with linear interpolation between measurement points. The chromophore map was superimposed on a 3-D scanned bust model in order to show the approximate relative position of the measurement. There is a region of increased ctTHb at the top left corner of the chromophore map, likely due to increased blood supply to the areola as compared to standard breast tissue, which has been observed in prior DOS studies.25

Graphic Jump Location
Fig. 9
F9 :

ctTHb measured on the right breast of a healthy volunteer. A quadrant of the nipple is located in the lower left corner of the image.

A new dDOS system was developed to increase access to patients, while adding new capabilities in measurement speed. Traditional network analyzer-based FD-DOS systems are large and costly, and miniaturization of analog RF circuits is technically challenging. The hardware-integrated dDOS system developed here utilizes direct digital signal generation and detection and wavelength multiplexing, with frequency sweep measurements rates up to 97 Hz, allowing rapid physiological changes to be detected well above the cardiac rate. Previous systems achieving high sampling rates have utilized a small number of discrete modulation frequencies and wavelengths,14,21 whereas the dDOS system sweeps a range of hundreds of modulation frequencies and utilizes six discrete wavelengths.

The dDOS system performance is comparable to previous systems. For example, Pham et al. reported an accuracy of ±5% and ±3% for μa and μs, respectively, whereas the dDOS system had an accuracy of ±5.3% and ±5.5% for the same parameters. The dDOS design has similarities to that of the modular digital FD diffuse optical monitor developed by Durduran et al., which utilized digital signal synthesis and a dual-channel ACD to sample and downconvert the signal from three laser wavelengths at a single modulation frequency;21 however, our design differs in that we use six multiplexed laser wavelengths modulated at hundreds of discrete frequencies in the RF range. Zimmerman et al. recently developed a FD NIR high speed imaging system that utilized two wavelengths, two modulation frequencies, and an array of detectors; with this system, a measurement rate of 90 Hz was achieved, compared to the 97-Hz measurement rate achieved in the dDOS system with six wavelengths and a sweep of modulation frequencies.14

In the future, the portability of the system can be further improved by implementing miniaturized fiber-coupled laser diodes and miniaturized current driver modules that can be housed in a portable device enclosure. This will help to more easily measure patients in a clinical setting for a range of applications, such as chemotherapy monitoring or rapid hemodynamic analysis. The speed of the system also opens the possibility of monitoring rapid physiological changes at the cardiac rate, which may be useful in monitoring hemodynamics related to conditions, such as peripheral artery disease.

We have developed an inexpensive clinical FD-DOS platform based on DDSs and ADCs implemented on custom hardware-integrated PCBs. The system has high stability and accuracy in phantom measurements and is capable of detecting rapid physiological fluctuations above the cardiac rate using six wavelengths and a sweep of modulation frequencies. The dDOS system may improve access to patients while helping to open new areas of investigation.

The authors have no conflicts of interest.

We would like to thank Dr. Eric Hazen and the Boston University Electronics Design Facility for their ongoing support with hardware and firmware development. The authors gratefully acknowledge funding from the American Cancer Society (Grant No. RSG-14-014-01-CCE) and the Department of Defense (Award No. W81XWH-15-1-0070).

Tromberg  B. J.  et al., “Non-invasive measurements of breast tissue optical properties using frequency-domain photon migration,” Philos. Trans. R. Soc. Lond. B. Biol. Sci.. 352, , 661 –668 (1997). 0962-8436 CrossRef
Tromberg  B. J.  et al., “Optical property measurements in turbid media using frequency domain photon migration,” Proc. SPIE. 1525, , 52 –58 (1991). 0277-786X CrossRef
McBride  T. O.  et al., “A parallel-detection frequency-domain near-infrared tomography system for hemoglobin imaging of the breast in vivo,” Rev. Sci. Instrum.. 72, , 1817 –1824 (2001). 0034-6748 CrossRef
Culver  J. P.  et al., “Three-dimensional diffuse optical tomography in the parallel plane transmission geometry: evaluation of a hybrid frequency domain/continuous wave clinical system for breast imaging,” Med. Phys.. 30, , 235 –247 (2003). 0094-2405 CrossRef
Pham  T. H.  et al., “Broad bandwidth frequency domain instrument for quantitative tissue optical spectroscopy,” Rev. Sci. Instrum.. 71, , 2500 –2513 (2000). 0034-6748 CrossRef
Netz  U. J., , Beuthan  J., and Hielscher  A. H., “Multipixel system for gigahertz frequency-domain optical imaging of finger joints,” Rev. Sci. Instrum.. 79, , 034301  (2008). 0034-6748 CrossRef
Yang  Y. S.  et al., “Low-cost frequency-domain photon migration instrument for tissue spectroscopy, oximetry, and imaging,” Opt. Eng.. 36, , 1021 –1569 (1997).CrossRef
Yu  G.  et al., “Frequency-domain multiplexing system for in vivo diffuse light measurements of rapid cerebral hemodynamics,” Appl. Opt.. 42, , 2931 –2939 (2003). 0003-6935 CrossRef
Zhao  Y.  et al., “Portable, parallel 9-wavelength near-infrared spectral tomography (NIRST) system for efficient characterization of breast cancer within the clinical oncology infusion suite,” Biomed. Opt. Express. 7, , 2186 –2201 (2016). 2156-7085 CrossRef
Pogue  B.  et al., “Instrumentation and design of a frequency-domain diffuse optical tomography imager for breast cancer detection,” Opt. Express. 1, , 391 –403 (1997). 1094-4087 CrossRef
Madsen  S. J.  et al., “Keywords: frequency-domain photon migration, optical properties, diode laser, network analyzer,” Proc. SPIE. 2389, , 257 –263 (1995). 0277-786X CrossRef
No  K.-S.  et al., “Design and testing of a miniature broadband frequency domain photon migration instrument,” J. Biomed. Opt.. 13, , 050509  (2009). 1083-3668 CrossRef
O’Sullivan  T. D.  et al., “Diffuse optical imaging using spatially and temporally modulated light,” J. Biomed. Opt.. 17, , 071311  (2012). 1083-3668 CrossRef
Zimmermann  B. B.  et al., “Frequency domain near-infrared multiwavelength imager design using high-speed, direct analog-to-digital conversion,” J. Biomed. Opt.. 21, , 016010  (2016). 1083-3668 CrossRef
Jung  J., , Istfan  R., and Roblyer  D., “Note: a simple broad bandwidth undersampling frequency-domain digital diffuse optical spectroscopy system,” Rev. Sci. Instrum.. 85, , 076108  (2014). 0034-6748 CrossRef
Roblyer  D.  et al., “Feasibility of direct digital sampling for diffuse optical frequency domain spectroscopy in tissue,” Meas. Sci. Technol.. 24, , 045501  (2013). 0957-0233 CrossRef
Jung  J., , Istfan  R., and Roblyer  D., “Note: a simple broad bandwidth undersampling frequency-domain digital diffuse optical spectroscopy system,” Rev. Sci. Instrum.. 85, , 076108  (2014). 0034-6748 CrossRef
Zijlstra  W. G., and Buursma  A., “Spectrophotometry of hemoglobin: absorption spectra of bovine oxyhemoglobin, deoxyhemoglobin, carboxyhemoglobin, and methemoglobin,” Comput. Biochem. Physiol.. 118, , 743 –749 (1997).CrossRef
Bland  J. M., and Altman  D. G., “Statistical methods for assessing agreement between two methods of clinical measurement,” Int. J. Nurs. Stud.. 47, , 931 –936 (2010).CrossRef
Shah  N.  et al., “Spatial variations in optical and physiological properties,” J. Biomed. Opt.. 9, , 534 –540 (2004). 1083-3668 CrossRef
Weigel  U. M.  et al., “A new, modular frequency domain diffuse optical monitor in the digital domain,” in  Digital Holography and Three-Dimensional Imaging ,  Optical Society of America  (2012).
Yong  J., , Foo  A., and Lim  C. S., “Pulse transit time based on piezoelectric,” J. Clin. Monit. Comput.. 20, , 185 –192 (2006).CrossRef
Lee  W.  et al., “Systolic time intervals derived from electrocardiographic gated intra- renal artery Doppler waveform associated with left ventricular systolic function,” Sci. Rep.. 6, , 29293  (2016). 2045-2322 CrossRef
Vyas  M.  et al., “Augmentation index and central aortic stiffness in middle-aged to elderly individuals,” Am. J. Hypertens.. 20, , 642 –647 (2007).CrossRef
Roblyer  D.  et al., “Optical imaging of breast cancer oxyhemoglobin flare correlates with neoadjuvant chemotherapy response one day after starting treatment,” Proc. Natl. Acad. Sci. U. S. A.. 108, , 14626 –14631 (2011).CrossRef
© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation

Alyssa Torjesen ; Raeef Istfan and Darren Roblyer
"Ultrafast wavelength multiplexed broad bandwidth digital diffuse optical spectroscopy for in vivo extraction of tissue optical properties", J. Biomed. Opt. 22(3), 036009 (Mar 10, 2017). ; http://dx.doi.org/10.1117/1.JBO.22.3.036009


Figures

Graphic Jump Location
Fig. 1
F1 :

Schematic overview of dDOS system. GUI, graphical user interface; ADC, analog-to-digital converter; FPGA, field programmable gate array; RAM, rapid access memory; DDS, direct digital synthesizer; APD, avalanche photodiode.

Graphic Jump Location
Fig. 3
F3 :

Noise floor when collecting 4096 samples per frequency step in (b) a single wavelength measurement or (a) a multiplexed measurement with a 10 MHz offset between wavelengths.

Graphic Jump Location
Fig. 4
F4 :

Amplitude (in arbitrary units, A.U.) and phase measured with dDOS system at each wavelength. Red lines are calibrated data and black lines are model fits.

Graphic Jump Location
Fig. 5
F5 :

(a and b) Scatterplots and (c and d) Bland–Altman plots of optical properties obtained with the dDOS system compared to those obtained with the benchtop system. Top and bottom dashed lines represent 1.96 standard deviations from the mean difference; center dashed lines represent the mean difference.

Graphic Jump Location
Fig. 6
F6 :

2-h drift test showing dDOS system stability. Measurements were taken approximately every minute and show variation from value at the first time point (t=0). (a) Amplitude drift at 50 MHz and 850 nm and (b) phase drift at 50 MHz and 850 nm.

Graphic Jump Location
Fig. 7
F7 :

Forearm cuff occlusion test. Cuff was occluded during time interval between dashed lines. Measurements were taken approximately every 8 s and a three-point moving average is plotted for each chromophore.

Graphic Jump Location
Fig. 9
F9 :

ctTHb measured on the right breast of a healthy volunteer. A quadrant of the nipple is located in the lower left corner of the image.

Graphic Jump Location
Fig. 8
F8 :

High-speed dDOS measurements were taken on the index finger. An ECG signal was obtained simultaneously with the optical measurements and temporally aligned with optical property and chromophore values. From top to bottom, plots represent the ECG, μa, μs, ctTHb, ctO2Hb, and ctHHb. The leftmost vertical dashed line represents the approximate time of the peak of the ECG R-wave, the middle dashed line represents the approximate time of the start of the upstroke of the optical signal following the arterial transit time, and the rightmost dashed line represents the approximate timing of the dicrotic notch based on visual inspection.

Tables

Table Grahic Jump Location
Table 1System stability during a 2-h drift test with measurements taken approximately every minute.
Table Grahic Jump Location
Table 2dDOS performance specifications.

References

Tromberg  B. J.  et al., “Non-invasive measurements of breast tissue optical properties using frequency-domain photon migration,” Philos. Trans. R. Soc. Lond. B. Biol. Sci.. 352, , 661 –668 (1997). 0962-8436 CrossRef
Tromberg  B. J.  et al., “Optical property measurements in turbid media using frequency domain photon migration,” Proc. SPIE. 1525, , 52 –58 (1991). 0277-786X CrossRef
McBride  T. O.  et al., “A parallel-detection frequency-domain near-infrared tomography system for hemoglobin imaging of the breast in vivo,” Rev. Sci. Instrum.. 72, , 1817 –1824 (2001). 0034-6748 CrossRef
Culver  J. P.  et al., “Three-dimensional diffuse optical tomography in the parallel plane transmission geometry: evaluation of a hybrid frequency domain/continuous wave clinical system for breast imaging,” Med. Phys.. 30, , 235 –247 (2003). 0094-2405 CrossRef
Pham  T. H.  et al., “Broad bandwidth frequency domain instrument for quantitative tissue optical spectroscopy,” Rev. Sci. Instrum.. 71, , 2500 –2513 (2000). 0034-6748 CrossRef
Netz  U. J., , Beuthan  J., and Hielscher  A. H., “Multipixel system for gigahertz frequency-domain optical imaging of finger joints,” Rev. Sci. Instrum.. 79, , 034301  (2008). 0034-6748 CrossRef
Yang  Y. S.  et al., “Low-cost frequency-domain photon migration instrument for tissue spectroscopy, oximetry, and imaging,” Opt. Eng.. 36, , 1021 –1569 (1997).CrossRef
Yu  G.  et al., “Frequency-domain multiplexing system for in vivo diffuse light measurements of rapid cerebral hemodynamics,” Appl. Opt.. 42, , 2931 –2939 (2003). 0003-6935 CrossRef
Zhao  Y.  et al., “Portable, parallel 9-wavelength near-infrared spectral tomography (NIRST) system for efficient characterization of breast cancer within the clinical oncology infusion suite,” Biomed. Opt. Express. 7, , 2186 –2201 (2016). 2156-7085 CrossRef
Pogue  B.  et al., “Instrumentation and design of a frequency-domain diffuse optical tomography imager for breast cancer detection,” Opt. Express. 1, , 391 –403 (1997). 1094-4087 CrossRef
Madsen  S. J.  et al., “Keywords: frequency-domain photon migration, optical properties, diode laser, network analyzer,” Proc. SPIE. 2389, , 257 –263 (1995). 0277-786X CrossRef
No  K.-S.  et al., “Design and testing of a miniature broadband frequency domain photon migration instrument,” J. Biomed. Opt.. 13, , 050509  (2009). 1083-3668 CrossRef
O’Sullivan  T. D.  et al., “Diffuse optical imaging using spatially and temporally modulated light,” J. Biomed. Opt.. 17, , 071311  (2012). 1083-3668 CrossRef
Zimmermann  B. B.  et al., “Frequency domain near-infrared multiwavelength imager design using high-speed, direct analog-to-digital conversion,” J. Biomed. Opt.. 21, , 016010  (2016). 1083-3668 CrossRef
Jung  J., , Istfan  R., and Roblyer  D., “Note: a simple broad bandwidth undersampling frequency-domain digital diffuse optical spectroscopy system,” Rev. Sci. Instrum.. 85, , 076108  (2014). 0034-6748 CrossRef
Roblyer  D.  et al., “Feasibility of direct digital sampling for diffuse optical frequency domain spectroscopy in tissue,” Meas. Sci. Technol.. 24, , 045501  (2013). 0957-0233 CrossRef
Jung  J., , Istfan  R., and Roblyer  D., “Note: a simple broad bandwidth undersampling frequency-domain digital diffuse optical spectroscopy system,” Rev. Sci. Instrum.. 85, , 076108  (2014). 0034-6748 CrossRef
Zijlstra  W. G., and Buursma  A., “Spectrophotometry of hemoglobin: absorption spectra of bovine oxyhemoglobin, deoxyhemoglobin, carboxyhemoglobin, and methemoglobin,” Comput. Biochem. Physiol.. 118, , 743 –749 (1997).CrossRef
Bland  J. M., and Altman  D. G., “Statistical methods for assessing agreement between two methods of clinical measurement,” Int. J. Nurs. Stud.. 47, , 931 –936 (2010).CrossRef
Shah  N.  et al., “Spatial variations in optical and physiological properties,” J. Biomed. Opt.. 9, , 534 –540 (2004). 1083-3668 CrossRef
Weigel  U. M.  et al., “A new, modular frequency domain diffuse optical monitor in the digital domain,” in  Digital Holography and Three-Dimensional Imaging ,  Optical Society of America  (2012).
Yong  J., , Foo  A., and Lim  C. S., “Pulse transit time based on piezoelectric,” J. Clin. Monit. Comput.. 20, , 185 –192 (2006).CrossRef
Lee  W.  et al., “Systolic time intervals derived from electrocardiographic gated intra- renal artery Doppler waveform associated with left ventricular systolic function,” Sci. Rep.. 6, , 29293  (2016). 2045-2322 CrossRef
Vyas  M.  et al., “Augmentation index and central aortic stiffness in middle-aged to elderly individuals,” Am. J. Hypertens.. 20, , 642 –647 (2007).CrossRef
Roblyer  D.  et al., “Optical imaging of breast cancer oxyhemoglobin flare correlates with neoadjuvant chemotherapy response one day after starting treatment,” Proc. Natl. Acad. Sci. U. S. A.. 108, , 14626 –14631 (2011).CrossRef

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging & repositioning the boxes below.

Related Book Chapters

Topic Collections

Advertisement
  • Don't have an account?
  • Subscribe to the SPIE Digital Library
  • Create a FREE account to sign up for Digital Library content alerts and gain access to institutional subscriptions remotely.
Access This Article
Sign in or Create a personal account to Buy this article ($20 for members, $25 for non-members).
Access This Proceeding
Sign in or Create a personal account to Buy this article ($15 for members, $18 for non-members).
Access This Chapter

Access to SPIE eBooks is limited to subscribing institutions and is not available as part of a personal subscription. Print or electronic versions of individual SPIE books may be purchased via SPIE.org.