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.

3D-printed gear system for antenna motion in an MR environment: initial phantom imaging experiments.

We have developed a microwave imaging device for breast cancer imaging that can be used concurrently inside an MR imaging system. The microwave measurement system is comprised of a horizontal array of 16 monopole antennas that can be moved vertically for full 3D coverage of the breast. All compatibility issues have been addressed. The motion is achieved using a novel 3D printed gearing device. Initial results demonstrate that the system is capable of accurately recovering the size, shape, location and properties of a 3D shape varying object. This is a critical step towards clinical microwave breast imaging in the MR.

Impact of Skin on Microwave Tomography in the Lossy Coupling Medium.

In microwave imaging, the effects of skin on recovering property distributions of tissue underneath the surface may be significant because it has high dielectric contrast with subcutaneous fat, which inevitably causes significant signal reflections. While the thickness of skin, especially relative to the wavelengths in use, would presumably have minor effects, it can introduce practical difficulties, for instance, in reflection-based imaging techniques, where the impact of the skin is large-often as high as two orders of magnitude greater than that of signals from underlying tumors in the breast imaging setting. However, in tomography cases utilizing transmission-based measurement data and lossy coupling materials, the situation is considerably different. Accurately implementing a skin layer for numerical modeling purposes is challenging because of the need to discretize the size and shape of the skin without increasing computational overhead substantially. In this paper, we assess the effects of the skin on field solutions in a realistic 3D model of a human breast. We demonstrate that the small changes in transmission field values introduced by including the skin cause minor differences in reconstructed images.

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.

Open-Ended Transmission Coaxial Probes for Sarcopenia Assessment.

We developed a handheld, side-by-side transmission-based probe for interrogating tissue to diagnose sarcopenia-a condition largely characterized by muscle loss and replacement by fat. While commercial microwave reflection-based probes exist, they can only be used in a lab for a variety of applications. The penetration depth of these probes is only in the order of 0.3 mm, which does not even traverse the skin layer, and minor motion of the coaxial feedlines can completely dismantle the calibration. Our device builds primarily on the transmission-based concept that allows for substantially greater signal penetration depth operating over a very broad bandwidth. Additional features were integrated to further improve the penetration, optimize the geometry for a more focused planar excitation, and juxtapose the coaxial apertures for more controlled interrogation. The larger coaxial apertures further increased the penetration depth while retaining the broadband performance. Three-dimensional printing technology made it possible for the apertures to be compressed into ellipses for interrogation in a near-planar geometry. Finally, fixed side-by-side positioning provided repeatable and reliable performance. The probes were also not susceptible to multipath signal corruption due to the close proximity of the transmitting and receiving apertures. The new concept worked from 100 MHz to over 8 GHz and could sense property changes as deep as 2-3 cm. While the signal changes due to deeper feature aberrations were more subtle than for signals emanating from the skin and subcutaneous fat layers, the large property contrast between muscle and fat is a sarcopenic indication that helps to distinguish even the deepest objects. This device has the potential to provide needed specificity information about the relevant underlying tissue.

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.

Expansion of the Nodal-Adjoint Method for Simple and Efficient Computation of the 2D Tomographic Imaging Jacobian Matrix.

This paper focuses on the construction of the Jacobian matrix required in tomographic reconstruction algorithms. In microwave tomography, computing the forward solutions during the iterative reconstruction process impacts the accuracy and computational efficiency. Towards this end, we have applied the discrete dipole approximation for the forward solutions with significant time savings. However, while we have discovered that the imaging problem configuration can dramatically impact the computation time required for the forward solver, it can be equally beneficial in constructing the Jacobian matrix calculated in iterative image reconstruction algorithms. Key to this implementation, we propose to use the same simulation grid for both the forward and imaging domain discretizations for the discrete dipole approximation solutions and report in detail the theoretical aspects for this localization. In this way, the computational cost of the nodal adjoint method decreases by several orders of magnitude. Our investigations show that this expansion is a significant enhancement compared to previous implementations and results in a rapid calculation of the Jacobian matrix with a high level of accuracy. The discrete dipole approximation and the newly efficient Jacobian matrices are effectively implemented to produce quantitative images of the simplified breast phantom from the microwave imaging system.

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.

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.

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.

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.

Two-step inversion with a logarithmic transformation for microwave breast imaging.

PURPOSE: The authors have developed a new two-step microwave tomographic image reconstruction process specifically designed to incorporate logarithmic transformed microwave imaging algorithms as a means of significantly improving spatial resolution and target property recovery. Log transform eliminates the need for a priori information, but spatial filtering often integrated as part of the regularization required to stabilize image recovery, generally smooths image features and reduces object definition. The new implementation begins with this smoothed image as the first step, but then utilizes it as the starting estimate for a second step which continues the iterative process with a standard weighted Euclidean distance regularization. The penalty term of the latter restricts the new image to a multi-dimensional location close to the original but allows the algorithm to optimize the image without excessive smoothing. METHODS: The overall approach is based on a Gauss-Newton iterative scheme which incorporates a log transformation as a way of making the reconstruction more linear. It has been shown to be robust and not require a priori information as a condition for convergence, but does produce somewhat smoothed images as a result of associated regularization. The new two-step process utilizes the previous technique to generate a smoothed initial estimate and then uses the same reconstruction process with a weighted Euclidean distance penalty term. A simple and repeatable method has been implemented to determine the weighting factor without significant computational burden. The reconstructions are assessed according to conventional parameter estimation metrics. RESULTS: We apply the approach to phantom experiments using large, high contrast canonical shapes followed by a set of images recovered from an actual patient exam. The image improvements are substantial in regards to improved property recovery and feature delineation without inducing unwanted artifacts. Analysis of the residual vector after the reconstruction process further emphasizes that the minimization criterion is efficient with minimal biases. CONCLUSIONS: The outcome is a novel synergism of an established stable reconstruction algorithm with a conventional regularization technique. It maintains the ability to recover high quality microwave tomographic images without the bias of a priori information while substantially improving image quality. The results are confirmed on both phantom experiments and patient exams.

A Coaxial Dielectric Probe Technique for Distinguishing Tooth Enamel from Dental Resin.

For purposes of biodosimetry in the event of a large scale radiation disaster, one major and very promising point-of contact device is assessing dose using tooth enamel. This technique utilizes the capabilities of electron paramagnetic resonance to measure free radicals and other unpaired electron species, and the fact that the deposition of energy from ionizing radiation produces free radicals in most materials. An important stipulation for this strategy is that the measurements, need to be performed on a central incisor that is basically intact, i.e. which has an area of enamel surface that is as large as the probing tip of the resonator that is without decay or restorative care that replaces the enamel. Therefore, an important consideration is how to quickly assess whether the tooth has sufficient enamel to be measured for dose and whether there is resin present on the tooth being measured and to be able to characterize the amount of surface that is impacted. While there is a relatively small commercially available dielectric probe which could be used in this context, it has several disadvantages for the intended use. Therefore, a smaller, 1.19mm diameter 50 ohm, open-ended, coaxial dielectric probe has been developed as an alternative. The performance of the custom probe was validated against measurement results of known standards. Measurements were taken of multiple teeth enamel and dental resin samples using both probes. While the probe contact with the teeth samples was imperfect and added to measurement variability, the inherent dielectric contrast between the enamel and resin was sufficient that the probe measurements could be used as a robust means of distinguishing the two material types. The smaller diameter probe produced markedly more definitive results in terms of distinguishing the two materials.

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.

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.

Comparison of no-prior and soft-prior regularization in biomedical microwave imaging.

Microwave imaging for medical applications is attractive because the range of dielectric properties of different soft tissues can be substantial. Breast cancer detection and monitoring of treatment response are areas where this technology could be important because of the contrast between normal and malignant tissue. Unfortunately, the technique is unable to achieve the high spatial resolution at depth in tissue which is available from other conventional modalities such as x-ray computed tomography (CT) or magnetic resonance imaging (MRI). We have incorporated a soft-prior regularization strategy within our microwave reconstruction algorithm and compared it with the images obtained with traditional no-prior (Levenberg-Marquardt) regularization. Initial simulation and phantom results show a significant improvement of the recovered electrical properties. Specifically, errors in the microwave property estimates were improved by as much as 95%. The effects of a false-inclusion region were also evaluated and the results show that a small residual property bias of 6% in permittivity and 15% in conductivity can occur that does not otherwise degrade the property recovery accuracy of inclusions that actually exist. The work sets the stage for integrating microwave imaging with MR for improved resolution and functional imaging of the breast in the future.

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.