3D parallel-detection microwave tomography for clinical breast imaging.

A biomedical microwave tomography system with 3D-imaging capabilities has been constructed and translated to the clinic. Updates to the hardware and reconfiguration of the electronic-network layouts in a more compartmentalized construct have streamlined system packaging. Upgrades to the data acquisition and microwave components have increased data-acquisition speeds and improved system performance. By incorporating analog-to-digital boards that accommodate the linear amplification and dynamic-range coverage our system requires, a complete set of data (for a fixed array position at a single frequency) is now acquired in 5.8 s. Replacement of key components (e.g., switches and power dividers) by devices with improved operational bandwidths has enhanced system response over a wider frequency range. High-integrity, low-power signals are routinely measured down to -130 dBm for frequencies ranging from 500 to 2300 MHz. Adequate inter-channel isolation has been maintained, and a dynamic range >110 dB has been achieved for the full operating frequency range (500-2900 MHz). For our primary band of interest, the associated measurement deviations are less than 0.33% and 0.5° for signal amplitude and phase values, respectively. A modified monopole antenna array (composed of two interwoven eight-element sub-arrays), in conjunction with an updated motion-control system capable of independently moving the sub-arrays to various in-plane and cross-plane positions within the illumination chamber, has been configured in the new design for full volumetric data acquisition. Signal-to-noise ratios (SNRs) are more than adequate for all transmit/receive antenna pairs over the full frequency range and for the variety of in-plane and cross-plane configurations. For proximal receivers, in-plane SNRs greater than 80 dB are observed up to 2900 MHz, while cross-plane SNRs greater than 80 dB are seen for 6 cm sub-array spacing (for frequencies up to 1500 MHz). We demonstrate accurate recovery of 3D dielectric property distributions for breast-like phantoms with tumor inclusions utilizing both the in-plane and new cross-plane data.

Comparison of fibroglandular tissue distributions for microwave tomographic breast images with complementary MR T2 weighted images.

We have recently demonstrated good correlation between the recovered permittivity from microwave imaging (MIS) and the recovered water content from near infrared imaging (NIR) for a common set of normal patients undergoing associated breast examinations. We have subsequently conducted a small sample of comparison breast examinations between microwave imaging and MR to assess possible correlation between the location and extent of the fibroglandular as seen on MR images with increased permittivity zones of the microwave images. From various physiological and MR breast studies, it has been shown that the fibroglandular regions are generally comprised of significantly higher levels of water than the more dominant adipose tissue. The initial results of this study are quite encouraging and demonstrate obvious correlations between the permittivity and MR-recovered fibroglandular regions for a set of patients with widely varying tissue type variations. In addition, they illustrate the value of extracting diagnostic information from multiple modalities especially where the amount of direct in vivo property measurements is limited or nonexistent.

Image accuracy improvements in microwave tomographic thermometry: phantom experience.

Evaluation of a laboratory-scale microwave imaging system for non-invasive temperature monitoring has previously been reported with good results in terms of both spatial and temperature resolution. However, a new formulation of the reconstruction algorithm in terms of the log-magnitude and phase of the electric fields has dramatically improved the ability of the system to track the temperature-dependent electrical conductivity distribution. This algorithmic enhancement was originally implemented as a way of improving overall imaging capability in cases of large, high contrast permittivity scatterers, but has also proved to be sensitive to subtle conductivity changes as required in thermal imaging. Additional refinements in the regularization procedure have strengthened the reliability and robustness of image convergence. Imaging experiments were performed for a single heated target consisting of a 5.1 cm diameter PVC tube located within 15 and 25 cm diameter monopole antenna arrays, respectively. The performance of both log-magnitude/phase and complex-valued reconstructions when subjected to four different regularization schemes has been compared based on this experimental data. The results demonstrate a significant accuracy improvement (to 0.2 degrees C as compared with 1.6 degrees C for the previously published approach) in tracking thermal changes in phantoms where electrical properties vary linearly with temperature over a range relevant to hyperthermia cancer therapy.

Microwave thermal imaging: initial in vivo experience with a single heating zone.

The deployment of hyperthermia as a routine adjuvant to radiation or chemotherapy is limited largely by the inability to devise treatment plans which can be monitored through temperature distribution feedback during therapy. A non-invasive microwave tomographic thermal imaging system is currently being developed which has previously exhibited excellent correlation between the recovered electrical conductivity of a heated zone and its actual temperature change during phantom studies. To extend the validation of this approach in vivo, the imaging system has been re-configured for small animal experiments to operate within the bore of a CT scanner for anatomical and thermometry registration. A series of 5-7 day old pigs have been imaged during hyperthermia with a monopole antenna array submerged in a saline tank where a small plastic tube surgically inserted the length of the abdomen has been used to create a zone of heated saline at pre-selected temperatures. Tomographic microwave data over the frequency range of 300-1000 MHz of the pig abdomen in the plane perpendicular to the torso is collected at regular intervals after the tube saline temperatures have settled to the desired settings. Images are reconstructed over a range of operating frequencies. The tube location is clearly visible and the recovered saline conductivity varies linearly with the controlled temperature values. Difference images utilizing the baseline state prior to heating reinforces the linear relationship between temperature and imaged saline conductivity. Demonstration of in vivo temperature recovery and correlation with an independent monitoring device is an important milestone prior to clinical integration of this non-invasive imaging system with a thermal therapy device.

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.

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.

Evaluation of an iterative reconstruction method for quantitative elastography.

This paper describes an inverse reconstruction technique based on a modified Newton Raphson iterative scheme and the finite element method, which has been developed for computing the spatial distribution of Young's modulus from within soft tissues. Computer simulations were conducted to determine the relative merits of reconstructing tissue elasticity using knowledge of (a) known displacement boundary conditions (DBC), and (b) known stress boundary conditions (SBC). The results demonstrated that computing Young's modulus using knowledge of SBC allows accurate quantification of Young's modulus. However, the quality of the images produced using this reconstruction approach was dependent on the Young's modulus distribution assumed at the start of the reconstruction procedure. Computing Young's modulus from known DBC provided relative estimates of tissue elasticity which, despite the disadvantage of not being able to accurately quantify Young's modulus, formed images that were generally superior in quality to those produced using the known SBC, and were not affected by the trial solution. The results of preliminary experiments on phantoms demonstrated that this reconstruction technique is capable in practice of improving the fidelity of tissue elasticity images, reducing the artefacts otherwise present in strain images, and recovering Young's modulus images that possess excellent spatial and contrast resolution.

The intensity dependence of lesion position shift during focused ultrasound surgery.

Knowledge of the spatial distribution of intensity loss from an ultrasonic beam is critical for predicting lesion formation in focused ultrasound (US) surgery (FUS). To date, most models have used linear propagation models to predict intensity profiles required to compute the temporally varying temperature distributions used to compute thermal dose contours. These are used to predict the extent of thermal damage. However, these simulations fail to describe adequately the abnormal lesion formation behaviour observed during ex vivo experiments in cases for which the transducer drive levels are varied over a wide range. In such experiments, the extent of thermal damage has been observed to move significantly closer to the transducer with increased transducer drive levels than would be predicted using linear-propagation models. The first set of simulations described herein use the KZK (Khokhlov-Zabolotskaya-Kuznetsov) nonlinear propagation model with the parabolic approximation for highly focused US waves to demonstrate that both the peak intensity and the lesion positions do, indeed, move closer to the transducer. This illustrates that, for accurate modelling of heating during FUS, nonlinear effects should be considered. Additionally, a first order approximation has been employed that attempts to account for the abnormal heat deposition distributions that accompany high transducer drive level FUS exposures where cavitation and boiling may be present. The results of these simulations are presented. It is suggested that this type of approach may be a useful tool in understanding thermal damage mechanisms.

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

Fixed-array microwave imaging with multisensor data acquisition can suffer from nonactive antenna element interactions which cause distortions in the measurements. In Part I of a two-part paper, we develop a nonactive antenna compensation model for incorporation in model-based near-field microwave image reconstruction methods. The model treats the nonactive members of the antenna array as impedance boundary conditions applied over a cylindrical surface of finite radius providing two parameters, the effective antenna radius and impedance factor, which can be determined empirically from measured data. Results show that the effective radius and impedance factor provide improved fits to experimental data in homogeneous phantoms where measurements are obtained with and without the presence of the nonactive antenna elements. Once deduced, these parameters are incorporated into the nonactive antenna compensation model and lead to systematic data-model match improvements in heterogeneous phantoms. While the improvements afforded by the nonactive antenna model are small on a per measurement basis, they are not insignificant. As shown in Part II, inclusion of this new model for nonactive antenna compensation produces significantly higher quality image reconstructions from measurements obtained with a fixed-array data acquisition system over the frequency band 500-900 MHz.

Nonactive antenna compensation for fixed-array microwave imaging: Part II--Imaging results.

Model-based imaging techniques utilizing microwave signal illumination rely heavily on the ability to accurately represent the wave propagation with a suitable numerical model. To date, the highest quality images from our prototype system have been achieved utilizing a single transmitter/single receiver measurement system where both antennas are manually repositioned to facilitate multiple illuminations of the imaging region, thus requiring long data acquisition times. In an effort to develop a system that can acquire data in a real time manner, a 32-channel network has been fabricated with all ports capable of being electronically selected for either transmit or receive mode. The presence of a complete array of antenna elements at data collection time perturbs the field distributions being measured, which can subsequently degrade the image reconstruction due to increased data-model mismatch. Incorporating the nonactive antenna-compensation model from Part I of this paper into our hybrid element near field image reconstruction algorithm is shown to restore image quality when fixed antenna-array data acquisition is used. Improvements are most dramatic for inclusions located in near proximity to the antenna array itself, although cases of improvement in the recovery of centered heterogeneities are also illustrated. Increases in the frequency of illumination are found to warrant an increased need for nonactive antenna compensation. Quantitative measures of recovered inclusion shape and position reveal a systematic improvement in image reconstruction quality when the nonactive antenna-compensation model is employed. Improvements in electrical property value recovery of localized heterogeneities are also observed. Image reconstructions in freshly excised breast tissue illustrate the applicability of the approach when used with our two-dimensional microwave imaging system.

A 3-D finite-element model for computation of temperature profiles and regions of thermal damage during focused ultrasound surgery exposures.

Although there have been numerous models implemented for modeling thermal diffusion effects during focused ultrasound surgery (FUS), most have limited themselves to representing simple situations for which analytical solutions and the use of cylindrical geometries sufficed. For modeling single lesion formation and the heating patterns from a single exposure, good results were achieved in comparison with experimental results for predicting lesion size, shape and location. However, these types of approaches are insufficient when considering the heating of multiple sites with FUS exposures when the time interval between exposures is short. In such cases, the heat dissipation patterns from initial exposures in the lesion array formation can play a significant role in the heating patterns for later exposures. Understanding the effects of adjacent lesion formation, such as this, requires a three-dimensional (3-D) representation of the bioheat equation. Thus, we have developed a 3-D finite-element representation for modeling the thermal diffusion effects during FUS exposures in clinically relevant tissue volumes. The strength of this approach over past methods is its ability to represent arbitrarily shaped 3-D situations. Initial simulations have allowed calculation of the temperature distribution as a function of time for adjacent FUS exposures in excised bovine liver, with the individually computed point temperatures comparing favorably with published measurements. In addition to modeling these temperature distributions, the model was implemented in conjunction with an algorithm for calculating the thermal dose as a way of predicting lesion shape. Although used extensively in conventional hyperthermia applications, this thermal dose criterion has only been applied in a limited number of simulations in FUS for comparison with experimental measurements. In this study, simulations were run for focal depths 2 and 3 cm below the surface of pig's liver, using multiple intensity levels and exposure times. The results also compare favorably to published in vitro experimental measurements, which bodes well for future application to more complex problems, such as the modeling of multiple lesion arrays within complex anatomical geometries.

Microwave imaging for tissue assessment: initial evaluation in multitarget tissue-equivalent phantoms.

A prototype microwave imaging system is evaluated for its ability to recover two-dimensional (2-D) electrical property distributions under transverse magnetic (TM) illumination using multitarget tissue equivalent phantoms. Experiments conducted in a surrounding lossy saline tank demonstrate that simultaneous recovery of both the real and imaginary components of the electrical property distribution is possible using absolute imaging procedures over a frequency range of 300-700 MHz. Further, image reconstructions of embedded tissue-equivalent targets are found to be quantitative not only with respect to geometrical factors such as object size and location but also electrical composition. Quantitative assessments based on full-width half-height criteria reveal that errors in diameter estimates of reconstructed targets are less than 10 mm in all cases, whereas, positioning errors are less than 1 mm in single object experiments but degrade to 4-10 mm when multiple targets are present. Recovery of actual electrical properties is found to be frequency dependent for the real and imaginary components with background values being typically within 10-20% of their correct size and embedded object having similar accuracies as a percentage of the electrical contrast, although errors as high as 50% can occur. The quantitative evaluation of imaging performance has revealed potential advantages in a two-tiered receiver antenna configuration whose measured field values are more sensitive to target region changes than the typical tomographic type of approach which uses reception sites around the full target region perimeter. This measurement strategy has important implications for both the image reconstruction algorithm where there is a premium on minimizing problem size without sacrificing image quality and the hardware system design which seeks to economize on the amount of measured data required for quantitative image reconstruction while maximizing its sensitivity to target perturbations.

An active microwave imaging system for reconstruction of 2-D electrical property distributions.

The goal of this work is to develop a microwave-based imaging system for hyperthermia treatment monitoring and assessment. Toward this end, a four transmit channel and four receive channel hardware device and concomitant image reconstruction algorithm have been realized. The hardware is designed to measure electric fields (i.e., amplitude and phase) at various locations in a phantom tank with and without the presence of various heterogeneities using standard heterodyning principles. Particular attention has been paid to designing a receiver with better than 115 dB of linear dynamic range which is necessary for imaging biological tissue which often has very high conductivity, especially for tissues with high water content. A calibration procedure has been developed to compensate for signal loss due to three-dimensional radiation in the measured data, since the reconstruction process is only two-dimensional at the present time. Results are shown which demonstrate the stability and accuracy of the measurement system, the extent to which the forward computational model agrees with the measured field distribution when the electrical properties are known, and image reconstructions of electrically unknown targets of varying diameter. In the latter case, images of both the reactive and resistive component of the electrical property distribution have been recoverable. Quantitative information on object location, size, and electrical properties results when the target is approximately one-half wavelength in size. Images of smaller objects lack the same level of quantitative information, but remain qualitatively correct.

A dual mesh scheme for finite element based reconstruction algorithms.

The finite element (FE) method has found several applications in emerging imaging modalities, especially microwave imaging which has been shown to be potentially useful in a number of areas including thermal estimation. In monitoring temperature distributions, the biological phenomena of temperature variations of tissue dielectric properties is exploited. By imaging these properties and their changes during such therapies as hyperthermia, temperature distributions can be deduced using difference imaging techniques. The authors focus on a microwave imaging problem where the hybrid element (HE) method is used in conjunction with a dual mesh scheme in an effort to image complex wavenumbers, k(2). The dual mesh scheme is introduced to improve the reconstructed images of tissue properties and is ideally suited for systems using FE methods as their computational base. Since the electric fields typically vary rapidly over a given body when irradiated by high-frequency electromagnetic sources, a dense mesh is needed for these fields to be accurately represented. Conversely, k(2) may be fairly constant over subregions of the body which would allow for a less dense sampling of this parameter in those regions. In the dual mesh system employed, the first mesh, which is uniformly dense, is used for calculating the electric fields over the body whereas the second mesh, which is nonuniform and less dense, is used for representing the k(2) distribution within the region of interest. The authors examine the 2-D TM polarization case for a pair of dielectric distributions on both a large and small problem to demonstrate the flexibility of the dual mesh method along with some of the difficulties associated with larger imaging problems. Results demonstrate the capabilities of the dual mesh concept in comparison to a single mesh approach for a variety of test cases, suggesting that the dual mesh method is critical for FE based image reconstruction where rapidly varying physical quantities are used to recover smoother property profiles, as can occur in microwave imaging of biological bodies.

Finite element simulation of Sigma 60 heating in the Utah phantom: computed and measured data compared.

An initial series of comparisons are made between finite element computations and laboratory measurements obtained during heterogeneous phantom heating with the Sigma 60 applicator. The phantom is a relatively complex, though still idealized, rendering of the pelvic area which has been used to study the deep heating characteristics of the Sigma 60 in this anatomy. Direct electric field measurements as well as inferred SAR through transient temperature analysis are plotted against computed results along 11 one-dimensional tracks through the phantom. Quantitative comparisons provided through the track-by-track analysis show generally good agreement between computation and measurement. The finite element method is found to predict well the jumps in the electric field when polarized perpendicularly to a muscle/fat interface. Visualizations of the complete three-dimensional distributions are also highlighted and correlate well with physical reasoning about the expected behaviour of the fields produced. Some discrepancies in the data persist and are discussed and analysed in depth. They underscore the difficulties that can arise in performing comparisons between measured and computed results and stress the need for careful and thorough investigations when attempting these types of model validation studies.