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.

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.

We propose a 3-D-printed breast phantom for use in preclinical experimental microwave imaging studies. The phantom is derived from an MRI of a human subject; thus, it is anthropomorphic, and its interior is very similar to an actual distribution of fibroglandular tissues. Adipose tissue in the breast is represented by the solid plastic (printed) regions of the phantom, while fibroglandular tissue is represented by liquid-filled voids in the plastic. The liquid is chosen to provide a biologically relevant dielectric contrast with the printed plastic. Such a phantom enables validation of microwave imaging techniques. We describe the procedure for generating the 3-D-printed breast phantom and present the measured dielectric properties of the 3-D-printed plastic over the frequency range 0.5–3.5 GHz. We also provide an example of a suitable liquid for filling the fibroglandular voids in the plastic.

This paper presents an algorithm for estimating the location of the breast surface from scattered ultrawideband (UWB) microwave signals recorded across an antenna array. Knowing the location of the breast surface can improve imaging performance if incorporated as a priori information into recently proposed microwave imaging algorithms. These techniques transmit low-power microwaves into the breast using an antenna array, which in turn measures the scattered microwave signals for the purpose of detecting anomalies or changes in the dielectric properties of breast tissue. Our proposed surface identification algorithm consists of three procedures, the first of which estimates M points on the breast surface given M channels of measured microwave backscatter data. The second procedure applies interpolation and extrapolation to these M points to generate N > M points that are approximately uniformly distributed over the breast surface, while the third procedure uses these N points to generate a 3-D estimated breast surface. Numerical as well as experimental tests indicate that the maximum absolute error in the estimated surface generated by the algorithm is on the order of several millimeters. An error analysis conducted for a basic microwave radar imaging algorithm (least-squares narrowband beamforming) indicates that this level of error is acceptable. A key advantage of the algorithm is that it uses the same measured signals that are used for UWB microwave imaging, thereby minimizing patient scan time and avoiding the need for additional hardware.