Nonactive antenna compensation for fixed-array microwave imaging--Part I: Model development

Low Cost, High Performance, 16-Channel Microwave Measurement System for Tomographic Applications

We have developed a multichannel software defined radio-based transceiver measurement system for use in general microwave tomographic applications. The unit is compact enough to fit conveniently underneath the current illumination tank of the Dartmouth microwave breast imaging system. The system includes 16 channels that can both transmit and receive and it operates from 500 MHz to 2.5 GHz while measuring signals down to −140 dBm. As is the case with multichannel systems, cross-channel leakage is an important specification and must be lower than the noise floors for each receiver. This design exploits the isolation inherent when the individual receivers for each channel are physically separate; however, these challenging specifications require more involved signal isolation techniques at both the system design level and the individual, shielded component level. We describe the isolation design techniques for the critical system elements and demonstrate specification compliance at both the component and system level.

Enhancing Tumor Detection in IR-UWB Breast Cancer System

An ultra-wideband (UWB) microwave system for breast cancer detection is presented. The proposed system includes monocycle pulse generator, antipodal Vivaldi antenna, breast model, and calibration algorithm for tumor detection. Firstly, our pulse generator employs transmission gate in glitch generator to achieve several advantages such as low power consumption and low ringing level. Secondly, the antipodal Vivaldi antenna is designed assuming FR4 dielectric substrate material, and developed antenna element (80 × 80 mm²) features a −10 dB return loss and bandwidth ranges from 2.3 GHz to more than 11 GHz. Thirdly, the phantom breast can be modeled as a layer of skin, fat, and then tumor is inserted in this layer. Finally, subtract and add algorithm (SAD) is used as a calibration algorithm in tumor detection system. The proposed system suggested that horizontal antenna position with 90° between transmitting and receiving antennas is localized as a suitable antenna position with different rotating location and a 0.5 cm near to phantom. The mean advantages of this localization and tracking position around breast is a high received power signal approximately around mv as a higher recognized signal in tumor detection. Using our proposed system we can detect tumor in 5 mm diameter.

Integration of microwave tomography with magnetic resonance for improved breast imaging

PURPOSE

Breast magnetic resonance imaging is highly sensitive but not very specific for the detection of breast cancer. Opportunities exist to supplement the image acquisition with a more specific modality provided the technical challenges of meeting space limitations inside the bore, restricted breast access, and electromagnetic compatibility requirements can be overcome. Magnetic resonance (MR) and microwave tomography (MT) are complementary and synergistic because the high resolution of MR is used to encode spatial priors on breast geometry and internal parenchymal features that have distinct electrical properties (i.e., fat vs fibroglandular tissue) for microwave tomography.

METHODS

The authors have overcome integration challenges associated with combining MT with MR to produce a new coregistered, multimodality breast imaging platform--magnetic resonance microwave tomography, including: substantial illumination tank size reduction specific to the confined MR bore diameter, minimization of metal content and composition, reduction of metal artifacts in the MR images, and suppression of unwanted MT multipath signals.

RESULTS

MR SNR exceeding 40 dB can be obtained. Proper filtering of MR signals reduces MT data degradation allowing MT SNR of 20 dB to be obtained, which is sufficient for image reconstruction. When MR spatial priors are incorporated into the recovery of MT property estimates, the errors between the recovered versus actual dielectric properties approach 5%.

CONCLUSIONS

The phantom and human subject exams presented here are the first demonstration of combining MT with MR to improve the accuracy of the reconstructed MT images.

Microwave imaging for neoadjuvant chemotherapy monitoring: initial clinical experience

INTRODUCTION

Microwave tomography recovers images of tissue dielectric properties, which appear to be specific for breast cancer, with low-cost technology that does not present an exposure risk, suggesting the modality may be a good candidate for monitoring neoadjuvant chemotherapy.

METHODS

Eight patients undergoing neoadjuvant chemotherapy for locally advanced breast cancer were imaged longitudinally five to eight times during the course of treatment. At the start of therapy, regions of interest (ROIs) were identified from contrast-enhanced magnetic resonance imaging studies. During subsequent microwave examinations, subjects were positioned with their breasts pendant in a coupling fluid and surrounded by an immersed antenna array. Microwave property values were extracted from the ROIs through an automated procedure and statistical analyses were performed to assess short term (30 days) and longer term (four to six months) dielectric property changes.

RESULTS

Two patient cases (one complete and one partial response) are presented in detail and demonstrate changes in microwave properties commensurate with the degree of treatment response observed pathologically. Normalized mean conductivity in ROIs from patients with complete pathological responses was significantly different from that of partial responders (P value = 0.004). In addition, the normalized conductivity measure also correlated well with complete pathological response at 30 days (P value = 0.002).

CONCLUSIONS

These preliminary findings suggest that both early and late conductivity property changes correlate well with overall treatment response to neoadjuvant therapy in locally advanced breast cancer. This result is consistent with earlier clinical outcomes that lesion conductivity is specific to differentiating breast cancer from benign lesions and normal tissue.

A study of matching fluid loss in a biomedical microwave tomography system

PURPOSE

Effective imaging of human tissue with microwave tomography systems requires a matching fluid to reduce the wave reflections at the tissue boundary. Further, in order to match the idealized mathematical model used for imaging with the complicated physical measurement environment, loss must be added to the matching fluid. Both too little and too much loss result in low-quality images, but due to the nonlinear nature of the imaging problem, the exact nature of loss-to-image quality cannot be predicted a priori. Possible optimal loss levels include a single, highly sensitive value, or a broad range of acceptable losses. Herein, the authors outline a process of determining an appropriate level of loss inside the matching fluid and attempt to determine the bounds for which the images are the highest quality.

METHODS

Our biomedical microwave tomography system is designed for 2D limb imaging, operating from 0.8 to 1.2 GHz. Our matching fluid consists of deionized water with various levels of loss introduced by the addition of table salt. Using two homogeneous tissue-mimicking phantoms, and eight different matching fluids of varying salt concentrations, the authors introduce quantitative image quality metrics based on L-norms, mean values, and standard deviations to test the tomography system and assess image quality. Images are generated with a balanced multiplicative regularized contrast source inversion algorithm. The authors further generate images of a human forearm which may be analyzed qualitatively.

RESULTS

The image metrics for the phantoms support the claim that the worst images occur at the extremes of high and low salt concentrations. Importantly, the image metrics show that there exists a broad range of salt concentrations that result in high-quality images, not a single optimal value. In particular, 2.5-4.5 g of table salt per liter of deionized water provide the best reconstruction quality for simple phantoms. The authors argue that qualitatively, the human forearm data provide the best images at approximately the same salt concentrations.

CONCLUSIONS

There exists a relatively large-range of matching fluid losses (i.e., salt concentrations) that provide similar image quality. In particular, it is not necessary to spend time highly optimizing the level of loss in the matching fluid.

Fast 3-d tomographic microwave imaging for breast cancer detection

Microwave breast imaging (using electromagnetic waves of frequencies around 1 GHz) has mostly remained at the research level for the past decade, gaining little clinical acceptance. The major hurdles limiting patient use are both at the hardware level (challenges in collecting accurate and noncorrupted data) and software level (often plagued by unrealistic reconstruction times in the tens of hours). In this paper we report improvements that address both issues. First, the hardware is able to measure signals down to levels compatible with sub-centimeter image resolution while keeping an exam time under 2 min. Second, the software overcomes the enormous time burden and produces similarly accurate images in less than 20 min. The combination of the new hardware and software allows us to produce and report here the first clinical 3-D microwave tomographic images of the breast. Two clinical examples are selected out of 400+ exams conducted at the Dartmouth Hitchcock Medical Center (Lebanon, NH). The first example demonstrates the potential usefulness of our system for breast cancer screening while the second example focuses on therapy monitoring.

Clinical microwave tomographic imaging of the calcaneus: a first-in-human case study of two subjects

We have acquired 2-D and 3-D microwave tomographic images of the calcaneus bones of two patients to assess correlation of the microwave properties with X-ray density measures. The two volunteers were selected because each had one leg immobilized for at least six weeks during recovery from a lower leg injury. A soft-prior regularization technique was incorporated with the microwave imaging to quantitatively assess the bulk dielectric properties within the bone region. Good correlation was observed between both permittivity and conductivity and the computed tomography-derived density measures. These results represent the first clinical examples of microwave images of the calcaneus and some of the first 3-D tomographic images of any anatomical site in the living human.

Microwave tomography in the context of complex breast cancer imaging

The notion of applying microwave imaging to breast cancer imaging has been studied at various levels by numerous scientists. The earliest appeal of this concept related to the presumably high property contrast between benign and malignant tissue that was unique to the breast. Subsequent published studies have shown that this assumption was overly simplistic and that the tissue property heterogeneity is considerable within the breast. As we have expanded the clinical use of our microwave tomographic system, we are now using this approach to monitor tumor progressions during neoadjuvant chemotherapy. In these cases, while we can still characterize and track the tumor progression, we have observed a new phenomenon. Very often these cancer patients exhibit skin thickening near the tumor site. Our images have reconstructed elevated dielectric properties along the breast surface associated with the accompanying edema. These observations further add to the complex nature of breast dielectric properties and the challenges for imaging them using microwave interrogation.

Importance of phase unwrapping for the reconstruction of microwave tomographic images

Microwave image reconstruction is typically based on a regularized least-square minimization of either the complex-valued field difference between recorded and modeled data or the logarithmic transformation of these field differences. Prior work has shown anecdotally that the latter outperforms the former in limited surveys of simulated and experimental phantom results. In this paper, we provide a theoretical explanation of these empirical findings by developing closed form solutions for the field and the inverted electromagnetic property parameters in one dimension which reveal the dependency of the estimated permittivity and conductivity on the absolute (unwrapped) phase of the measured signal at the receivers relative to the source transmission. The analysis predicts the poor performance of complex-valued field minimization as target size and/or frequency and electromagnetic contrast increase. Such poor performance is avoided by logarithmic transformation and preservation of absolute measured signal phase. Two-dimensional experiments based on both synthetic and clinical data are used to confirm these findings. Robustness of the logarithmic transformation to variation in the initial guess of the reconstructed target properties is also shown. The results are generalizable to three dimensions and indicate that the minimization form with logarithmic transformation offers image reconstruction performance characteristics that are much more desirable for medial microwave imaging applications relative to minimizing differences in complex-valued field quantities.

Importance of phase unwrapping for the reconstruction of microwave tomographic images

Microwave image reconstruction is typically based on a regularized least-square minimization of either the complex-valued field difference between recorded and modeled data or the logarithmic transformation of these field differences. Prior work has shown anecdotally that the latter outperforms the former in limited surveys of simulated and experimental phantom results. In this paper, we provide a theoretical explanation of these empirical findings by developing closed form solutions for the field and the inverted electromagnetic property parameters in one dimension which reveal the dependency of the estimated permittivity and conductivity on the absolute (unwrapped) phase of the measured signal at the receivers relative to the source transmission. The analysis predicts the poor performance of complex-valued field minimization as target size and/or frequency and electromagnetic contrast increase. Such poor performance is avoided by logarithmic transformation and preservation of absolute measured signal phase. Two-dimensional experiments based on both synthetic and clinical data are used to confirm these findings. Robustness of the logarithmic transformation to variation in the initial guess of the reconstructed target properties is also shown. The results are generalizable to three dimensions and indicate that the minimization form with logarithmic transformation offers image reconstruction performance characteristics that are much more desirable for medial microwave imaging applications relative to minimizing differences in complex-valued field quantities.

Initial clinical experience with microwave breast imaging in women with normal mammography

RATIONALE AND OBJECTIVES

We have developed a microwave tomography system for experimental breast imaging.

MATERIALS AND METHODS

In this article, we illustrate a strategy for optimizing the coupling liquid for the antenna array based on in vivo measurement data. We present representative phantom experiments to illustrate the imaging system's ability to recover accurate property distributions over the range of dielectric properties expected to be encountered clinically. To demonstrate clinical feasibility and assess the microwave properties of the normal breast in vivo, we summarize our initial experience with microwave breast exams of 43 women with negative mammography according to the Breast Imaging Reporting and Data System (BI-RADS 1).

RESULTS

The clinical results show a high degree of bilateral symmetry in the whole breast average microwave properties. Focal assessments of microwave properties are associated with breast tissue composition evaluated through radiographic density categorization verified through magnetic resonance image correlation in selected cases. Specifically, both whole-breast average and local microwave properties increase with increasing radiographic density, in which the latter exhibits a more substantial rise.

CONCLUSION

These findings support our hypothesis that water content variations in the breast play an influential role in dictating the overall dielectric property distributions and indicate that the microwave properties in the breast are more heterogeneous than previously believed based on ex vivo property measurements reported in the literature.

Microwave image reconstruction from 3-D fields coupled to 2-D parameter estimation

An efficient Gauss-Newton iterative imaging technique utilizing a three-dimensional (3-D) field solution coupled to a two-dimensional (2-D) parameter estimation scheme (3-D/2-D) is presented for microwave tomographic imaging in medical applications. While electromagnetic wave propagation is described fully by a 3-D vector field, a 3-D scalar model has been applied to improve the efficiency of the iterative reconstruction process with apparently limited reduction in accuracy. In addition, the image recovery has been restricted to 2-D but is generalizable to three dimensions. Image artifacts related primarily to 3-D effects are reduced when compared with results from an entirely two-dimensional inversion (2-D/2-D). Important advances in terms of improving algorithmic efficiency include use of a block solver for computing the field solutions and application of the dual mesh scheme and adjoint approach for Jacobian construction. Methods which enhance the image quality such as the log-magnitude/unwrapped phase minimization were also applied. Results obtained from synthetic measurement data show that the new 3-D/2-D algorithm consistently outperforms its 2-D/2-D counterpart in terms of reducing the effective imaging slice thickness in both permittivity and conductivity images over a range of inclusion sizes and background medium contrasts.

Quantification of 3-D field effects during 2-D microwave imaging

Two-dimensional (2-D) approaches to microwave imaging have dominated the research landscape primarily due to the moderate levels of measurement data, data-acquisition time, and computational costs required. Three-dimensional (3-D) approaches have been investigated in simulation, phantom, and animal experiments. While 3-D approaches are certainly important in terms of the potential to improve image quality, their associated costs are significant at this time. In addition, benchmarks are needed to evaluate these new generation systems as more 3-D methods begin to appear. In this paper, we present a systematic series of experiments which assess the capability of our 2-D system to image classical 3-D geometries. We demonstrate where current methods suffer from 3-D effects but also identify situations where they remain quite useful. Comparisons between reconstructions utilizing phantom measurements and simulated 3-D data are also shown to validate the results. These findings suggest that for certain biomedical applications, 2-D approaches remain quite attractive.

Time-domain reconstruction for thermoacoustic tomography in a spherical geometry

Reconstruction-based microwave-induced thermoacoustic tomography in a spherical configuration is presented. Thermoacoustic waves from biological tissue samples excited by microwave pulses are measured by a wide-band unfocused ultrasonic transducer, which is set on a spherical surface enclosing the sample. Sufficient data are acquired from different directions to reconstruct the microwave absorption distribution. An exact reconstruction solution is derived and approximated to a modified backprojection algorithm. Experiments demonstrate that the reconstructed images agree well with the original samples. The spatial resolution of the system reaches 0.5 mm.

Pre-scaled two-parameter Gauss-Newton image reconstruction to reduce property recovery imbalance

Gauss-Newton image reconstruction in microwave imaging can be formulated in terms of a single complex quantity, the wave number squared (k2), with the understanding that the relative permittivity and conductivity images can be extracted afterwards through a simple constitutive relationship. However, this approach ignores the fact that the magnitude of the average real and imaginary components can be considerably out of balance depending on the operating frequency and tissue characteristics which can inadvertently imbalance the process in favour of one parameter over the other. In an effort to achieve property recovery which is balanced, we introduce a pre-scaling procedure at the property update stage of the reconstruction. Utilization of this concept in conjunction with our two-step regularization process for both simulation and phantom experiments demonstrates that the penalty term weighting parameters for the optimal mean-squared property errors for the two recovered distributions (relative permittivity and conductivity) together with that yielding the lowest least-squared electric field error coincide only when the scaling is applied. The scheme provides a means for simultaneous optimization of the two permittivity and conductivity images.

A two-stage microwave image reconstruction procedure for improved internal feature extraction

We have developed a two-stage Gauss-Newton reconstruction process with an automatic procedure for determining the regularization parameter. The combination is utilized by our microwave imaging system and has facilitated recovery of quantitatively improved images. The first stage employs a Levenberg-Marquardt regularization along with a spatial filtering technique for a few iterations to produce an intermediate image. In effect, the first set of iterative image reconstruction steps synthesizes a priori information from the measurement data versus actually requiring physical prior information on the interrogated object. Because of the interaction of the Levenberg-Marquardt regularization and spatial filtering at each iteration, the intermediate image produced from the first reconstruction stage represents an improvement in terms of the least squared error over the initial uniform guess; however, it has not completely converged in a least squared sense. The second stage involves using this distribution as a priori information in an iteratively regularized Gauss-Newton reconstruction with a weighted Euclidean distance penalty term. The penalized term restricts the final image to a vicinity (determined by the scale of the weighting parameter) about the intermediate image while allowing more flexibility in extracting internal object structures. The second stage makes use of an empirical Bayesian/random effects model that enables an optimal determination of the weighting parameter of the penalized term. The new approach demonstrates quantifiably improved images in simulation, phantom and in vivo experiments with particularly striking improvements with respect to the recovery of heterogeneities internal to large, high contrast scatterers such as encountered when imaging the human breast in a water-coupled configuration.

Microwave image reconstruction utilizing log-magnitude and unwrapped phase to improve high-contrast object recovery

Reconstructing images of large high-contrast objects with microwave methods has proved difficult. Successful images have generally been obtained by using a priori information to constrain the image reconstruction to recover the correct electromagnetic property distribution. In these situations, the measured electric field phases as a function of receiver position around the periphery of the imaging field-of-view vary rapidly often undergoing changes of greater than pi radians especially when the object contrast and illumination frequency increase. In this paper, we introduce a modified form of a Maxwell equation model-based image reconstruction algorithm which directly incorporates log-magnitude and phase of the measured electric field data. By doing so, measured phase variation can be unwrapped and distributed over more than one Rieman sheet in the complex plane. Simulation studies and microwave imaging experiments demonstrate that significant image quality enhancements occur with this approach for large high-contrast objects. Simple strategies for visualizing and unwrapping phase values as a function of the transmitter and receiver positions within our microwave imaging array are described. Metrics of the degree of phase variation expressed in terms of the amount and extent of phase wrapping are defined and found to be figures-of-merit which estimate when it is critical to deploy the new image reconstruction approach. In these cases, the new algorithm recovers high-quality images without resorting to the use of a priori information on object contrast and/or size as previously required.