Simultaneous near-infrared diffusive light and ultrasound imaging

Optical Breast Imaging: A Review of Physical Principles, Technologies, and Clinical Applications

Optical imaging involves the propagation of light through tissue. Current optical breast imaging technologies, including diffuse optical spectroscopy, diffuse optical tomography, and photoacoustic imaging, capitalize on the selective absorption of light in the near-infrared spectrum by deoxygenated and oxygenated hemoglobin. They provide information on the morphological and functional characteristics of different tissues based on their varied interactions with light, including physiologic information on lesion vascular content and anatomic information on tissue vascularity. Fluorescent contrast agents, such as indocyanine green, are used to visualize specific tissues, molecules, or proteins depending on how and where the agent accumulates. In this review, we describe the physical principles, spectrum of technologies, and clinical applications of the most common optical systems currently being used or developed for breast imaging. Most notably, US co-registered photoacoustic imaging and US-guided diffuse optical tomography have demonstrated efficacy in differentiating benign from malignant breast masses, thereby improving the specificity of diagnostic imaging. Diffuse optical tomography and diffuse optical spectroscopy have shown promise in assessing treatment response to preoperative systemic therapy, and photoacoustic imaging and diffuse optical tomography may help predict tumor phenotype. Lastly, fluorescent imaging using indocyanine green dye performs comparably to radioisotope mapping of sentinel lymph nodes and appears to improve the outcomes of autologous tissue flap breast reconstruction.

Evaluation of a pipeline for simulation, reconstruction, and classification in ultrasound-aided diffuse optical tomography of breast tumors

SIGNIFICANCE

Diffuse optical tomography is an ill-posed problem. Combination with ultrasound can improve the results of diffuse optical tomography applied to the diagnosis of breast cancer and allow for classification of lesions.

AIM

To provide a simulation pipeline for the assessment of reconstruction and classification methods for diffuse optical tomography with concurrent ultrasound information.

APPROACH

A set of breast digital phantoms with benign and malignant lesions was simulated building on the software VICTRE. Acoustic and optical properties were assigned to the phantoms for the generation of B-mode images and optical data. A reconstruction algorithm based on a two-region nonlinear fitting and incorporating the ultrasound information was tested. Machine learning classification methods were applied to the reconstructed values to discriminate lesions into benign and malignant after reconstruction.

RESULTS

The approach allowed us to generate realistic US and optical data and to test a two-region reconstruction method for a large number of realistic simulations. When information is extracted from ultrasound images, at least 75% of lesions are correctly classified. With ideal two-region separation, the accuracy is higher than 80%.

CONCLUSIONS

A pipeline for the generation of realistic ultrasound and diffuse optics data was implemented. Machine learning methods applied to a optical reconstruction with a nonlinear optical model and morphological information permit to discriminate malignant lesions from benign ones.

Prospective assessment of adjunctive ultrasound-guided diffuse optical tomography in women undergoing breast biopsy: Impact on BI-RADS assessments

PURPOSE

To assess the impact of adjunctive ultrasound guided diffuse optical tomography (US-guided DOT) on BI-RADS assessment in women undergoing US-guided breast biopsy.

METHOD

This prospective study enrolled women referred for US-guided breast biopsy between 3/5/2019 and 3/19/2020. Participants underwent US-guided DOT immediately before biopsy. The US-guided DOT acquisition generated average maximum total hemoglobin (HbT) spatial maps and quantitative HbT values. Four radiologists blinded to histopathology assessed conventional imaging (CI) to assign a CI BI-RADS assessment and then integrated DOT information in assigning a CI&DOT BI-RADS assessment. HbT was compared between benign and malignant lesions using an ANOVA test and Tukey's test. Benign biopsies were tabulated, deeming BI-RADS ≥ 4A as positive. Reader agreement was assessed.

RESULTS

Among 61 included women (mean age 48 years), biopsy demonstrated 15 (24.6%) malignant and 46 (75.4%) benign lesions. Mean HbT was 55.3 ± 22.6 µM in benign lesions versus 85.4 ± 15.6 µM in cancers (p < .001). HbT threshold of 78.5 µM achieved sensitivity 80% (12/15) and specificity 89% (41/46) for malignancy. Across readers and patients, 197 pairs of CI BI-RADS and CI&DOT BI-RADS assessments were assigned. Adjunctive US-guided DOT achieved a net decrease in 23.5% (31/132) of suspicious (CI BI-RADS ≥ 4A) assessments of benign lesions (34 correct downgrades and 3 incorrect upgrades). 38.3% (31/81) of 4A assessments were appropriately downgraded. No cancer was downgraded to a non-actionable assessment. Interreader agreement analysis demonstrated kappa = 0.48-0.53 for CI BI-RADS and kappa = 0.28-0.44 for CI&DOT BI-RADS.

CONCLUSIONS

Integration of US-guided DOT information achieved a 23.5% reduction in suspicious BI-RADS assessments for benign lesions. Larger studies are warranted, with attention to improved reader agreement.

Deep learning-based method to accurately estimate breast tissue optical properties in the presence of the chest wall

SIGNIFICANCE

In general, image reconstruction methods used in diffuse optical tomography (DOT) are based on diffusion approximation, and they consider the breast tissue as a homogenous, semi-infinite medium. However, the semi-infinite medium assumption used in DOT reconstruction is not valid when the chest wall is underneath the breast tissue.

AIM

We aim to reduce the chest wall's effect on the estimated average optical properties of breast tissue and obtain accurate forward model for DOT reconstruction.

APPROACH

We propose a deep learning-based neural network approach where a convolution neural network (CNN) is trained to simultaneously obtain accurate optical property values for both the breast tissue and the chest wall.

RESULTS

The CNN model shows great promise in reducing errors in estimating the optical properties of the breast tissue in the presence of a shallow chest wall. For patient data, the CNN model predicted the breast tissue optical absorption coefficient, which was independent of chest wall depth.

CONCLUSIONS

Our proposed method can be readily used in DOT and diffuse spectroscopy measurements to improve the accuracy of estimated tissue optical properties.

Enhanced diffuse optical tomographic reconstruction using concurrent ultrasound information

Multimodal imaging is an active branch of research as it has the potential to improve common medical imaging techniques. Diffuse optical tomography (DOT) is an example of a low resolution, functional imaging modality that typically has very low resolution due to the ill-posedness of its underlying inverse problem. Combining the functional information of DOT with a high resolution structural imaging modality has been studied widely. In particular, the combination of DOT with ultrasound (US) could serve as a useful tool for clinicians for the formulation of accurate diagnosis of breast lesions. In this paper, we propose a novel method for US-guided DOT reconstruction using a portable time-domain measurement system. B-mode US imaging is used to retrieve morphological information on the probed tissues by means of a semi-automatical segmentation procedure based on active contour fitting. A two-dimensional to three-dimensional extrapolation procedure, based on the concept of distance transform, is then applied to generate a three-dimensional edge-weighting prior for the regularization of DOT. The reconstruction procedure has been tested on experimental data obtained on specifically designed dual-modality silicon phantoms. Results show a substantial quantification improvement upon the application of the implemented technique. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 2'.

Ultrasound segmentation-guided edge artifact reduction in diffuse optical tomography using connected component analysis

Ultrasound (US)-guided diffuse optical tomography (DOT) has demonstrated potential value for breast cancer diagnosis and treatment response assessment. However, in clinical use, the chest wall, poor probe-tissue contact, and tissue heterogeneity can all cause image artifacts. These image artifacts, appearing commonly as hot spots in the non-lesion regions (edge artifacts), can decrease the reconstruction accuracy and cause misinterpretation of lesion images. Here we introduce an iterative, connected component analysis-based image artifact reduction algorithm. A convolutional neural network (CNN) is used to segment co-registered US images to extract the lesion location and size to guide the artifact reduction. We demonstrate its performance using Monte Carlo simulations on VICTRE digital breast phantoms and breast patient images. In simulated tissue mismatch models, this algorithm successfully reduces edge artifacts without significantly changing the reconstructed target absorption coefficients. With clinical data it improves the optical contrast between malignant and benign groups, from 1.55 without artifact reduction to 1.91 with artifact reduction. The proposed algorithm has a broad range of applications in other modality-guided DOT imaging.

The LUCA device: a multi-modal platform combining diffuse optics and ultrasound imaging for thyroid cancer screening

We present the LUCA device, a multi-modal platform combining eight-wavelength near infrared time resolved spectroscopy, sixteen-channel diffuse correlation spectroscopy and a clinical ultrasound in a single device. By simultaneously measuring the tissue hemodynamics and performing ultrasound imaging, this platform aims to tackle the low specificity and sensitivity of the current thyroid cancer diagnosis techniques, improving the screening of thyroid nodules. Here, we show a detailed description of the device, components and modules. Furthermore, we show the device tests performed through well established protocols for phantom validation, and the performance assessment for in vivo. The characterization tests demonstrate that LUCA device is capable of performing high quality measurements, with a precision in determining in vivo tissue optical and dynamic properties of better than 3%, and a reproducibility of better than 10% after ultrasound-guided probe repositioning, even with low photon count-rates, making it suitable for a wide variety of clinical applications.

Effect and correction of optode coupling errors in breast imaging using diffuse optical tomography

In diffuse optical tomography (DOT) and spectroscopy (DOS) using handheld probes, tissue curvature can cause bad fiber-to-tissue contact. Understanding and minimizing image artifacts caused by these coupling errors would significantly improve DOT and DOS image quality. In this work, we utilized Monte Carlo simulations and experiments with gelatin-Intralipid phantoms to systematically study the influence of source or detector (optode) coupling errors. Optode coupling errors can increase the amplitude and decrease the phase of the measured diffuse reflectance, creating artifacts in the reconstructed absorption maps, such as hot spots on the edges. We propose an outlier removal algorithm that can correct these image artifacts, and we demonstrate its performance using simulations, phantom experiments, and breast patient data acquired with bad probe contact due to a dense or small breast. Further, we designed and implemented a new resistance-type thin-film force sensor array that provides real-time optode coupling feedback and guides the outlier removal to minimize optode coupling errors. Our approaches and study results have significant implications for reducing image artifacts arising from handheld probes, which are commonly used with mobile and wearable DOT and DOS devices.

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.

Reducing image artifact in diffuse optical tomography by iterative perturbation correction based on multiwavelength measurements

Ultrasound (US) guided diffuse optical tomography has demonstrated great potential for breast cancer diagnosis, treatment monitoring, and chemotherapy response prediction. Optical measurements of four different wavelengths are used to reconstruct unknown optical absorption maps, which are then used to calculate the hemoglobin concentration distribution of the US visible lesion. Reconstructed absorption maps are prone to image artifacts from outliers in measurement data from tissue heterogeneity, bad coupling between tissue and light guides, and motion by patient or operator. We propose an automated iterative perturbation correction algorithm to reduce image artifacts based on the structural similarity index (SSIM) of absorption maps of four optical wavelengths. The initial image is estimated from the truncated pseudoinverse solution. The SSIM was calculated for each wavelength to assess its similarity with other wavelengths. An absorption map is repeatedly reconstructed and projected back into measurement space to quantify projection error. Outlier measurements with highest projection errors are iteratively removed until all wavelength images are structurally similar with SSIM values greater than a threshold. Clinical data demonstrate statistically significant improvement in image artifact reduction.

An Automated Preprocessing Method for Diffuse Optical Tomography to Improve Breast Cancer Diagnosis

The ultrasound-guided diffuse optical tomography is a noninvasive imaging technique for breast cancer diagnosis and treatment monitoring. The technique uses a handheld probe capable of providing measurements of multiple wavelengths in a few seconds. These measurements are used to estimate optical absorptions of lesions and calculate the total hemoglobin concentration. Any measurement errors caused by low signal to noise ratio data and/or movements during data acquisition would reduce the accuracy of reconstructed total hemoglobin concentration. In this article, we introduce an automated preprocessing method that combines data collected from multiple sets of lesion measurements of 4 optical wavelengths to detect and correct outliers in the perturbation. Two new measures of correlation between each pair of wavelength measurements and a wavelength consistency index of all reconstructed absorption maps are introduced. For phantom and patients' data without evidence of measurement errors, the correlation coefficient between each pair of wavelength measurements was above 0.6. However, for patients with measurement errors, the correlation coefficient was much lower. After applying the correction method to 18 patients' data with measurement errors, the correlation has improved and the wavelength consistency index is in the same range as the cases without wavelength-dependent measurement errors. The results show an improvement in classification of malignant and benign lesions.

Compact ultrasound-guided diffuse optical tomography system for breast cancer imaging

Near-infrared diffuse optical tomography (DOT) has demonstrated a great potential as an adjunct modality for differentiation of malignant and benign breast lesions and for monitoring treatment response in patients with locally advanced breast cancers. The path toward commercialization of DOT techniques depends upon the improvement of robustness and user-friendliness of this technique in hardware and software. In this study, we introduce our recently developed ultrasound-guided DOT system, which has been improved in system compactness, robustness, and user-friendliness by custom-designed electronics, automated data preprocessing, and implementation of a new two-step reconstruction algorithm. The system performance has been tested with several sets of solid and blood phantoms and the results show accuracy in reconstructed absorption coefficients as well as blood oxygen saturation. A clinical example of a breast cancer patient, who was undergoing neoadjuvant chemotherapy, is given to demonstrate the system performance.

Spread spectrum time-resolved diffuse optical measurement system for enhanced sensitivity in detecting human brain activity

Diffuse optical spectroscopy (DOS) and imaging methods have been widely applied to noninvasive detection of brain activity. We have designed and implemented a low cost, portable, real-time one-channel time-resolved DOS system for neuroscience studies. Phantom experiments were carried out to test the performance of the system. We further conducted preliminary human experiments and demonstrated that enhanced sensitivity in detecting neural activity in the cortex could be achieved by the use of late arriving photons.

Improving breast cancer diagnosis by reducing chest wall effect in diffuse optical tomography

We have developed the ultrasound (US)-guided diffuse optical tomography technique to assist US diagnosis of breast cancer and to predict neoadjuvant chemotherapy response of patients with breast cancer. The technique was implemented using a hand-held hybrid probe consisting of a coregistered US transducer and optical source and detector fibers which couple the light illumination from laser diodes and photon detection to the photomultiplier tube detectors. With the US guidance, diffused light measurements were made at the breast lesion site and the normal contralateral reference site which was used to estimate the background tissue optical properties for imaging reconstruction. However, background optical properties were affected by the chest wall underneath the breast tissue. We have analyzed data from 297 female patients, and results have shown statistically significant correlation between the fitted optical properties ( ? a and ? s ? ) and the chest wall depth. After subtracting the background ? a at each wavelength, the difference of computed total hemoglobin (tHb) between malignant and benign lesion groups has improved. For early stage malignant lesions, the area-under-the-receiver operator characteristic curve (AUC) has improved from 88.5% to 91.5%. For all malignant lesions, the AUC has improved from 85.3% to 88.1%. Statistical test has revealed the significant difference of the AUC improvements after subtracting background tHb values.

Ultrasound-Guided Diffuse Optical Tomography for Predicting and Monitoring Neoadjuvant Chemotherapy of Breast Cancers: Recent Progress

In this manuscript, we review the current progress of utilizing ultrasound-guided diffuse optical tomography (US-guided DOT) for predicting and monitoring neoadjuvant chemotherapy (NAC) outcomes of breast cancer patients. We also report the recent advance on optical tomography systems toward portable and robust clinical use at multiple clinical sites. The first patient who has been closely monitored before NAC, at day 2, day 8, end of first three cycles of NAC, and before surgery is given as an example to demonstrate the potential of US-guided DOT technique.

Ultrasound-Guided Diffuse Optical Tomography for Predicting and Monitoring Neoadjuvant Chemotherapy of Breast Cancers: Recent Progress

In this manuscript, we review the current progress of utilizing ultrasound-guided diffuse optical tomography (US-guided DOT) for predicting and monitoring neoadjuvant chemotherapy (NAC) outcomes of breast cancer patients. We also report the recent advance on optical tomography systems toward portable and robust clinical use at multiple clinical sites. The first patient who has been closely monitored before NAC, at day 2, day 8, end of first three cycles of NAC, and before surgery is given as an example to demonstrate the potential of US-guided DOT technique.

Calibration and optimization of 3D digital breast tomosynthesis guided near infrared spectral tomography

Calibration of a three-dimensional multimodal digital breast tomosynthesis (DBT) x-ray and non-fiber based near infrared spectral tomography (NIRST) system is challenging but essential for clinical studies. Phantom imaging results yielded linear contrast recovery of total hemoglobin (HbT) concentration for cylindrical inclusions of 15 mm, 10 mm and 7 mm with a 3.5% decrease in the HbT estimate for each 1 cm increase in inclusion depth. A clinical exam of a patient's breast containing both benign and malignant lesions was successfully imaged, with greater HbT was found in the malignancy relative to the benign abnormality and fibroglandular regions (11 μM vs. 9.5 μM). Tools developed improved imaging system characterization and optimization of signal quality, which will ultimately improve patient selection and subsequent clinical trial results.

Light-emitting diode-based multiwavelength diffuse optical tomography system guided by ultrasound

Laser diodes are widely used in diffuse optical tomography (DOT) systems but are typically expensive and fragile, while light-emitting diodes (LEDs) are cheaper and are also available in the near-infrared (NIR) range with adequate output power for imaging deeply seated targets. In this study, we introduce a new low-cost DOT system using LEDs of four wavelengths in the NIR spectrum as light sources. The LEDs were modulated at 20 kHz to avoid ambient light. The LEDs were distributed on a hand-held probe and a printed circuit board was mounted at the back of the probe to separately provide switching and driving current to each LED. Ten optical fibers were used to couple the reflected light to 10 parallel photomultiplier tube detectors. A commercial ultrasound system provided simultaneous images of target location and size to guide the image reconstruction. A frequency-domain (FD) laser-diode-based system with ultrasound guidance was also used to compare the results obtained from those of the LED-based system. Results of absorbers embedded in intralipid and inhomogeneous tissue phantoms have demonstrated that the LED-based system provides a comparable quantification accuracy of targets to the FD system and has the potential to image deep targets such as breast lesions.

Two-step reconstruction method using global optimization and conjugate gradient for ultrasound-guided diffuse optical tomography

Ultrasound-guided diffuse optical tomography (DOT) is a promising method for characterizing malignant and benign lesions in the female breast. We introduce a new two-step algorithm for DOT inversion in which the optical parameters are estimated with the global optimization method, genetic algorithm. The estimation result is applied as an initial guess to the conjugate gradient (CG) optimization method to obtain the absorption and scattering distributions simultaneously. Simulations and phantom experiments have shown that the maximum absorption and reduced scattering coefficients are reconstructed with less than 10% and 25% errors, respectively. This is in contrast with the CG method alone, which generates about 20% error for the absorption coefficient and does not accurately recover the scattering distribution. A new measure of scattering contrast has been introduced to characterize benign and malignant breast lesions. The results of 16 clinical cases reconstructed with the two-step method demonstrates that, on average, the absorption coefficient and scattering contrast of malignant lesions are about 1.8 and 3.32 times higher than the benign cases, respectively.

A digital x-ray tomosynthesis coupled near infrared spectral tomography system for dual-modality breast imaging

A Near Infrared Spectral Tomography (NIRST) system has been developed and integrated into a commercial Digital Breast Tomosynthesis (DBT) scanner to allow structural and functional imaging of breast in vivo. The NIRST instrument uses an 8-wavelength continuous wave (CW) laser-based scanning source assembly and a 75-element silicon photodiode solid-state detector panel to produce dense spectral and spatial projection data from which spectrally constrained 3D tomographic images of tissue chromophores are produced. Integration of the optical imaging system into the DBT scanner allows direct co-registration of the optical and DBT images, while also facilitating the synergistic use of x-ray contrast as anatomical priors in optical image reconstruction. Currently, the total scan time for a combined NIRST-DBT exam is ~50s with data collection from 8 wavelengths in the optical scan requiring ~42s to complete. The system was tested in breast simulating phantoms constructed using intralipid and blood in an agarose matrix with a 3 cm x 2 cm cylindrical inclusion at 1 cm depth from the surface. Diffuse image reconstruction of total hemoglobin (HbT) concentration resulted in accurate recovery of the lateral size and position of the inclusion to within 6% and 8%, respectively. Use of DBT structural priors in the NIRST reconstruction process improved the quantitative accuracy of the HbT recovery, and led to linear changes in imaged versus actual contrast, underscoring the advantages of dual-modality optical imaging approaches. The quantitative accuracy of the system can be further improved with independent measurements of scattering properties through integration of frequency or time domain data.

Mono- and multimodal registration of optical breast images

Optical breast imaging offers the possibility of noninvasive, low cost, and high sensitivity imaging of breast cancers. Poor spatial resolution and a lack of anatomical landmarks in optical images of the breast make interpretation difficult and motivate registration and fusion of these data with subsequent optical images and other breast imaging modalities. Methods used for registration and fusion of optical breast images are reviewed. Imaging concerns relevant to the registration problem are first highlighted, followed by a focus on both monomodal and multimodal registration of optical breast imaging. Where relevant, methods pertaining to other imaging modalities or imaged anatomies are presented. The multimodal registration discussion concerns digital x-ray mammography, ultrasound, magnetic resonance imaging, and positron emission tomography.

FPGA-based reconfigurable processor for ultrafast interlaced ultrasound and photoacoustic imaging

In this paper, we report, to the best of our knowledge, a unique field-programmable gate array (FPGA)-based reconfigurable processor for real-time interlaced co-registered ultrasound and photoacoustic imaging and its application in imaging tumor dynamic response. The FPGA is used to control, acquire, store, delay-and-sum, and transfer the data for real-time co-registered imaging. The FPGA controls the ultrasound transmission and ultrasound and photoacoustic data acquisition process of a customized 16-channel module that contains all of the necessary analog and digital circuits. The 16-channel module is one of multiple modules plugged into a motherboard; their beamformed outputs are made available for a digital signal processor (DSP) to access using an external memory interface (EMIF). The FPGA performs a key role through ultrafast reconfiguration and adaptation of its structure to allow real-time switching between the two imaging modes, including transmission control, laser synchronization, internal memory structure, beamforming, and EMIF structure and memory size. It performs another role by parallel accessing of internal memories and multi-thread processing to reduce the transfer of data and the processing load on the DSP. Furthermore, because the laser will be pulsing even during ultrasound pulse-echo acquisition, the FPGA ensures that the laser pulses are far enough from the pulse-echo acquisitions by appropriate time-division multiplexing (TDM). A co-registered ultrasound and photoacoustic imaging system consisting of four FPGA modules (64-channels) is constructed, and its performance is demonstrated using phantom targets and in vivo mouse tumor models.

Investigation of a diffuse optical measurements-assisted quantitative photoacoustic tomographic method in reflection geometry

Photoacoustic tomography provides the distribution of absorbed optical energy density, which is the product of optical absorption coefficient and optical fluence distribution. We report the experimental investigation of a novel fitting procedure that quantitatively determines the optical absorption coefficient of chromophores. The experimental setup consisted of a hybrid system of a 64-channel photoacoustic imaging system with a frequency-domain diffused optical measurement system. The fitting procedure included a complete photoacoustic forward model and an analytical solution of a target chromophore using the diffusion approximation. The fitting procedure combines the information from the photoacoustic image and the background information from the diffuse optical measurements to minimize the photoacoustic measurements and forward model data and recover the target absorption coefficient quantitatively. 1-cm-cube phantom absorbers of high and low contrasts were imaged at depths of up to 3.0 cm. The fitted absorption coefficient results were at least 80% of their true values. The sensitivities of this fitting procedure to target location, target radius, and background optical properties were also investigated. We found that this fitting procedure was most sensitive to the accurate determination of the target radius and depth. Blood sample in a thin tube of radius 0.58 mm, simulating a blood vessel, was also studied. The photoacoustic images and fitted absorption coefficients are presented. These results demonstrate the clinical potential of this fitting procedure to quantitatively characterize small lesions in breast imaging.

Digital optical tomography system for dynamic breast imaging

Diffuse optical tomography has shown promising results as a tool for breast cancer screening and monitoring response to chemotherapy. Dynamic imaging of the transient response of the breast to an external stimulus, such as pressure or a respiratory maneuver, can provide additional information that can be used to detect tumors. We present a new digital continuous-wave optical tomography system designed to simultaneously image both breasts at fast frame rates and with a large number of sources and detectors. The system uses a master-slave digital signal processor-based detection architecture to achieve a dynamic range of 160 dB and a frame rate of 1.7 Hz with 32 sources, 64 detectors, and 4 wavelengths per breast. Included is a preliminary study of one healthy patient and two breast cancer patients showing the ability to identify an invasive carcinoma based on the hemodynamic response to a breath hold.

Artifact reduction method in ultrasound-guided diffuse optical tomography using exogenous contrast agents

In diffuse optical tomography (DOT), a typical perturbation approach requires two sets of measurements obtained at the lesion breast (lesion or target site) and a contra-lateral location of the normal breast (reference site) for image reconstruction. For patients who have a small amount of breast tissue, the chest-wall underneath the breast tissue at both sites affects the imaging results. In this group of patients, the perturbation, which is the difference between measurements obtained at the lesion and reference sites, may include the information of background mismatch which can generate artifacts or affect the reconstructed quantitative absorption coefficient of the lesion. Also, for patients who have a single breast due to prior surgery, the contra-lateral reference is not available. To improve the DOT performance or overcome its limitation, we introduced a new method based on an exogenous contrast agent and demonstrate its performance using animal models. Co-registered ultrasound was used to guide the lesion localization. The results have shown that artifacts caused by background mismatch can be reduced significantly by using this new method.

A multichannel time-domain scanning fluorescence mammograph: performance assessment and first in vivo results

We present a scanning time-domain fluorescence mammograph capable to image the distribution of a fluorescent contrast agent within a female breast, slightly compressed between two parallel glass plates, with high sensitivity. Fluorescence of the contrast agent is excited using a near infrared picosecond diode laser module. Four additional picosecond diode lasers with emission wavelengths between 660 and 1066 nm allow to measure the intrinsic optical properties of the breast tissue. By synchronously moving a source fiber and seven detection fiber bundles across the breast, distributions of times of flight of photons are recorded simultaneously for selected source-detector combinations in transmission and reflection geometry either at the fluorescence wavelength or at the selected laser wavelengths. To evaluate the performance of the mammograph, we used breastlike rectangular phantoms comprising fluorescent and absorbing objects using the fluorescent dye Omocyanine as contrast agent excited at 735 nm. We compare two-dimensional imaging of the phantom based on transmission and reflection data. Furthermore, we developed an improved tomosynthesis algorithm which permits three-dimensional reconstruction of fluorescence and absorption properties of lesions with good spatial resolution. For illustration, we present fluorescence mammograms of one patient recorded 30 min after administration of the contrast agent indocyanine green showing the carcinoma at high contrast originating from fluorescence of the extravasated dye, excited at 780 nm.

Optical mammography: Diffuse optical imaging of breast cancer

Existing imaging modalities for breast cancer screening, diagnosis and therapy monitoring, namely X-ray mammography and magnetic resonance imaging, have been proven to have limitations. Diffuse optical imaging is a set of non-invasive imaging modalities that use near-infrared light, which can be an alternative, if not replacement, to those existing modalities. This review covers the background knowledge, recent clinical outcome, and future outlook of this newly emerging medical imaging modality.

Diffuse Optics for Tissue Monitoring and Tomography

This review describes the diffusion model for light transport in tissues and the medical applications of diffuse light. Diffuse optics is particularly useful for measurement of tissue hemodynamics, wherein quantitative assessment of oxy- and deoxy-hemoglobin concentrations and blood flow are desired. The theoretical basis for near-infrared or diffuse optical spectroscopy (NIRS or DOS, respectively) is developed, and the basic elements of diffuse optical tomography (DOT) are outlined. We also discuss diffuse correlation spectroscopy (DCS), a technique whereby temporal correlation functions of diffusing light are transported through tissue and are used to measure blood flow. Essential instrumentation is described, and representative brain and breast functional imaging and monitoring results illustrate the workings of these new tissue diagnostics.

Automated coregistered imaging using a hand-held probe-based optical imager

Near-infrared optical imaging holds a promise as a noninvasive technology toward cancer diagnostics and other tissue imaging applications. In recent years, hand-held based imagers are of great interest toward the clinical translation of the technology. However hand-held imagers developed to date are typically designed to obtain surface images and not tomography information due to lack of coregistration facilities. Herein, a recently developed hand-held probe-based optical imager in our Optical Imaging Laboratory has been implemented with novel coregistration facilities toward real-time and tomographic imaging of tissue phantoms. Continuous-wave fluorescence-enhanced optical imaging studies were performed using an intensified charge coupled device camera based imaging system in order to demonstrate the feasibility of automated coregistered imaging of flat phantom surfaces, using a flexible probe that can also contour to curvatures. Three-dimensional fluorescence tomographic reconstructions were also demonstrated using coregistered frequency-domain measurements obtained using the hand-held based optical imager. It was also observed from preliminary studies on cubical phantoms that multiple coregistered scans differentiated deeper targets (approximately 3 cm) from artifacts that were not feasible from a single coregistered scan, demonstrating the possibility of improved target depth detectability in the future.

Fluorescence imaging of vascular endothelial growth factor in tumors for mice embedded in a turbid medium

We demonstrate the feasibility of fluorescence imaging of deeply seated tumors using mice injected with an angiogenesis tracer, a vascular endothelial growth factor conjugated with the infrared dye cyanine 7 (VEGF/Cy7). Our optical-only imaging reconstruction method separately estimates the target depth, and then applies this information to reconstruct functional information such as fluorophore concentration. Fluorescence targets with concentrations as low as sub-25 nM are well reconstructed at depths up to 2 cm in both homogeneous and heterogeneous media with this technique.

Effect of the chest wall on breast lesion reconstruction

The chest wall underneath the breast tissue affects near-infrared (NIR) diffusive waves measured with reflection geometry. With the assistance of a co-registered ultrasound, the depth and the tilting angle of the chest wall can be determined and are used to model the breast as a two-layer medium. Finite element method (FEM) is suitable for modeling complex boundary conditions and is adapted to model the breast tissue and chest wall. Four parameters of bulk absorption and reduced scattering coefficients of these two layers are estimated and used for imaging reconstruction. Using a two-layer model, we have systematically investigated the effect of the chest wall on breast lesion reconstruction. Results have shown that chest-wall depth, titling angle, and difference between optical properties of two layers of lesion and reference sites affect the lesion reconstruction differently. Our analysis will be valuable and informative to researchers who are using reflectance geometry for breast imaging. The analysis can also provide guidelines for imaging operators to minimize image artifacts and to produce the best reconstruction results.

Ultrasound-modulated fluorescence based on a fluorophore-quencher-labeled microbubble system

Ultrasound-modulated fluorescence from a fluorophore-quencher-labeled microbubble system driven by a single ultrasound pulse was theoretically quantified by solving a modified Herring equation (for bubble oscillation), a two-energy-level rate equation (for fluorophore excitation), and a diffusion equation (for light propagation in tissue). The efficiency of quenching caused by fluorescence resonance energy transfer (FRET) between the fluorophore and the quencher was modulated when the microbubble oscillates in size driven by the ultrasound pulse. Both intensity- and lifetime-based imaging methods are discussed in three different illumination modes of the excitation light: continuous wave (DC), frequency domain (FD), and time domain (TD). Results show that microbubble expansion opens a time period during which the quenching efficiency is dramatically reduced so that the emitted fluorescence strength and fluorophore lifetime are significantly increased. The modulation efficiency may even reach 100%. In addition, an important finding in this study is that in TD illumination mode, the modulated fluorescence photons may be temporally separated from the unmodulated photons, which makes the modulation efficiency limited only by thermal noise of the measurement system.

Noninvasive monitoring of breast cancer during neoadjuvant chemotherapy using optical tomography with ultrasound localization

The purposes of this study were 1) to investigate the feasibility of using optical tomography in the near-infrared (NIR) spectrum combined with ultrasound (US) localization (NIR/US) in monitoring tumor vascular changes and assessing tumor pathological response during chemotherapy and 2) to compare the accuracy of NIR/US with magnetic resonance imaging (MRI) in predicting residual cancer after neoadjuvant chemotherapy. Eleven female patients were studied during treatments with a combined imager consisting of a commercially available US system coupled to an NIR imager. Contrast-enhanced MRI was performed before treatment and surgery. Tumor vascular content was assessed based on total hemoglobin concentration and volume obtained from NIR data. A percentage blood volume index (%BVI) was calculated as the percentage ratio of the product of total hemoglobin concentration and volume normalized to pretreatment values. At treatment completion, pathologic assessment revealed three response groups: complete or near-complete responders (A), partial responders (B), and nonresponders (C). The mean %BVIs of groups A, B, and C at the treatment completion were 29.1 +/- 6.9%, 46.3 +/- 3.7%, and 86.8 +/- 30.1%, respectively (differences statistically significant, P < .04). At the end of cycle 2, the %BVI of group A was noticeably lower than that of the partial (P = .091) and nonresponder groups (P = .075). Both NIR/US and MRI were equally effective in distinguishing different response groups in this pilot study. Our initial findings indicate that NIR/US using %BVI can be used during chemotherapy to repeatedly monitor tumor vascular changes. NIR/US also may evaluate pathologic response during treatment allowing for tailoring therapies to response.

Progress in multimodality imaging: truly simultaneous ultrasound and magnetic resonance imaging

Multimodality medical imaging takes advantage of the strengths of different imaging modalities to provide a more complete picture of the anatomy under investigation. Many complementary modalities have been combined to form such systems and some are gaining use clinically. One combination that has not been developed, in large part due to technical difficulties, is a combined magnetic resonance (MR) and ultrasound (US) imaging system. Such a system offers the potential to combine the strengths of these modalities in a wide range of diagnostic and therapeutic applications. The goal of this study was to evaluate the feasibility of performing simultaneous multimodality US and MR imaging. An US imaging system capable of operation in a clinical MR imager was developed, and methods to perform simultaneous imaging were investigated. Simultaneous imaging was feasible without any mutual interference by either filtering the transmitted and received US signal, or by synchronizing data acquisition between the two imaging systems. Spatial registration between the two modalities was achieved by using a reference phantom with implanted glass beads in orthogonal planes. Excellent agreement was observed between spatial measurements of an object made with both modalities, and the feasibility of using this system in vivo was demonstrated in a rabbit model. Simultaneous US and MR imaging is achievable, and can provide complementary information about an object under investigation. This demonstration of technical feasibility and the development of a prototype system open up the potential to investigate the promising clinical applications of this combined technology.

Optical tomography with ultrasound localization for breast cancer diagnosis and treatment monitoring

Optical tomography with ultrasound (US) localization uses coregistered ultrasound images to guide optical imaging reconstruction. To simultaneously acquire US images and optical measurements, the authors used a hand-held probe consisting of a commercial US transducer and near-infrared optical imaging sensors of multiple wavelengths. A novel image scheme was used to map the ultrasound-visible lesions for optical imaging reconstruction. As a result, the problem of intense light scattering caused by breast tissue was overcome and reliable tumor hemoglobin concentration and blood oxygen saturation distributions from a group of patients were obtained. These functional parameters are valuable for aiding US diagnosis and for assessing chemotherapy response.

Controlled Monte Carlo method for light propagation in tissue of semi-infinite geometry

The controlled Monte Carlo method is generalized to model photon migration in turbid media of arbitrary geometries. Its implementation for the reflection geometry is exemplified in this paper. The most probable diffuse direction of light is used as the local attractive vector that serves as the basis of biased sampling of scattering angles. Consequently, path-length resolved photon trajectories can be generated with a significantly improved efficiency. We report a more than 29 times reduction in simulation time for early arriving photons in a typical configuration.

Optimal probing of optical contrast of breast lesions of different size located at different depths by US localization

We report a frequency domain optical tomography system utilizing three RF modulation frequencies, which are optimized for probing breast lesions of different size located at different depths. A real-time co-registered ultrasound scanner is used to provide on-site estimation of lesion size and location. Based on the lesion information, an optimal light modulation frequency can be selected, which may yield more accurate estimates of lesion angiogenesis and hypoxia. Phantom experiments have demonstrated that a high modulation frequency, such as 350Mhz, is preferable for probing small lesions closer to the surface while a low modulation frequency, such as 50Mhz, is desirable for imaging deeper and larger lesions. A clinical example of a large invasive carcinoma is presented to demonstrate the application of this novel technique.

Separately reconstructing the structural and functional parameters of a fluorescent inclusion embedded in a turbid medium

We report a novel imaging technique for fluorescence diffuse optical tomography (FDOT). Unlike conventional FDOT, this technique separates the imaging procedure into two steps to respectively reconstruct the structural information (such as the center position and the radius), and the functional information (such as the fluorophore concentration and/or lifetime) of a fluorescing target embedded in a turbid medium. The structural parameters of the target were estimated from the amplitude ratio and phase difference of fluorescence signals received at different detectors, because the amplitude ratio and phase difference were found independent of, or weakly related to, the functional parameters. Based on the estimated structural parameters, a dual-zone mesh technique was utilized to reconstruct the fluorophore concentration. Results of simulations and phantom experiments showed that the structural parameters could be accurately recovered, without knowing the functional information, and that the reconstruction accuracy of the functional parameter was greater than 80%.

Analytical solution for light propagation in a two-layer tissue structure with a tilted interface for breast imaging

Reflectance measurement of breast tissue is influenced by the underlying chest wall, which is often tilted as seen by the detection probe. We develop an analytical solution of light propagation in a two-layer tissue structure with tilted interface and refractive index difference between the layers. We validate the analytical solution with Monte Carlo simulations and phantom experiments, and a good agreement is seen. The influence of varying the tilting angle of the interface on the reflectance is discussed for two types of layered structures. Further, we apply the developed analytical solution to obtain the optical properties of breast tissue and chest wall from clinical data. Inverse calculation using the developed solution applied to the data obtained from Monte Carlo simulations shows that the optical properties of both layers are obtained with higher accuracy as compared to using a simple two-layer model ignoring the interface tilt. This is expected to improve the accuracy in estimating the optical properties of breast tissue, thus enhancing the accuracy of optical tomography of breast tumors.

Benign versus malignant breast masses: optical differentiation with US-guided optical imaging reconstruction

PURPOSE

To investigate prospectively the feasibility of using optical tomography with ultrasonographic (US) localization to differentiate malignant from benign breast masses and to compare optical tomography with color Doppler US.

MATERIALS AND METHODS

The study was approved by the local internal review board committee and by the Human Subjects Research Review Board of Army Medical Research and Materiel Command. Signed informed consent was obtained, and the study was HIPAA compliant. Between May 2003 and March 2004, 65 consecutive women (mean age, 51 years; age range, 24-80 years) with 81 breast lesions underwent US-guided biopsy and were scanned with a combined imager. The hand-held probe, which consisted of a centrally located US transducer surrounded by near-infrared sensors, was used to simultaneously acquire coregistered US images and optical data. The lesion location obtained at US was used to guide optical imaging reconstruction. Light absorption was measured at two wavelengths. From these measurements, tumor angiogenesis was assessed on the basis of calculated total hemoglobin concentration. A Student t distribution was used to calculate the statistical significance of mean maximum and mean average hemoglobin concentrations obtained in malignant and benign lesion groups, and P < .001 was considered to indicate a statistically significant difference.

RESULTS

Biopsy results revealed eight early stage invasive carcinomas (malignant group) and 73 benign lesions (benign group). The mean maximum and mean average hemoglobin concentrations in the malignant group were 122 micromol/L +/- 26.8 (+/- standard deviation) and 88 micromol/L +/- 24.5, respectively. The mean maximum and mean average hemoglobin concentrations in the benign group were 55 micromol/L +/- 24.8 and 38 micromol/L +/- 17.4, respectively. Both the maximum and average total hemoglobin concentrations were significantly higher in the malignant group compared with the benign group (P < .001). When a maximum hemoglobin concentration of 95 micromol/L was used as the threshold value, the sensitivity, specificity, positive predictive value, and negative predictive value of optical tomography were 100%, 96%, 73%, and 100%, respectively, and the sensitivity, specificity, positive predictive value, and negative predictive value of color Doppler US were 63%, 69%, 19%, and 94%, respectively.

CONCLUSION

Findings indicate that optical tomography with US localization is feasible for differentiating benign and early stage malignant breast lesions.

Optical tomography with ultrasound localization: initial clinical results and technical challenges

Optical tomography with ultrasound localization utilizes co-registered ultrasound lesion structure information to guide optical imaging reconstruction. A hand-held probe consisting of a commercial ultrasound transducer and near infrared optical imaging sensors was used to simultaneously acquire ultrasound images and optical measurements. A dual-mesh scheme was used to map the ultrasound-visible lesions to finer-grid lesion regions and coarser-grid background regions for optical imaging reconstruction. As a result, optical imaging reconstruction was well-conditioned for inverse mapping of lesion hemoglobin concentration and blood oxygen saturation. Initial clinical results have shown that early stage invasive cancers may be distinguished by a two-fold greater total hemoglobin concentration compared with fibroadenomas and other benign lesions. Initial results of advanced cancers have shown that the hemoglobin distribution is highly distorted and heterogeneous and the distorted distributions correlate with histological microvessel density counts and could be used to assess chemotherapy response.

Utilizing optical tomography with ultrasound localization to image heterogeneous hemoglobin distribution in large breast cancers

PURPOSE

Angiogenesis in advanced breast cancers is highly distorted and heterogeneous. Noninvasive imaging that can monitor angiogenesis may be invaluable initially for diagnosis and then for assessing tumor response to treatment. By combining ultrasound (US) and near-infrared (NIR) optical imaging, a reliable new technique has emerged for localizing and characterizing tumor angiogenesis within the breast.

METHODS

This new technique employs a commercial US transducer coupled with an array of NIR optical fibers mounted on a hand-held probe. The US image is used for lesion localization and for guiding optical imaging reconstruction. Optical sensors are used for imaging tumor total hemoglobin distribution, which is directly related to tumor angiogenesis.

RESULTS

Six large breast carcinomas were studied and microvessel density count was then performed on tissue samples obtained from these cancers. Two patients had locally advanced breast cancers and received neoadjuvant chemotherapy for 3 months. In one patient, before chemotherapy, the total hemoglobin distribution showed a high concentration at the cancer periphery; the distribution was later confined to the core area after 3 months of treatment. In another patient, as treatment progressed, the maximum hemoglobin concentration decreased from 255.3, to 147.5, to 76.9 micromol/l with an associated reduction in spatial extension. The other four patients had cancers of 2.0 to 3.0 cm in size and were imaged either at the time of core biopsy or definitive surgery. The histologic microvessel density counts from these tumor samples correlate to hemoglobin distributions with a correlation coefficient of 0.64 (P < .05).

CONCLUSION

These initial results suggest that this new imaging technique may have great potential in imaging the heterogeneous vascular distribution of larger breast cancers in vivo and in monitoring treatment-related changes in angiogenesis during chemotherapy.

Estimation of chest-wall-induced diffused wave distortion with the assistance of ultrasound

The chest-wall layer underneath breast tissue consists of muscles and bones, which induce distortion in near-infrared diffused waves measured at distant source--detector pairs when reflection geometry is used. A priori information on chest-wall depth obtained from coregistered real-time ultrasound can be used to assist in the removal of distant measurements. We applied Monte Carlo simulation to a simple two-layer model consisting of breast tissue and a chest wall to investigate chest-wall-induced distortion. The Monte Carlo method indicates that, when more than 50% of the received photons travel through the breast tissue layer before being detected, the detected signal may be useful for image reconstruction. The results of phantom experiments obtained from the two-layer model further validate the distortion problem and demonstrate imaging improvement after distant measurements have been filtered out. Clinical examples have shown similar imaging improvements on reconstructed absorption maps. Clinical data obtained from 20 patients with the chest-wall depths of less than 2 cm from the skin surface suggest that the cutoff distances of distorted measurements are largely related to the chest-wall depth and are relatively independent of the optical properties of tissue.

Ultrasound-guided optical tomographic imaging of malignant and benign breast lesions: initial clinical results of 19 cases

The diagnosis of solid benign and malignant tumors presents a unique challenge to all noninvasive imaging modalities. Ultrasound is used in conjunction with mammography to differentiate simple cysts from solid lesions. However, the overlapping appearances of benign and malignant lesions make ultrasound less useful in differentiating solid lesions, resulting in a large number of benign biopsies. Optical tomography using near-infrared diffused light has great potential for imaging functional parameters of 1) tumor hemoglobin concentration, 2) oxygen saturation, and 3) metabolism, as well as other tumor distinguishing characteristics. These parameters can differentiate benign from malignant lesions. However, optical tomography, when used alone, suffers from low spatial resolution and target localization uncertainty due to intensive light scattering. Our aim is to combine diffused light imaging with ultrasound in a novel way for the detection and diagnosis of solid lesions. Initial findings of two early-stage invasive carcinomas, one combined fibroadenoma and fibrocystic change with scattered foci of lobular neoplasia/lobular carcinoma in situ, and 16 benign lesions are reported in this paper. The invasive cancer cases reveal about two-fold greater total hemoglobin concentration (mean 119 micromol) than benign cases (mean 67 micromol), and suggest that the discrimination of benign and malignant breast lesions might be enhanced by this type of achievable optical quantification with ultrasound localization. Furthermore, the small invasive cancers are well localized and have wavelength-dependent appearance in optical absorption maps, whereas the benign lesions appear diffused and relatively wavelength-independent.

Portable near-infrared diffusive light imager for breast cancer detection

We present a frequency-domain near-infrared optical tomography system designed for breast cancer detection, in conjunction with conventional ultrasound. It features fast optical switching, three-wavelength excitations, and avalanche photodiode as detectors. Laser diodes at 660, 780, and 830 nm are used as light sources and their outputs are distributed sequentially to one of nine source fibers. An equivalent 130-dB isolation between electrical signals from different source channels is achieved with the optical switches of very low crosstalk. Ten detection channels, each of which includes a silicon avalanche photodiode, detect diffusive photon density waves simultaneously. The dynamic range of an avalanche photodiode is about 20 to 30 dB higher than that of a photomultiplier tube, thus eliminating the need for multistep system gain control. The entire system is compact in size (<0.051 m(3)) and fast in data acquisition (less than 2 sec for a complete scan). Calibration and the clinical experiment results are presented in the paper.

Dual-mesh optical tomography reconstruction method with a depth correction that uses a priori ultrasound information

A dual-mesh reconstruction method with a depth correction for near-infrared diffused wave imaging with ultrasound localization is demonstrated by use of phantoms and clinical cancer cases. Column normalization is applied to the weight matrix obtained from the Born approximation to correct the depth-dependent problem in the reconstructed absorption maps as well as in the total hemoglobin concentration maps. With the depth correction, more uniform absorption maps for target layers at different depths are obtained from the phantoms, and the correlation between the reconstructed hemoglobin concentration maps of deeply located, large cancers and the histological microvessel density counts are dramatically improved.