PubMed Health⌕ Search

Biomedical subjects

P M Meaney

Publications and source records attributed to P M Meaney.

10 recordsLinked to original sources

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.

Algorithms↗

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.

Algorithms↗

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.

Animals↗

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.

Computer Simulation↗

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.

Diagnostic Imaging↗

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.

Algorithms↗

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.

Algorithms↗

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.

Algorithms↗

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.

Algorithms↗

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.

Biophysical Phenomena↗