Acoustoelastic full-waveform inversion for transcranial ultrasound computed tomography

Marty, Patrick; Boehm, Christian; Fichtner, Andreas

doi:10.1117/12.2581029

Cited by 20 publications

(17 citation statements)

References 33 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The operations summarized are those of the largest data sizes in Figures 2 and 3; single precision a • x + y (1D SAXPY), one-dimensional central differencing using a 3-point stencil (1D CD), two-dimensional element wise function (2D EW), two-dimensional Laplacian using a 5-point stencil (2D LA) and twodimensional Laplacian using a 9-point stencil (2D LA9p). man-made structures [12,13], and in medical tomography when imaging the human body [14,15].…”

Section: Accelerating Existing Codes: Elastic Wave Propagation On The...mentioning

confidence: 99%

Seamless GPU acceleration for C++ based physics with the Metal Shading Language on Apple's M series unified chips

Gebraad¹,

Fichtner²

2022

Preprint

Self Cite

View full text Add to dashboard Cite

The M series of chips produced by Apple have proven a capable and power-efficient alternative to mainstream Intel and AMD x86 processors for everyday tasks. Additionally, the unified design integrating the central processing and graphics processing unit, have allowed these M series chips to excel at many tasks with heavy graphical requirements without the need for a discrete graphical processing unit (GPU), and in some cases even outperforming discrete GPUs.In this work, we show how the M series chips can be leveraged using the Metal Shading Language (MSL) to accelerate typical array operations in C++. More importantly, we show how the usage of MSL avoids the typical complexity of CUDA or OpenACC memory management, by allowing the central processing unit (CPU) and GPU to work in unified memory. We demonstrate how performant the M series chips are on standard onedimensional and two-dimensional array operations such as array addition, SAXPY and finite difference stencils, with respect to serial and OpenMP accelerated CPU code. The reduced complexity of implementing MSL also allows us to accelerate an existing elastic wave equation solver (originally based on OpenMP accelerated C++) using MSL, with minimal effort, while retaining all CPU and OpenMP functionality.The resulting performance gain of simulating the wave equation is near an order of magnitude for specific settings. This gain attained from using MSL is similar to other GPU-accelerated wave-propagation codes with respect to their CPU variants, but does not come at much increased programming complexity that prohibits the typical scientific programmer to leverage these accelerators. This result shows how unified processing units can be a valuable tool to seismologists and computational scientists in general, lowering the bar to writing performant codes that leverage modern GPUs.

show abstract

Section: Accelerating Existing Codes: Elastic Wave Propagation On The...mentioning

confidence: 99%

Seamless GPU acceleration for C++ based physics with the Metal Shading Language on Apple's M series unified chips

Gebraad¹,

Fichtner²

2022

Preprint

Self Cite

View full text Add to dashboard Cite

show abstract

“…In this paper, we show that the conditions for elastodynamic Green's function retrieval are significantly better when an elastic volume can be accessed from two sides, as is the case in particular laboratory experiments [12], non-destructive testing [13,14], brain imaging [15,16], transcranial ultrasound focusing [17,18], transcranial photoacoustics [19,20], or by using auxiliary downhole receivers in seismic data acquisition [21,22]. Although the underlying representations of our work could be extended to account for lateral variations [11], the presence of a free surface [23,24] and intrinsic attenuation [25,26], we restrict ourselves to a layered lossless medium for simplicity.…”

Section: Introductionmentioning

confidence: 97%

Marchenko Green's Function Retrieval in Layered Elastic Media from Two-Sided Reflection and Transmission Data

Neut¹,

Brackenhoff²,

Meles³

et al. 2022

Preprint

View full text Add to dashboard Cite

By solving a Marchenko equation, Green's functions at an arbitrary (inner) depth level inside an unknown elastic layered medium can be retrieved from single-sided reflection data, which are to be collected at the top of the medium. So far, an exact solution could only be obtained if the medium obeys stringent monotonicity conditions and if all forward-scattered (non-converted and converted) transmissions between the acquisition level and the inner depth level are known a-priori. We introduce an alternative Marchenko equation by revising the window operators that are applied in its derivation. We also introduce an auxiliary equation for transmission data, which are to be collected at the bottom of the medium, and a coupled equation, which is based on both reflection and transmission data. We show that the joint system of the Marchenko equation, the auxiliary equation and the coupled equation can be succesfully inverted when broadband reflection and transmission data are available. This results in a novel methodology for elastodynamic Green's function retrieval from two-sided data. Apart from these data, our approach requires P- and S-wave transmission times between the inner depth level and the top of the medium, as well as two angle-dependent amplitude scaling factors, which can be estimated from the data by enforcing energy conservation.

show abstract

“…Such possibilities in real NDT experiments are presented in [14][15][16] assuming the scalar wave model. Some publications are concerned with solving inverse problems of tomographic ultrasonic diagnostics in a vector model, which supports transverse, surface, and other types of waves propagating in the medium in addition to the longitudinal wave [17][18][19].…”

Section: Introductionmentioning

confidence: 99%

On mathematical problems of two-coefficient inverse problems of ultrasonic tomography

Goncharsky,

Romanov,

Seryozhnikov

2024

Inverse Problems

View full text Add to dashboard Cite

This paper proves the theorem of uniqueness for the solution of a coefficient inverse problem for the wave equation in R^2 with two unknown coefficients: speed of sound and absorption. The original nonlinear coefficient inverse problem is reduced to an equivalent system of two uniquely solvable linear integral equations of the first kind with respect to the sound speed and absorption coefficients. Estimates are made, substantiating the multistage method for two unknown coefficients. These estimates show that given sufficiently low frequencies and small inhomogeneities, the residual functional for the nonlinear inverse problem approaches a convex one. This solution method for nonlinear coefficient inverse problems is not linked to the limit approach as frequency tends to zero, but assumes solving the inverse problem using sufficiently low, but not zero, frequencies at the first stage. For small inhomogeneities that are typical, for instance, for medical tasks, carrying out real experiments at such frequencies does not present major difficulties. The capabilities of the method are demonstrated on a model inverse problem with unknown sound speed and absorption coefficients. The method effectively solves the nonlinear problem with parameter values typical for tomographic diagnostics of soft tissues in medicine. A resolution of approximately 2 mm was achieved using an average sounding pulse wavelength of 5 mm.

show abstract

Acoustoelastic full-waveform inversion for transcranial ultrasound computed tomography

Cited by 20 publications

References 33 publications

Seamless GPU acceleration for C++ based physics with the Metal Shading Language on Apple's M series unified chips

Seamless GPU acceleration for C++ based physics with the Metal Shading Language on Apple's M series unified chips

Marchenko Green's Function Retrieval in Layered Elastic Media from Two-Sided Reflection and Transmission Data

On mathematical problems of two-coefficient inverse problems of ultrasonic tomography

Contact Info

Product

Resources

About