Spectral tomography with diffuse near-infrared light: inclusion of broadband frequency domain spectral data

Single-stage approach for estimating optical parameters in spectral quantitative photoacoustic tomography

In quantitative photoacoustic tomography, the optical parameters of a target, most importantly the concentrations of chromophores such as deoxygenated and oxygenated hemoglobin, are estimated from photoacoustic data measured on the boundary of the target. In this work, a numerical approximation of a forward model for spectral quantitative photoacoustic tomography is constructed by utilizing the diffusion approximation for light propagation, the acoustic wave equation for ultrasound propagation, and spectral models of optical absorption and scattering to describe the wavelength dependence of the optical parameters. The related inverse problem is approached in the framework of Bayesian inverse problems. Concentrations of four chromophores (deoxygenated and oxygenated hemoglobin, water, and lipid), two scattering parameters (reference scattering and scattering power), and the Grüneisen parameter are estimated in a single-stage from photoacoustic data. The methodology is evaluated using numerical simulations in different full-view and limited-view imaging settings. The results show that, utilizing spectral data and models, the spectral optical parameters and the Grüneisen parameter can be simultaneously estimated. Furthermore, the approach can also be utilized in limited-view imaging situations.

Multiwavelength photo-magnetic imaging algorithm improved for direct chromophore concentration recovery using spectral constraints

Multiwavelength photo-magnetic imaging (PMI) is a novel combination of diffuse optics and magnetic resonance imaging, to the best of our knowledge, that yields tissue chromophore concentration maps with high resolution and quantitative accuracy. Here, we present the first experimental results, to the best of our knowledge, obtained using a spectrally constrained PMI image reconstruction method, where chromophore concentration maps are directly recovered, unlike the conventional two-step approach that requires an intermediate step of reconstructing wavelength-dependent absorption coefficient maps. The imposition of the prior spectral information into the PMI inverse problem improves the reconstructed image quality and allows recovery of highly quantitative concentration maps, which are crucial for effective cancer detection and characterization. The obtained results demonstrate the higher performance of the direct reconstruction method. Indeed, the reconstructed concentration maps are not only of higher quality but also obtained approximately 2 times faster than the conventional method.

Resolving tissue chromophore concentration at MRI resolution using multi-wavelength photo-magnetic imaging

Photo-magnetic imaging (PMI) is an emerging optical imaging modality that showed great performance on providing absorption maps with high resolution and quantitative accuracy. As a multi-modality technology, PMI warms up the imaged object using a near infrared laser while temperature variation is measured using magnetic resonance imaging. By probing tissue at multiple wavelengths, concentration of the main tissue chromophores such as oxy- and deoxy-hemoglobin, lipid, and water are obtained then used to derive functional parameters such as total hemoglobin concentration and relative oxygen saturation. In this paper, we present a multi-wavelength PMI system that was custom-built to host five different laser wavelengths. After recovering the high-resolution absorption maps, a least-squared minimization process was used to resolve the different chromophore concentration. The performance of the system was experimentally tested on a phantom with two different dyes. Their concentrations were successfully assessed with high spatial resolution and average accuracy of nearly 80%.

Method for Quantitative Broadband Diffuse Optical Spectroscopy of Tumor-Like Inclusions

A hybrid reflectance-based diffuse optical imaging (DOI) technique combining discrete wavelength frequency-domain (FD) near-infrared spectroscopy (NIRS) with broadband continuous wave NIRS measurements was developed to quantify the broadband optical properties of deep tumor-like inclusions. This method was developed to more accurately measure the broadband optical properties of human tumors using a compact handheld imaging probe and without requiring a priori spectral constraints. We simulated the reconstruction of absorption and scattering spectra (650–1000 nm) of human breast tumors in a homogeneous background at depths of 0 to 10 mm. The hybrid DOI technique demonstrated enhanced performance in reconstruction of optical absorption with a mean accuracy over all 71 wavelengths of 8.39% versus 32.26% for a 10 mm deep tumor with the topographic DOI method. The new hybrid technique was also tested and validated on two heterogeneous tissue-simulating phantoms with inclusion depths of 2, 7, and 9 mm. The mean optical absorption accuracy over all wavelengths was similarly improved up to 5x for the hybrid DOI method versus topographic DOI for the deepest inclusions.

Frequency domain diffuse optical spectroscopy with a near-infrared tunable vertical cavity surface emitting laser

We present an approach for performing frequency domain diffuse optical spectroscopy (fd-DOS) utilizing a near-infrared tunable vertical cavity surface emitting laser (VCSEL) that enables high spectral resolution optical sensing in a miniature format. The tunable VCSEL, designed specifically for deep tissue imaging and sensing, utilizes an electrothermally tunable microelectromechanical systems topside mirror to tune the laser cavity resonance. At room temperature, the laser is tunable across 14nm from 769 to 782nm with single mode CW output and a peak output power of 1.3mW. We show that the tunable VCSEL is suitable for use in fd-DOS by measuring the optical properties of a tissue-simulating phantom over the tunable range. Optical properties were recovered within 0.0006mm^-1 (absorption) and 0.09mm^-1 (reduced scattering) compared to a broadband fd-DOS reference system. Our results indicate that tunable VCSELs may be an attractive choice to enable high spectral resolution optical sensing in a wearable format.

Automatic and robust calibration of optical detector arrays for biomedical diffuse optical spectroscopy

The design and testing of a new, fully automated, calibration approach is described. The process was used to calibrate an image-guided diffuse optical spectroscopy system with 16 photomultiplier tubes (PMTs), but can be extended to any large array of optical detectors and associated imaging geometry. The design goals were accomplished by developing a routine for robust automated calibration of the multi-detector array within 45 minutes. Our process was able to characterize individual detectors to a median norm of the residuals of 0.03 V for amplitude and 4.4 degrees in phase and achieved less than 5% variation between all the detectors at the 95% confidence interval for equivalent measurements. Repeatability of the calibrated data from the imaging system was found to be within 0.05 V for amplitude and 0.2 degrees for phase, and was used to evaluate tissue-simulating phantoms in two separate imaging geometries. Spectroscopic imaging of total hemoglobin concentration was recovered to within 5% of the true value in both cases. Future work will focus on streamlining the technology for use in a clinical setting with expectations of achieving accurate quantification of suspicious lesions in the breast.

Time-domain diffuse optical tomography processing by using the Mellin-Laplace transform

We investigate the use of the Mellin-Laplace transform for reconstructing optical parameters from time-resolved optical tomography in diffusive media. We present here its definition, its mathematical properties, and its sensitivity to variations of optical properties. The method was validated on two-dimensional reconstructions from simulation in the reflection geometry. We conclude that reconstructions based on the Mellin-Laplace transform are more robust to noise than the methods using first moments.

Diffuse Optical Monitoring of the Neoadjuvant Breast Cancer Therapy

Recent advances in the use of diffuse optical techniques for monitoring the hemodynamic, metabolic and physiological signatures of the neoadjuvant breast cancer therapy effectiveness is critically reviewed. An extensive discussion of the state-of-theart diffuse optical mammography is presented alongside a discussion of the current approaches to breast cancer therapies. Overall, the diffuse optics field is growing rapidly with a great deal of promise to fill an important niche in the current approaches to monitor, predict and personalize neoadjuvant breast cancer therapies.

Implementation of edge-preserving regularization for frequency-domain diffuse optical tomography

In this study, we first propose the use of edge-preserving regularization in optimizing an ill-conditioned problem in the reconstruction procedure for diffuse optical tomography to prevent unwanted edge smoothing, which usually degrades the attributes of images for distinguishing tumors from background tissues when using Tikhonov regularization. In the edge-preserving regularization method presented here, a potential function with edge-preserving properties is introduced as a regularized term in an objective function. With the minimization of this proposed objective function, an iterative method to solve this optimization problem is presented in which half-quadratic regularization is introduced to simplify the minimization task. Both numerical and experimental data are employed to justify the proposed technique. The reconstruction results indicate that edge-preserving regularization provides a superior performance over Tikhonov regularization.

The utility of direct-current as compared to frequency domain measurements in spectrally-constrained diffuse optical tomography toward cancer imaging

This work investigates, by means of analytical and simulation studies, the performance of spectrally-constrained image reconstruction in Continuous-Wave or Direct-Current (DC) and Frequency Domain (FD) near-infrared optical tomography. A recent analytic approach for estimating the accuracy of target recovery and the level of background artifact for optical tomography at single wavelength, based on the analysis of parametric reconstruction uncertainty level (PRUL), is extended to spectrally-constrained optical tomography. The analytical model is implemented to rank three sets of wavelengths that had been used as spectral prior in an independent experimental study. Subsequent simulation appraises the recovery of oxygenated hemoglobin (HbO), deoxygenated hemoglobin (Hb), water (H2O), scattering amplitude (A), and scattering power (b) using DC-only, DC-excluded FD, and DC-included FD, based on the three sets of wavelengths as the spectral prior. The simulation results support the analytic ranking of the performance of the three sets of spectral priors, and generally agree with the performance outcome of DC-only versus that of DC-excluded FD and DC-included FD. Specifically, this study indicate that: 1) the rank of overall quality of chromophore recovery is Hb, H2O, and HbO from the highest to lowest; and in the scattering part the A is always better recovered than b. This outcome does suggest that the DC-only information gives rise to unique solution to the image reconstruction routine under the given spectral prior. 2) DC-information is not-redundant in FD-reconstruction, as the artifact levels of DC-included FD reconstruction are always lower than those of DC-excluded FD. 3) The artifact level as represented by the noise-to-contrast-ratio is almost always the lowest in DC-only, leading to generally better resolution of multiple targets of identical contrasts over the background than in FD. However, the FD could outperform DC in the recovery of scattering properties including both A and b when the spectral prior is less optimal, implying the benefit of phase-information in scattering recovery in the context of spectrally-constrained optical tomography.

Remote positioning optical breast magnetic resonance coil for slice-selection during image-guided near-infrared spectroscopy of breast cancer

The design and testing of a pneumatic optical positioning interface produced with the goal of improving fiber positioning in magnetic resonance (MR)-guided diffuse spectral imaging of breast cancer is presented. The system was created for vertical positioning of optical fibers inside the MR bore during a patient exam to target suspicious lesions with MR scans for reference and collect multiple planes of optical data. The interface includes new fiber plates for mechanical and optical coupling to the breast, and was tested in phantoms and human imaging. Reconstructions with data taken in the new interface show acceptable linearity over different absorber concentrations (residual norm = 0.067), and exhibit good contrast recovery at different imaging planes, which is consistent with previous work. An example of human breast imaging through the new interface is shown and a discussion of how it compares to other patient interfaces for breast imaging is presented. Design goals of increasing the available degrees of freedom for fiber positioning while maintaining good patient-fiber contact and comfort were accomplished. This interface allows improved volumetric imaging with interactive and accurate slice selection to quantify targeted suspicious lesions.

Hyperspectral image reconstruction for diffuse optical tomography

We explore the development and performance of algorithms for hyperspectral diffuse optical tomography (DOT) for which data from hundreds of wavelengths are collected and used to determine the concentration distribution of chromophores in the medium under investigation. An efficient method is detailed for forming the images using iterative algorithms applied to a linearized Born approximation model assuming the scattering coefficient is spatially constant and known. The L-surface framework is employed to select optimal regularization parameters for the inverse problem. We report image reconstructions using 126 wavelengths with estimation error in simulations as low as 0.05 and mean square error of experimental data of 0.18 and 0.29 for ink and dye concentrations, respectively, an improvement over reconstructions using fewer specifically chosen wavelengths.

MR water quantitative priors improves the accuracy of optical breast imaging

Magnetic resonance (MR) guided optical breast imaging is a promising modality to improve the specificity of breast imaging, because it provides high-resolution quantitative maps of total hemoglobin, oxygen saturation, water content, and optical scattering. These properties have been shown to distinguish malignant from benign lesions. However, the optical detection hardware required for deep tissue imaging has poor spectral sensitivity which limits accurate water quantification; this reduces the accuracy of hemoglobin quantification. We present a methodology to improve optical quantification by utilizing the ability of Dixon MR imaging to quantitatively estimate water and fat; this technique effectively reduces optical crosstalk between water and oxyhemoglobin. The techniques described in this paper reduce hemoglobin quantification error by as much as 38%, as shown in a numerical phantom, and an experimental phantom. Error is reduced by as much 20% when imperfect MR water quantification is given. These techniques may also increase contrast between diseased and normal tissue, as shown in breast tissue in vivo. It is also shown that using these techniques may permit fewer wavelengths to be used with similar quantitative accuracy, enabling higher temporal resolution. In addition, it is shown that these techniques can improve the ability of MRI to quantify water in the presence of bias in the Dixon water/fat separation.

In vivo quantitative imaging of normal and cancerous breast tissue using broadband diffuse optical tomography

PURPOSE

A NIR tomography system that combines frequency domain (FD) and continuous wave (CW) measurements was used to image normal and malignant breast tissues.

METHODS

FD acquisitions were confined to wavelengths less than 850 nm because of detector limitations, whereas light from longer wavelengths (up to 948 nm) was measured in CW mode with CCD-coupled spectrometer detection. The two data sets were combined and processed in a single spectrally constrained reconstruction to map concentrations of hemoglobin, water, and lipid, as well as scattering parameters in the breast.

RESULTS

Chromophore concentrations were imaged in the breasts of nine asymptomatic volunteers to evaluate their intrasubject and intersubject variability. Normal subject data showed physiologically expected trends. Images from three cancer patients indicate that the added CW data is critical to recovering the expected increases in water and decreases in lipid content within malignancies. Contrasts of 1.5 to twofold in hemoglobin and water values were found in cancers.

CONCLUSIONS

In vivo breast imaging with instrumentation that combines FD and CW NIR data acquisition in a single spectral reconstruction produces more accurate hemoglobin, water, and lipid results relative to FD data alone.

Near-infrared tomography of breast cancer hemoglobin, water, lipid, and scattering using combined frequency domain and cw measurement

In this study, near-IR tomography was implemented in the wavelength range from 661 to 948 nm to characterize breast tumors in vivo. Frequency-domain measurements provide amplitude and phase transmitted at wavelengths below 850 nm, where photomultiplier tube detection is efficient. Continuous-wave detection at additional longer wavelengths (903, 912, and 948 nm) was collected using a CCD-based spectrometer. Phantom validation experiments showed improved accuracy in hemoglobin and water concentrations using this technique. Three women with malignant breast tumors were studied. The addition of cw data at longer wavelengths increased the recovered contrast of water in the tumor region relative to surrounding tissue and allowed quantification of lipid.

Diffuse optical imaging

Diffuse optical imaging is a medical imaging technique that is beginning to move from the laboratory to the hospital. It is a natural extension of near-infrared spectroscopy (NIRS), which is now used in certain niche applications clinically and particularly for physiological and psychological research. Optical imaging uses sophisticated image reconstruction techniques to generate images from multiple NIRS measurements. The two main clinical applications--functional brain imaging and imaging for breast cancer--are reviewed in some detail, followed by a discussion of other issues such as imaging small animals and multimodality imaging. We aim to review the state of the art of optical imaging.

Near-infrared spectral imaging of the female breast for quantitative oximetry in optical mammography

We present a hybrid continuous-wave, frequency-domain instrument for near-infrared spectral imaging of the female breast based on a tandem, planar scanning of one illumination optical fiber and one collection optical fiber configured in a transmission geometry. The spatial sampling rate of 25 points/cm(2) is increased to 400 points/cm(2) by postprocessing the data with a 2D cubic spline interpolation. We then apply a previously developed spatial second-derivative algorithm to an edge-corrected intensity image (N-image) to enhance the visibility and resolution of optical inhomogeneities in breast tissue such as blood vessels and tumors. The spectral data at each image pixel consist of 515-point spectra over the 650-900 nm wavelength range, thus featuring a spectral density of two data points per nanometer. We process the measured spectra with a paired-wavelength spectral analysis method to quantify the oxygen saturation of detected optical inhomogeneities, under the assumption that they feature a locally higher hemoglobin concentration. Our initial measurements on two healthy human subjects have generated high-resolution optical mammograms displaying a network of blood vessels with values of hemoglobin saturation typically falling within the 60%-95% range, which is physiologically reasonable. This approach to spectral imaging and oximetry of the breast has the potential to efficiently exploit the high intrinsic contrast provided by hemoglobin in breast tissue and to contribute a useful tool in the detection, diagnosis, and monitoring of breast pathologies.

Broadband frequency-domain near-infrared spectral tomography using a mode-locked Ti:sapphire laser

Frequency-domain near-infrared (NIR) diffuse spectral tomography with a mode-locked Ti:sapphire laser is presented, providing tunable multiwavelength quantitative spectroscopy with maximal power for thick tissue imaging. The system was developed to show that intrinsically high stability can be achieved with many wavelengths in the NIR range, using a mode-locked signal of 80 MHz with heterodyned lock-in detection. The effect of cumulative noise from multiple wavelengths of data on the reconstruction process was studied, and it was shown that inclusion of more wavelengths can reduce skew in the noise distribution. This normalization of the data variance then minimizes errors in estimation of chromophore concentrations. Simulations and tissue phantom experiments were used to quantify this improvement in image accuracy for recovery of tissue hemoglobin and oxygen saturation.

Wavelength band optimization in spectral near-infrared optical tomography improves accuracy while reducing data acquisition and computational burden

Multispectral near-infrared (NIR) tomographic imaging has the potential to provide information about molecules absorbing light in tissue, as well as subcellular structures scattering light, based on transmission measurements. However, the choice of possible wavelengths used is crucial for the accurate separation of these parameters, as well as for diminishing crosstalk between the contributing chromophores. While multispectral systems are often restricted by the wavelengths of laser diodes available, continuous-wave broadband systems exist that have the advantage of providing broadband NIR spectroscopy data, albeit without the benefit of the temporal data. In this work, the use of large spectral NIR datasets is analyzed, and an objective function to find optimal spectral ranges (windows) is examined. The optimally identified wavelength bands derived from this method are tested using both simulations and experimental data. It is found that the proposed method achieves images as qualitatively accurate as using the full spectrum, but improves crosstalk between parameters. Additionally, the judicious use of these spectral windows reduces the amount of data needed for full spectral tomographic imaging by 50%, therefore increasing computation time dramatically.