一种应用于超声成像的稳健宽带自适应波束形成方法
【技术领域】
[0001] 本发明涉及医学超声诊断及超声无损检测等领域，尤其涉及到一种应用于超声成 像的稳健宽带自适应波束形成方法。
【背景技术】
[0002] 传统的医学超声诊断及超声无损检测系统中普遍使用了延时叠加的波束形成方 法，这种方法简单可靠，易于实现，但性能却比较一般。为提高超声成像的质量，最小方差自 适应波束形成方法受到了广泛的关注，然而这种方法面临如下的两个主要问题：（1)经典 的最小方差波束形成方法是根据窄带信号的模型所提出的，而超声成像领域的回波信号一 般属于宽带信号的范畴，所以直接将最小方差波束形成方法应用于超声成像中并未能充分 展现出原方法的优势；(2)对于非均匀介质的成像过程中往往存在着相位畸变，导致最小 方差波束形成方法中的协方差矩阵和导向矢量的估计存在一定误差，继而影响了该方法的 稳健性，难以在较为复杂的应用环境中实现理想的性能。
[0003] 目前存在的稳健自适应波束形成方法，大都围绕在协方差矩阵的重构及导向矢量 的估计上，如在最小方差波束形成方法中，常使用对角加载技术来提高协方差矩阵估计的 稳健性，这种方法简单易行，而对角加载因子这一参数的选取则对性能有着重要的影响。采 用前后向相结合的空间平滑方式对超声回波信号去相干，无需使用对角加载即可在一定程 度上保证协方差矩阵估计的稳健性，同时也能够提高最终成像的对比度。但上述方法都集 中在时域内处理，并未充分考虑到超声回波信号的宽带性特点，也未能对相位畸变条件下 失配的导向矢量进行有效的估计，导致最终的性能提升较为有限。
【发明内容】
[0004] 本发明针对现有技术的不足，对最小方差自适应波束形成方法进行改进，充分利 用超声回波信号的宽带性特点，在频域内对真实环境下的导向矢量进行重构，使其能在相 位畸变的条件下具备较为稳健的性能，并同时提高非均匀介质下的超声成像质量。
[0005] 为实现上述目的，本发明提供了一种应用于超声成像的稳健宽带自适应波束形成 方法，应用于等间距的M个阵元的超声换能器阵列中，该方法主要使用在对接收的回波数 据进行波束形成处理的阶段，包括以下的步骤：
[0006] 将所得到的超声回波数据由时域转换至频域，然后进行频带的划分；
[0007] 计算所述频带划分后的各子频带的协方差矩阵；
[0008] 根据得到的所述协方差矩阵，利用特征空间分解方法将假定的导向矢量投影在信 号子空间中；
[0009] 通过凸优化方法估计真实的导向矢量；
[0010] 利用最小方差方法的基本原理，得到各子频带空间的波束形成处理结果；
[0011] 对所述各子频带空间的波束形成处理结果进行叠加，并转换到时域，得到用于最 终成像的输出结果。
[0012] 进一步地，所述将所得到的超声回波数据由时域转换至频域，然后进行频带的划 分步骤，包括：
[0013] 对B模式超声成像下每一根扫描线相应的各个阵元所接收到回波数据进行短时 傅立叶变换，将变换后所得到的频域数据划分为K个子频带。
[0014] 更进一步地，所述将所得到的超声回波数据由时域转换至频域，然后进行频带的 划分步骤包括：
[0015] 根据发射的脉冲宽度预先设计出所述短时傅立叶变换的窗函数的宽度，然后进行 回波数据空间的划分；
[0016] 对所述回波数据空间划分后的各段数据进行"补零"处理，然后由时域转换至频 域；
[0017] 对所述频域数据进行频带的划分。
[0018] 进一步地，所述计算所述频带划分后的各子频带的协方差矩阵步骤，包括：
[0019] 对各子频带的数据进行空间平滑处理，其中以L个阵元为一组来计算协方差矩 阵，共分为M-L+1组，计算的公式如下：
[0021] 其中，:&(?,)是协方差矩阵，是频域上的回波数据，的共辄转 置，i=l，2，…，K，K为子频带的个数。
[0022] 更进一步地，所述利用特征空间分解方法将假定的导向矢量投影在信号子空间中 的步骤，包括：
[0023] 按下式对得到的协方差矩阵在特征空间进行分解：
[0024] R=UAU^ -UsAsU/+UpApUf
[0025] 其中A=diag[入!，入2, ? ? ?，入_…，入J，且特征值（入！多入2多…入num多…彡入L)， U= [1?,?,... ,Unum,...uj是与之相对应的特征矢量，选取参数P介于0. 1~0. 5,满足 入num彡员入1;
[0026] 将整个特征空间划分为信号子空间Us及噪声子空间Up两部分。
[0027] 利用得到的信号子空间Us，计算假定的导向矢量&在信号子空间内的投影：
[0028] a二
[0029] 其中，Us为信号子空间，USH为Us的共辄转置，假定的导向矢量&中的元素均取为 "1"。
[0030] 更进一步地，所述通过凸优化方法估计真实的导向矢量，包括：
[0031] 利用凸优化方法估计真实情况下的导向矢量a，计算公式如下：：
[0034] 其中，免1是协方差矩阵的逆，导向矢量S是假定导向矢量在信号子空间内的投影， 参数e是所设定的误差范围，一般介于
[0035] 更进一步地，所述利用最小方差方法的基本原理，得到各子频带空间的波束形成 处理结果，包括：
[0036] 根据得到的协方差矩阵及真实情况下的导向矢量，利用最小方差波束形成的基本 原理计算出一组最优权向量；
[0037] 具体地，按如下公式进行计算：
[0039] 其中wMV(to是最优向量，R .V(辫）是得到的协方差矩阵的逆，导向矢量a是估计的 真实情况下的结果；
[0040] 根据得到的最优权向量以及如下公式计算出各子频带数据的波束形成的输出结 果：
[0042] 其中，<V⑷）是Wmv(CO1)的共辄转置，Y(CO1)是第i个子频带的波束形成的输出 结果，Xd(W1)是频域上的回波数据。
[0043] 进一步地，所述对所述各子频带数据的波束形成处理结果进行叠加，得到频域空 间的波束形成输出结果步骤，包括：
[0044] 对每一段数据中的各子频带空间的波束形成的处理结果进行累加，得到频域空间 的波束形成输出结果，计算的公式如下：
[0046] 其中，Y(co)是频域空间的波束形成输出结果，Y(CO1)是第i个子频带的波束形成 的输出结果。
[0047] 本发明的优点在于：
[0048] 本发明通过短时傅立叶变换将时域的回波数据转换到各个子频带中处理，充分利 用了超声回波信号的宽带特点，能够取得更为理想的性能；同时在各个频带内将假定的导 向矢量投影在信号子空间内，并采用凸优化方法进行重构，提高了算法在相位畸变条件下 的稳健性。
【附图说明】
[0049] 图1为本发明实施例提供的流程示意图；
[0050] 图2为本发明实施例提供的对原始时域上的回波数据进行时延预处理的示意图；
[0051] 图3为本发明实施例提供的得到最终输出结果的原理框图；
[0052] 图4为本发明实施例提供的各子频带空间进行波束形成的处理过程。
【具体实施方式】
[0053] 下面结合附图和实施例，对本发明的技术方案做进一步的详细描述。
[0054] 图1为本发明实施例提供的一种自适应波束形成方法流程示意图，本发明实施例 应用于B模式下的超声成像。等间距的M个阵元的超声换能器阵列工作在B模式成像线性 扫描方式下，对所接收到的回波数据进行波束形成处理，每一根扫描线对应着M个阵元的 接收信号（即M为阵元数目），并且每个阵元的采样数目为N，在B模式超声成像下每一根 扫描线对应的所有阵元接收的回波数据可以用M*N的矩阵来描述。需要说明的是，这些时 域上的回波数据已经经过时延的预处理，图2给出了时延预处理的示意图，以下步骤中即 以预处理得到的X。作为默认的处理对象，同时需要指出的是由于非均匀介质的影响，回波 数据并不能保证得到完美的"对齐"，而本发明也无需对回波数据的时延预处理结果做出额 外的假设。图3是本发明实施例提供的得到最终输出结果的原理框图，由图1和图3所示， 可知本发明实施例的的波束形成处理方法包括步骤101-103 :
[0055] 步骤101，将所得到的超声回波数据由时域转换至频域，然后进行频带的划分；
[0056] 具体的，在B模式超声成像下每一根扫描线中的所有阵元接收的回波数据可以用 M*N的矩阵来描述，其中M为阵元数目，N为采样数目，需要说明的是，该回波数据已经完成 了时延预处理，对该回波数据进行短时傅立叶变换，并将变换后所得到的频域数据划分为K 个子频带，为在各个子空间内的波束形成完成必要的准备。
[0057] 下面进一步详细介绍所述步骤101的处理方法：
[0058] 步骤1011，根据发射的脉冲宽度预先设计出所述短时傅立叶变换的窗函数的宽 度，然后进行回波数据空间的划分；
[0059] 具体地，根据发射的脉冲宽度设计出合适的窗函数的宽度N1，将整个数据空间划 分为Ng段；其中，N为采样数目，Ng=N/N1并取其整数部分，即每一段对应的数据空间大 小为MXN1。需要进一步说明的是，为保证最终成像的纵向分辨率，Nl-般不会超过发射脉 冲的宽度。
[0060] 步骤1012,对所述回波数据空间划分后的各段数据进行"补零"处理，然后由时域 转换至频域；
[0061] 具体地，在每个阵元所对应的时域上的回波数据后添加（Nf-Nl)个"0"，然后对每 一段数据进行傅里叶变换处理，即将时域上的回波数据转换到频域，经过短时傅里叶变换 后的回波数据将有Ng段；
[0062] 在处理过程中，为了短时傅立叶变换的便捷计算，Nf-般取2的指数幂，并满足 Nf>2Nl，从而每一段数据空间的大小变为M*Nf，即变换后得到的每一段频域数据都是M*Nf 的矩阵；
[0063] 步骤1013,对所述频域数据进行频带的划分；
[0064] 具体地，上述步骤1012共得到Ng段M*Nf的频域数据，对得到的频域数据，进行频 带划分，将每一段频带划分为K个子频带，其中，Nl〈K〈Nf/2。
[0065] 步骤102,对所述频带划分后的各子频带空间进行波束形成处理，得到所述各子频 带空间的波束形成处理结果；
[0066] 具体地，如图4所示，在所划分的各个子频带中，在获得各子频带的协方差矩阵 后，先通过特征空间分解的方法将假定的导向矢量投影在信号子空间内，然后通过凸优化 方法估计实际情形下的导向矢量，继而利用最小方差方法的基本原理，得到对各个子频带 空间波束形成输出结果。
[0067] 下面以第i(i= 1，2，…，K)个子频带为例，详细介绍步骤102的处理过程：
[0068] 步骤1021