CN101149433A - 一种使用反向Golay码译码的超声运动检测方法 - Google Patents

Abstract

本发明一种使用反向Golay码译码的超声运动检测方法，对接收电路的接受信号进行译码时采用一种由两块译码模块分别对两条数据路径进行处理的方法，使用原始正交的Golay对译码模块对第一条数据路径进行译码并通过叠加模块叠加；使用反向的Golay对译码模块对第二条数据路径进行译码并通过叠加模块进行叠加；最后通过抵消模块对第一条数据路径的结果和第二条数据路径结果进行处理，对第二条数据路径的叠加结果通过延迟模块做一合适的时间延迟。本发明中的第一条数据路径将使成像系统在改善SNR和增强穿透力的同时，去除伪像；第二条数据路径将探测和增强B型成像中的运动显示，增强的运动区域可被一个清晰的静态图像覆盖。

Description

一种使用反向Golay码译码的超声运动检测方法

技术领域

本发明涉及超声信号处理、超声编码激励技术，尤其是涉及一种使用反向Golay码译码的超声运动检测方法。

背景技术

众所周知，在超声成像中存在着分辨率和穿透性间折中的问题，增加系统的穿透性有两种方法：(1)、增加发射超声波的振幅。(2)、增加发射超声波的持续时间。然而声波的振幅不能无限的增加，在目前的超声系统中它已经很接近M工了，那唯一可行的办法就是增加超声波的持续时间，但是简单的增加脉冲的长度会降低系统的轴向分辨率，造成图像模糊。目前编码激励技术能很好的解决这个问题，编码激励技术在雷达中已经成功的应用了50多年了，在雷达系统中的编码激励技术能在不增加峰值发射功率的前提下，显著提高平均发射功率，从而提高系统的信噪比(SNR)。然后对于超声系统，超声波在人体内传播的时候有很强的频率衰减，产生非线性效应，使编码脉冲的波形扭曲，减低了编码激励的性能，所以编码激励技术只是最近5，6年才成功的应用于超声系统中。编码激励技术成像的框图见图1：

编码激励成像系统与传统的脉冲回波成像系统的不同之处在于：

(1)发射电路采用编码发射激励，必要时还需要对发射编码进行调制；

(2)接收电路中需要对回波信号做脉冲译码；

编码激励系统的基本工作原理如图2所示。图2(a)所示的传统脉冲发射系统使用单脉冲进行激励，图2(b)所示的编码激励系统，使用一长串编码脉冲进行激励。

Newhouse  在1974年提出了白噪编码的超声成像和多谱勒测量系统。在此后的30多年时间内，包括M序列伪随机码、Barker码、Golay码、Chirp和伪Chrip码等各种编码方法被用于超声编码激励技术的研究。

目前在超声编码激励系统中使用的比较多的编码是Chrip码、Golay码、Barker码。其中Chrip码是唯一在频率上编码的，而Golay码也是唯一需要多次发射的。下面分别简要的介绍下Chrip码和Golay码。

(1)Chrip码的编码激励技术

线性调频(Chrip)信号早雷达系统中的应用很广泛。由于Chrip信号有很好的自相关特性，很适合应用于超声编码发射。考虑到超声探头的带通特性，编码发射中使用的Chrip信号为：

$S\left(t\right)=\cos [2\mathrm{\π}(f_0-\frac{B}{2})t+\frac{\mathrm{\α}}{2}{t}^2]$

其中：f₀为探头的中心频率；B为探头的带宽；α为线性调频的斜率。但是Chrip码的编码激励的缺点很明显，它需要专门设计发射电路而且电路很复杂，不容易实现。目前超声系统中采用较多的是伪Chrip码，伪Chrip用一串0、1脉冲激励发射换能器，电路结构相对简单。伪Chrip信号定义如下：

$e\left(t\right)=\mathrm{sign}\left\{\cos \right[2\mathrm{\π}(f_0-\frac{B}{2})t+\frac{\mathrm{\α}}{2}{t}^2]\}$

其中，f₀为探头的中心频率；B为探头的带宽；α为线性调频的斜率，sign为符号函数。

(2)Golay码编码激励技术

在所有的编码激励技术中，Golay码编码激励技术在理论上是唯一一个可以完全消除旁瓣的激励技术，脉冲译码过程如图3所示：

Golay码又称Golay互补序列对，其定义为一对由两种元素构成的等长、有限序列。且在任何给定的间隔下，一个序列中的相同元素对的个数等于另一个序列中相异元素对的个数。在很多的实际的编码中，Golay码是唯一具有δ函数自相关这个优良特性的，自相关的和仍然为一个δ函数。数学上，一对Golay编码的脉冲可以被定义为：

[A₁，A₂]＝[G₁*S，G₂*S]

它由一对Golay序列[G₁，G₂]和传输脉冲S做卷积得到。由于Golay码的自相关的和是一个δ函数，即：

G₁_G₁+G₂_G₂＝δ

在接收信号的时候，使用一对正交的Golay对[G₁，G₂]对信号进行译码，得到：

G₁_(A₁*Q)+G₂_(A₂*Q)＝[(G₁_G₁+G₂_G₂)*S]*Q＝S*Q    (1)

，其中Q定义传播媒介。从上式我们可以看到除了增加SNR外，译码的Golay脉冲实际上和直接传送的基序脉冲的作用是一样的。

理论上，Golay互补序列对的编码发射可以在保持主瓣宽度不变的情况下，完全消除旁瓣。但Golay码编码激励技术存在以下两个问题：

(1)两次发射间组织或者血流的运动会严重影响其性能，使其往往达不到理想的效果。

(2)采用Golay码会使得图像的帧频降低一半。

在上面的论述中看到了Golay编码激励技术存在的两个问题。因为在现代的超声系统中都采用多焦区成像，这种技术本身就会使得帧频降低，所以使用Golay码编码激励造成的帧频降低可以隐含的解决。如果只使用单焦区成像，这个问题也是无论如何都解决不了的。所以主要研究两次发射间组织或者血流的运动对Golay编码激励技术的影响。

目前，GE公司取得了关于使用Golay编码激励技术在波束形成方面的一些专利，他们声称当Golay编码激励技术应用于对比谐波成像、组织谐波成像和B模式血流成像时，它能够提高超声系统的信噪比(SNR)并且获得更好的穿透性。然而，在GE的专利中，当他们使用成对的Golay码激励时，并没有解决组织或血流运动带来的问题。他们对这个问题采取的唯一的对策是：“对每一条波束，在两次Golay编码序列间使用两次连续的发射以降低组织运动的影响。据我们的分析可知，若仅使用交错发射的方法仍会产生严重的伪像。同时，在临床上的血流运动的有效范围内，伪信号可能比实际的信号更大。

发明内容

本发明的目的是提供一种去除运动伪影的使用反向Golay码译码的超声运动检测方法。

为了达到上述目的，本发明采用如下技术方案：一种使用反向Golay码译码的超声运动检测方法，在Golay编码激励的超声成像系统中，在对接收电路的接受信号进行译码时采用一种由两块译码模块分别对两条数据路径进行处理的方法，其方法为：

(1)、使用原始正交的Golay对译码模块对第一条数据路径进行译码，并通过叠加模块对译码结果进行叠加，获得叠加结果；

(2)、使用反向的Golay对译码模块对第二条数据路径进行译码，并通过叠加模块对译码结果进行叠加，获得叠加结果；

(3)、最后通过抵消模块对第一条数据路径的叠加结果和第二条数据路径的叠加结果进行处理获得没有运动伪影的组织信号，对第二条数据路径的叠加结果通过延迟模块做一合适的时间延迟获得运动伪影信号。

本发明的有益效果是：(1)、去除运动伪影；(2)、提出了运动检测这个新应用。第一条数据路径将使成像系统在改善SNR和增强穿透力的同时，去除伪像。第二条条数据路径帮助探测和增强B型成像中的运动显示，增强的运动区域可被一个清晰的静态图像覆盖。

附图说明

图1是编码激励技术成像框图；

图2是编码激励系统的基本工作原理图；

图3是Golay编码脉冲译码过程图；

图4是对本发明通过功能模块进行阐述的功能结构图；

图5是没有使用Golay编码激励时，静止靶线信号图；

图6是使用Golay编码激励时，静止靶线信号图；

图7是靶线移动2个采样点，接收到的靶线信号图；

图8是移动2个采样点时，得到的运动伪影信号图；

图9是移动2个采样点时，通过论文方法恢复得靶线；

图10中的表格为计算出在采样频率为36MHz的情况下，在不同深度的一些阀值速度。

具体实施方式

本发明一种使用反向Golay码译码的超声运动检测方法见图4、图5、图6、图7、图8、图9所示：一种使用反向Golay码译码的超声运动检测方法，在Golay编码激励的超声成像系统中(包括发射电路和接收电路，为已有电路)，Golay编码激励的超声成像系统为现有技术，在对接收电路的接受信号进行译码时采用一种由两块译码模块分别对两条数据路径进行处理的方法，其方法为：

(1)、使用原始正交的Golay对译码模块对第一条数据路径进行译码，并通过叠加模块对译码结果进行叠加(就是把译码结果加在一起)，获得叠加结果；见图4中，此步骤1由缓冲器1功能单元、缓冲器2功能单元、解码1功能单元完成。

(2)、使用反向的Golay对译码模块对第二条数据路径进行译码(就是把译码结果加在一起)，并通过叠加模块对译码结果进行叠加，获得叠加结果；见图4中，此步骤2由缓冲器1功能单元、缓冲器2功能单元、解码2功能单元完成。

(3)、最后通过抵消模块对第一条数据路径的叠加结果和第二条数据路径的叠加结果进行处理获得没有运动伪影的组织信号(该抵消模块的目的是把第一条数据路径的叠加结果减去第二条数据路径的叠加结果)，对第二条数据路径的叠加结果通过延迟模块做一合适的时间延迟获得运动伪影信号，其它静态靶线信号被压制了。此步骤3由解码1功能单元、解码2功能单元完成，把第一条数据路径的叠加结果减去第二条数据路径的叠加结果送静止组织增强处理，对第二条数据路径的叠加结果通过延迟模块做一合适的时间延迟获得运动伪影信号(即将解码1功能单元得到的结果送血流运动增强处理，再送往B模式组织血流重叠)。。

在上述方法中，原始正交的Golay对译码模块的作用就是将第一条数据路径进行译码，具体译码方法见下面详述；叠加模块的作用就是将第一条数据路径进行译码后将译码结果加在一起，具体叠加方法见下面详述；反向的Golay对译码模块的作用就是将第二条数据路径进行译码，具体译码方法见下面详述；叠加模块的作用就是将第二条数据路径进行译码后将译码结果加在一起，具体叠加方法见下面详述；抵消模块的作用就是将第一条数据路径的叠加结果减去第二条数据路径的叠加结果，具体抵消方法见下面详述；延迟模块的作用就是将第二条数据路径的结果作一合适的时间延迟，具体延迟方法见下面详述。使用反向Golay码译码的超声运动检测方法具体详述如下：

一对Golay码可以写成下面(2)式中这样的奇偶冲激对的形式，

G₁(t)＝δ(t)+δ(t-N_g)

G₂(t)＝δ(t)-δ(t-N_g)(2)

其中N_g+1是Golay码的长度，在系统实现时，可以把N_g大小设置为脉冲长度一样，把Golay码与双极脉冲发射信号做卷积，就得到了Golay编码序列，并把它存储在传输序列存储器中。Golay编码序列与换能器激励响应做卷积，就产生了下面(3)中的激励脉冲，

A₁(t)＝G₁*S＝S(t)+S(t-N_g)

A₂(t)＝G₂*S＝S(t)-S(t-N_g)(3)

式子(3)中的S(t)被设计为优化超声脉冲形状和频谱能量的形式，例如以下的一个余弦调制的高斯函数

$S\left(t\right)=e^{-{\left(\frac{t}{w}\right)}^2}\cos \left(2\mathrm{\π}{F}_{c}t\right)---\left(4\right)$

式子(3)中Golay编码激励脉冲A₁和A₂按顺序的发射到媒介Q(一序列靶线)，形式如

$Q\left(t\right)=\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}\mathrm{\δ}(t-k),$

得到信号如下面(5)式的形式：

$A_1*Q=\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-k)+\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-N_g-k)$

$A_2*Q=\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-k)-\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-N_g-k)---\left(5\right)$

在接收信号时，将有两条数据路径，第一条使用正交的Golay对译码，得到G₁_(A₁*Q)和G₂_(A₂*Q)的形式如下(6)和(7)所示，其中_表示互相关运算：

$G_1\⊗(A_1*Q)=[\mathrm{\δ}\left(t\right)+\mathrm{\δ}(t-N_g)]\⊗[\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-k)+\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-N_g-k)]$

$=\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-k)+\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-N_g-k)+\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t+N_g-k)+\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-k)---\left(6\right)$

$=\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t+N_g-k)+2\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-k)+\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-N_g-k)$

$G_2\⊗(A_2*Q)=-\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t+N_g-k)+2\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-k)-\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-N_g-k)---\left(7\right)$

很明显的，把等式(6)和(7)相加结果将是：

$G_1\⊗(A_1*Q)+G_2\⊗(A_2*Q)=4\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-k)---\left(8\right)$

(8)式意味着在理想的条件下(例如组织不动)，使用Golay编码激励理论上可以实现12dB信号强度的增加。

当在组织和血流运动的情况下，可以定义为第一条靶线和第二条靶线间的时间改变为τ这样的形式如：

$Q_1\left(t\right)=\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}\mathrm{\δ}(t-k)$

$Q_2\left(t\right)=\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}\mathrm{\δ}(t-k+\mathrm{\τ})$

第一条数据路径的译码过程将如下所示：

$G_1\⊗(A_1*Q)=\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t+N_g-k)+2\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-k)+\underset{k=k}{\overset{K}{\mathrm{\Σ}}}S(t-N_g-k)$

$G_2\⊗(A_2*Q)=-\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t+N_g-k+\mathrm{\τ})+2\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-k+\mathrm{\τ})-\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-N_g-k+\mathrm{\τ})---\left(9\right)$

在一个很小的时间改变的时候，宽带信号S(t)可以写为：

S(t+τ)＝-c(t)S(t)                      (10)

其中，系数c依赖于脉冲形状S和时间改变τ。在最坏的情况下，即由于产生的时间改变接近180度的相位差的时候，它等同于采样点移动的数目将是：

$n_{\mathrm{\τ}}=\frac{F_s}{2F_c}---\left(11\right)$

其中F_s和F_c分别表示采样频率和中心频率。

通过把上式应用于等式(9)，我们得到等式(12)，如下：

$G_1\⊗(A_1*Q)+G_2\⊗(A_2*Q)=\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}(1+c_1)S(t+N_g-k)+2\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}(1-c_2)S(t-k)+\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}(1+c_3)S(t-N_g-k)---\left(12\right)$

(12)式的第一项是由于做卷积和互相关运算的时候产生的镜像序列，但从T_s＝k₀开始使用接收信号，k₀是第一个靶线的深度。在实际的系统中，深度T_s是ROI窗口的深度。(12)式中的第二项是衰减的靶线信号。在180度相位改变这种最坏的情况下，c₂的值将接近1，以至于(12)式中第三项的伪影信号可能比实际的靶线信号还要大。

在第二条数据路径上，我们使用相反(正交)Golay码译码，

$G_2\⊗(A_1*Q_1)=[\mathrm{\δ}\left(t\right)-\mathrm{\δ}(t-N_g)]\⊗[\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-k)+\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-N_g-k)]$

$=\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-k)+\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-N_g-k)-\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t+N_g-k)-\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-k)---\left(13\right)$

$=-\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t+N_g-k)+\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-N_g-k)$

$G_1\⊗(A_2*Q_2)=\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t+N_g-k+\mathrm{\τ})-\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}S(t-N_g-k+\mathrm{\τ})---\left(14\right)$

可以清楚的看到，如果不考虑运动的情况，两个式子的和将是0。然而通过运用式子(10)，得到第二条数据路径的结果，

$G_2\⊗(A_1*Q_1)+G_1\⊗(A_2*Q_2)=-\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}(1+c_1)S(t+N_g-k)+\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}(1+c_3)S(t-N_g-k)---\left(15\right)$

同样(15)式的第一项是由互相关运算产生的镜像序列，可以被忽略掉。

式子(15)表示使用反向Golay码译码可以去除时间改变为k₀+N_g的靶线信号，这将产生运动检测这个新的应用。

如果定义深度Z和声速V_c，则回程声波时延将为：

$d_t=\frac{2z}{V_c}$

并且这个时间对应着采样点移动的数目n为：

$n=\frac{V_f{d}_t}{V_c}{F}_s=\frac{2V_{f}Z{F}_s}{V_c^2}$

在这里V_f表示运动靶线的速度(例如血流速度)，F_s表示采样频率。这个式子等价于

$V_f=\frac{n{V}_c^2}{2Z{F}_s}---\left(16\right)$

如果设置式子(16)中的n＝1，就获得的了阈值速度，它意味着任何运动靶线速度大于阈值数值的任何的时候，都能从第一个数据路径中检测出伪影信号(使用原始Golay译码)，在第二个数据路径中检测出流动信号(使用反向Golay译码)。

表格(见图10)计算出在采样频率为36MHz的情况下，在不同深度的一些阈值速度。

从式子(16)可以清楚的看到，可以不只一次的交错发射第一个和第二个Golay编码脉冲来扩展检测运动速度的范围。例如发射第一个脉冲两次在一个连续的位置，同样的在同一位置发射两次第二个脉冲。这个将使可以检测的速度范围扩大两倍，例如可以检测出原速度一半的这样的慢速。在数学上我们有：

$V_f=\frac{n{V}_c^2}{2\mathrm{mZ}{F}_s}---\left(17\right)$

，m为交替发射脉冲的数目。

这样就得到了两条数据路径的数据：在第一条数据路径上使用正交的Golay对译码，对于静止的靶线我们得到了式子(8)，对于运动的靶线我们得到了式子(12)；在第二条数据路径上使用反向的Golay码译码，对于静止的靶线得到了0，对于运动的靶线，我们得到了式子(15)。

因此，通过把第一条数据路径的结果减去第二条数据路径的结果，得到了更高SNR并且没有运动伪影的静态信号：

$[G_1\⊗(A_1*Q_1)+G_2\⊗(A_2*Q_2)]-[G_2\⊗(A_1*Q_1)+G_1\⊗(A_2*Q_2)]$

$=2\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}(1+c_1)S(t+N_g-k)+2\underset{k=k_0}{\overset{K}{\mathrm{\Σ}}}(1-c_2)S(t-k)---\left(18\right)$

同样，(18)式的第一项可以被忽略掉，从T_s＝k₀开始观察信号。等式(18)给出了一个受损的靶线信号，但去除了由运动产生的伪影信号。注意到受损的靶线信号依赖于脉冲的形状和运动产生的时间改变。这将导致靶线信号轴向分辨率的降低。从模拟的结果得到，(18)式在运动速度从21cm/s到125cm/s(即在采样频率为36MHz时，采样点由1到6时的平均信号峰值比传统的方法高6.7dB。

因此，将两条译码路径的结果用于新的应用：(1)式子(15)表示的第二条路径仅仅包含运动信号，静态信号被压制了。(2)把第一条数据路径的结果减去第二条数据路径的结果，得到没有运动伪影的静态信号。

以上方法(1)可做一个运动检测器，定性的显示沿着波束方向的B型血流模型，它比传统的彩色血流显示具有高帧频和更好的轴向分辨率。方法(2)将用于具有高SNR和好的穿透性，并且没有运动伪影的B模式成像。

这个方法可以解决Goaly编码激励技术存在的第一个问题。它不但可以去除运动产生的伪影信号，还可以提高超声系统的SNR和穿透性。具体的效果和优点可以在计算机仿真和结果中看到。测试信号S来自超声系统，中心频率F_s＝3.5MHz，采样频F_c＝40MHz。靶线信号T的采样频率为40MHz，它是在1000个采样点和2000个采样点处的两个冲激信号。这样把T和S做卷积得到了RX-Target1信号，峰值为3.6，它将作为参考信号，如图5所示。

在前部分的论述中知道Golay编码激励对于静态的靶线，除了增加SNR外没有什么区别，通过计算机仿真实验得到图6。

可以计算得到，静态靶线信号的强度增加了12dB。然而，由于靶线的运动，接收到的实际信号不会像图3那么理想，将会产生运动伪影。前面等式(9)给出了在运动存在时，接收的靶线信号的形式，其中式子中的第三项即为运动产生的伪影信号，把靶线移动2个采样点，得到图7。

从实验中还看到当靶线移动6个采样点时，图7中最后两个信号比实际的靶线信号还要大。通过计算知道6个采样点为最差的情况，即相当于180度的相位旋转。

在第二条数据路径上，对接收到的编码脉冲做的是反向译码，得到的仅仅是由于运动而产生的那部分信号。当靶线移动2个采样点时得到了图8。

等式(13)是通过使用Golay反向译码去除运动伪影信号的结果。从图7和图8中也可以很清楚的看到，图7的最后两项和图8是相同的，通过计算机仿真得到图9。

从图9中很清楚的看到，完全消除了运动产生的伪影信号，并且提高了SNR。过大量的实验数据证明，运用这种方法可以平均增加7.3dB的信号强度。这样就克服了Golay编码激励对于运动敏感这个缺点，使其更适合超声影像系统。

Claims (1)

1.一种使用反向Golay码译码的超声运动检测方法，其特征在于：在Golay编码激励的超声成像系统中，在对接收电路的接受信号进行译码时采用一种由两块译码模块分别对两条数据路径进行处理的方法，其方法为：

(1)、使用原始正交的Golay对译码模块对第一条数据路径进行译码，并通过叠加模块对译码结果进行叠加，获得叠加结果；

(2)、使用反向的Golay对译码模块对第二条数据路径进行译码，并通过叠加模块对译码结果进行叠加，获得叠加结果；

(3)、最后通过抵消模块对第一条数据路径的叠加结果和第二条数据路径的叠加结果进行处理获得没有运动伪影的组织信号，对第二条数据路径的叠加结果通过延迟模块做一合适的时间延迟获得运动伪影信号。

Images