Objective: Develop a dense algorithm for calculating the speed-of-sound shift between consecutive acoustic acquisitions as a noninvasive means to evaluating temperature change during thermal ablation. 
Methods: An algorithm for dense speed-of-sound shift imaging (DSI) was developed to simultaneously incorporate information from the entire field of view using a combination of dense optical flow and inverse problem regularization, thus speeding up the calculation and introducing spatial agreement between pixels natively. Thermal ablation monitoring consisted of two main steps: pixel shift tracking using Farneback optical flow, and mathematical modeling of the relationship between the pixel displacement and temperature change as an inverse problem to find the speed-of-sound shift. A calibration constant translates from speed-of-sound shift to temperature change. The method performance was tested in ex-vivo samples and compared to standard thermal strain imaging (TSI) methods. 
Main results: Thermal ablation at a frequency of 2 MHz was applied to an agarose phantom that created a speed-of-sound shift measured by an L12-5 imaging transducer. A focal spot was reconstructed by solving the inverse problem. Next, a thermocouple measured the temperature rise during thermal ablation of ex-vivo chicken breast to calibrate the setup. Temperature changes between 3-15ºC was measured with high thermometry precision of less than 2ºC error for temperature changes as low as 8ºC. The DSI method outperformed standard TSI in both spatial coherence and runtime in HIFU-induced hyperthermia
Significance: Dense ultrasonic speed-of-sound shift imaging can successfully monitor the speed-of-sound shift introduced by thermal ablation. This technique is faster and more robust than current methods, and therefore can be used as a noninvasive, real time and cost-effective thermometry method, with high clinical applicability.