Purpose: Breast density measurement has the potential to play an important role in individualized breast cancer risk assessment and prevention decisions. Routine evaluation of breast density will require the availability of a low-cost, nonionizing, three-dimensional ͑3-D͒ tomographic imaging modality that exploits a strong properties contrast between dense fibroglandular tissue and less dense adipose tissue. The purpose of this computational study is to investigate the performance of 3-D tomography using low-power microwaves to reconstruct the spatial distribution of breast tissue dielectric properties and to evaluate the modality for application to breast density characterization. Methods: State-of-the-art 3-D numerical breast phantoms that are realistic in both structural and dielectric properties are employed. The test phantoms include one sample from each of four classes of mammographic breast density. Since the properties of these phantoms are known exactly, these testbeds serve as a rigorous benchmark for the imaging results. The distorted Born iterative imaging method is applied to simulated array measurements of the numerical phantoms. The forward solver in the imaging algorithm employs the finite-difference time-domain method of solving the timedomain Maxwell's equations, and the dielectric profiles are estimated using an integral equation form of the Helmholtz wave equation. A multiple-frequency, bound-constrained, vector field inverse scattering solution is implemented that enables practical inversion of the large-scale 3-D problem. Knowledge of the frequency-dependent characteristic of breast tissues at microwave frequencies is exploited to obtain a parametric reconstruction of the dispersive dielectric profile of the interior of the breast. Imaging is performed on a high-resolution voxel basis and the solution is bounded by a known range of dielectric properties of the constituent breast tissues. The imaging method is validated using a breast phantom with a single, high-contrast interior scattering target in an otherwise homogeneous interior. The method is then used to image a set of realistic numerical breast phantoms of varied fibroglandular tissue density. Results: Imaging results are presented for each numerical phantom and show robustness of the method relative to tissue density. In each case, the distribution of fibroglandular tissues is well represented in the resulting images. The resolution of the images at the frequencies employed is wider than the feature dimensions of the normal tissue structures, resulting in a smearing of their reconstruction. Conclusions:The results of this study support the utility of 3-D microwave tomography for imaging the distribution of normal tissues in the breast, specifically, dense fibroglandular tissue versus less dense adipose tissue, and suggest that further investigation of its use for volumetric evaluation of breast density is warranted.

No abstract

Breast imaging via microwave tomography involves estimating the distribution of dielectric properties within the patient's breast on a discrete mesh. The number of unknowns in the discrete mesh can be very large for three-dimensional imaging, and this results in computational challenges. We propose a new approach where the discrete mesh is replaced with a relatively small number of smooth basis functions. The dimension of the tomography problem is reduced by estimating the coefficients of the basis functions instead of the dielectric properties at each element in the discrete mesh. The basis functions are constructed using knowledge of the location of the breast surface. The number of functions used in the basis can be varied to balance resolution and computational complexity. The reduced dimension of the inverse problem enables application of a computationally efficient, multiple-frequency inverse scattering algorithm in 3-D. The efficacy of the proposed approach is verified using two 3-D anatomically realistic numerical breast phantoms. It is shown for the case of single-frequency microwave tomography that the imaging accuracy is comparable to that obtained when the original discrete mesh is used, despite the reduction of the dimension of the inverse problem. Results are also shown for a multiple-frequency algorithm where it is computationally challenging to use the original discrete mesh.

The detection of early-stage tumors in the breast by microwave imaging is challenged by both the moderate endogenous dielectric contrast between healthy and malignant glandular tissues and the spatial resolution available from illumination at microwave frequencies. The high endogenous dielectric contrast between adipose and fibroglandular tissue structures increases the difficulty of tumor detection due to the high dynamic range of the contrast function to be imaged and the low level of signal scattered from a tumor relative to the clutter scattered by normal tissue structures. Microwave inverse scattering techniques, used to estimate the complete spatial profile of the dielectric properties within the breast, have the potential to reconstruct both normal and cancerous tissue structures. However, the ill-posedness of the associated inverse problem often limits the frequency of microwave illumination to the UHF band within which early-stage cancers have subwavelength dimensions. In this computational study, we examine the reconstruction of small, compact tumors in three-dimensional numerical breast phantoms by a multiple-frequency inverse scattering solution. Computer models are also employed to investigate the use of exogenous contrast agents for enhancing tumor detection. Simulated array measurements are acquired before and after the introduction of the assumed contrast effects for two specific agents currently under consideration for breast imaging: microbubbles and carbon nanotubes. Differential images of the applied contrast demonstrate the potential of the approach for detecting the preferential uptake of contrast agents by malignant tissues.