结合有限元和非刚体运动估计的左心室应变分析  

Strain Analysis of Left Ventricle by Integrating Finite Element Model and Non-rigid Motion Estimation

作　　者：刘复昌[1] 尤建洁[1] 郭亮[2] 王平安[3] 韦志辉[1] 夏德深[1] 

机构地区：[1]南京理工大学计算机科学与技术学院,南京210094 [2]苏州大学附属第一医院放射科,苏州215006 [3]香港中文大学计算机科学与工程学系

出　　处：《计算机辅助设计与图形学学报》2007年第8期1028-1033,共6页Journal of Computer-Aided Design & Computer Graphics

基　　金：江苏省自然科学基金(BK2006704-2);香港特区政府研究资助局资助项目(CUHK/4180/01E;CUHK1/00C)

摘　　要：提出一种直观地分析左心室在心脏收缩期的应变情况的方法.首先利用心脏核磁共振图像数据对左心室进行三维实体建模,再进行有限元单元划分;然后,使用非刚体运动估计算法——单位法线及高斯曲率混合算法,在左心室收缩期相邻时刻的有限元模型表面上寻找最相似单元结点作为对应单元结点,计算各对应结点间的位移,并用α剪枝平均对位移过大的可疑结点进行处理;最后,利用弹塑性力学中应变与位移间的几何关系求解左心室表面各单元结点的应变,获得位移和应变的彩色云图.实验结果表明,该方法能直观、有效地反映左心室表面在整个心脏收缩期的形变情况和应变分布.A strain analysis approach of the left ventricle （LV） is proposed. In this approach, we can intuitively analyze the strain of the LV during systole. First, we build a 3D solid model by utilizing the cardiac magnetic resonance （MR） images, and then construct a 3D finite element model. Second, we calculate the correspondence points of element nodes on the finite element model surface by utilizing hybrid unit normal/Gaussian curvature （UNGCK） algorithm. UNGCK algorithm is a non-rigid motion estimation technique. Lastly, we use α-trimmed means to correct the displacement of the suspicious node. We estimate the strain of the LV by using the knowledge of elasticity and plasticity, and then acquire color cloud pictures of the deformation and strain distribution. Our experiments indicate that the strain analysis approach of the LV can intuitively and effectively reflect the strain distribution and deformation trends of surfaces of the LV during systole.

关 键 词：有限元模型 单位法线及高斯曲率混合算法 左心室应变分析 α剪枝平均 心脏核磁共振图像 

分 类 号：O313.3[理学—一般力学与力学基础]