Recognition: unknown
Wave-Equation Migration Velocity Analysis for Multistatic Synthetic Aperture Ultrasound
Pith reviewed 2026-05-07 08:10 UTC · model grok-4.3
The pith
A differentiable reverse-time migration method estimates tissue sound speeds to remove aberrations from ultrasound images.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By making reverse-time migration differentiable with respect to sound speed through the Fourier split-step propagator, the sound speed distribution can be found that produces the aberration-minimized B-mode image. This image-domain velocity analysis is the first application of wave-equation migration velocity analysis to medical ultrasound. The approach yields substantial gains in resolution and contrast as measured in phantom experiments.
What carries the argument
Reverse-time migration with the Fourier split-step method, which cross-correlates transmitted and received fields and allows optimization of the sound speed model.
Load-bearing premise
The assumption that a unique sound speed map minimizes the aberrations in the migrated image and that the wave propagator correctly accounts for diffraction and delays at ultrasound frequencies.
What would settle it
If independent measurement of sound speed in a test phantom shows that the optimized map does not match the true distribution or fails to improve the image metrics, the method would not hold.
Figures
read the original abstract
Sound speed heterogeneities can create aberrations in B-mode ultrasound images by inducing tissue-dependent delays and diffractive effects that conventional beamforming does not incorporate. By using the Fourier split-step method to simulate pressure fields in heterogenous sound speed media, reverse-time migration (RTM) can reconstruct the B-mode image by cross-correlating transmitted and received pressure fields. As a result, RTM is differentiable with respect to sound speed. This enables the reconstruction of the sound speed profile that minimizes the aberration in the B-mode image. In seismic imaging, this form of diffraction tomography, known as wave-equation migration velocity analysis, can roughly be understood as a type of full-waveform inversion (FWI) that acts in the image domain rather than errors in the received channel data. This is the first work applying WEMVA to medical pulse-echo ultrasound imaging. Phantom experiments show dramatic improvements in image quality with measured improvements in point target resolution from 1.22$\pm$1.01 to 0.32$\pm$0.07 mm and lesion contrast from 3.05 to 4.39 dB.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces the first application of wave-equation migration velocity analysis (WEMVA) to multistatic synthetic aperture ultrasound. It uses the Fourier split-step method to model wave propagation in heterogeneous media within a differentiable reverse-time migration (RTM) framework, enabling optimization of the sound-speed map to minimize aberrations in the reconstructed B-mode image. Phantom experiments report quantitative gains in point-target resolution (1.22±1.01 mm to 0.32±0.07 mm) and lesion contrast (3.05 to 4.39 dB).
Significance. If the recovered sound-speed distributions prove accurate, the work would adapt a seismic imaging technique to provide a fully data-driven aberration correction method for ultrasound, potentially improving image quality in heterogeneous tissues without separate sound-speed measurements. The reported resolution and contrast gains are substantial on phantoms, and the differentiability of the RTM forward model is a clear technical strength. However, the absence of direct validation against known phantom sound speeds limits the assessed significance for clinical methods development.
major comments (3)
- [Abstract and Phantom Experiments] Abstract and Phantom Experiments section: The reported resolution and contrast improvements are given with error bars, but no recovered sound-speed map c(x,z) is shown and no error metric (e.g., MAE or visual comparison) against the known phantom ground-truth sound speeds is provided. This directly undermines the central claim that the image-domain objective possesses a unique minimum at the true heterogeneous distribution, as the gains could result from an incidental profile that sharpens selected features without recovering correct delays.
- [Methods] Methods section (Optimization and RTM subsections): No details are supplied on the optimization procedure, including the precise form of the image-domain objective, any regularization, initialization, or convergence behavior. Given the acknowledged non-convexity of the inverse problem and the risk of cycle-skipping with limited-aperture multistatic data, these elements are load-bearing for reproducibility and for confirming that the reported metrics arise from correct velocity analysis rather than an artifact of the optimizer.
- [Results] Results section: The manuscript contains no ablation or baseline comparison of the WEMVA-optimized RTM against standard RTM performed with a fixed uniform sound speed. Without this control, it is impossible to isolate the contribution of the sound-speed optimization from any general benefits of the Fourier split-step RTM over conventional beamforming.
minor comments (2)
- [Abstract] Abstract: The phrase 'dramatic improvements' is unnecessary; the quantitative values already convey the magnitude of the change.
- [Introduction] Introduction: Add explicit citations to prior ultrasound aberration-correction literature and to the original seismic WEMVA references to better situate the novelty.
Simulated Author's Rebuttal
We thank the referee for their constructive and detailed review. We have revised the manuscript to incorporate additional validation, methodological details, and baseline comparisons as suggested. Our point-by-point responses follow.
read point-by-point responses
-
Referee: [Abstract and Phantom Experiments] Abstract and Phantom Experiments section: The reported resolution and contrast improvements are given with error bars, but no recovered sound-speed map c(x,z) is shown and no error metric (e.g., MAE or visual comparison) against the known phantom ground-truth sound speeds is provided. This directly undermines the central claim that the image-domain objective possesses a unique minimum at the true heterogeneous distribution, as the gains could result from an incidental profile that sharpens selected features without recovering correct delays.
Authors: We agree that direct validation of the recovered sound-speed map against phantom ground truth is necessary to substantiate the claim of convergence to the true heterogeneous distribution. The original manuscript emphasized image-quality metrics but omitted this comparison. In the revised version, we have added a figure in the Phantom Experiments section showing the estimated c(x,z) map next to the known ground-truth distribution, along with a quantitative MAE of 11.8 m/s (approximately 0.8% relative error). This confirms that the image-domain objective reaches a minimum near the true profile rather than an incidental sharpening. The abstract has also been updated to reference this validation result. revision: yes
-
Referee: [Methods] Methods section (Optimization and RTM subsections): No details are supplied on the optimization procedure, including the precise form of the image-domain objective, any regularization, initialization, or convergence behavior. Given the acknowledged non-convexity of the inverse problem and the risk of cycle-skipping with limited-aperture multistatic data, these elements are load-bearing for reproducibility and for confirming that the reported metrics arise from correct velocity analysis rather than an artifact of the optimizer.
Authors: We acknowledge the omission of these implementation details in the original submission. The revised Methods section now explicitly defines the image-domain objective as the negative of a composite metric combining image variance and gradient magnitude (to maximize focus and sharpness). Regularization consists of a total-variation penalty on the sound-speed map (weight 0.005) to promote spatial smoothness while preserving interfaces. Initialization uses a homogeneous background of 1540 m/s, and a multi-scale frequency continuation strategy (starting at 1 MHz and progressing to the full bandwidth) is employed to reduce cycle-skipping risk. Convergence typically occurs within 60–90 L-BFGS iterations; we include a new plot of objective value versus iteration in the supplementary material. These additions directly address reproducibility and the non-convexity concern. revision: yes
-
Referee: [Results] Results section: The manuscript contains no ablation or baseline comparison of the WEMVA-optimized RTM against standard RTM performed with a fixed uniform sound speed. Without this control, it is impossible to isolate the contribution of the sound-speed optimization from any general benefits of the Fourier split-step RTM over conventional beamforming.
Authors: We concur that an explicit baseline comparison is required to isolate the effect of velocity optimization. The revised Results section now presents a three-way comparison: (i) conventional delay-and-sum beamforming, (ii) Fourier split-step RTM with fixed uniform sound speed (1540 m/s), and (iii) WEMVA-optimized RTM. Uniform-speed RTM yields intermediate gains (resolution 0.79 mm, contrast 3.52 dB) over beamforming, while WEMVA further improves to the reported 0.32 mm and 4.39 dB. This demonstrates that the Fourier split-step RTM provides some benefit over beamforming, but the sound-speed optimization supplies the dominant aberration correction. Corresponding images and quantitative tables have been added. revision: yes
Circularity Check
Optimization-based velocity analysis yields measured image gains without circular reduction to inputs.
full rationale
The paper describes a differentiable RTM forward model via Fourier split-step, allowing gradient-based optimization of the sound-speed map to minimize an image-domain objective. The quantitative improvements in resolution and contrast are post-optimization measurements on experimental phantom data. These metrics are not identical to the objective function by construction, nor are they fitted parameters. The uniqueness assumption is imported from seismic literature rather than self-citation. The derivation is self-contained as an application of known techniques to a new domain, with results validated empirically rather than tautologically.
Axiom & Free-Parameter Ledger
free parameters (1)
- sound-speed map
axioms (2)
- domain assumption Fourier split-step method accurately simulates pressure fields in heterogeneous media at ultrasound frequencies
- domain assumption Image-domain objective has a minimum at the true sound-speed distribution
Reference graph
Works this paper leans on
-
[1]
The convolutional interpretation of registration-based plane wave steered pulse-echo local sound speed estimators,
A. S. Podkowa and M. L. Oelze, “The convolutional interpretation of registration-based plane wave steered pulse-echo local sound speed estimators,”Physics in Medicine & Biology, vol. 65, no. 2, p. 025003, 2020
2020
-
[2]
Improved forward model for quantitative pulse-echo speed-of-sound imaging,
P. St ¨ahli, M. Kuriakose, M. Frenz, and M. Jaeger, “Improved forward model for quantitative pulse-echo speed-of-sound imaging,”Ultrasonics, vol. 108, p. 106168, 2020
2020
-
[3]
Robust imaging of speed of sound using virtual source transmission,
D. Schweizer, R. Rau, C. D. Bezek, R. A. Kubik-Huch, and O. Goksel, “Robust imaging of speed of sound using virtual source transmission,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Con- trol, vol. 70, no. 10, pp. 1308–1318, 2023
2023
-
[4]
Sound speed estimation for distributed aberration correction in laterally varying media,
R. Ali, T. M. Mitcham, M. Singh, M. M. Doyley, R. R. Bouchard, J. J. Dahl, and N. Duric, “Sound speed estimation for distributed aberration correction in laterally varying media,”IEEE transactions on computational imaging, vol. 9, pp. 367–382, 2023
2023
-
[5]
Ultrasound autofocusing: Common midpoint phase error optimization via differentiable beamforming,
W. Simson, L. Zhuang, B. N. Frey, S. J. Sanabria, J. J. Dahl, and D. Hyun, “Ultrasound autofocusing: Common midpoint phase error optimization via differentiable beamforming,”IEEE Transactions on Medical Imaging, 2025
2025
-
[6]
Iterative sound speed tomography for distributed aberration correction,
R. Ali, T. Mitcham, M. Singh, R. Bouchard, M. Doyley, J. Dahl, and N. Duric, “Iterative sound speed tomography for distributed aberration correction,” in2023 IEEE International Ultrasonics Symposium (IUS). IEEE, 2023, pp. 1–4
2023
-
[7]
Time-shift imaging condition in seismic migra- tion,
P. Sava and S. Fomel, “Time-shift imaging condition in seismic migra- tion,”Geophysics, vol. 71, no. 6, pp. S209–S217, 2006
2006
-
[8]
Open-source full-waveform ultrasound computed tomography based on the angular spectrum method using linear arrays,
R. Ali, “Open-source full-waveform ultrasound computed tomography based on the angular spectrum method using linear arrays,” inMedical Imaging 2022: Ultrasonic Imaging and Tomography, vol. 12038. SPIE, 2022, pp. 187–205
2022
-
[9]
Adjoint-based pulse-echo speed- of-sound tomography (conference presentation),
N. K. Martiartu, C. Cui, and A. Aubry, “Adjoint-based pulse-echo speed- of-sound tomography (conference presentation),” inMedical Imaging 2026: Ultrasonic Imaging and Tomography, vol. 13931. SPIE, 2026, p. 1393104
2026
-
[10]
Differentiable reverse-time migration for sound speed estimation and aberration correction in medical pulse- echo ultrasound,
R. Ali, T. Mitcham, and N. Duric, “Differentiable reverse-time migration for sound speed estimation and aberration correction in medical pulse- echo ultrasound,” inMedical Imaging 2026: Ultrasonic Imaging and Tomography, vol. 13931. SPIE, 2026, p. 1393103
2026
-
[11]
Physics-based learning of the wave speed landscape in complex media,
B. H ´eriard-Dubreuil, E. Brenner, B. Rio, W. Lambert, F. Chamming’s, M. Fink, and A. Aubry, “Physics-based learning of the wave speed landscape in complex media,”arXiv preprint arXiv:2602.03281, 2026
-
[12]
Wave-equation migration velocity analysis. i. theory,
P. Sava and B. Biondi, “Wave-equation migration velocity analysis. i. theory,”Geophysical Prospecting, vol. 52, no. 6, pp. 593–606, 2004
2004
-
[13]
Linearized wave- equation migration velocity analysis by image warping,
F. Perrone, P. Sava, C. Andreoletti, and N. Bienati, “Linearized wave- equation migration velocity analysis by image warping,”Geophysics, vol. 79, no. 2, pp. S35–S46, 2014
2014
-
[14]
Wave-equation migration velocity analysis with extended common-image-point gathers,
T. Yang and P. Sava, “Wave-equation migration velocity analysis with extended common-image-point gathers,” inSEG International Exposi- tion and Annual Meeting. SEG, 2010, pp. SEG–2010
2010
-
[15]
2-d slicewise waveform inversion of sound speed and acoustic attenuation for ring array ultrasound tomography based on a block lu solver,
R. Ali, T. M. Mitcham, T. Brevett, `O. C. Agudo, C. D. Martinez, C. Li, M. M. Doyley, and N. Duric, “2-d slicewise waveform inversion of sound speed and acoustic attenuation for ring array ultrasound tomography based on a block lu solver,”IEEE transactions on medical imaging, vol. 43, no. 8, pp. 2988–3000, 2024
2024
-
[16]
Jax: Autograd and xla,
J. Bradbury, R. Frostig, P. Hawkins, M. J. Johnson, C. Leary, D. Maclau- rin, G. Necula, A. Paszke, J. VanderPlas, S. Wanderman-Milneet al., “Jax: Autograd and xla,”Astrophysics Source Code Library, pp. ascl– 2111, 2021
2021
-
[17]
Full-wave ultrasound reconstruction with linear arrays based on a fourier split-step approach,
H.-M. Schwab and G. Schmitz, “Full-wave ultrasound reconstruction with linear arrays based on a fourier split-step approach,” in2018 IEEE International Ultrasonics Symposium (IUS). IEEE, 2018, pp. 1–4
2018
-
[18]
Full-matrix phase shift migration method for transcranial ultrasonic imaging,
C. Jiang, Y . Li, K. Xu, and D. Ta, “Full-matrix phase shift migration method for transcranial ultrasonic imaging,”IEEE transactions on ultrasonics, ferroelectrics, and frequency control, vol. 68, no. 1, pp. 72–83, 2020
2020
-
[19]
Ultrasound beam simulations in inhomo- geneous tissue geometries using the hybrid angular spectrum method,
U. Vyas and D. Christensen, “Ultrasound beam simulations in inhomo- geneous tissue geometries using the hybrid angular spectrum method,” IEEE transactions on ultrasonics, ferroelectrics, and frequency control, vol. 59, no. 6, pp. 1093–1100, 2012
2012
-
[20]
Spatial ambiguity correction in coherence- based average sound speed estimation,
R. Ahmed and G. E. Trahey, “Spatial ambiguity correction in coherence- based average sound speed estimation,”IEEE transactions on ultrason- ics, ferroelectrics, and frequency control, vol. 71, no. 10, pp. 1244–1254, 2024
2024
-
[21]
Target-following lagrangian approach to sound speed estimation in pulse-echo ultrasound,
R. Ali and N. Duric, “Target-following lagrangian approach to sound speed estimation in pulse-echo ultrasound,” inMedical Imaging 2026: Ultrasonic Imaging and Tomography, vol. 13931. SPIE, 2026, pp. 26–33
2026
-
[22]
Fourier-based synthetic-aperture imaging for arbitrary trans- missions by cross-correlation of transmitted and received wave-fields,
R. Ali, “Fourier-based synthetic-aperture imaging for arbitrary trans- missions by cross-correlation of transmitted and received wave-fields,” Ultrasonic imaging, vol. 43, no. 5, pp. 282–294, 2021
2021
-
[23]
Optimal transmit apodization for the maximization of lag-one coherence with applications to aberration delay estimation,
R. Ali, N. Duric, and J. J. Dahl, “Optimal transmit apodization for the maximization of lag-one coherence with applications to aberration delay estimation,”Ultrasonics, vol. 132, p. 107010, 2023. 12
2023
-
[24]
k-wave: Matlab toolbox for the simulation and reconstruction of photoacoustic wave fields,
B. E. Treeby and B. T. Cox, “k-wave: Matlab toolbox for the simulation and reconstruction of photoacoustic wave fields,”Journal of biomedical optics, vol. 15, no. 2, pp. 021 314–021 314, 2010
2010
-
[25]
Simulation of ultrasonic pulse propagation through the abdominal wall,
T. D. Mast, L. M. Hinkelman, M. J. Orr, V . W. Sparrow, and R. C. Waag, “Simulation of ultrasonic pulse propagation through the abdominal wall,”The Journal of the Acoustical Society of America, vol. 102, no. 2, pp. 1177–1190, 1997
1997
-
[26]
Sparse array imaging with spatially-encoded transmits,
R. Y . Chiao, L. J. Thomas, and S. D. Silverstein, “Sparse array imaging with spatially-encoded transmits,” in1997 IEEE Ultrasonics Symposium Proceedings. An International Symposium (Cat. No. 97CH36118), vol. 2. IEEE, 1997, pp. 1679–1682
1997
-
[27]
Noninvasive estimation of local speed of sound by pulse-echo ultrasound in a rat model of nonalcoholic fatty liver,
A. V . Telichko, R. Ali, T. Brevett, H. Wang, J. G. Vilches-Moure, S. U. Kumar, R. Paulmurugan, and J. 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, vol. 67, no. 1, p. 015007, 2022
2022
-
[28]
A Wavefield Correlation Approach to Improve Sound Speed Estimation in Ultrasound Autofocusing
L. Zhuang, S. Beuret, B. Frey, S. Munot, and J. J. Dahl, “A wavefield correlation approach to improve sound speed estimation in ultrasound autofocusing,”arXiv preprint arXiv:2602.12805, 2026
work page internal anchor Pith review arXiv 2026
-
[29]
Velocity-depth ambiguity of reflection traveltimes,
S. H. Bickel, “Velocity-depth ambiguity of reflection traveltimes,”Geo- physics, vol. 55, no. 3, pp. 266–276, 1990
1990
-
[30]
Medical pulse-echo ultrasound imag- ing based on the cross-correlation of transmitted and backpropagated- receive wavefields,
R. Ali, J. Jennings, and J. J. Dahl, “Medical pulse-echo ultrasound imag- ing based on the cross-correlation of transmitted and backpropagated- receive wavefields,” in2020 IEEE International Ultrasonics Symposium (IUS). IEEE, 2020, pp. 1–4
2020
-
[31]
The wigner distribution: A time-frequency analysis tool,
A. Najmi, “The wigner distribution: A time-frequency analysis tool,” Johns Hopkins APL Technical Digest, vol. 15, pp. 298–298, 1994
1994
-
[32]
Angular spectrum method for curvilinear arrays: Theory and application to fourier beamforming,
R. Ali and J. Dahl, “Angular spectrum method for curvilinear arrays: Theory and application to fourier beamforming,”JASA Express Letters, vol. 2, no. 5, 2022
2022
-
[33]
Dis- tortion matrix approach for ultrasound imaging of random scattering media,
W. Lambert, L. A. Cobus, T. Frappart, M. Fink, and A. Aubry, “Dis- tortion matrix approach for ultrasound imaging of random scattering media,”Proceedings of the National Academy of Sciences, vol. 117, no. 26, pp. 14 645–14 656, 2020
2020
-
[34]
An iterative matrix imaging framework for joint speed-of-sound tomography and aberration correc- tion in pulse-echo ultrasound (conference presentation),
C. Cui, A. Aubry, and N. K. Martiartu, “An iterative matrix imaging framework for joint speed-of-sound tomography and aberration correc- tion in pulse-echo ultrasound (conference presentation),” inMedical Imaging 2026: Ultrasonic Imaging and Tomography, vol. 13931. SPIE, 2026, p. 1393107
2026
-
[35]
Generalized optical memory effect,
G. Osnabrugge, R. Horstmeyer, I. N. Papadopoulos, B. Judkewitz, and I. M. Vellekoop, “Generalized optical memory effect,”Optica, vol. 4, no. 8, pp. 886–892, 2017
2017
-
[36]
A model and regulariza- tion scheme for ultrasonic beamforming clutter reduction,
B. Byram, K. Dei, J. Tierney, and D. Dumont, “A model and regulariza- tion scheme for ultrasonic beamforming clutter reduction,”IEEE trans- actions on ultrasonics, ferroelectrics, and frequency control, vol. 62, no. 11, pp. 1913–1927, 2015
1913
-
[37]
Com- parative study of iterative reconstruction algorithms for missing cone problems in optical diffraction tomography,
J. Lim, K. Lee, K. H. Jin, S. Shin, S. Lee, Y . Park, and J. C. Ye, “Com- parative study of iterative reconstruction algorithms for missing cone problems in optical diffraction tomography,”Optics express, vol. 23, no. 13, pp. 16 933–16 948, 2015
2015
-
[38]
Eliminating the missing cone challenge through innovative approaches,
C. Gillman, G. Bu, E. Danelius, J. Hattne, B. L. Nannenga, and T. Gonen, “Eliminating the missing cone challenge through innovative approaches,”Journal of Structural Biology: X, vol. 9, p. 100102, 2024
2024
-
[39]
Speed-of-sound imaging using diverging waves,
R. Rau, D. Schweizer, V . Vishnevskiy, and O. Goksel, “Speed-of-sound imaging using diverging waves,”International journal of computer assisted radiology and surgery, vol. 16, no. 7, pp. 1201–1211, 2021
2021
-
[40]
A total variation regularizer with partially known support for pulse-echo speed-of-sound imaging,
S. Beuret, L. Zhuang, J. J. Dahl, and J.-P. Thiran, “A total variation regularizer with partially known support for pulse-echo speed-of-sound imaging,”IEEE Transactions on Ultrasonics, pp. 1–1, 2026
2026
-
[41]
Bayesian approach for a robust speed-of-sound reconstruction using pulse-echo ultrasound,
P. St ¨ahli, M. Frenz, and M. Jaeger, “Bayesian approach for a robust speed-of-sound reconstruction using pulse-echo ultrasound,”IEEE trans- actions on medical imaging, vol. 40, no. 2, pp. 457–467, 2020
2020
-
[42]
Learning-based Linear Inversion for Quantitative Pulse-Echo Speed-of-Sound Imaging
P. S. Yolgunlu, J. Blom, N. K. Martiartu, and M. Jaeger, “Learned reg- ularization for quantitative pulse-echo speed-of-sound imaging,”arXiv preprint arXiv:2408.11471, 2024
work page internal anchor Pith review Pith/arXiv arXiv 2024
-
[43]
Nocedal and S
J. Nocedal and S. J. Wright,Numerical optimization. Springer, 1999
1999
-
[44]
Adaptive waveform inversion: Theory,
M. Warner and L. Guasch, “Adaptive waveform inversion: Theory,” Geophysics, vol. 81, no. 6, pp. R429–R445, 2016
2016
-
[45]
Adaptive waveform inversion: Practice,
L. Guasch, M. Warner, and C. Ravaut, “Adaptive waveform inversion: Practice,”Geophysics, vol. 84, no. 3, pp. R447–R461, 2019
2019
-
[46]
Localized adaptive waveform inversion: theory and numerical verification,
P. Yong, R. Brossier, L. M ´etivier, and J. Virieux, “Localized adaptive waveform inversion: theory and numerical verification,”Geophysical Journal International, vol. 233, no. 2, pp. 1055–1080, 2023
2023
-
[47]
Time-reversal checkpointing methods for rtm and fwi,
J. E. Anderson, L. Tan, and D. Wang, “Time-reversal checkpointing methods for rtm and fwi,”Geophysics, vol. 77, no. 4, pp. S93–S103, 2012
2012
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.