一种瞬时弹性成像快速时移估计方法

摘要

本发明公开了一种瞬时弹性成像快速时移估计方法，该方法包括以下步骤：接收不同时刻的N个超声回波信号s₁，…，s_N；根据预先设定的窗长和重叠率，将所述N个超声回波信号s₁，…，s_N分为若干个信号片段，所述N个超声回波信号s₁，…，s_N上相同深度的信号片段分别表示为s_1k，…，s_Nk，计算所述信号片段s_2k相对所述超声回波信号s₁的时移d_2k，并以所述时移d_2k作为先验值计算所述信号片段s_3k相对所述超声回波信号s₁的时移d_3k，以此类推，直到以d_(n-1)k作为先验值计算所述信号片段s_nk和所述超声回波信号s₁的时移d_nk。本发明可以有效地缩小时移估计的范围，提高计算速度。

主权项

一种瞬时弹性成像快速时移估计方法，其特征在于，包括以下步骤：接收不同时刻的N个超声回波信号s₁、…、s_N；根据预先设定的窗长和重叠率，将所述N个超声回波信号s₁、…、s_N分为若干个信号片段，所述N个超声回波信号s₁、…、s_N上相同深度的信号片段分别表示为s_1k、…、s_Nk，所述s_nk＝s_n(k:M1+k‑1)，其中信号片段s_nk为超声回波信号信号s_n的第k至第M1+k‑1个采样点，n＝1、…、N，M1为窗长；计算所述信号片段s_2k相对所述超声回波信号s₁的时移d_2k，并以所述时移d_2k作为先验值计算所述信号片段s_3k相对所述超声回波信号s₁的时移d_3k，以此类推，直到以d_(n‑1)k作为先验值计算所述信号片段s_nk和所述超声回波信号s₁的时移d_nk，其中n＝2、…、N；所述计算所述信号片段s_2k相对所述超声回波信号s₁的时移d_2k步骤包括：分别计算信号片段s_2k与超声回波信号s₁(k+i:M+k+i‑1)的绝对差和SAD值sad(i)，i＝‑L、…、L，其中L是一个正整数，为时移估计的搜索范围，即$sad\left(i\right)=\underset{q=k}{\overset{M1+k-1}{\Sigma }}|s_2\left(q\right)-s_1(q+i)|,i=-L,...,L$计算d_2k＝arg(min(sad(i))，其中i＝‑L、…、L，即找出sad(i)的最小值出现的位置d，所述信号片段s_2k相对s₁的时移d_2k＝d；所述以d_(n‑1)k作为先验值计算所述信号片段s_nk和所述超声回波信号s₁的时移d_nk步骤包括：分别计算s_nk与s₁(k+i:M+k+i‑1)的SAD值sad(i)，i＝d_(n‑1)k‑L、…、d_(n‑1)k+L，即$sad\left(i\right)=\underset{q=k}{\overset{M1+k-1}{\Sigma }}|s_n\left(q\right)-s_1(q+i)|,i=d_{(n-1)k}-L,...,d_{(n-1)k}+L$计算d_nk＝arg(min(sad(i))，其中i＝d_(n‑1)k‑L、…、d_(n‑1)k+L，即找出sad(i)的最小值出现的位置d，所述信号片段s_nk相对s₁的时移d_nk＝d。