Acceleration-Based Kalman Tracking for Super-Resolution Ultrasound Imaging in vivo

Biao Huang,Jipeng Yan,Megan Morris,Victoria Sinnett,Navita Somaiah,Meng-Xing Tang

2023-04-03

Abstract:Super-resolution ultrasound can image microvascular structure and flow at sub-wave-diffraction resolution based on localising and tracking microbubbles. Currently, tracking microbubbles accurately under limited imaging frame rates and high microbubble concentrations remains a challenge, especially under the effect of cardiac pulsatility and in highly curved vessels. In this study, an acceleration-incorporated microbubble motion model is introduced into a Kalman tracking framework. The tracking performance was evaluated using simulated microvasculature with different microbubble motion parameters and acquisition frame rates, and in vivo human breast tumour ultrasound datasets. The simulation results show that the acceleration-based method outperformed the non-acceleration-based method at different levels of acceleration and acquisition frame rates and achieved significant improvement in true positive rate (up to 10.03%), false negative rate (up to 28.61%) and correctly pairing fraction (up to 170.14%). The proposed method can also reduce errors in vasculature reconstruction via the acceleration-based nonlinear interpolation, compared with linear interpolation (up to 19 um). The tracking results from temporally downsampled low frame rate in vivo datasets from human breast tumours show that the proposed method has better microbubble tracking performance than the baseline method, if using results from the initial high frame data as reference. Finally, the acceleration estimated from tracking results also provides a spatial speed gradient map that may contain extra valuable diagnostic information.

Image and Video Processing

What problem does this paper attempt to address?

This paper attempts to solve the problem of how to accurately track microbubbles in super - resolution ultrasound imaging under the conditions of limited frame rate and high microbubble concentration, especially under the influence of cardiac pulsatility and highly curved blood vessels. Specifically, the paper proposes an acceleration - based Kalman tracking framework to improve the accuracy of microbubble tracking.

### Main Problems and Solutions

1. **Challenges of Microbubble Tracking at Low Frame Rates**
- **Problem**: Current methods are difficult to accurately track microbubbles at low frame rates and high microbubble concentrations, especially in the presence of cardiac pulsatility and highly curved blood vessels.
- **Solution**: A microbubble motion model including acceleration was introduced and integrated into the Kalman tracking framework. By considering the acceleration change of microbubbles between adjacent frames, the motion trajectory of microbubbles can be predicted and tracked more accurately.

2. **Microbubble Motion in Curved Blood Vessels**
- **Problem**: The linear motion model assumes that microbubbles move at a constant speed between adjacent frames, which is not valid in curved blood vessels, resulting in trajectory reconstruction errors.
- **Solution**: An acceleration - based non - linear interpolation method was proposed for reconstructing microbubble trajectories. This method can better adapt to the microbubble motion in curved blood vessels and provide more accurate trajectory reconstruction.

3. **Extraction of Diagnostic Information in Clinical Applications**
- **Problem**: Existing super - resolution ultrasound imaging techniques do not fully utilize the acceleration information in microbubble motion.
- **Solution**: By estimating the acceleration in microbubble motion, a spatial velocity gradient map was generated, providing additional diagnostic information, which may be helpful in distinguishing between benign and malignant tumors.

### Formula Summary

- **Kalman Filtering State - Transition Model**
\[
S_{k|k - 1}=F S_{k - 1|k - 1}+Q
\]
\[
\begin{bmatrix}
x(k)\\
v_x(k)\\
a_x(k)\\
y(k)\\
v_y(k)\\
a_y(k)
\end{bmatrix}
=
\begin{bmatrix}
1&\Delta t&\frac{\Delta t^2}{2}&0&0&0\\
0&1&\Delta t&0&0&0\\
0&0&1&0&0&0\\
0&0&0&1&\Delta t&\frac{\Delta t^2}{2}\\
0&0&0&0&1&\Delta t\\
0&0&0&0&0&1
\end{bmatrix}
\times
\begin{bmatrix}
x(k - 1)\\
v_x(k - 1)\\
a_x(k - 1)\\
y(k - 1)\\
v_y(k - 1)\\
a_y(k - 1)
\end{bmatrix}
+
\begin{bmatrix}
\frac{\Delta t^4}{4}&\frac{\Delta t^3}{2}&\frac{\Delta t^2}{2}&0&0&0\\
\frac{\Delta t^3}{2}&\Delta t^2&\Delta t&0&0&0\\
\frac{\Delta t^2}{2}&\Delta t&1&0&0&0\\
0&0&0&\frac{\Delta t^4}{4}&\frac{\Delta t^3}{2}&\frac{\Delta t^2}{2}\\
0&0&0&\frac{\Delta t^3}{2}&\Delta t^2&\Delta t\\
0&
\end{bmatrix}