[发明专利]血管单位时间血流量与血流速度的计算方法

权利要求书

1.一种血管单位时间血流量的计算方法，该方法包括：

确定血管的感兴趣区域；

选择起始时间为造影剂充盈之前，计算每帧造影中感兴趣区域内总灰度值，并根据灰度值拟合出灰度随时间变化的灰度拟合曲线；

确定一预定时间区间内的最大灰度值曲线或最小灰度值曲线；该步骤进一步包括：确定第一时间点，以及以第一时间点为中心的一预定时间区间内灰度拟合曲线的最大值和最小值；所述最大灰度值曲线是以预定时间区间内灰度拟合曲线的最大值为纵坐标所做的曲线；所述最小灰度值曲线是以预定时间区间内灰度拟合曲线的最小值为纵坐标所做的曲线；

计算预定时间区间内最大灰度值曲线或最小灰度值曲线与灰度拟合曲线所包围的区域面积值；

基于该区域面积值，获取该面积值所对应的单位时间血流量。

2.一种血管血流速度的计算方法，该方法包括：

确定血管的感兴趣区域；

选择起始时间为造影剂充盈之前，计算每帧造影中感兴趣区域内总灰度值，并根据灰度值拟合出灰度随时间变化的灰度拟合曲线；

确定一预定时间区间内的最大灰度值曲线或最小灰度值曲线；该步骤进一步包括：确定第一时间点，以及以第一时间点为中心的一预定时间区间内灰度拟合曲线的最大值和最小值；所述最大灰度值曲线是以预定时间区间内灰度拟合曲线的最大值为纵坐标所做的曲线；所述最小灰度值曲线是以预定时间区间内灰度拟合曲线的最小值为纵坐标所做的曲线；

计算预定时间区间内最大灰度值曲线或最小灰度值曲线与灰度拟合曲线所包围的区域面积值；

基于该区域面积值，获取该面积值所对应的单位时间血流量；

基于所述单位时间血流量以及血管管腔面积，获得所述血管的血流速度。

3.根据权利要求1或2所述的方法，其特征在于，所述感兴趣区域包括注入造影剂的主支血管及其分支。

4.根据权利要求1或2所述的方法，其特征在于，通过目标图像跟踪检测感兴趣区域在不同心跳时刻下的位置的变化，从而获得最佳的感兴趣区域。

5.根据权利要求1或2所述的方法，其特征在于，接受血管的X射线造影影像序列，选择感兴趣区域。

6.根据权利要求1或2所述的方法，其特征在于，

当所述灰度拟合曲线变化趋势为下降时，所述第一时间点为灰度拟合曲线中灰度值下降最快的时间点；

当所述灰度拟合曲线变化趋势为上升时，所述第一时间点为灰度拟合曲线中灰度值上升最快的时间点。

7.根据权利要求6所述的方法，其特征在于，

当灰度拟合曲线变化趋势为下降时，计算灰度拟合曲线上每点的斜率，获取斜率为负值，且斜率绝对值最大的点，所述点即为灰度值下降最快的时间点；

当灰度拟合曲线变化趋势为上升时，计算灰度拟合曲线上每点的斜率，获取斜率为正值，所述斜率的值最大的点即为灰度值上升最快的时间点。

8.根据权利要求6所述的方法，其特征在于，所述区域面积值的计算过程进一步包括：

当灰度拟合曲线变化趋势为下降时，获取所述灰度拟合曲线中的第一时间点，在所述以第一时间点为中心的预定时间区间内对灰度拟合曲线进行积分，计算该预定时间区间内最大灰度值曲线与所述灰度拟合曲线所包围的区域面积值；

当灰度拟合曲线变化趋势为上升时，获取所述灰度拟合曲线中的第一时间点，在所述以第一时间点为中心的预定时间区间内对灰度拟合曲线进行积分，计算该预定时间区间内最小灰度值曲线与所述灰度拟合曲线所包围的区域面积值。

9.根据权利要求1或2中任一所述的方法，其特征在于，所述预定时间区间为整数个心动周期，所述整数大于等于1。

10.根据权利要求9所述的方法，其特征在于，所述预定时间区间为一个心动周期，包括第一时间点前后的各二分之一心动周期；

其中，所述第一时间点前后二分之一心动周期时间区间为造影剂开始充盈血管感兴趣区域后，在没有完全充盈感兴趣区域前的这一时间段。

11.根据权利要求1或2中任一所述的方法，其特征在于，通过查对应表，可获得该区域面积值所对应的单位时间血流量，所述对应表为不同面积值与单位时间不同血流量之间的对应表格。