[发明专利]一种超声心动图粒子图像测速速度场修正方法

摘要

一种超声心动图粒子图像测速速度场修正方法，它涉及速度场修正方法。本发明解决了解决现有超声医学图像进行噪声抑制的方法中空域滤波法会引起图像模糊和细节丢失、小波变换法对超声图像的乘性噪声效果不佳及各向异性扩散降噪法算法难度大的技术问题。方法：1、读入t时刻的超声心动图的初始速度场；2、计算各速度矢量对应的协方差矩阵；3、对各速度矢量进行各向异性扩散；4、各速度矢量进行中值滤波；5、各速度矢量移流计算；6、判断是否是最后的循环；如果是，则执行步骤9的结束；如果不是，则执行步骤7的计算下一循环的初始速度场；步骤8、令循环次数加1后，返回至步骤2；本方法能保持速度场细节，可对各种粒子图像测速场进行处理。

主权项

1.一种超声心动图粒子图像测速速度场修正方法，其特征在于超声心动图粒子图像测速速度场修正方法按以下步骤进行：步骤一：用超声仪器检测超声心动图，然后读入用粒子图像测速法测量的t时刻的超声心动图的初始速度场该初始速度场u₁由m×n个网格组成，每个网格中有一个速度矢量其中1≤i≤m，1≤j≤n；用k表示循环次数，K为循环次数的最大值，K为3～5；初始化k＝1；步骤二：计算初始速度场中2≤i≤m-1且2≤j≤n-1的各个速度矢量对应的协方差矩阵S_c，其中，Sc=ψxxψxyψyxψyy,]]>其ψxx=∂2ux∂x2,]]>ψyy=∂2uy∂x2,]]>ψxy=∂2ux∂x∂y,]]>ψyx=∂2uy∂y∂x;]]>u_x为速度矢量在x轴方向上的分量，u_y为速度矢量在y轴方向上的分量；步骤三：以步骤二中得到的协方差矩阵作为扩散核，对初始速度场中2≤i≤m-1且2≤j≤n-1的各个速度矢量分别按u2k(i,j)=Σ-1≤α≤1-1≤β≤1u1k(i+α,j+β)exp[-(α,β)T·Sc-1·(α,β)]Σ-1≤α≤1-1≤β≤1exp[-(α,β)T·Sc-1·(α,β)]]]>进行各向异性扩散，得到速度场α、β为整数；S_c指步骤二中计算得到的对应的扩散核，即协方差矩阵；为S_c的逆阵；(α，β)^T为向量(α，β)的转置；步骤四、对步骤三得到的速度场中各速度矢量实施中值滤波，中值滤波的范围是当前计算点的3×3邻域，得到速度场步骤五：将步骤四的速度场中各速度矢量进行移流计算，反推出Δt时刻前该速度矢量所在的位置以及在该位置的速度值再用这个速度矢量值替代得到末态速度矢量从而得到末态速度场其中Δt为时间间隔，Δt取1/300秒～1/30000秒；步骤六：判断k＝K是否成立，如果否，执行步骤七；如果是，执行步骤九；步骤七：按计算第k+1次循环的初始速度场ε为比例系数，取0.05～0.15；步骤八；令k加1，执行步骤二；步骤九：结束。