Computational modeling of microwave ablation with thermal accelerants.

PURPOSE: To develop a computational model of microwave ablation (MWA) with a thermal accelerant gel and apply the model toward interpreting experimental observations in ex vivo bovine and in vivo porcine liver. METHODS: A 3D coupled electromagnetic-heat transfer model was implemented to characterize thermal profiles within ex vivo bovine and in vivo porcine liver tissue during MWA with the HeatSYNC thermal accelerant. Measured temperature dependent dielectric and thermal properties of the HeatSYNC gel were applied within the model. Simulated extents of MWA zones and transient temperature profiles were compared against experimental measurements in ex vivo bovine liver. Model predictions of thermal profiles under in vivo conditions in porcine liver were used to analyze thermal ablations observed in prior experiments in porcine liver in vivo. RESULTS: Measured electrical conductivity of the HeatSYNC gel was ∼83% higher compared to liver at room temperature, with positive linear temperature dependency, indicating increased microwave absorption within HeatSYNC gel compared to tissue. In ex vivo bovine liver, model predicted ablation zone extents of (31.5 × 36) mm with the HeatSYNC, compared to (32.9 ± 2.6 × 40.2 ± 2.3) mm in experiments (volume differences 4 ± 4.1 cm3). Computational models under in vivo conditions in porcine liver suggest approximating the HeatSYNC gel spreading within liver tissue during ablations as a plausible explanation for larger ablation zones observed in prior in vivo studies. CONCLUSION: Computational models of MWA with thermal accelerants provide insight into the impact of accelerant on MWA, and with further development, could predict ablations with a variety of gel injection sites.

Discrete Dipole Approximation-Based Microwave Tomography for Fast Breast Cancer Imaging.

This paper describes a fast microwave tomography reconstruction algorithm based on the two-dimensional discrete dipole approximation. Synthetic data from a finite-element based solver and experimental data from a microwave imaging system are used to reconstruct images and to validate the algorithm. The microwave measurement system consists of 16 monopole antennas immersed in a tank filled with lossy coupling liquid and a vector network analyzer. The low-profile antennas and lossy nature of system make the discrete dipole approximation an ideal forward solver in the image reconstructions. The results show that the algorithm can readily reconstruct a 2D plane of a cylindrical phantom. The proposed forward solver combined with the nodal adjoint method for computing the Jacobian matrix enables the algorithm to reconstruct an image within 6 seconds. This implementation provides a significant time savings and reduced memory requirements and is a dramatic improvement over previous implementations.

Comments on "Investigation of Histology Region Dielectric Measurements of Heterogeneous Tissues".

In the paper, "Investigation of histology region in dielectric measurements of heterogeneous tissues," by Porter and O'Halloran, the authors utilize a flexible phantom in a layered material dielectric property analysis to quantify the effective sensing volume of a coaxial dielectric probe. Ostensibly, this test has been used by others to characterize the region for which percent variation in the material composition in front of the probe corresponds to percent variation in the computed effective dielectric properties. By employing a compressible material, the authors fail to isolate features that are attributable solely to the probe, itself, and inadvertently incorporate confounding characteristics associated with the compressible nature of the material. The net effect is to exaggerate the probe's sensing volume which undermines conclusions drawn from the subsequent tissue dielectric property studies.

Electrical Characterization of Glycerin: Water Mixtures: Implications for Use as a Coupling Medium in Microwave Tomography.

We examine the broadband behavior of complex electrical properties of glycerin and water mixtures over the frequency range of 0.1 - 25.0 GHz, especially as they relate to using these liquids as coupling media for microwave tomographic imaging. Their combination is unique in that they are mutually miscible over the full range of concentrations which allows them to be tailored to dielectric property matching for biological tissues. While the resultant mixture properties are partially driven by differences in the inherent low frequency permittivity of each constituent, relaxation frequency shifts play a disproportionately larger role in increasing the permittivity dispersion while also dramatically increasing the effective conductivity over the frequency range of 1 to 3 GHz. For the full range of mixture ratios, the relaxation frequency shifts from 17.5 GHz for 0% glycerin to less than 0.1 GHz for 100% glycerin. Of particular interest is the fact that the conductivity stays above 1.0 S/m over the 1-3 GHz range for glycerin mixture ratios (70-90% glycerin) we use for microwave breast tomography. The high level of attenuation is critical for suppressing unwanted multipath signals. This paper presents a full characterization of these liquids along with a discussion of their benefits and limitations in the context of microwave tomography.

Open-Ended Coaxial Dielectric Probe Effective Penetration Depth Determination.

We have performed a series of experiments which demonstrate the effect of open-ended coaxial diameter on the depth of penetration. We used a two layer configuration of a liquid and movable cylindrical piece of either Teflon or acrylic. The technique accurately demonstrates the depth in a sample for which a given probe diameter provides a reasonable measure of the bulk dielectric properties for a heterogeneous volume. In addition we have developed a technique for determining the effective depth for a given probe diameter size. Using a set of simulations mimicking four 50 Ω coaxial cable diameters, we demonstrate that the penetration depth in both water and saline has a clear dependence on probe diameter but is remarkably uniform over frequency and with respect to the intervening liquid permittivity. Two different 50 Ω commercial probes were similarly tested and confirm these observations. This result has significant implications to a range of dielectric measurements, most notably in the area of tissue property studies.

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.

Surface scanning through a cylindrical tank of coupling fluid for clinical microwave breast imaging exams.

PURPOSE: Microwave tomographic image quality can be improved significantly with prior knowledge of the breast surface geometry. The authors have developed a novel laser scanning system capable of accurately recovering surface renderings of breast-shaped phantoms immersed within a cylindrical tank of coupling fluid which resides completely external to the tank (and the aqueous environment) and overcomes the challenges associated with the optical distortions caused by refraction from the air, tank wall, and liquid bath interfaces. METHODS: The scanner utilizes two laser line generators and a small CCD camera mounted concentrically on a rotating gantry about the microwave imaging tank. Various calibration methods were considered for optimizing the accuracy of the scanner in the presence of the optical distortions including traditional ray tracing and image registration approaches. In this paper, the authors describe the construction and operation of the laser scanner, compare the efficacy of several calibration methods-including analytical ray tracing and piecewise linear, polynomial, locally weighted mean, and thin-plate-spline (TPS) image registrations-and report outcomes from preliminary phantom experiments. RESULTS: The results show that errors in calibrating camera angles and position prevented analytical ray tracing from achieving submillimeter accuracy in the surface renderings obtained from our scanner configuration. Conversely, calibration by image registration reliably attained mean surface errors of less than 0.5 mm depending on the geometric complexity of the object scanned. While each of the image registration approaches outperformed the ray tracing strategy, the authors found global polynomial methods produced the best compromise between average surface error and scanner robustness. CONCLUSIONS: The laser scanning system provides a fast and accurate method of three dimensional surface capture in the aqueous environment commonly found in microwave breast imaging. Optical distortions imposed by the imaging tank and coupling bath diminished the effectiveness of the ray tracing approach; however, calibration through image registration techniques reliably produced scans of submillimeter accuracy. Tests of the system with breast-shaped phantoms demonstrated the successful implementation of the scanner for the intended application.

Motion-based microwave tomographic measurement device for three-dimensional coverage in a magnetic resonance system.

PURPOSE: We have developed a fully 3D data acquisition system for microwave breast imaging which can operate simultaneously inside a magnetic resonance imaging (MRI). MRI is used regularly for breast imaging to distinguish tumors from normal tissue. It generally has poor specificity unless used with a gadolinium contrast agent. Microwave imaging could fill this need because of the good endogenous tumor:normal tissue property contrast, especially in light of safety concerns for gadolinium. The antenna array consists of 16 monopole antennas positioned in a horizontal circle surrounding the breast which can then be moved vertically for 3D coverage of the breast. The tank system materials were chosen to minimize artifacts in the MR image within the specific shared imaging zone. The support rods are stainless steel, albeit positioned sufficiently far from the imaging target to have little effect. The mechanical motion parts are all 3D printed plastic. Unlike many conventional antennas, the monopoles consist of just the center conductor and insulator of the coaxial cable, making it one of the least possible metallic structures. METHODS: Data were acquired both inside and outside of the MR bore to confirm that the MR bore did not have adverse effects on the microwave imaging process. The imaging tank was filled with a mixture of glycerin and water to both provide a reasonable property match to the phantom and to highly attenuate the fields which also acted to suppress multi-path signals. Microwave images were reconstructed using our Gauss-Newton scheme combined with a log transformation for a more linear convergence. MR images were also acquired to assess the effects of the microwave tank structures on the imaging. RESULTS: The microwave measurement data were acquired in log magnitude and phase format at 200 MHz increments from 700-1900 MHz. Each antenna acted sequentially as a transmitter while the complement of 15 acted as a receiver. The single frequency images were reconstructed using a Gauss-Newton iterative technique with a standard log transformation to linearize the process. The data showed that the signal strengths were between 7-10 dB lower for the case when the array was inside the MRI versus when not. Notwithstanding, the image quality was still high because of the significant signal to noise ratio. The reconstructed images in both situations demonstrated good 3D object recovery of the vertically size and shaped varying object. The MR images were not adversely affected by the presence of antennas or feed structures. CONCLUSIONS: We have demonstrated that our technique can recover high-quality images of a 3D varying object within an MRI system. Compatibility issues have been addressed for both the microwave and MRI systems. The reduced SNR for the case operating in the MRI did not adversely affect the images. To the best of our knowledge, this is the first example of a microwave imaging system operating in an MRI with full 3D volumetric capability.

Viable Three-Dimensional Medical Microwave Tomography: Theory and Numerical Experiments.

Three-dimensional microwave tomography represents a potentially very important advance over 2D techniques because it eliminates associated approximations which may lead to more accurate images. However, with the significant increase in problem size, computational efficiency is critical to making 3D microwave imaging viable in practice. In this paper, we present two 3D image reconstruction methods utilizing 3D scalar and vector field modeling strategies, respectively. Finite element (FE) and finite-difference time-domain (FDTD) algorithms are used to model the electromagnetic field interactions in human tissue in 3D. Image reconstruction techniques previously developed for the 2D problem, such as the dual-mesh scheme, iterative block solver, and adjoint Jacobian method are extended directly to 3D reconstructions. Speed improvements achieved by setting an initial field distribution and utilizing an alternating-direction implicit (ADI) FDTD are explored for 3D vector field modeling. The proposed algorithms are tested with simulated data and correctly recovered the position, size and electrical properties of the target. The adjoint formulation and the FDTD method utilizing initial field estimates are found to be significantly more effective in reducing the computation time. Finally, these results also demonstrate that cross-plane measurements are critical for reconstructing 3D profiles of the target.

Magnetic resonance microwave absorption imaging: feasibility of signal detection.

PURPOSE: Magnetic resonance (MR) technique was used to detect small displacements induced by localized absorption of pulsed 434 MHz microwave power as a potential method for tumor detection. METHODS: Phase contrast subtraction was used to separate the phase change due to motion from thermoelastic expansion from other contributions to phase variation such as the bulk temperature rise of the medium and phase offsets from the MR scanner itself. A simple set of experiments was performed where the motion was constrained to be one dimensional which provided controls on the data acquisition and motion extraction procedures. Specifically, the MR-detected motion signal was isolated by altering the direction of the microwave-induced motion and sampling the response with motion encoding gradients in all three directions when the microwave power was turned on and turned off. RESULTS: Successful signal detection, as evidenced by the recording of a systematic alternating (zigzag) phase pattern, occurred only when the motion encoding was in parallel with either the vertical or horizontal direction of the microwave-induced motion on both 10 and 4 mm spatial scales. CONCLUSIONS: These results demonstrate, for the first time, that motion associated with thermoelastic expansion from the absorption of pulsed microwave power can be detected with MR.

The correlation of in vivo and ex vivo tissue dielectric properties to validate electromagnetic breast imaging: initial clinical experience.

Electromagnetic (EM) breast imaging provides low-cost, safe and potentially a more specific modality for cancer detection than conventional imaging systems. A primary difficulty in validating these EM imaging modalities is that the true dielectric property values of the particular breast being imaged are not readily available on an individual subject basis. Here, we describe our initial experience in seeking to correlate tomographic EM imaging studies with discrete point spectroscopy measurements of the dielectric properties of breast tissue. The protocol we have developed involves measurement of in vivo tissue properties during partial and full mastectomy procedures in the operating room (OR) followed by ex vivo tissue property recordings in the same locations in the excised tissue specimens in the pathology laboratory immediately after resection. We have successfully applied all of the elements of this validation protocol in a series of six women with cancer diagnoses. Conductivity and permittivity gauged from ex vivo samples over the frequency range 100 Hz-8.5 GHz are found to be similar to those reported in the literature. A decrease in both conductivity and permittivity is observed when these properties are gauged from ex vivo samples instead of in vivo. We present these results in addition to a case study demonstrating how discrete point spectroscopy measurements of the tissue can be correlated and used to validate EM imaging studies.

Application of Two-Dimensional Discrete Dipole Approximation in Simulating Electric Field of a Microwave Breast Imaging System.

The two-dimensional electric field distribution of the microwave imaging system is numerically simulated for a simplified breast tumour model. The proposed two-dimensional discrete dipole approximation (DDA) has the potential to improve computational speed compared to other numerical methods while retaining comparable accuracy. We have modeled the field distributions in COMSOL Multiphysics as baseline results to benchmark the DDA simulations. We have also investigated the adequate sampling size and the effect of inclusion size and property contrast on solution accuracy. In this way, we can utilize the 2D DDA as an alternative, fast and reliable forward solver for microwave tomography. From a mathematical perspective, the derivation of the 2D DDA and its application to microwave imaging is new and not previously implemented. The simulation results and the measurements show that the 2D DDA is a well-grounded forward solver for the specified microwave breast imaging system.

3-D Microwave Tomography Using the Soft Prior Regularization Technique: Evaluation in Anatomically Realistic MRI-Derived Numerical Breast Phantoms.

OBJECTIVE: Fusion of magnetic resonance imaging (MRI) breast images with microwave tomography is accomplished through a soft prior technique, which incorporates spatial information (from MRI), i.e., accurate boundary location of different regions of interest, into the regularization process of the microwave image reconstruction algorithm. METHODS: Numerical experiments were completed on a set of three-dimensional (3-D) breast geometries derived from MR breast data with different parenchymal densities, as well as a simulated tumor to evaluate the performance over a range of breast shapes, sizes, and property distributions. RESULTS: When the soft prior regularization technique was applied, both permittivity and conductivity relative root mean square error values decreased by more than 87% across all breast densities, except in two cases where the error decrease was only 55% and 78%. In addition, the incorporation of structural priors increased contrast between tumor and fibroglandular tissue by 59% in permittivity and 192% in conductivity. CONCLUSION: This study confirmed that the soft prior algorithm is robust in 3-D and can function successfully across a range of complex geometries and tissue property distributions. SIGNIFICANCE: This study demonstrates that our microwave tomography is capable of recovering accurate tissue property distributions when spatial information from MRI is incorporated through soft prior regularization.

Microwave thermal imaging of scanned focused ultrasound heating: phantom results.

We are developing a microwave tomographic imaging system capable of monitoring thermal distributions based on the temperature dependence of the recovered dielectric properties. The system has been coupled to a high intensity focused ultrasound (HIFU) therapy device which can be mechanically steered under computer control to generate arbitrarily shaped heating zones. Their integration takes advantage of the focusing capability of ultrasound for the therapy delivery and the isolation of the microwave imaging signal from the power deposition source to allow simultaneous treatment monitoring. We present several sets of phantom experiments involving different types of heating patterns that demonstrate the quality of both the spatial and temporal thermal imaging performance. This combined approach is adaptable to multiple anatomical sites and may have the potential to be developed into a viable alternative to current clinical temperature monitoring devices for HIFU, such magnetic resonance (MR) imaging.

Addressing Multipath Signal Corruption in Microwave Tomography and the Influence on System Design and Algorithm Development.

In developing a microwave tomography system, we started by examining the fundamental signal measurement challenges-i.e., how to interrogate the target while suppressing unwanted multi-path signals. Beginning with a lossy coupling bath to suppress unwanted surface waves, we have developed a robust and reliable system that is both simple and low profile. However, beyond the basic measurement configuration, the lossy coupling medium concept has also informed our choice of array antenna and imaging algorithms. The synergism of these concepts has produced a novel concept which is embodied in a system that has been successfully translated to the clinic.

Effects of the Plastic of the Realistic GeePS-L2S-Breast Phantom.

A breast phantom developed at the Supelec Institute was interrogated to study its suitability for microwave tomography measurements. A microwave measurement system based on 16 monopole antennas and a vector network analyzer was used to study how the S-parameters are influenced by insertion of the phantom. The phantom is a 3D-printed structure consisting of plastic shells that can be filled with tissue mimicking liquids. The phantom was filled with different liquids and tested with the measurement system to determine whether the plastic has any effects on the recovered images or not. Measurements of the phantom when it is filled with the same liquid as the surrounding coupling medium are of particular interest. In this case, the phantom plastic has a substantial effects on the measurements which ultimately detracts from the desired images.

Image Registration for Microwave Tomography of the Breast Using Priors From Nonsimultaneous Previous Magnetic Resonance Images.

Microwave imaging is a low-cost imaging method that has shown promise for breast imaging and, in particular, neoadjuvant chemotherapy monitoring. The early studies of microwave imaging in the therapy monitoring setting are encouraging. For the neoadjuvant therapy application, it would be desirable to achieve the most accurate possible characterization of the tissue properties. One method to achieve increased resolution and specificity in microwave imaging reconstruction is the use of a soft prior regularization. The objective of this study is to develop a method to use magnetic resonance (MR) images, taken in a different imaging configuration, as this soft prior. To enable the use of the MR images as a soft prior, it is necessary to register the MR images to the microwave imaging space. Registration fiducials were placed around the breast that are visible in both the MRI and with an optical scanner integrated into the microwave system. Utilizing these common registration locations, numerical algorithms have been developed to warp the original breast MR images into a geometry closely resembling that in which the breast is pendant in the microwave system.

Log transformation benefits parameter estimation in microwave tomographic imaging.

Microwave tomographic imaging falls under a broad category of nonlinear parameter estimation methods when a Gauss-Newton iterative reconstruction technique is used. A fundamental requirement in using these approaches is evaluating the appropriateness of the regression model. While there have been numerous investigations of regularization techniques to improve overall image quality, few, if any, studies have explored the underlying statistical properties of the model itself. The ordinary least squares (OLS) approach is used most often, but there are other options such as the weighted least squares (WLS), maximum likelihood (ML), and maximum a posteriori (MAP) that may be more appropriate. In addition, a number of variance stabilizing transformations can be applied to make the inversion intrinsically more linear. In this paper, a statistical analysis is performed of the properties of the residual errors from the reconstructed images utilizing actual measured data and it is demonstrated that the OLS algorithm with a log transformation (OLSlog) is clearly advantageous relative to the more commonly used OLS approach by itself. In addition, several high contrast imaging experiments are performed, which demonstrate that different subsets of data are emphasized in each method and may contribute to the overall image quality differences.

3-Point support mechanical steering system for high intensity focused ultrasound.

We have developed a simple approach for mechanically scanning a focused bowl ultrasound (US) transducer for either hyperthermia or tissue ablation therapies called the 3-point support (3PS) mechanical steering technique. The scanning involves translation of the required 3D motion of the ultrasound transducer to the more manageable linear movement of three support rods. It is a cost-effective alternative, especially compared with electronic scanning and other previous implementations of mechanically scanned systems. The 3PS approach is particularly well suited for integration with our microwave breast imaging technique--the combination of which could be an effective, low-cost thermal therapy/monitoring approach. The results show that the US focus can be moved laterally in a spiral pattern 3 cm below the surface in a gel phantom and that similar patterns can be moved to multiple locations within the phantom volume in succession. The feasibility of simultaneously acquiring microwave thermal images is also demonstrated.

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.