技術特征：

1.一種基于可變?yōu)V波器長度的超聲彈性成像方法，其特征在于，包括如下步驟：

S1：根據(jù)變形前后兩幀超聲射頻信號，利用數(shù)值優(yōu)化算法計算出位移場；

S2：利用濾波器長度為固定值的SG差分濾波器對位移場進行差分運算，得出初應變場；

S3：對初應變場進行差分運算求出含有異物邊界信息的應變梯度場；

S4：通過映射函數(shù)對應變梯度場的絕對值進行線性變換，得出與異物邊界信息相關的濾波器長度分布場；

S5：按照濾波器長度分布場，對步驟S1中求出的位移場進行SG差分濾波器差分運算，得出最終的應變場。

2.根據(jù)權(quán)利要求1所述的基于可變?yōu)V波器長度的超聲彈性成像方法，其特征在于，步驟S1中的數(shù)值優(yōu)化算法具體為：

假設圖像窗口在運動前后滿足仿射變換關系，運動前窗口中的任意一點P(x,y)對應于運動后窗口中點Q(x′,y′)，兩者的坐標滿足：

$\left\{\begin{array}{c}{x}^{\′}=x+\mathrm{\ξ}\left(x,y,V_{rs}\right)=x+u+\frac{\∂u}{\∂x}x+\frac{\∂u}{\∂y}y\\ y^{\′}=y+\mathrm{\η}\left(x,y,V_{rs}\right)=y+v+\frac{\∂v}{\∂x}x+\frac{\∂v}{\∂y}y\end{array}\right.$

式中，向量為窗口的移動和變形，u，v為剛體位移，為窗口的應變，即窗口的一階變形，將目標函數(shù)Ob(V_rs)最小化即可同時得到向量V_rs中的6個參數(shù)值，其中，目標函數(shù)

$Ob\left(V_{rs}\right)=\underset{y=1}{\overset{L}{\Σ}}\underset{x=1}{\overset{N}{\Σ}}\[P(x,y)-Q(x+\mathrm{\ξ}(x,y,V_{rs}),y+\mathrm{\η}(x,y,V_{rs}\left)\right)\].$

3.根據(jù)權(quán)利要求1所述的基于可變?yōu)V波器長度的超聲彈性成像方法，其特征在于，步驟S2中采用1階或2階多項式擬合的SG差分濾波器進行應變計算，假設SG差分濾波器的長度為2M+1，則1階或2階多項式擬合的SG差分濾波器的輸入輸出關系為：

4.根據(jù)權(quán)利要求1所述的基于可變?yōu)V波器長度的超聲彈性成像方法，其特征在于，步驟S4中的映射函數(shù)為：FL_var＝round{ζ[1-(1-σ)Φ]}，F(xiàn)L_var為濾波器長度，round()為四舍五入函數(shù)，式中，ζ為一常數(shù)，σ為調(diào)整因子，Φ為應變梯度歸一化函數(shù)，式中，MB為應變梯度場，MB_max為應變梯度場的最大值，MB_min為應變梯度場的最小值。

5.根據(jù)權(quán)利要求4所述的基于可變?yōu)V波器長度的超聲彈性成像方法，其特征在于，調(diào)整因子σ取值為0.05～0.95。