Author to whom any correspondence should be addressed.
Phys Med Biol. 2020 Jan 17;65(2):025003. doi: 10.1088/1361-6560/ab6071.
Pulse-echo reconstruction of sound speed has long been considered a difficult problem within the domain of quantitative biomedical ultrasound. However, recent results (Jaeger 2015 Ultrasound Med. Biol. 41 235-50; Jaeger and Frenz 2015 Ultrasonics 62 299-304; Jaeger et al 2015 Phys. Med. Biol. 60 4497-515) have demonstrated that pulse-echo reconstructions of sound speed are achievable by exploiting correlations in post-beamformed data from steered, plane-wave excitations in the presence of diffuse scatterers. Despite these recent advances, a coherent theoretical imaging framework for describing the approach and results is lacking in the literature. In this work, the problem of sound speed reconstruction using steered, plane-wave excitations is reformulated as a truncated convolutional problem, and the theoretical implications of this reformulation are explored. Additionally, a matrix-free algorithm is proposed that leverages the computational and storage advantages of the fast Fourier transform (FFT) while simultaneously avoiding FFT wraparound artifacts. In particular, the storage constraints of the approach are reduced down from [Formula: see text] to [Formula: see text] over full matrix reconstruction, making this approach a better candidate for large reconstructions on clinical machines. This algorithm was then tested in the open source simulation package k-Wave to assess its robustness to modeling error and resolution reduction was demonstrated under full-wave propagation conditions relative to ideal straight-ray simulations. The method was also validated in a phantom experiment.
声速的脉冲回波重建长期以来一直被认为是定量生物医学超声领域的一个难题。然而,最近的研究结果(Jaeger 2015 Ultrasound Med. Biol. 41 235-50;Jaeger 和 Frenz 2015 Ultrasonics 62 299-304;Jaeger 等人 2015 Phys. Med. Biol. 60 4497-515)表明,通过利用定向平面波激励后波束形成数据中的相关性,可以实现声速的脉冲回波重建,前提是存在漫散射体。尽管有这些最新进展,但文献中缺乏用于描述该方法和结果的相干理论成像框架。在这项工作中,使用定向平面波激励进行声速重建的问题被重新表述为截断卷积问题,并探讨了这种重新表述的理论意义。此外,还提出了一种无矩阵算法,该算法利用快速傅里叶变换(FFT)的计算和存储优势,同时避免了 FFT 缠绕伪影。特别是,该方法的存储限制从[公式:见文本]降低到[公式:见文本],适用于完整矩阵重建,这使得该方法成为临床机器上进行大型重建的更好候选方案。然后,在开源仿真包 k-Wave 中对该算法进行了测试,以评估其对建模误差的鲁棒性,并证明了在全波传播条件下相对于理想直线射线模拟的分辨率降低。该方法还在体模实验中得到了验证。