基于ＤＦＴ的全方向Ｍ型超声心动图运动曲线的频域滤波.docVIP

基于ＤＦＴ的全方向Ｍ型超声心动图运动曲线的频域滤波 　　摘 要：介绍一种基于离散傅里叶变换(DFT)的全方向M型超声心动图运动曲线频域滤波技术。全方向M型心动图是利用B超获取的心脏各结构的序列图像，沿各个方向提取的心脏某结构部位灰度点位置变化在时间轴上展开得到的运动轨迹，其边缘线即全方向M型心动图边界代表了了该灰度点的运动方程。依据其特点及要求，通过对不同区域的频域滤波效果比较，提出了分段滤波方法，实现了消除噪声干扰及减弱阶梯效应的影响的目的，从而保证了对全方向M型心动图的动态信息进一步分析的可靠性。?? 　　关键词：全方向M型超声心动图；离散傅里叶变换(DFT)；分段滤波；曲线频域滤波?? 　　中图分类号：TN919.8文献标识码：B 　　文章编号：1004-373X(2008)07-179-03?オ? 　　 　　Frequency Domain Filtering of Motion Curve in Omni-directional M-mode ?? 　　Echocardiography Images Based on DFT?? 　　ZHONG Kehua,LIN Qiang?? 　　(Fuzhou University,Fuzhou,350002,China)?オ? 　　Abstract：The paper presents a method of motion curve frequency domain filtering in Omni-directional M-mode echocardiography images based on Discrete Fourier Transform (DFT).According to the particularity and requirement of images edge,by comparing the frequency range result of different region,this article presents the method of region filtering,which can erase the interfering noise and weaken the influence of quantization effect.Therefore,it can offer reliability for Omni-directional M-mode echocardiography images to be further analyzed.?? 　　Keywords：Omni-directional M-mode echocardiography;Discrete Fourier Transform (DFT);region filtering;curve frequency domain filtering?? 　　 　　1 引 言?? 　　全方向M型超声心动图(图1)是由二维超声序列图像中采样线上的灰度(位置)信息随扫描展开而重建出的波型图[1]。全方向M型超声成像技术简便易行，能比较清晰显示局部心血管组织结构快速运动变化的细微状态。提出全方向M型心动图的运动曲线，有利于准确定量检测局部运动的幅度、速度及其与心动周期有关的室壁厚度、心血管腔内径，在心脏局部和整体功能的定量评估方面有重要的应用价值。 　　 　　但由于“数字化”的因素，图像由各个像素点组成，所得到的曲线并不是连续变化的曲线。如图2所示，对边界图2(a)进行微分求速度，理想情况下得到的速度应为图2(c)??曲线，但实际得到的速度为图2(b)曲线，其中存在许多突变点。同样，提取出的边界线由于“数字化”及分辨率大小的影响，阶梯效应十分显著(如图2所示)，当时间间隔取得很短时，所得到的速度存在许多突变点，时间间隔取得越短，一阶微分得到的速度则越大，甚至趋于无穷大，在物理学上相当于心脏某部位所受到的力为无穷大，而这对于心肌运动来说是不可能的，再加上噪声的干扰，产生了许多错误信息，对有用的信息造成很大的干扰，这样所提取出的数字化的边界就无法利用了。?? 　　因此，为了能够利用数字化边界对心脏某部位的运动进行更深入的分析，必须尽量克服边界中的阶梯效应与噪声干扰造成的影响。?? 　　本文针对运动曲线的傅里叶频谱分析处理，在保证边界曲线能较准确地代表波形的运动轨迹的基础上，通过分段滤波方法合理地滤除运动曲线的高频分量从而达到抑制量化噪声及一些现有医疗分析中未能利用的细节信息，有利于提高运动曲线微分得到的速度场加速度场检测数据的稳定性。?? 　　 　　2 滤波器类型的选择与截止频率的确定?? 　　一维运动曲线经过离散傅里叶变换(