图像数值分析在人工机械心脏瓣膜运动特征中的应用
- 格式:pdf
- 大小:959.34 KB
- 文档页数:6
第15卷第11期2015年4月 1671—1815(2015)11—0031—06 科学技术与工程
Science Technology and Engineering Vo1.15 No.11 Apr.2015
⑥2015 Sci.Tech.Engrg.
医药卫生
图像数值分析在人工机械心脏瓣膜 运动特征中的应用
周 丰 崔园园 陈 颖 游天雪 黄楠 (西南交通大学材料先进技术教育部重点实验室,材料科学与工程学院,成都610031)
摘要动态图像测量技术是一种非接触式、动态、精确的检测运动物体特性的分析方式。人工心脏瓣膜体外动态开闭性能 测试,依靠高速图像摄影技术捕捉瓣膜的细微运动,快速图像和精准的时间记录,将瓣膜运动转化为图像信息,分析瓣膜及各 组件运动时的相对位移和速度变化,有效地解析瓣膜运动状况,为评价瓣膜性能提供参考。在人工心脏瓣膜脉动流平台上, 选择双叶型机械心脏瓣膜在正常心动条件下,采用高速摄影技术通过标记点跟踪瓣膜瓣叶运动,记录瓣膜的运动轨迹及运动 状态;通过图像数值分析瓣膜开启和关闭的运动情况,将瓣膜运动过程离散为多步静止状态的图像,离散后的图像数据导入 Matlab数值分析,解析瓣膜中瓣叶的转动速度及角加速度等运动状态。以此改善瓣膜结构与流场之间的相互作用,改进瓣膜 使用性能,结合瓣膜临床病症特征,解决瓣膜类相关并发症。 关键词 人工机械心脏瓣膜 高速摄像 图像数值分析 Matlab 中图法分类号PO18.11; 文献标志码A
图像测量技术是以光学为基础,将光电技术、计 算机技术、图像处理等现代科学技术融为一体的检 测技术。实时采集物体与景物的相对运动图像,通 过图像分析转化为物体的运动状态值。特点是短时 间内完成对高速目标的快速、多次采样,再将目标过 程放慢,得到清晰、准确的运动过程的动态或静态图 像。具有实时目标捕获、快速记录、清晰直观等 优点。 高速摄像系统将高速运动目标通过自身的发光 或反射光,光经物镜成像后,落在光电成像器件的感 光面上,光电元件会对成像面上的目标像快速响应, 根据感光面上目标像的光能量分布,在各采样点及 像素点产生响应的电荷包,完成图像的光电转换,即 光信号转化为图像信息。带有图像信息的电信号迅 速转移存储到临时寄存器中。信号经信号处理(去 噪)后传输至电脑硬盘,由电脑对图像进行读取输 出,并由Matlab分析软件对获取图像信息进行分析 处理。整个高速成像系统过程包括光学成像、光电 成像、信号传输、控制、图像存储、处理等部分。 采用高速摄像对人工机械心脏瓣膜测量分析试 验,血液流动为非牛顿流体高速流场,对于此类速度 测试分析方法,包括粒子影像速度法(particle image velocimetry,PIV)¨ ,多普勒超声法(Doppler ultra一 2014年12月05日收到 第一作者简介:周丰,男。E-mail:34030207@qq.corn。 sonic method)L2]等检测评价瓣膜性能。然而,这些 方法受实验条件限制,需要昂贵的实验设备和准备 时间长,不易应用于瓣膜前期研发试验。 在人工心脏瓣膜脉动流试验平台(TH 200) 上,测试瓣膜在75 ̄/min心动周期脉动条件下瓣 膜运动,图i表示瓣膜在一个循环周期内经过瓣膜 的血流流量(Y轴)随时间( 轴)变化波形曲线。
/ 一 £
图1瓣膜流量波形 Fig.1 Flow waveforrn through valve
循环过程中四个运动行程阶段 .4]: 开启保持阶段——瓣膜从完全开启到开始关闭 之间的时间段。 关闭阶段——瓣膜从80%最大开口面积位置 到瓣膜完全关闭时问段,也是瓣膜不稳定扇动阶段。 开启阶段——瓣膜从开始开启到完全开启之间 时间段。 关闭保持阶段——瓣膜完全关闭到下一阶段瓣 膜开启之间时间段。 32 科学技术与工程 15卷 开启阶段和关闭阶段是瓣膜的运动阶段,其它 两个阶段为非运动阶段。
1 实验方法 在人工心脏瓣膜体外脉动流试验平台 上,测 试机械心脏瓣膜体外动态开闭过程,瓣膜的运动图 像经观察窗投影到平面镜折射后,补光灯增强背景 光,由GigaView高速摄像仪(South SVSI,USA)平行 对焦录制瓣膜开闭运动状态,将光信号转化为图像 信息存储于摄影仪高速数据缓存区,通过数据线传 输到计算机硬盘存储。 高速摄像测量是将实验物体反射的光信号转化 为电信号,由摄像机快速、多次采样,收集存储,并以 常规速度播放或对采样图片逐个分析。高速摄像可 以达到1O ~10 帧/s的记录速度,以高频率实时记 录瓣膜整体运动动态,将动态图像分为多个静态连 贯的图像按照精确的间隔时间播放。 将采集的图像数据通过像素点转换为瓣膜运动 位移项,图像二值化处理 J,将图像像素取1或0离 散值来表示,用中值滤波或均值滤波来减少图片噪 音,将取亮点所在图像范围域内,自动拾取坐标位置 点,手动修正少数错误位置点坐标,将数值导人,分 析亮点运动轨迹。根据轨迹就能算出瓣叶的运动状 况,如位移及转动角度、角速度、角加速等。 图像处理初始化后,图像数值化处理通过Mat— lab软件分析,Matlab是一个交互式系统,基于数字 单元无维度阵列。瓣膜运动过程图像,转化为参数 图表…,通过Matlab数据处理瓣膜标记点运动,分 析瓣叶的瞬时转动变化,求解瓣叶转动的开启速度 数值曲线,分析瓣膜的运动规律,包括瓣膜的线位移 速度,角速度等瓣膜运动特征,评价和改进瓣膜结构 设计。 图像测量系统包括硬件和软件两个部分,硬件 系统包括瓣膜运动测试平台,高速摄像仪、计算机、 补光灯。软件系统功能包括控制图像(测试台)、采 集和转化处理图像、计算分析图像。 实验现场如图2所示,采集图像通过直角反光 镜折射到高速摄影仪,用摄影灯对准实物照亮,增强 光线对比度。系统整体包括三个部分:①疲劳实验 平台,用于模拟在体内血流往复运动瓣膜模式;②高 速摄影仪,瓣膜运动过程经过光线折射,光电转换信 息将数值图像存储于高速硬盘,通过光纤数据线传 递给Pc计算机储存;③计算机分析,将图像离散的
一系列连续图像,根据瓣膜上标记点位置参数变化 分析瓣膜运动,并计算表征出瓣膜运动状态参数。 瓣膜安装到疲劳测试仪活塞处,瓣膜作为阻流
体安装在圆柱直管内作垂直往复运动,其运动速率 和频率由下部分控制器控制,选择瓣膜运动频率等 于平均心率75 ̄,Jmin(0.8 s/次),待瓣膜运动达到 稳定后,高速摄影仪通过平面镜反射后开始记录瓣 膜垂直投影的影像,实时录制瓣膜运动。
脉动流实验 平台
镜中图像 图2 瓣膜实验系统示意图 Fig.2 The valvular experimental system
试验平台模拟瓣膜真实的运动状态,瓣膜固 定在管道内夹持,由马达驱动和垂直往复运动。 控制台设备瓣膜管腔用有机玻璃制造,可直接观 察流体和瓣膜移动,并通过下端透明窗口观察瓣 膜开闭运动情况。中介流体为纯净水,其密度、雷 诺数(Re)和斯特劳哈尔数(St)应接近于血液特 性。随着活塞的垂直往复运动,瓣膜受流体作用 而相应的开闭。 双叶型机械心脏瓣膜最大开启角和闭合角各自 是85。和28。,瓣膜周期运动时,瓣叶相对于流体冲 击推动反向转动,当瓣膜向下运动时,瓣膜关闭;瓣 膜向上运动时,瓣膜随之迅速开启。
2 实验过程 瓣膜固定于实验装置活塞中,在电机的驱动下, 瓣膜在管腔内相对流体沿轴线方向作竖直上下周期 往复运动,模拟生理条件下血流对瓣膜的作用,流体 流经瓣膜促使其开闭。同时管腔上下端相通,使得 瓣膜在开闭过程中,瓣膜上下端始终保持恒定的压 力差。 瓣膜开闭运动由高速摄影仪捕捉记录,并存储 于高速计算机。标记在瓣环瓣叶边缘上的记录标记 点,有利于观察瓣膜位置运动变化情况,使用Matlab 统计图像中点位置坐标,将瓣膜标记点的位置变化 转化为瓣膜同瓣叶位移的相对变化及瓣叶开闭角 度、角速度变化情况等。瓣膜动态体外模拟实验采 集分析图像过程如下: 2.1点坐标位置标定 为了测试瓣膜运动,选择瓣膜特征位置标记瓣 11期 周丰,等:图像数值分析在人工机械心脏瓣膜运动特征中的应用1 35 叶运动,使瓣叶开启时曲线性波动。 3)标记点坐标位置变化造成瓣叶角度波动。 因标记点坐标位置是由软件程序自动提取,虽然误 差较小,但实验过程中,瓣叶随瓣膜的运动,会造成 瓣膜影像的变化,因而使得计算提取瓣膜坐标值出 现波动。因此通过数据曲线平滑拟合后更能真实反 映瓣膜开闭状态的行为。
3结果及讨论 通过高速摄影仪和数值图像处理技术对瓣膜开 闭行为进行研究,试验采用高速摄像技术实现非接 触测试分析机械瓣膜运动,为评价心瓣开闭性能提 供了有效的分析方法。 在研究中,数据获取和数据进程靠高速图像、非 接触测量方式直接观察人工机械瓣膜运动(包括开 启角,开启和关闭速度,有效开口面积)。测量方式 特点如下: (1)无损耗:人工瓣膜在流动腔内模拟心脏瓣 膜的开闭运动过程,由于采用无接触测量获取数据, 瓣膜在此状态下运动影像转换为高速图像,瓣膜不 会在实验条件下损坏,在瓣膜上的标记容易被区分, 不影响瓣膜运动; (2)准确可靠:检测方法能准确捕捉瓣膜操作 的特征参数和相关运动过程,且试验具有可重复性 和对外界因素无影响; (3)经济性和效率:瓣膜测试效率高,过程稳定 可靠,能有效重复试验操作。
参考文献 1 Chandran K B,Cabell G N,Khalighi B,et a1.Laser anemometry
measurements of pulsatile flow past ao ̄ic valve prostheses.Joumal of biomechanics,1983;16(10):865—873 2 Satomura S.Ultrasonic Doppler method for the inspection of cardiac functions.The Journal of the Acoustical Society of America,1957; 29(11):1181—1185 3 Simon H A,Leo H L,Carberry J,et a1.Comparison of the hinge flow fields of two bileanet mechanical heart valyes under aortic and mitral conditions.Annals of Biomedical Engineering,2004;32 (12):1607—1617 4陈琛.人工心瓣疲劳试验原理初步研究.成都:电子科技大学. 1991 Chen Chen.A preliminary study of artificial heart valves fatigue test principle.Chengdu:University of Electronic Science and Technology of China.1991 5周丰,崔园园,武亮亮,等.新型人工曲面机械双叶瓣体外脉 动流实验.功能材料,2014;45(20):20090--20094 Zhou Feng,Cui Yuanyuan,Wu Liangliang.A new artificial mechani— cal curved surface bi—leaflet cardiac valve:in vitro result.Journal of Functional Materials,2014;45(2O):20090--20094 6赵永志,彭国华.一种有效的图像二值化方法.科学技术与工 程,2007;7(1):139—141 Zhao Yongzhi,Peng Guohus.An efficient approach to image binariza- tion.Science Technology and Engineering,2007;7(1):139—14】 7 Gonzalez R C,Woods R E,Eddins S L.Digital image processing using Matlab.India:Pearson Education,2004 8 Berne R M.Cardiovascular physiology.Annual Review of Physiolo- gY,1981;43(1):357—058 9 Nobili M.Morbiducci U,Ponzini R.et a1.Numerical simulation of the dynamics of a bileaflet prosthetic heart valve using a fluid--struc-- ture interaction approach.Journal of Biomechanics,2008;41(1 1): 2539--2550 10 Dapeng C,Jianhai Z.A comparative study of the opening and clos— ing process of two types of mechanical heart valves using ALE finite element method.Applied Mathematics and Mechanics,1996;17 (4):299—3O8