[发明专利]三维心肌形变应变计算方法

权利要求书

1.一种三维心肌形变应变计算方法，其特征在于：在完成加标记的心脏核磁共振影像序列分割处理的基础上，首先进行数据后处理，对心肌按时间、空间坐标划分多个局部区间，拟合局部位移场，然后对位移场进行计算，再通过牛顿迭代法求解心肌质点的运动和心肌形变参数，最后使用非线性局部插补技术计算心肌应变参数。

2.根据权利要求1所述的三维心肌形变应变计算方法，其特征在于：含有短轴方向和长轴方向的影像序列作分割数据后处理，即来自对加标记的心脏核磁共振影像Tagged MRI序列的分割结果，其短轴方向的tag线为网格型，间隔为ds；其长轴方向的tag线为平行线，间隔为dl；影像序列的n个采样时间点T＝{t₀，t₁，...t_n-1)，其中t₀为舒张末期作为参考时刻，短轴影像面间隔为dz；长轴影像面间隔为dθ，在对Tagged MRI序列进行分割获得各层图像tag线的网格交点和tag线与心肌内外轮廓的交点，经分割处理后得到所有交点的编码i、交点参考时刻坐标p_i和t_j时刻的位移分量Δx(t_j)及Δy(t_j)的短轴分割处理结果，其中编码i＝1，2，...m_s；p_i为x_i0＝x(t₀)，y_i0＝y(t₀)，z_i0＝z(t₀)；Δx(t_j)＝x(t_j)-x(t₀)，Δy(t_j)＝y(t_j)-y(t₀)，j＝1，2...n-1；以及包括所有交点的编码i、交点参考时刻极坐标q_i和t_j时刻的位移Δz(t_j)的长轴分割结果，其中编码i＝1，2，...m_l·q_i为θ_i0＝θ(t₀)，r_i0＝r(t₀)，z_i0＝z(t₀)；θ(t₀)＝arctg(y(t₀)/x(t₀))，r(t0)=x2(t0)+y2(t0),]]>Δz(t_j)＝z(t_j)-z(t₀)，j＝1，2...n-1。