pith. sign in

arxiv: 2408.11471 · v2 · submitted 2024-08-21 · ⚛️ physics.med-ph

Learning-based Linear Inversion for Quantitative Pulse-Echo Speed-of-Sound Imaging

Pith reviewed 2026-05-23 21:34 UTC · model grok-4.3

classification ⚛️ physics.med-ph
keywords speed of soundultrasound tomographyecho modelinear inversionlearning-basedbias reductionquantitative imagingliver imaging
0
0 comments X

The pith

Training a linear operator on random tissue models reduces biases in speed-of-sound maps from ultrasound echoes.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that a linear operator trained to minimize speed-of-sound errors on average across many random tissue models yields maps with strongly reduced biases relative to those obtained with spatial gradient regularization. A sympathetic reader would care because such biases have been observed to depend on tissue layer geometry and can affect quantitative accuracy in applications such as liver imaging. The learned operator can be applied directly to echo-shift measurements or used to correct maps already produced by gradient regularization, with the direct route showing slightly better performance while the correction route is computationally lighter. Residual biases persist because the underlying physical mapping is only partially linear.

Core claim

By training a linear operator to minimize SoS errors on average over a large number of random tissue models that sample the distribution of geometries and SoS values expected in vivo, biases in the resulting speed-of-sound maps are strongly reduced. The operator can be applied either directly to echo-shift data or to maps estimated with gradient regularization. Residual biases that remain are the result of a partial non-linearity in the actual physical problem. The approach transfers from simulation to real ultrasound data, as confirmed by experimental phantom results.

What carries the argument

The trained linear operator that maps echo-shift data to SoS maps by minimizing average reconstruction errors over simulated tissue models.

If this is right

  • Biases are strongly reduced in an extensive simulation study on liver imaging.
  • Direct application to echo-shift data slightly outperforms post-processing of gradient-regularized maps.
  • The method remains effective when transferred to real ultrasound data in phantom experiments.
  • Residual biases are attributable to partial non-linearity in the physical forward model.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If real anatomy deviates systematically from the random models used in training, the operator may introduce new biases on clinical scans.
  • The same training strategy could be tested on other linearized inverse problems in pulse-echo imaging where geometry-dependent regularization errors appear.
  • Incorporating a small number of nonlinear forward simulations into the training loss might further shrink the residual bias floor.

Load-bearing premise

The distribution of random tissue models used for training accurately samples the geometries and speed-of-sound values expected in living tissue.

What would settle it

Applying the trained operator to echo-shift data from tissue geometries or speed-of-sound values outside the training distribution and observing that bias reduction disappears would show the sampling assumption does not hold.

Figures

Figures reproduced from arXiv: 2408.11471 by (2) University of Twente, (3) Ecole polytechnique f\'ed\'erale de Lausanne), Jules Blom (2), Michael Jaeger (1) ((1) University of Bern, Naiara Korta Martiartu (3), Parisa Salemi Yolgunlu (1).

Figure 1
Figure 1. Figure 1: Demonstration of optimum regularization parameters and minimum loss performance for GR on 10000 test [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Demonstration of liver-SoS biases in GR-SoS maps, for three example tissue models (I) to (III). (a) Target [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Loss performance of LR when trained for correcting GR-SD maps that were obtained with the linear model. [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Distribution of liver-SoS biases when using LR for correcting GR-SD maps that were obtained with the linear [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: LR-SoS maps when LR is optimized on N = 50000 training examples for correcting GR-SD maps that were obtained with the linear model. LR-SoS maps (b) are shown vs. the corresponding target SoS maps (a) for the three tissues models (I) to (III) [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Loss performance of LR when trained for correcting GR-SD maps that were obtained with the full-wave [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Distribution of liver-SoS biases when using LR for correcting GR-SD maps that were obtained with the [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: LR-SoS maps when LR is optimized on N = 50000 training examples for correcting GR-SD maps that were obtained with the full-wave model. LR-SoS maps (b) are shown vs. the corresponding target SoS maps (a) for the three tissues models (I) to (III). width is only slightly reduced above N = 1500. The STD values are 9.0 m/s (N = 50), 6.5 m/s (N = 1500), 6.1 m/s (N = 10000), and 6.1 m/s (N = 50000). The larger di… view at source ↗
Figure 9
Figure 9. Figure 9: Loss performance of LR when trained for correcting GR-SD maps that were obtained with the full-wave [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Distribution of liver-SoS biases when using LR for correcting GR-SD maps that were obtained with the [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: LR-SoS maps when LR is optimized on N = 50000 training examples for correcting GR-SD maps that were obtained with the full-wave model and high-resolution target SD maps. LR-SoS maps (b) are shown vs. the corresponding high-resolution target SoS maps (a) for the three tissues models (I) to (III). ones observed for SD-based LR, the values for N = 10000 and N = 50000 are smaller. This indicates that ES-based… view at source ↗
Figure 12
Figure 12. Figure 12: Loss performance of LR when trained to invert SD from ES maps that were obtained with the full-wave [PITH_FULL_IMAGE:figures/full_fig_p012_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Distribution of liver-SoS biases of LR when trained to invert SD from ES maps that were obtained with the [PITH_FULL_IMAGE:figures/full_fig_p012_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: LR-SoS maps when LR is optimized on N = 50000 training examples to directly derive SD from ES maps that were obtained with the full-wave model. LR-SoS maps (b) are shown vs. the corresponding target SoS maps (a) for the three tissues models (I) to (III). process is taken into account (full-wave model). For ES following the linear model, the LR approach reconstructed liver SoS values with a precision of 1.… view at source ↗
Figure 15
Figure 15. Figure 15: Experimental results of GR and LR approach in three physical phantoms with different layer geometry (I) to [PITH_FULL_IMAGE:figures/full_fig_p014_15.png] view at source ↗
read the original abstract

Computed ultrasound tomography in echo mode generates maps of tissue speed of sound (SoS) from the shift of echoes when detected under varying steering angles. It solves a linearized inverse problem that requires regularization to complement the echo shift data with a priori constraints. Spatial gradient regularization has been used to enforce smooth solutions, but SoS estimates were found to be biased depending on tissue layer geometry. Here, we propose to train a linear operator to minimize SoS errors on average over a large number of random tissue models that sample the distribution of geometries and SoS values expected in vivo. In an extensive simulation study on liver imaging, we demonstrate that biases are strongly reduced, with residual biases being the result of a partial non-linearity in the actual physical problem. This approach can either be applied directly to echo-shift data or to the SoS maps estimated with gradient regularization, where the former shows slightly better performance, but the latter is computationally more efficient. Experimental phantom results confirm the transferability of our results to real ultrasound data.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes training a linear operator on a large ensemble of random tissue models to minimize average speed-of-sound (SoS) estimation errors for pulse-echo ultrasound tomography. This data-driven regularization is shown in liver-imaging simulations to strongly reduce geometry-dependent biases compared with spatial-gradient regularization; residual errors are attributed to partial nonlinearity of the forward problem. The learned operator can be applied either directly to echo-shift measurements or as a post-correction to gradient-regularized maps, with the former yielding slightly better accuracy. Phantom experiments are presented to demonstrate transferability from simulation to real data.

Significance. If the training distribution accurately represents in-vivo liver anatomy and SoS statistics, the approach offers a computationally attractive route to bias-reduced quantitative SoS maps without requiring nonlinear solvers or additional hardware. The simulation study and phantom validation constitute concrete evidence of bias reduction; the explicit attribution of residuals to nonlinearity is a useful diagnostic insight.

major comments (2)
  1. [Abstract] Abstract: the assertion that the random tissue models 'sample the distribution of geometries and SoS values expected in vivo' is load-bearing for the generalization claim, yet no quantitative validation (e.g., comparison of first- and second-order moments, layer-thickness histograms, or spatial correlation lengths) against real liver data is supplied.
  2. [Simulation study (inferred from abstract)] The central claim that residual biases arise solely from partial nonlinearity would be strengthened by an explicit quantification (e.g., a table or figure showing the magnitude of the nonlinearity residual versus the learned-operator residual) rather than a qualitative statement.
minor comments (2)
  1. The abstract would benefit from stating the number of training models, the SoS range, and the precise error metric (bias, RMSE, etc.) used to train and evaluate the operator.
  2. Clarify whether the learned operator is a single matrix or a set of angle-dependent matrices, and how its application scales with the number of steering angles.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive comments. We respond to each major comment below, indicating planned revisions where appropriate.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the assertion that the random tissue models 'sample the distribution of geometries and SoS values expected in vivo' is load-bearing for the generalization claim, yet no quantitative validation (e.g., comparison of first- and second-order moments, layer-thickness histograms, or spatial correlation lengths) against real liver data is supplied.

    Authors: We agree the abstract phrasing is strong. The random models were constructed from literature-reported ranges for liver SoS (typically 1480-1620 m/s) and abdominal layer thicknesses, but the manuscript contains no explicit statistical comparison (moments, histograms, or correlation lengths) to in-vivo datasets. In revision we will change the abstract wording to 'random tissue models based on literature-reported distributions of geometries and SoS values' and add a short methods paragraph describing the chosen parameter ranges and their sources. revision: yes

  2. Referee: [Simulation study (inferred from abstract)] The central claim that residual biases arise solely from partial nonlinearity would be strengthened by an explicit quantification (e.g., a table or figure showing the magnitude of the nonlinearity residual versus the learned-operator residual) rather than a qualitative statement.

    Authors: The manuscript states that residuals remain because the learned operator is linear while the underlying echo-shift physics is partially nonlinear. We will strengthen this by adding a quantitative comparison: on the held-out simulation test set we will report the L2 norm of the difference between the linear forward operator and a nonlinear reference simulation, and directly compare its magnitude to the residual SoS error after the learned operator is applied. This will appear as a new panel in an existing figure or a short supplementary table. revision: yes

Circularity Check

0 steps flagged

No significant circularity: operator trained on independent random simulations

full rationale

The derivation trains a linear operator by minimizing average SoS error over an ensemble of independently generated random tissue models whose distribution is chosen to approximate in vivo statistics. This training set is external to any particular target measurement, so the learned operator is not equivalent to a fit or redefinition of the input echo-shift data. No equations reduce a claimed prediction to the input by construction, no uniqueness theorems are imported from prior self-work, and no ansatz is smuggled via citation. Phantom experiments supply an external check on transferability. The central claim therefore remains self-contained against the provided benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that simulated random tissue models capture in-vivo statistics and that the linear approximation is sufficiently accurate for the learned operator to transfer.

free parameters (1)
  • Tissue model distribution parameters
    The statistics used to generate the large number of random tissue models for training are chosen to represent expected in-vivo geometries and SoS values.
axioms (1)
  • domain assumption The echo-shift inverse problem admits a useful linear approximation
    The paper states it solves a linearized inverse problem and trains within that linear framework.

pith-pipeline@v0.9.0 · 5756 in / 1235 out tokens · 18756 ms · 2026-05-23T21:34:00.317696+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Wave-Equation Migration Velocity Analysis for Multistatic Synthetic Aperture Ultrasound

    physics.med-ph 2026-04 conditional novelty 8.0

    Wave-equation migration velocity analysis applied to medical ultrasound corrects sound speed aberrations, improving point target resolution from 1.22 mm to 0.32 mm and lesion contrast from 3.05 dB to 4.39 dB in phanto...

Reference graph

Works this paper leans on

28 extracted references · 28 canonical work pages · cited by 1 Pith paper

  1. [1]

    In vivo breast sound-speed imaging with ultrasound tomography

    Cuiping Li, Nebojsa Duric, Peter Littrup, and Lianjie Huang. In vivo breast sound-speed imaging with ultrasound tomography. Ultrasound in medicine & biology , 35(10):1615–1628, 2009

  2. [2]

    Non-linear inverse scattering: High resolution quantitative breast tissue tomography

    James Wiskin, David T Borup, SA Johnson, and M Berggren. Non-linear inverse scattering: High resolution quantitative breast tissue tomography. The Journal of the Acoustical Society of America , 131(5):3802–3813, 2012

  3. [3]

    Frequency domain ultrasound waveform tomography: breast imaging using a ring transducer

    GY Sandhu, Cuiping Li, Olivier Roy, S Schmidt, and Neb Duric. Frequency domain ultrasound waveform tomography: breast imaging using a ring transducer. Physics in Medicine & Biology , 60(14):5381–5398, 2015

  4. [4]

    L. Ruby, S. J. Sanabria, K. Martini, K. J. Dedes, D. V orburger, E. Oezkan, T. Frauenfelder, O. Goksel, and M. B. Rominger. Breast Cancer Assessment With Pulse-Echo Speed of Sound Ultrasound From Intrinsic Tissue Reflections: Proof-of-Concept. Investigative Radiology, 54:419–427, 2019

  5. [5]

    Detection of breast cancer with ultrasound tomography: First results with the Computed Ultrasound Risk Evaluation (CURE) prototype

    Nebojsa Duric, Peter Littrup, Lou Poulo, Alex Babkin, Roman Pevzner, Earle Holsapple, Olsi Rama, and Carri Glide. Detection of breast cancer with ultrasound tomography: First results with the Computed Ultrasound Risk Evaluation (CURE) prototype. Medical Physics, 34(2):773–785, 2007

  6. [6]

    Evaluation of sound speed for detection of liver fibrosis: prospective comparison with transient dynamic elastography and histology

    Bita Boozari, Andrej Potthoff, Ingmar Mederacke, Andreas Hahn, Ansgar Reising, Kinan Rifai, Heiner Wedemeyer, Matthias Bahr, Stefan Kubicka, Michael Manns, et al. Evaluation of sound speed for detection of liver fibrosis: prospective comparison with transient dynamic elastography and histology. Journal of Ultrasound in Medicine , 29(11):1581–1588, 2010

  7. [7]

    Robust sound speed estimation for ultrasound-based hepatic steatosis assessment

    Marion Imbault, Alex Faccinetto, Bruno-Félix Osmanski, Antoine Tissier, Thomas Deffieux, Jean-Luc Gennisson, Valérie Vilgrain, and Mickaël Tanter. Robust sound speed estimation for ultrasound-based hepatic steatosis assessment. Physics in Medicine & Biology , 62(9):3582–3598, 2017

  8. [8]

    Ultrasonic adaptive sound speed estimation for the diagnosis and quantification of hepatic steatosis: a pilot study

    Marco Dioguardi Burgio, Marion Imbault, Maxime Ronot, Alex Faccinetto, Bernard E Van Beers, Pierre- Emmanuel Rautou, Laurent Castera, Jean-Luc Gennisson, Mickael Tanter, and Valérie Vilgrain. Ultrasonic adaptive sound speed estimation for the diagnosis and quantification of hepatic steatosis: a pilot study. Ultraschall in der Medizin-European Journal of U...

  9. [9]

    Noninvasive estimation of local speed of sound by pulse-echo ultrasound in a rat model of nonalcoholic fatty liver

    Arsenii V Telichko, Rehman Ali, Thurston Brevett, Huaijun Wang, Jose G Vilches-Moure, Sukumar U Kumar, Ramasamy Paulmurugan, and Jeremy J Dahl. Noninvasive estimation of local speed of sound by pulse-echo ultrasound in a rat model of nonalcoholic fatty liver. Physics in Medicine & Biology , 67(1):015007, 2022

  10. [10]

    Stähli, C

    P. Stähli, C. Becchetti, N. Korta Martiartu, A. Berzigotti, M. Frenz, and M. Jaeger. First-in-human diagnostic study of hepatic steatosis with computed ultrasound tomography in echo mode. Communications Medicine, 3(1):176, 2023

  11. [11]

    Measurement of velocity of propagation from ultrasonic pulse-echo data

    DE Robinson, F Chen, and LS Wilson. Measurement of velocity of propagation from ultrasonic pulse-echo data. Ultrasound in medicine & biology , 8(4):413–420, 1982

  12. [12]

    An evaluation of an in vivo local sound speed estimation technique by the crossed beam method

    Masafumi Kondo, Kinya Takamizawa, Makoto Hirama, Kiyoshi Okazaki, Kazuhiro Iinuma, and Yasuaki Takehara. An evaluation of an in vivo local sound speed estimation technique by the crossed beam method. Ultrasound in medicine & biology, 16(1):65–72, 1990

  13. [13]

    Nonlinear simultaneous reconstruction of inhomogeneous compressibility and mass density distributions in unidirectional pulse-echo ultrasound imaging

    Markus C Hesse, Leili Salehi, and Georg Schmitz. Nonlinear simultaneous reconstruction of inhomogeneous compressibility and mass density distributions in unidirectional pulse-echo ultrasound imaging. Physics in Medicine & Biology, 58(17):6163, 2013

  14. [14]

    Local speed of sound estimation in tissue using pulse-echo ultrasound: Model-based approach

    Marko Jakovljevic, Scott Hsieh, Rehman Ali, Gustavo Chau Loo Kung, Dongwoon Hyun, and Jeremy J Dahl. Local speed of sound estimation in tissue using pulse-echo ultrasound: Model-based approach. The Journal of the Acoustical Society of America , 144(1):254–266, 2018

  15. [15]

    Spatial domain reconstruction for imaging speed-of-sound with pulse-echo ultrasound: simulation and in vivo study

    Sergio J Sanabria, Ece Ozkan, Marga Rominger, and Orcun Goksel. Spatial domain reconstruction for imaging speed-of-sound with pulse-echo ultrasound: simulation and in vivo study. Physics in Medicine & Biology , 63(21):215015, 2018

  16. [16]

    The convolutional interpretation of registration-based plane wave steered pulse-echo local sound speed estimators

    Anthony S Podkowa and Michael L Oelze. The convolutional interpretation of registration-based plane wave steered pulse-echo local sound speed estimators. Physics in Medicine & Biology , 65(2):025003, 2020

  17. [17]

    Micha Feigin, Daniel Freedman, and Brian W. Anthony. A deep learning framework for single-sided sound speed inversion in medical ultrasound. IEEE Transactions on Biomedical Engineering , 67(4):1142–1151, 2020. 15

  18. [18]

    Refraction-based speed of sound estimation in layered media: an angular approach

    Baptiste Hériard-Dubreuil, Adrien Besson, Frédéric Wintzenrieth, Claude Cohen-Bacrie, and Jean-Philippe Thiran. Refraction-based speed of sound estimation in layered media: an angular approach. IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, 70(6):486–497, 2023

  19. [19]

    Simson, Magdalini Paschali, Vasiliki Sideri-Lampretsa, Nassir Navab, and Jeremy J

    Walter A. Simson, Magdalini Paschali, Vasiliki Sideri-Lampretsa, Nassir Navab, and Jeremy J. Dahl. Investigating pulse-echo sound speed estimation in breast ultrasound with deep learning. Ultrasonics, 137:107179, 2024

  20. [20]

    Windowed Radon Transform for Robust Speed-of-Sound Imaging With Pulse-Echo Ultrasound.IEEE Transactions on Medical Imaging, 43(4):1579–1593, 2024

    Samuel Beuret, Baptiste Hériard-Dubreuil, Naiara Korta Martiartu, Michael Jaeger, and Jean-Philippe Thiran. Windowed Radon Transform for Robust Speed-of-Sound Imaging With Pulse-Echo Ultrasound.IEEE Transactions on Medical Imaging, 43(4):1579–1593, 2024

  21. [21]

    Computed Ultrasound Tomography in Echo Mode for Imaging Speed of Sound Using Pulse-Echo Sonography: Proof of Principle

    Michael Jaeger, Gerrit Held, Sara Peeters, Stefan Preisser, Michael Grünig, and Martin Frenz. Computed Ultrasound Tomography in Echo Mode for Imaging Speed of Sound Using Pulse-Echo Sonography: Proof of Principle. Ultrasound in Medicine & Biology , 41(1):235–250, 2015

  22. [22]

    Improved forward model for quantitative pulse-echo speed-of-sound imaging

    Patrick Stähli, Maju Kuriakose, Martin Frenz, and Michel Jaeger. Improved forward model for quantitative pulse-echo speed-of-sound imaging. Ultrasonics, 108:106168, 2020

  23. [23]

    Pulse-echo speed-of-sound imaging using convex probes

    Michael Jaeger, Patrick Stähli, Naiara Korta Martiartu, Parisa Salemi Yolgunlu, Thomas Frappart, Christophe Fraschini, and Martin Frenz. Pulse-echo speed-of-sound imaging using convex probes. Physics in Medicine & Biology, 67(21):215016, 2022

  24. [24]

    Bayesian Approach for a Robust Speed-of-Sound Reconstruction Using Pulse-Echo Ultrasound

    Patrick Stähli, Martin Frenz, and Michael Jaeger. Bayesian Approach for a Robust Speed-of-Sound Reconstruction Using Pulse-Echo Ultrasound. IEEE Transactions on Medical Imaging, 40(2):457–467, 2021

  25. [25]

    Excluding Echo Shift Noise in Real-Time Pulse-Echo Speed-of-Sound Imaging

    Parisa Salemi Yolgunlu, Naiara Korta Martiartu, Urs Richard Gerber, Martin Frenz, and Michael Jaeger. Excluding Echo Shift Noise in Real-Time Pulse-Echo Speed-of-Sound Imaging. Sensors, 23(12), 2023

  26. [26]

    Learnt correction for regularization-related biases in pulse-echo speed-of-sound imaging

    Michael Jaeger, Urs Richard Gerber, Parisa Salemi Yolgunlu, Naiara Korta Martiartu, and Martin Frenz. Learnt correction for regularization-related biases in pulse-echo speed-of-sound imaging. In 2023 IEEE International Ultrasonics Symposium (IUS), pages 1–4, 2023

  27. [27]

    Ultrasound beam simulations in inhomogeneous tissue geometries using the hybrid angular spectrum method

    Urvi Vyas and Douglas Christensen. Ultrasound beam simulations in inhomogeneous tissue geometries using the hybrid angular spectrum method. IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, 59(6):1093–1100, 2012

  28. [28]

    Tanter, J

    Gabriel Montaldo, M. Tanter, J. Bercoff, Nicolas Benech, and Mathias Fink. Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography.IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, 56(3):489–506, 2009. 16