动态超声彩色血流图像的计算机仿真.pdfVIP

第44卷 第 1期 复 旦 学报 (自然科学版) Vo1．44No．1 2005年 2月 JournalofFudanUniversity(NaturalScience) Feb．2005 文章编号 ：0427—7104(2005)01—0026—05 动态超声彩色血流图像的计算机仿真 方 昕，汪源源，王威琪 (复旦大学 电子工程系，上海 200433) 摘 要 ：计算机仿真超声彩色血流图像，可以为各种超声彩色血流 图像处理方法 的性能 比较提供有效的信号 源 ．在仿真某一时刻静态超声彩色血流 图像 的基础上 ，根据给定的一个心动周期 的最大血流速度曲线 ，计算 出 血管内的血流速度分布随时间的变化 ，对一个心动周期划分 24时相 ，仿真 出动态超声彩色血流图像 ． 关键词 ：超声 ；彩色血流成像 ；动态 图像 ；声脉冲响应 中图分类号：TB559 文献标识码 ：A 目前 ，B型超声成像技术、超声多普勒血流检测技术和超声彩色血流成像技术是医学超声应用的几个 重要方面．B型超声成像技术能够实现软组织解剖的动态显示…，超声多普勒技术能够无损检测血管中 的血流速度[，超声彩色血流成像技术通过图像直观、动态地显示出血管内的血流速度及其分布[． 超声彩色血流图像计算机仿真方法的建立 ，可以为各种超声彩色血流图像处理方法的性能比较提供 方便而有效的信号源．目前的仿真方法只能得到某个时刻静态的彩色血流图像 ，对于观察血流速度变化 的意义并不大．本文根据最大频率与血液速度的关系，通过一个心动周期 内的最大频率曲线计算出血流 速度及其在血管内的分布随时间的变化 ，进而仿真出24帧 ／J动周期的彩色血流图像，再用动态显示的方 法得到动态超声彩色血流图像． 1 原理和方法 1．1 声脉冲响应方法计算血流的回波信号 将体内看成是有许多随机分布点散射体的三维空间，其中红血球被视为运动的点散射体 ，其他组织被 视为静止的点散射体．换能器接收的超声回波信号看成是由每一个散射体的回波信号累加而成，而每一 个散射体的回波信号由声脉响应方法[4]计算出． 定义空间某一点处的声脉冲响应 ()： f dS； (1) 式中 ，是散射体的空间位置向量 ，C是超声波在体内的传播速度，为冲激函数 ，dS为接收回波信号换能 器探头的增量面积元，，是散射体的空间位置，到换能器增量面积元 dS的距离． 换能器得到的输出电压 V(，，t)=五f (t)×h(，，t)×h，(，，t)； (2) 式中h(，，t)和 h，(，，t)分别是换能器发送和接收的声脉冲响应 (本文中设二者相同，设为 h(，，t))， (t)是发射的脉冲(一般取高斯脉冲)，走为散射体的散射系数． 将(2)式看作单个散射体的回波信号，则总的回波信号只要将所有散射体的回波信号叠加即可， 收稿 日期 ：2003—11-10 基金项 目：教育部优秀年轻教师基金跟踪计划资助项 目(001704) 作者简介：方 ~(1981--)，男，硕士研究生；通讯联系人汪源源教授 ，博士生导师；Email：yywang@Iudan．edu．ca 第 1期 方 昕等：动态超声彩色血流图像的计算机仿真 27 V(r，t)=Pe(t)×∑kih(，t)×hr(，，t)； (3) l 1 式 中N是散射体总数． 仿真中需要将回波信号在空间和时间上离散化 ．选用如图 1所示的仿真空间模型和圆形换能器 ，将 换能器焦平面处某一径向上的散射体按一定间隔(Ar)均匀排列．对这些散射体 ，(2)式的离散化计算公 式为 ： (gAr，At)=hg(gAr，nAt)×hg(gA