[发明专利]改进的直径速度环方法测量升主动脉脉搏波速度的方法

权利要求书

1.改进的直径速度环方法测量升主动脉脉搏波速度PWV的方法，其特征在于，包括以下步骤：

步骤S10、利用心电图的R波波峰分别对升主动脉的M模式超声图像、D模式超声图像分割心动周期；

步骤S20、从分割好心动周期的M模式超声图像中提取升主动脉的直径波形；

步骤S30、再从分割好心动周期的D模式超声图像中提取升主动脉的血流速度波形；

步骤S40、分别获取直径波形和血流速度波形中的上升点，并利用RANSAC计算PWV值；

所述步骤S20的具体过程为：

步骤S21、在分割好心动周期的M模式超声图像中，将相邻的心动周期图像叠加，得到叠加图像；

步骤S22、通过阈值法得到步骤S21中叠加图像的血管上下壁，并对血管上下壁分别使用不同阈值处理图像；

步骤S23、对步骤S22中处理后的图像分别通过阈值法、重心法得到血管上下壁边界，并使用梯度方法滤除通过阈值法得到的血管上下壁边界的异常点；

步骤S24、调节通过阈值法得到的血管上下壁边界和通过重心法得到的血管上下壁边界的比例系数α，然后进行叠加；

步骤S25、最后使用样条平滑的方法对步骤S24中叠加图的血管上下壁进行插值平滑，血管上下壁相减得到升主动脉的直径波形；

所述步骤S30的具体过程为：

步骤S31、通过阈值法得到分割好心动周期的M模式超声图像中的最大包络；

步骤S32、再采用梯度方法滤除最大包络异常的离群点；

步骤S33、最后通过样条平滑得到一个心动周期随时间变化的血流速度波形；

所述步骤S40的具体过程为：

步骤S41、分别在直径波形和血流速度波形中采用逐点向下的线性回归，找到线性回归决定系数R²小于0.985的点分别作为直径波形的上升点、血流速度波形的上升点；

步骤S42、直径波形和速度波形按照线性回归查找的上升点对齐得到直径速度环，使用RANSAC方法拟合确定直径速度环中的无反射期中直线部分，该直线部分斜率的一半即为PWV值。

2.根据权利要求1所述的改进的直径速度环方法测量升主动脉脉搏波速度PWV的方法，其特征在于，所述步骤S31中阈值法中对一个心动周期的D模式超声图像的左1/8的像素采用4倍于中间部分的阈值，图像右2/5的像素采用3倍于中间部分的阈值。

3.根据权利要求1所述的改进的直径速度环方法测量升主动脉脉搏波速度PWV的方法，其特征在于，在步骤S20和步骤S30中分别提取若干个直径波形、血流速度波形，计算直径波形与血流速度波形所有的匹配组合，并计算出所有组合的PWV值，将所有组合的PWV值排序并滤除前后25%的可能异常值，计算剩下值的平均值作为最终的PWV值。