多普勒信号频谱计算方法与装置制造方法及图纸

本发明专利技术公开了一种用于计算多普勒信号频谱的方法与装置。所述方法包括估计步骤，用于对多普勒信号进行参数估计，所述参数用于表征多普勒信号的强度、频率中心、以及频率离散度；平滑步骤，用于对估计的参数进行平滑处理；以及计算步骤，用于利用所述参数求解预定的频谱函数，所述频谱函数至少是用于表征多普勒信号强度、频率中心、以及频率离散度的参数的函数。其中所述装置包括估计模块，平滑模块，以及计算模块。按照本发明专利技术技术方案的多普勒信号频谱计算方法与装置，获得的声谱图具有很高的时间和频率分辨率，而且无需进行复杂的运算。可以利用很低的成本实现较高质量的声谱图。

【技术实现步骤摘要】

本专利技术涉及一种多普勒信号频谱计算方法与装置，尤其是涉及 一种在超声诊断系统中用于计算多普勒血流或者组织的信号频谱的 方法与装置。
技术介绍
频谱多普勒技术被广泛用于人体内血流(或者组织运动)速度 的无损检测和测量。以超声多普勒血流信号检测为例，经过发射换 能器将电信号转变为超声信号，然后，入射到人体内。超声信号遇 到人体细胞(例如组织细胞或者血细胞)会发生散射，部分散射信 号到达接收换能器。接收换能器再将散射的声信号转变为电信号， 接收到的微弱电信号经过低噪声放大、正交解调、低通滤波，得到 两路音频正交多普勒信号。由于组织和血管壁回波具有比血细胞回 波大得多的散射幅度，因此往往要使获得的正交多普勒信号先通过 一个高通滤波器(又称壁滤波器)，以滤除极低频率的组织和管壁 回波信号。对获得的多普勒信号进行频谱分析，就可获得该信号的功率i普。对不同时刻的多普勒信号求功率谱就可获得多普勒信号的 声谱图，该声语图可送到显示器进行显示。同时，采集到的正交多 普勒信号经过方向分离，分别获得对应于正向血流和反向血流的音 频多普勒信号，该信号分别送到一对立体声扬声器的左、右声道输 出。另外，可以在声谱图上提取平均频率曲线、最大频率曲线，并 从中计算出一些重要的临床诊断参数。现有超声诊断系统通常对多普勒信号进行傅立叶变换来计算多 普勒信号的功率谱。在频谱计算之前，先对采集到的一段多普勒信号用窗函数进行加权处理，对加窗后的信号进行固定长度(如256点) 的FFT计算，然后对计算得到的系数取模的平方。当信号的长度小于 FFT计算的点数的时候，对加窗后的信号进行补零，使得补零后的信 号长度为FFT计算所需要的点凄t (如256点)。经过上述计算获得的声 谱图的时间和频率分辨率完全取决于窗函数的种类和长度。而窗函 数的种类通常是固定的，因此窗函数的长度也就决定了声谱图的分 辨率。为了保证显示的声谱图有较好的时间分辨率，窗函数的长度 通常与多普勒信号的采样率(即脉冲多普勒信号的脉冲重复频率或 者连续波多普勒信号的釆样率)和相邻谱线更新的时间间隔相关。 比如，信号采样率增加，则窗函数长度增加；相邻谱线更新的时间 间隔增加，则窗函数长度增加。但是当待检测目标的运动速度较低 (例如信号采样率下降到1kHz以下)，且要求的谱线更新时间较短(例 如语线更新间隔为4ms)的时候，为了满足较好的时间分辨率，上述 传统方法所需的窗函数很短，从而使频谱的频率分辨率下降很多， 并导致很大的速度估计误差。另外，也可以利用Autoregressive模型(即 AR模型)对超声多普勒信号进行建模，通过对AR模型的参数进行Z变 换(也可以用FFT实现快速Z变换)可获得信号的频谱。由于AR模型 是纯极点模型，在极点对应的频率上会产生一个非常大的谱峰，因 此直接利用AR模型的参数计算的功率谱具有非常大的动态范围，且 随信号带宽的变化功率谱曲线会产生非常大的改变，声谱图显示形 态显得不自然。文献和结合上述两种方法提出了一种利用AR模型 对采集到的多普勒信号进行建模，然后利用该AR模型对采集到的信 号进行预测估计(即实现信号外插值)，将采集到的信号与预测估 计的信号拼接后进行加窗和FFT计算处理。该方法一定程度上解决了 时间和频率分辨率的问题，声谱图的显示形态也与传统的基于FFT的 方法一样。但是，AR模型参数估计和信号的预测将使得实现成本大 大提高，且AR模型参数的估计偏差可能导致预测信号和拼接信号具 有不同的频率信息，并最终导致声谱图显示的畸变。
技术实现思路
本专利技术针对现有技术存在的这些缺点，提供一种多普勒信号频 谱计算方法与装置，其直接根据估计出的多普勒信号的平均频率、功率和带宽，利用解析函数才莫型估计多普勒信号的功率谱。为了实 现这一目的，本专利技术所采取的才支术方案如下。按照本专利技术实施例的第 一 方面，提供一种用于计算多普勒信号频谱的方法，包括以下步骤估计步骤，用于对多普勒信号进行参 数估计，所述参数用于表征多普勒信号的强度、频率中心、以及频 率离散度；以及计算步骤，用于利用所述参数求解预定的频谱函数， 所述频谱函数至少是用于表征多普勒信号强度、频率中心、以及频 率离散度的参数的函数。可选的是，按照本专利技术实施例第一方面的用于计算多普勒信号 频谱的方法，还包括平滑步骤，用于对估计的参数进行平滑处理。 优选的是，所述平滑处理为线性或非线性滤波。其中所述参数包括 多普勒信号的功率、平均频率、以及带宽。所述平滑处理包括对估 计的功率、平均频率、以及带宽的其中一个或多个进行平滑处理。优选的是，在所述估计步骤中，利用自相关法进行所述参数估 计，或者利用一阶AR模型进行所述参数估计。其中所述预定的频谱函数是频率的连续函数。所述频率的连续 函数可以是高斯函数，其中所述高斯函数的幅度取决于表征多普勒 信号强度的参数，所述高斯函数的中心取决于表征多普勒信号的频 率中心的参数，所述高斯函数的宽度取决于表征多普勒信号的频率 离散度的参数。按照本专利技术实施例的第二方面，提供一种用于计算多普勒信号 频语的装置，包括以下模块估计模块，用于对多普勒信号进行参 数估计，所述参数用于表征多普勒信号的强度、频率中心、以及频 率离散度；以及计算模块，用于利用所述参数求解预定的频谱函数,所述频谱函数至少是用于表征多普勒信号强度、频率中心、以及频 率离散度的参数的函数。可选的是，按照本专利技术实施例第二方面的用于计算多普勒信号频谱的装置，还包括平滑模块，用于对估计的参数进行平滑处理。 优选的是，所述平滑模块为线性滤波器或非线性滤波器。其中所述 参数包括多普勒信号的功率、平均频率、以及带宽。所述平滑模块 对估计的功率、平均频率、以及带宽的其中一个或多个进行平滑处 理。优选的是，所述估计模块利用自相关法进行所述参数估计，或 者利用 一 阶AR模型进行所述参数估计。按照本专利技术技术方案的多普勒信号频谱计算方法与装置，获得 的声谱图具有^f艮高的时间和频率分辨率，而且无需进行复杂的运算， 例如，大于二阶的AR模型参数估计或者256点的FFT计算等等。可以 利用很低的成本实现较高质量的声i普图。附图说明图1是典型的超声诊断成像系统中多普勒处理单元组成框图； 图2是按照本专利技术实施例的多普勒信号频谱计算方法流程图； 图3是按照本专利技术实施例的计算多普勒信号频谱的装置的结构 示意图4是采用现有技术获得的声i普图5是按照本专利技术的方法或装置获得的声谱图。具体实施例方式1.计算多普勒信号频谱的方法超声诊断成像系统中多普勒处理单元的典型组成如图1所示。 射频超声回波信号经过波束合成100后，进行正交解调102以获得I(in-phase)、 Q (quadrature)正交多普勒信号。当系统同时工作在两 维成像模式和频谱多普勒模式的时候，多普勒信号存在间隙，需经 过间隙填充104以获得连续的频谱显示和声音输出。填充后的多普勒 信号经过高通滤波器106去除大量大幅度的组织和管壁低频信号。高 通滤波后的信号经过频语分析108获得多普勒信号的声谱图，利用声 谱图提取最大频率、平均频率等曲线，并从中计算出一些重要的临 床诊断参数110。这些参数经过扫描转换112后，其中声谱图和最大 频率曲线在显示器114上实时显示。高通滤波后的信号还本文档来自技高网...

【技术保护点】
一种用于计算多普勒信号频谱的方法，其特征在于，包括：　估计步骤，用于对多普勒信号进行参数估计，所述参数用于表征多普勒信号的强度、频率中心、以及频率离散度；以及　计算步骤，用于利用所述参数求解预定的频谱函数，所述频谱函数至少是用于表征多普勒信号强度、频率中心、以及频率离散度的参数的函数。

【技术特征摘要】

【专利技术属性】
技术研发人员：张羽，
申请(专利权)人：深圳迈瑞生物医疗电子股份有限公司，
类型：发明
国别省市：94[中国|深圳]