1.一种冠脉血管的运动参数获取方法,其特征在于,包括:
获取冠脉血管的冠脉造影血管序列;
基于所述冠脉造影血管序列,获取所述冠脉血管对应的基频心脏搏动信号的周期,其中,所述基频心脏搏动信号为所述冠脉血管对应的心脏搏动信号中频率最小的心脏搏动信号,所述基频心脏搏动信号的周期用于表征所述冠脉血管的运动周期;
基于所述冠脉造影血管序列,建立三维血管树结构,所述三维血管树结构包含所述冠脉血管在不同时刻的空间位置;
基于所述冠脉血管在不同时刻的空间位置、以及所述基频心脏搏动信号的周期,通过第一映射关系,确定第二映射关系,其中,所述第一映射关系表征了冠脉血管在不同时刻的空间位置、基频心脏搏动信号的周期、冠脉血管对应的光流场的能量与心脏搏动信号对应的谐波系数之间的映射关系,所述第二映射关系为所述第一映射关系中的冠脉血管在不同时刻的空间位置和基频心脏搏动信号的周期为已知的情况下的第一映射关系;
基于所述第二映射关系,确定第一参数值,其中,所述第一参数值为最小的光流场的能量对应的心脏搏动信号的谐波系数;
基于所述第一参数值,确定所述冠脉血管的运动参数。
2.根据权利要求1所述的方法,其特征在于,所述第一映射关系包括第三映射关系和第四映射关系,所述第三映射关系为所述冠脉血管在不同时刻的空间位置、所述基频心脏搏动信号的周期、冠脉血管的每相邻两个时刻的空间位置的偏移量与心脏搏动信号对应的谐波系数之间的映射关系,所述第四映射关系为所述冠脉血管在不同时刻的空间位置、所述基频心脏搏动信号的周期、冠脉血管对应的光流场的能量与冠脉血管的相邻两个时刻的空间位置的偏移量之间的映射关系;
所述基于所述冠脉血管在不同时刻的空间位置、以及所述基频心脏搏动信号的周期,通过第一映射关系,确定第二映射关系,包括:根据所述冠脉血管在不同时刻的空间位置、所述基频心脏搏动信号的周期和所述第三映射关系,确定第五映射关系,所述第五映射关系为在所述冠脉血管在不同时刻的空间位置和所述基频心脏搏动信号的周期已知的情况下的第三映射关系;
根据所述冠脉血管在不同时刻的空间位置、所述基频心脏搏动信号的周期和所述第四映射关系,确定第六映射关系,所述第六映射关系为在所述冠脉血管在不同时刻的空间位置和所述基频心脏搏动信号的周期已知的情况下的第四映射关系;
根据各所述偏移量、所述第五映射关系和所述第六映射关系,确定所述第二映射关系。
3.根据权利要求2所述的方法,其特征在于,所述基于所述第一参数值,确定所述冠脉血管的运动参数,包括:将所述第一参数值输入至所述第五映射关系,得到所述冠脉血管的每相邻两个时刻的空间位置的偏移量;
基于所述冠脉血管的每相邻两个时刻的空间位置的偏移量,计算所述冠脉血管的三维运动的幅度和相位,将所述三维运动的幅度和相位确定为所述冠脉血管的运动参数。
4.根据权利要求2所述的方法,其特征在于,所述第三映射关系通过以下公式表示:其中,(x,y,z)表示冠脉血管的空间位置,t表示时刻,u(x,y,z,t)、v(x,y,z,t)、o(x,y,z,t)分别表示冠脉血管在t时刻的空间位置和在t时刻相邻的时刻的空间位置之间在x、y、z三个方向上的偏移量,am,bm,cm,dm,em,fm表示第m个心脏搏动信号对应的谐波系数,M表示冠脉血管对应的心脏搏动信号的个数,角频率ωm=2πm/T,T表示基频心脏搏动信号的周期,π表示圆周率。
5.根据权利要求2所述的方法,其特征在于,所述第四映射关系通过以下公式表示:E(u,v,o)=E1(u,v,o)+λE2(u,v,o)
其中,E(u,v,o)表示光流场的能量,I(x,y,z,t)表示三维血管树结构中t时刻处于(x,y,z)的空间位置的像素值,(x,y,z)表示空间位置,t表示时刻,u、v、o分别表示冠脉血管相邻两个时刻的空间位置之间在x、y、z三个方向上的偏移量,T表示基频心脏搏动信号的周期,λ表示用于平衡E1(u,v,o)和E2(u,v,o)的能量项参数,D2表示光流场中的坐标范围,▽表示一阶差分运算。
6.根据权利要求1至5中任一项所述的方法,其特征在于,所述基于所述第二映射关系,确定第一参数值,包括:基于所述第二映射关系,利用欧拉‑拉格朗日方程建立用于计算所述最小的光流场的能量对应的谐波系数的参数值的线性方程;
利用迭代的雅克比方法求解所述线性方程,得到所述第一参数值。
7.根据权利要求1至5中任一项所述的方法,其特征在于,所述基于所述冠脉造影血管序列,建立三维血管树结构,包括:对所述冠脉造影血管序列进行去噪;
对去噪后的所述冠脉造影血管序列进行图像增强;
基于图像增强后的所述冠脉造影血管序列,确定中心线图像,其中,所述中心线图像用于表征冠脉血管树的结构;
对所述中心线图像进行运动补偿,对运动补偿后的所述中心线图像中的冠脉血管树的分叉点和冠脉血管段进行标记;
基于标记后的所述中心线图像,建立所述三维血管树结构。
8.一种冠脉血管的运动参数获取装置,其特征在于,包括:
第一获取模块,用于获取冠脉血管的冠脉造影血管序列;
第二获取模块,用于基于所述冠脉造影血管序列,获取所述冠脉血管对应的基频心脏搏动信号的周期,其中,所述基频心脏搏动信号为所述冠脉血管对应的心脏搏动信号中频率最小的心脏搏动信号,所述基频心脏搏动信号的周期用于表征所述冠脉血管的运动周期;
第一处理模块,用于基于所述冠脉造影血管序列,建立三维血管树结构,所述三维血管树结构包含所述冠脉血管在不同时刻的空间位置;
第二处理模块,用于基于所述冠脉血管在不同时刻的空间位置、以及所述基频心脏搏动信号的周期,通过第一映射关系,确定第二映射关系,其中,所述第一映射关系表征了冠脉血管在不同时刻的空间位置、基频心脏搏动信号的周期、冠脉血管对应的光流场的能量与心脏搏动信号对应的谐波系数之间的映射关系,所述第二映射关系为所述第一映射关系中的冠脉血管在不同时刻的空间位置和基频心脏搏动信号的周期为已知的情况下的第一映射关系;
第三处理模块,用于基于所述第二映射关系,确定第一参数值,其中,所述第一参数值为最小的光流场的能量对应的心脏搏动信号的谐波系数;
参数获取模块,用于基于所述第一参数值,确定所述冠脉血管的运动参数。
9.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质中存储有指令,当所述指令在终端设备上运行时,使得所述终端设备执行如权利要求1至7中任一项所述的冠脉血管的运动参数获取方法的步骤。
10.一种电子设备,包括存储器、处理器及存储在所述存储器上并在所述处理器上运行的程序,其特征在于,所述处理器执行所述程序时实现如权利要求1至7中任一项所述的冠脉血管的运动参数获取方法的步骤。