[发明专利]一种利用造影微泡测量血流速度的方法

说明书

技术领域

本发明涉及一种利用造影微泡测量血流速度的方法，属于超声图像处理技术领域。

背景技术

血管中的血流速度是人体健康的重要参数和指标，可以反映很多疾病发病的中间环节或前期征兆。因此，准确地测量血管中的血流速度，有助于医生对病情进行更好的分析，具有重要的临床价值与意义。

随着超声造影剂的发展，二次谐波成像技术得到了很好的运用。充有惰性气体的油脂、聚合物或者生化蛋白制成的微泡具有无创性。利用造影微泡能产生强可视化二次谐波信号而制成的谐波或次谐波成像仪器进行探测，可以得到效果理想的目标结构。目前已有报道采用超声造影微泡技术测量血流速度，该方法是将造影微泡注入到血管中；用超声显像仪记录微泡运动的视频并保存；对运动视频图像进行分析，得到微泡在不同帧超声图像中的位置：即将视频分成一帧一帧的图像，然后在每帧图像中，通过人工判断每一帧中的微泡的位置并记录下该位置的坐标；最终，根据这些位置信息以及帧率计算出微泡运动的速度。这种通过人工判断获取微泡位置然后计算速度的做法，局限性非常大：当图像帧数比较多时，一帧一帧去找，非常繁琐；而且因为是通过肉眼判断和确定运动微泡的位置，主观因素很大，不同的人来做得出的结果有明显差异，可重复性差。所以说设计一种利用超声造影微泡测量血流速度的自动化方法是十分必要的。

本发明针对人工测量方法的繁琐与重复性差的缺陷，提出了一种通过超声微泡测量血流速度的自动化方法，同时针对测量在血管中注入造影微泡不止一个的特点，可区分出多个造影微泡并分别对其进行自动跟踪，方法具有较高的准确性。至今为止，未见有关于利用超声微泡测量血流速度的自动化方法的报道出现。

发明内容

本发明所要解决的技术问题在于提供一种能够自动跟踪超声视频图像中的运动微泡，从而计算血流速度的方法，由该方法获得的实验结果跟手动跟踪测量获得的结果非常接近，可以取代手动测量。

本发明的技术方案，包含以下步骤：

（1）向血管中注入造影微泡，记录微泡运动的超声视频图像；

（2）对各帧超声图像进行去噪预处理，然后用运动目标分割算法分割出运动微泡；在初始帧中选定待跟踪微泡的模板，并用微泡的的初始位置初始化kalman滤波器；

（3）计算微泡模板的统计特征直方图，在后一帧图像中用kalman滤波器预测微泡的位置，在预测位置邻域内，依次选取跟微泡模板同样大小的区域，计算其LBP特征直方图，并用卡方法进行直方图匹配，匹配结果最佳的位置即为微泡在该帧图像中的位置；利用相邻两帧图像中目标微泡的位置信息及视频帧率计算微泡的速度，也即得到血流的速度；

（4）根据步骤（3）中最佳匹配点的位置（也就是下一次预测中的初始位置）更新kalman滤波器，依次读取余下的图像帧，转入步骤（2），直至处理完所有图像。

上述步骤（3）中，所述邻域为模板面积16倍大小。

下面进一步详细介绍本发明方法。

本发明上述步骤（2）中的预处理方法是经过本发明实践分析后选取的、能很好适应本发明中所处理的超声图像的各项异性扩散滤波(anisotropic diffusion)方法。各向异性扩散首先是由Perona和Malik提出的，其原理如下：

∂I∂t=div[c(|▿|)·▿I]]]>

I(t＝0)＝I₀

其中div是散度算子，c(·)是系数因子。上述方程模拟了热扩散方程，方程中扩散系数是位置的函数，随着位置的变化而变化，不同的方向亦是如此。由于我们的目的是去除噪声，保留边缘信息，所以在边缘处方向的法线方向，扩散系数应该是一个小量，而在其它区域和边缘方向则就是一个较大量。扩散系数的通常形式为: