[发明专利]定量超声系统中基于经验模态分解的组织微观结构检测方法

说明书

技术领域

本发明涉及一种利用超声背散射技术对表征组织微观结构的平均散射子间距进行检测的方法。

背景技术

医学超声广泛应用于人体组织器官成像中，虽然该技术在临床应用中已经具有相当成熟的应用基础，并能提供实时的高精度人体器官图像，但其不能提供由于组织钙化、变异等导致的微观结构变化。

目前一些等效散射子模型如二元混合模型、Faran圆柱模型以及弱散射子模型等常用来表征组织微观结构变化。在这些模型中，组织的检测转化为对超声衰减、背散射系数、散射子密度、散射子直径和平均散射子间距(MSS)等的检测。其中MSS是一种最常用的方法，该方法主要利用了生物组织在声学特性分布上具有一定的准周期性，这种生物周期性来自组织结构中的规则散射元和弥散散射元，会随着组织健康状况改变而发生变化，通过利用MSS测量该周期性，可以对生物组织微观结构进行测量，并对如血管动脉粥样硬化、骨质疏松等健康状况做出诊断。

在平均散射子间距估计中目前常用的方法主要有倒谱方法、功率谱方法和瞬态相关图等方法。这些方法在一定程度上能够有效检测平均散射子间距，但受背散射超声信号序列短、衰减大和噪声大的影响，检测结果不稳定。一些自相关谱分析方法、周期性小波变换方法、奇异谱分析方法等提高了检测鲁棒性，但计算量大，在去除噪声和弥散散射信号干扰方面有待提高。

发明内容

为了解决定量超声系统中利用现有技术方法无法有效消除超声背散射信号中的噪声和弥散散射信号干扰的问题，本发明提供了一种定量超声系统中基于经验模态分解(EMD)的组织微观结构检测方法。

本发明的目的是通过以下技术方案实现的：

一种定量超声系统中基于经验模态分解的组织微观结构检测方法，包括如下步骤：

步骤一、获取生物组织超声背散射信号；

步骤二、截取DOI(Domain ofInterest)区域的背散射信号作为输入数据；

步骤三、对输入数据进行EMD分解，得到三阶本征模态函数和一个残差项；

步骤四、执行能量谱筛选：提取三阶本征模态函数能量谱，并与原信号能量谱进行相似度分析，确定包含组织微结构信息的最优模态函数项(BIMF)；

步骤五、对确定出的BIMF进行估计，计算组织散射子间距信息。

所述步骤一中，超声背散射信号获取有仿真和实验两种获取方式。其中实验获取方式中，超声背散射信号收发由2MHz高频单阵元聚焦探头实现，探头聚焦距离为25mm，探测方式为单一探头背散射作用方式。

所述步骤二中，截取DOI区域信号时，采用矩形窗函数。仿真数据截取信号序列时间长度为4μs，实验数据截取信号序列长度为8.125μs。

所述步骤三中，基于经验模态分解的组织微观结构检测方法具体通过以下迭代进行检测：

1)输入DOI区域背散射信号x(t)(t＝1，2，...，N，N为采样点数)为原始函数；

2)确定输入的原始函数极值点，采用三次样条函数确定输入信号上下包络x_L(t)x_U(t)；

3)计算上下包络均值定义h_i(t)为原始函数与该包络均值的差值函数h_i(t)＝x_i(t)-m_i(t)；(i表示第i次分解，分解结束得到第i个本征模态函数)

4)迭代计算判断差值函数h₁(t)是否满足IMF条件，判断决策过程采用标准偏差(Standard Deviation，SD)准则：