-
Benchmark problems for transcranial ultrasound simulation: Intercomparison of compressional wave models
Authors:
Jean-Francois Aubry,
Oscar Bates,
Christian Boehm,
Kim Butts Pauly,
Douglas Christensen,
Carlos Cueto,
Pierre Gelat,
Lluis Guasch,
Jiri Jaros,
Yun Jing,
Rebecca Jones,
Ningrui Li,
Patrick Marty,
Hazael Montanaro,
Esra Neufeld,
Samuel Pichardo,
Gianmarco Pinton,
Aki Pulkkinen,
Antonio Stanziola,
Axel Thielscher,
Bradley Treeby,
Elwin van 't Wout
Abstract:
Computational models of acoustic wave propagation are frequently used in transcranial ultrasound therapy, for example, to calculate the intracranial pressure field or to calculate phase delays to correct for skull distortions. To allow intercomparison between the different modeling tools and techniques used by the community, an international working group was convened to formulate a set of numeric…
▽ More
Computational models of acoustic wave propagation are frequently used in transcranial ultrasound therapy, for example, to calculate the intracranial pressure field or to calculate phase delays to correct for skull distortions. To allow intercomparison between the different modeling tools and techniques used by the community, an international working group was convened to formulate a set of numerical benchmarks. Here, these benchmarks are presented, along with intercomparison results. Nine different benchmarks of increasing geometric complexity are defined. These include a single-layer planar bone immersed in water, a multi-layer bone, and a whole skull. Two transducer configurations are considered (a focused bowl and a plane piston), giving a total of 18 permutations of the benchmarks. Eleven different modeling tools are used to compute the benchmark results. The models span a wide range of numerical techniques, including the finite-difference time-domain method, angular-spectrum method, pseudospectral method, boundary-element method, and spectral-element method. Good agreement is found between the models, particularly for the position, size, and magnitude of the acoustic focus within the skull. When comparing results for each model with every other model in a cross comparison, the median values for each benchmark for the difference in focal pressure and position are less than 10\% and 1 mm, respectively. The benchmark definitions, model results, and intercomparison codes are freely available to facilitate further comparisons.
△ Less
Submitted 9 February, 2022;
originally announced February 2022.
-
Computationally efficient full-waveform inversion of the brain using frequency-adaptive grids and lossy compression
Authors:
Letizia Protopapa,
Carlos Cueto
Abstract:
A tomographic technique called full-waveform inversion has recently shown promise as a fast, affordable, and safe modality to image the brain using ultrasound. However, its high computational cost and memory footprint currently limit its clinical applicability. Here, we address these challenges through a frequency-adaptive discretisation of the imaging domain and lossy compression techniques. Beca…
▽ More
A tomographic technique called full-waveform inversion has recently shown promise as a fast, affordable, and safe modality to image the brain using ultrasound. However, its high computational cost and memory footprint currently limit its clinical applicability. Here, we address these challenges through a frequency-adaptive discretisation of the imaging domain and lossy compression techniques. Because full-waveform inversion relies on the adjoint-state method, every iteration involves solving the wave equation over a discretised spatiotemporal grid and storing the numerical solution to calculate gradient updates. The computational cost depends on the grid size, which is controlled by the maximum frequency being modelled. Since the propagated frequency typically varies during the reconstruction, we reduce reconstruction time and memory use by allowing the grid size to change throughout the inversion. Moreover, we combine this approach with multiple lossy compression techniques that exploit the sparsity of the wavefield to further reduce its memory footprint. We explore applying these techniques in the spatial, wavelet, and wave atom domains. Numerical experiments using a human-head model show that our methods lead to a 30% reduction in reconstruction time and up to three orders of magnitude less memory, while negligibly affecting the accuracy of the reconstructions.
△ Less
Submitted 8 November, 2021;
originally announced November 2021.
-
Stride: a flexible platform for high-performance ultrasound computed tomography
Authors:
Carlos Cueto,
Oscar Bates,
George Strong,
Javier Cudeiro,
Fabio Luporini,
Oscar Calderon Agudo,
Gerard Gorman,
Lluis Guasch,
Meng-Xing Tang
Abstract:
Advanced ultrasound computed tomography techniques like full-waveform inversion are mathematically challenging and orders of magnitude more computationally expensive than conventional ultrasound imaging methods. This computational and algorithmic complexity, and a lack of open-source libraries in this field, represent a barrier preventing the generalised adoption of these techniques, slowing the p…
▽ More
Advanced ultrasound computed tomography techniques like full-waveform inversion are mathematically challenging and orders of magnitude more computationally expensive than conventional ultrasound imaging methods. This computational and algorithmic complexity, and a lack of open-source libraries in this field, represent a barrier preventing the generalised adoption of these techniques, slowing the pace of research and hindering reproducibility. Consequently, we have developed Stride, an open-source Python library for the solution of large-scale ultrasound tomography problems. On one hand, Stride provides high-level interfaces and tools for expressing the types of optimisation problems encountered in medical ultrasound tomography. On the other, these high-level abstractions seamlessly integrate with high-performance wave-equation solvers and with scalable parallelisation routines. The wave-equation solvers are generated automatically using Devito, a domain specific language, and the parallelisation routines are provided through the custom actor-based library Mosaic. Through a series of examples, we show how Stride can handle realistic tomographic problems, in 2D and 3D, providing intuitive and flexible interfaces that scale from a local multi-processing environment to a multi-node high-performance cluster.
△ Less
Submitted 18 May, 2022; v1 submitted 7 October, 2021;
originally announced October 2021.