A new, freely available third party MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields is described. The toolbox, named k-Wave, is designed to make realistic photoacoustic modeling simple and fast. The forward simulations are based on a k-space pseudo-spectral time domain solution to coupled first-order acoustic equations for homogeneous or heterogeneous media in one, two, and three dimensions. The simulation functions can additionally be used as a flexible time reversal image reconstruction algorithm for an arbitrarily shaped measurement surface. A one-step image reconstruction algorithm for a planar detector geometry based on the fast Fourier transform (FFT) is also included. The architecture and use of the toolbox are described, and several novel modeling examples are given. First, the use of data interpolation is shown to considerably improve time reversal reconstructions when the measurement surface has only a sparse array of detector points. Second, by comparison with one-step, FFT-based reconstruction, time reversal is shown to be sufficiently general that it can also be used for finite-sized planar measurement surfaces. Last, the optimization of computational speed is demonstrated through parallel execution using a graphics processing unit.

/pdf/k-wave-matlab-toolbox-for-the-simulation-and-reconstruction-3uo74cqwwh.pdf

k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields

A closed-loop control system is disclosed for use with an electrosurgical generator that generates electrosurgical energy. The closed loop control system includes a user interface for allowing a user to select at least one pre-surgical parameter, such as the type of surgical instrument operatively connected to the generator, the type of tissue and the desired surgical effect. A sensor module is also included for continually sensing at least one of electrical and physical properties proximate a surgical site and generating at least one signal relating thereto. The system also includes a control module for continually receiving the at least one selected pre-surgical parameter from the user interface and each of the signals from the sensor module, and processing each of the signals in accordance with the at least one pre-surgical parameter using at least one of a microprocessor, computer algorithm and a mapping. The control module generates at least one corresponding control signal relating to each signal from the sensor module, and relays the control signal to the electrosurgical generator for controlling the same. A method is also disclosed for performing an electrosurgical procedure at a surgical site on a patient. The method includes the steps of: applying at least one electrical pulse to the surgical site; continually sensing electrical and physical properties proximate the surgical site; and varying pulse parameters of individual pulses of the at least one pulse in accordance with the continually-sensed properties.

Method and system for controlling output of RF medical generator

Numerical Techniques in Electromagnetics is designed to show the reader how to pose, numerically analyze, and solve electromagnetic (EM) problems. It gives them the ability to expand their problem-solving skills using a variety of available numerical methods. Topics covered include fundamental concepts in EM; numerical methods; finite difference methods; variational methods, including moment methods and finite element methods; transmission-line matrix or modeling (TLM); and Monte Carlo methods. The simplicity of presentation of topics throughout the book makes this an ideal text for teaching or self-study by senior undergraduates, graduate students, and practicing engineers.

Numerical Techniques in Electromagnetics, Second Edition

The CUDA programming model provides a straightforward means of describing inherently parallel computations, and NVIDIA's Tesla GPU architecture delivers high computational throughput on massively parallel problems. This article surveys experiences gained in applying CUDA to a diverse set of problems and the parallel speedups over sequential codes running on traditional CPU architectures attained by executing key computations on the GPU.

Parallel Computing Experiences with CUDA

After 20 years of innovation in techniques that specifically image the biomechanical properties of tissue, the evolution of elastographic imaging can be viewed from its infancy, through a proliferation of approaches to the problem to incorporation on research and then clinical imaging platforms. Ultimately this activity has culminated in clinical trials and improved care for patients. This remarkable progression represents a leading example of translational research that begins with fundamentals of science and engineering and progresses to needed improvements in diagnostic and monitoring capabilities applied to major categories of disease, surgery and interventional procedures. This review summarizes the fundamental principles, the timeline of developments in major categories of elastographic imaging, and concludes with recent results from clinical trials and forward-looking issues.

https://iopscience.iop.org/article/10.1088/0031-9155/56/1/R01/pdf

Imaging the elastic properties of tissue: the 20 year perspective.

An apparatus and method for rapid real time imaging with wavefield energy by inverse scattering using a C.P.U programmed to process data derived from wavefield energy that has been transmitted and scattered by an object so as to reconstruct a wavefield image of the object. Electronic signals are propagated and are transduced into wavefield energy waves which in turn are propagated toward the object. Detector means detect the wavefield energy waves scattered by the object. The detected wavefield energy waves are then electronically processed and input into a high-speed digital computer which may comprise a C.P.U. and/or a C.P.U in combination with an array or parallel processor. Data is also prepared and input to the computer representing the incident field and the computer then reconstructs a high-quality image of the object having high spacial resolution and including actual internal viscous and elastic properties of the object through the use of new inverse scattering techniques used in the data processing steps. The media in which the object is embedded may be fluid or solid, homogeneous, or layered (such as stratigraphic layering, or ocean velocity layers, or layering of composites in nondestructive imaging applications), or may consist of porous material (either sedimentary deposits or composites in nondestructive testing).

Apparatus and method for imaging with wavefields using inverse scattering techniques

A transmission wave field imaging method, comprising the transmission of an incident wave field into an object, the incident wave field propagating into the object and, at least, partially scattering. Also includes the measuring of a wave field transmitted, at least in part, through an object to obtain a measured wave field, the measured wave field based, in part, on the incident wave field and the object. Additionally, the processing of the measured wave field utilizing a parabolic approximation reconstruction algorithm to generate an image data set representing at least one image of the object.

Apparatus and method for imaging objects with wavefields

Although there are acceptable methods for calculating whole body electromagnetic absorption, no completely acceptable method for calculating the local specific absorption rate (SAR) at points within the body has been developed. Frequency domain methods, such as the method of moments (MoM) have achieved some success; however, MoM requires computer storage on the order of (3N) 2 and computation time on the order of (3N) 3 where N is the number of cells. The finite-difference time-domain (FDTD) method has been employed extensively in calculating the scattering of metallic objects, and recently is seeing some use in calculating the interaction of EM fields with complex, lossy dielectric bodies. Since the FDTD method has storage and time requirements proportional to N, it presents an attractive alternative to calculating SAR distribution in large bodies. This paper describes the FDTD method and evaluates it by comparing its results to analytic solutions in two and three dimensions. The utility of the FDTD method is demonstrated by a 3D scan of the human torso. The results obtained demonstrate that the FDTD method is capable of calculating internal SAR distribution with acceptable accuracy. With the availability of supercomputers, such as the CRAY II, the calculation of SAR distribution in a man model of 50 000 cells (1.27 cm per cell) appears to be feasible.

Use of the Finite-Difference Time-Domain Method in Calculating EM Absorption in Human Tissues

Berenger's perfectly matched layer (PML) absorbing boundary condition for electromagnetic (EM) waves is derived to absorb 2-D and 3-D acoustic waves in finite difference time domain (FDTD) simulation of acoustic wave propagation and scattering. A PML medium suitable for acoustic waves is constructed. Plane wave propagation in the PML medium is solved for both 2-D and 3-D cases and explicit FDTD boundary conditions are derived. The equations show that a matched PML medium is a perfect simulation of free space in that a plane wave does not change its direction of propagation or its speed when it propagates from free space into a matched PML medium. FDTD simulation of a pulsed point source propagating in two dimensions is carried out to test the performance of the PML boundary for acoustic waves. Results show that an eight layer PML boundary condition reduces the reflected error 40 dB over Mur's second order boundary condition.

Formulation and validation of Berenger's PML absorbing boundary for the FDTD simulation of acoustic scattering

Recent published results in inverse scattering generally show the difficulty in dealing with moderate to high contrast inhomogeneities when employing linearized or iteratively linearized algorithms (e.g., distorted Born iterative method). This paper presents a fully nonlinear algorithm utilizing full wave field data, that results in ultrasound computed tomographic images from a laboratory breast scanner, and shows several such unique images from volunteer subjects. The forward problem, data collection process and inverse scattering algorithm used are discussed. A functional that represents the "best fit" between predicted and measured data is minimized, and therefore requires a very fast forward problem solver, Jacobian calculation, and gradient estimation, all of which are described. The data collection device is described. The algorithm and device yield quantitative estimates of human breast tissue in vivo. Several high resolution images, measuring ∼150 by 150 wavelengths, obtained from the 2D inverse scattering algorithms, using data collected from a first prototype, are shown and discussed. The quantitative values are compared with previous published work.

/pdf/non-linear-inverse-scattering-high-resolution-quantitative-2xx0a87zor.pdf

David T. Borup

Papers

Apparatus and method for imaging with wavefields using inverse scattering techniques

Apparatus and method for imaging objects with wavefields

Use of the Finite-Difference Time-Domain Method in Calculating EM Absorption in Human Tissues

Formulation and validation of Berenger's PML absorbing boundary for the FDTD simulation of acoustic scattering

Non-linear inverse scattering: high resolution quantitative breast tissue tomography.