一种带运动补偿的合成孔径超声成像方法转让专利

申请号 : CN201410419588.4

文献号 : CN105403891B

文献日 :

基本信息:

PDF:

法律信息:

相似专利:

发明人 : 理华孟晓辉肖灵

申请人 : 中国科学院声学研究所

摘要 :

本发明涉及一种带运动补偿的合成孔径超声成像方法,应用于带有N个阵元的合成孔径超声换能器阵列中,各个阵元之间等间距,且各个阵元均能发射信号与接收信号;包括:采用发射、接收都是单焦点,且发射焦点与接收焦点相同的线性扫描方式得到第一低分辨率扫描线;采用偏移阵元接收的方式得到运动估计扫描线;利用运动估计扫描线进行估计,所得到的估计值用于对第一低分辨率扫描线做补偿,得到第二低分辨率扫描线;将虚拟源开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与虚拟源对应第二低分辨率扫描线上的相应采样点相加,得到高分辨率扫描线,进而得到高分辨率图像;其中,所述虚拟源为合成孔径超声换能器阵列中各个阵元的发射焦点。

权利要求 :

1.一种带运动补偿的合成孔径超声成像方法,应用于带有N个阵元的合成孔径超声换能器阵列中,各个阵元之间等间距,且各个阵元均能发射信号与接收信号;该方法包括:步骤1)、采用发射、接收都是单焦点,且发射焦点与接收焦点相同的线性扫描方式得到第一低分辨率扫描线;采用偏移阵元接收的方式得到运动估计扫描线;

步骤2)、利用步骤1)得到的运动估计扫描线进行估计,所得到的估计值用于对步骤1)所得到的第一低分辨率扫描线做补偿,从而得到第二低分辨率扫描线;

步骤3)、将虚拟源开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与虚拟源对应第二低分辨率扫描线上的相应采样点相加,得到高分辨率扫描线,进而得到高分辨率图像;其中,所述虚拟源为合成孔径超声换能器阵列中各个阵元的发射焦点;

所述步骤1)包括:

采用合成孔径超声换能器阵列中的第1到第M个阵元发射信号,然后采用第i到第M+i-1个阵元发射信号,以此类推,直到第N-M+1到第N个阵元发射信号,共N-M+1次发射信号,其中,将第i次发射的阶段记为Ti,1≤i≤N-M+1;在接收信号时,则会依据阵列中发射信号的阵元的位置来选取合适的接收阵元;具体包括:当i=1,即第T1次发射时,用所述阵列中的第1到第M个阵元发射信号,发射焦点为VS1;

由所述阵列中的第1到第M+1个阵元接收信号,每个阵元的采样点数为P,阵列中的阵元在接收信号时需经历正常接收与右移接收两个阶段;

正常接收:对第1到第M阵元所接收的信号进行波束形成,接收焦点也为VS1,得到第1条第一低分辨率扫描线L1,该扫描线有P个采样点;

右移接收:对第2到第M+1阵元所接收的信号进行波束形成,接收焦点也为VR1,VR1=VS1+Δ,得到第1条右移焦点扫描线LR1,该扫描线有P个采样点;

当1

正常接收:对第i到第M+i-1阵元所接收的信号进行波束形成,接收焦点也为VSi,得到第i条第一低分辨率扫描线Li,该扫描线有P个采样点;

右移接收:对第i+1到第M+i阵元所接收的信号进行波束形成,接收焦点也为VRi,VRi=VSi+Δ,得到第i条右移焦点扫描线LRi,该扫描线有P个采样点;

当i=N-M+1时,即第TN-M+1次发射时,用所述阵列中的第N-M+1到第N个阵元发射信号,发射焦点为VSN-M+1,由所述阵列中的第N-M到第N个阵元接收信号,每个阵元的采样点数为P,阵列中的阵元在接收信号时需经历左移接收与正常接收两个阶段;

左移接收:对第N-M到第N-1阵元所接收的信号进行波束形成,接收焦点也为VLN-M+1,VLN-M+1=VSN-M+1-Δ,得到第N-M条左移焦点扫描线LL(N-M+1),该扫描线有P个采样点;

正常接收:对第N-M+1到第N阵元所接收的信号进行波束形成,接收焦点也为VSN-M+1,得到第N-M+1条第一低分辨率扫描线LN-M+1,该扫描线有P个采样点;

依次重复上述步骤,直至得到N-M+1条第一低分辨率扫描线;其中,第一低分辨率扫描线L1,L2,……,LN-M+1构成一幅低分辨率图像LRI;所得到的左移焦点聚焦扫描线LL2,LL3,……,LL(N-M+1)用于与右移焦点聚焦扫描线LR1,LR2,……,LR(N-M)做运动估计处理。

2.根据权利要求1所述的带运动补偿的合成孔径超声成像方法,其特征在于,所述步骤

2)进一步包括:

步骤2-1)、计算第一条左移焦点聚焦扫描线LL2与第一条右移焦点聚焦扫描线LR1之间的时延,得到第一时延值D1;

步骤2-2)、计算第i条左移焦点聚焦扫描线LL(i+1)与第i条右移焦点聚焦扫描线LR(i)之间的时延,得到第i个时延值Di;

步骤2-3)、第一条第一低分辨率扫描线L1不做补偿,直接作为第一条第二低分辨率扫描线,记作CL1;

步骤2-4)、对于第二条第一低分辨率扫描线L2,利用第一时延值D1进行延时处理,得到补偿后的第二条第二低分辨率扫描线CL2;

步骤2-5)、利用第1至第i时延值之和D1+D2+…+Di对第i+1条第一低分辨率扫描线Li+1进行延时处理,得到补偿后的第i+1条第二低分辨率扫描线CLn+1,直至i的值达到N-M,得到补偿后的第N-M+1条第二低分辨率扫描线CLN-M+1;

其中,CL1、CL2……CLN-M+1称为第二低分辨率扫描线。

3.根据权利要求2所述的带运动补偿的合成孔径超声成像方法,其特征在于,所述步骤

3)进一步包括:

步骤3-1)、对第一条第二低分辨率扫描线进行第二次波束形成,即将VS1开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与第一条第二低分辨率扫描线上相应采样点相加,得到第一条高分辨率扫描线H1;

步骤3-2)、对第二条第二低分辨率扫描线进行第二次波束形成,即将VS2开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与第二条第二低分辨率扫描线上相应采样点相加,得到第二条高分辨率扫描线H2;

步骤3-3)、重复上述处理,直到对第N-M+1条第二低分辨率扫描线进行第二次波束形成,即将VSN-M+1开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与第N-M+1条第二低分辨率扫描线上相应采样点相加,得到第N-M+1条高分辨率扫描线HN-M+1。

4.根据权利要求3所述的带运动补偿的合成孔径超声成像方法,其特征在于,将虚拟源开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与虚拟源对应第二低分辨率扫描线上相应采样点相加,得到高分辨率扫描线的具体计算公式为:其中, 表示第二低分辨率扫描线 上延时td(k)对应的值,w为加权函

数,其值随深度z和扫描线的水平位置xk而改变;K是深度z的函数,K(z)表示虚拟源开角内的低分辨率扫描线数量;

K(z)的计算公式如下:

其中,zv为虚拟源深度,α为虚拟源的开角,L(z)为深度z处的声场宽度,Δ为虚拟源间距;

延时td(k)的计算公式为:

其中,zv为虚拟源深度, 为(x,z)点到线 上的虚拟源 的距离;+号表示(x,z)点深度大于zv,-号表示其深度小于zv,c为声速。

说明书 :

一种带运动补偿的合成孔径超声成像方法

技术领域

[0001] 本发明涉及合成孔径成像领域,特别涉及一种带运动补偿的合成孔径超声成像方法。

背景技术

[0002] 合成发射孔径通过单阵元发射,全孔径接收,经波束形成得到低分辨率图像(LRI),再通过将所有发射阵元得到的LRI相干相加得到高质量图像(HRI)。由于HRI对发射和接收都实现了动态聚焦,所以其较传统的超声成像方法就有更高的分辨率。但是,这种合成发射孔径成像方法面临三个主要问题:1)单阵元发射造成的穿透性差和回波信号信噪比低;2)低分辨率图像间的运动会降低其相干性,从而影响成像质量;3)系统实现过程中需要存储和传输大量的射频回波(RF)线,系统复杂度很高。
[0003] 对于合成发射孔径成像方法的上述三个问题,目前有多种改进方法。合成孔径序列波束方法就是一种比较容易实现的方法。合成孔径序列波束方法由于不需要存储和传输大量的单阵元接收的RF线数据,其实现复杂度较合成发射孔径方法大大降低。此外,由于是多阵元发射,所以也不存在合成发射孔径方法单阵元发射造成的穿透性差和回波信号信噪比低的问题。但是,合成孔径序列波束方法作为一种合成孔径方法,其中仍然存在运动影响成像质量的问题。
[0004] 传统的运动补偿算法需要增加多次发射次数,然后利用增加的接收回波进行运动估计,再进行补偿。这样的方法的缺点是降低了系统成像的扫描帧频,会造成图像不连贯,高速运动物体成像不佳等问题。

发明内容

[0005] 本发明的目的在于解决现有技术在存在轴向位移时成像质量较差的问题,从而提供一种无需增加发射次数却能得到高分辨率图像的方法.
[0006] 为了实现上述目的,本发明提供了一种带运动补偿的合成孔径超声成像方法,应用于带有N个阵元的合成孔径超声换能器阵列中,各个阵元之间等间距,且各个阵元均能发射信号与接收信号;该方法包括:
[0007] 步骤1)、采用发射、接收都是单焦点,且发射焦点与接收焦点相同的线性扫描 方式得到第一低分辨率扫描线;采用偏移阵元接收的方式得到运动估计扫描线;
[0008] 步骤2)、利用步骤1)得到的运动估计扫描线进行估计,所得到的估计值用于对步骤1)所得到的第一低分辨率扫描线做补偿,从而得到第二低分辨率扫描线;
[0009] 步骤3)、将虚拟源开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与虚拟源对应第二低分辨率扫描线上的相应采样点相加,得到高分辨率扫描线,进而得到高分辨率图像;其中,所述虚拟源为合成孔径超声换能器阵列中各个阵元的发射焦点。
[0010] 上述技术方案中,所述步骤1)包括:
[0011] 采用合成孔径超声换能器阵列中的第1到第M个阵元发射信号,然后采用第i到第M+i-1个阵元发射信号,以此类推,直到第N-M+1到第N个阵元发射信号,共N-M+1次发射信号,其中,将第i次发射的阶段记为Ti,1≤i≤N-M+1;在接收信号时,则会依据阵列中发射信号的阵元的位置来选取合适的接收阵元;具体包括:
[0012] 当i=1,即第T1次发射时,用所述阵列中的第1到第M个阵元发射信号,发射焦点为VS1;由所述阵列中的第1到第M+1个阵元接收信号,每个阵元的采样点数为P,阵列中的阵元在接收信号时需经历正常接收与右移接收两个阶段;
[0013] 正常接收:对第1到第M阵元所接收的信号进行波束形成,接收焦点也为VS1,得到第1条第一低分辨率扫描线L1,该扫描线有P个采样点;
[0014] 右移接收:对第2到第M+1阵元所接收的信号进行波束形成,接收焦点也为VR1,VR1=VS1+Δ,得到第1条右移焦点扫描线LR1,该扫描线有P个采样点;
[0015] 当1
[0016] 左移接收:对第i-1到第M+i-2阵元所接收的信号进行波束形成,接收焦点也为VLi,VLi=VSi-Δ,得到第i-1条左移焦点扫描线LLi,该扫描线有P个采样点;
[0017] 正常接收:对第i到第M+i-1阵元所接收的信号进行波束形成,接收焦点也为VSi,得到第i条第一低分辨率扫描线Li,该扫描线有P个采样点;
[0018] 右移接收:对第i+1到第M+i阵元所接收的信号进行波束形成,接收焦点也为 VRi,VRi=VSi+Δ,得到第i条右移焦点扫描线LRi,该扫描线有P个采样点;
[0019] 当i=N-M+1时,即第TN-M+1次发射时,用所述阵列中的第N-M+1到第N个阵元发射信号,发射焦点为VSN-M+1,由所述阵列中的第N-M到第N个阵元接收信号,每个阵元的采样点数为P,阵列中的阵元在接收信号时需经历左移接收与正常接收两个阶段;
[0020] 左移接收:对第N-M到第N-1阵元所接收的信号进行波束形成,接收焦点也为VLN-M+1,VLN-M+1=VSN-M+1-Δ,得到第N-M条左移焦点扫描线LL(N-M+1),该扫描线有P个采样点;
[0021] 正常接收:对第N-M+1到第N阵元所接收的信号进行波束形成,接收焦点也为VSN-M+1,得到第N-M+1条第一低分辨率扫描线LN-M+1,该扫描线有P个采样点;
[0022] 依次重复上述步骤,直至得到N-M+1条第一低分辨率扫描线;其中,第一低分辨率扫描线L1,L2,……,LN-M+1构成一幅低分辨率图像LRI;所得到的左移焦点聚焦扫描线LL2,LL3,……,LL(N-M+1)用于与右移焦点聚焦扫描线LR1,LR2,……,LR(N-M)做运动估计处理。
[0023] 上述技术方案中,所述步骤2)进一步包括:
[0024] 步骤2-1)、计算第一条左移焦点聚焦扫描线LL2与第一条右移焦点聚焦扫描线LR1之间的时延,得到第一时延值D1;
[0025] 步骤2-2)、计算第i条左移焦点聚焦扫描线LL(i+1)与第i条右移焦点聚焦扫描线LR(i)之间的时延,得到第i个时延值Di;
[0026] 步骤2-3)、第一条第一低分辨率扫描线L1不做补偿,直接作为第一条第二低分辨率扫描线,记作CL1;
[0027] 步骤2-4)、对于第二条第一低分辨率扫描线L2,利用第一时延值D1进行延时处理,得到补偿后的第二条第二低分辨率扫描线CL2;
[0028] 步骤2-5)、利用第1至第i时延值之和D1+D2+…+Di对第i+1条第一低分辨率扫描线Li+1进行延时处理,得到补偿后的第i+1条第二低分辨率扫描线CLn+1,直至i的值达到N-M,得到补偿后的第N-M+1条第二低分辨率扫描线CLN-M+1;
[0029] 其中,CL1、CL2……CLN-M+1称为第二低分辨率扫描线。
[0030] 上述技术方案中,所述步骤3)进一步包括:
[0031] 步骤3-1)、对第一条第二低分辨率扫描线进行第二次波束形成,即将VS1开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与第一条第二低分辨率扫描线上相应采样点相加,得到第一条高分辨率扫描线H1;
[0032] 步骤3-2)、对第二条第二低分辨率扫描线进行第二次波束形成,即将VS2开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与第二条第二低分辨率扫描线上相应采样点相加,得到第二条高分辨率扫描线H2;
[0033] 步骤3-3)、重复上述处理,直到对第N-M+1条第二低分辨率扫描线进行第二次波束形成,即将VSN-M+1开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与第N-M+1条第二低分辨率扫描线上相应采样点相加,得到第N-M+1条高分辨率扫描线HN-M+1。
[0034] 上述技术方案中,将虚拟源开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与虚拟源对应第二低分辨率扫描线上相应采样点相加,得到高分辨率扫描线的具体计算公式为:
[0035]
[0036] 其中, 表示第二低分辨率扫描线 上延时td(k)对应的值,w为加权函数,其值随深度z和扫描线的水平位置xk而改变;K是深度z的函数,K(z)表示虚拟源开角内的低分辨率扫描线数量;
[0037] K(z)的计算公式如下:
[0038]
[0039] 其中,zv为虚拟源深度,α为虚拟源的开角,L(z)为深度z处的声场宽度,Δ为虚拟源间距;
[0040] 延时td(k)的计算公式为:
[0041]
[0042] 其中,zv为虚拟源深度, 为(x,z)点到线 上的虚拟源 的距离;+号表示(x,z)点深度大于zv,-号表示其深度小于zv,c为声速。
[0043] 本发明的优点在于:
[0044] 本发明通过对第一阶段波束形成得到的低分辨率扫描线间运动进行估计和补偿,在第二阶段波速形成得到高分辨率扫描线,有效地解决了现有技术在存在轴向位移时成像质量较差的问题。同时,由于发射次数没有增加,在充分计算能力的保证下,成像的帧频不会降低。

附图说明

[0045] 图1为本发明合成孔径超声成像运动补偿方法的流程图;
[0046] 图2为本发明合成孔径超声成像运动补偿方法高分辨率扫描线形成原理图;
[0047] 图3为本发明合成孔径超声成像运动补偿方法第一低分辨率扫描线及第二低分辨率扫描线形成原理图。

具体实施方式

[0048] 现结合附图对本发明作进一步的描述。
[0049] 本发明方法所适用的合成孔径超声换能器阵列共有N个阵元,这些阵元形成一等间距阵列,阵列中阵元间距记为Δ,各个阵元均可发射信号与接收信号。
[0050] 以上述合成孔径超声换能器阵列为基础,对本发明的方法做进一步说明。
[0051] 参考图1-图3,本发明的带运动补偿的合成孔径超声成像方法包括:
[0052] 步骤101,采用发射、接收都是单焦点,且发射焦点与接收焦点相同的线性扫描方式得到第一低分辨率扫描线;采用偏移阵元接收的方式得到运动估计扫描线。
[0053] 所述合成孔径超声换能器阵列中的阵元在发射信号时需选取阵列中的部分阵元做信号发射操作,包括:首先采用合成孔径超声换能器阵列中的第1到第M(M为有效孔径,即表示合成孔径超声换能器阵列中一次操作可以使用的阵元,M
[0054] 步骤1011,当i=1,即第T1次发射时,用所述阵列中的第1到第M个阵元发射信号,发射焦点为VS1;由所述阵列中的第1到第M+1个阵元接收信号(接收阵元的数目较发射阵元的数目多一个),每个阵元的采样点数为P,阵列中的阵元在接收 信号时需经历正常接收与右移接收两个阶段;
[0055] 正常接收:对第1到第M阵元所接收的信号进行波束形成,接收焦点也为VS1,得到第1条第一低分辨率扫描线L1,该扫描线有P个采样点;
[0056] 右移接收:对第2到第M+1阵元所接收的信号进行波束形成,接收焦点也为VR1,VR1=VS1+Δ,得到第1条右移焦点扫描线LR1,该扫描线有P个采样点。
[0057] 步骤1012,当1
[0058] 左移接收:对第i-1到第M+i-2阵元所接收的信号进行波束形成,接收焦点也为VLi,VLi=VSi-Δ,得到第i-1条左移焦点扫描线LLi,该扫描线有P个采样点;
[0059] 正常接收:对第i到第M+i-1阵元所接收的信号进行波束形成,接收焦点也为VSi,得到第i条第一低分辨率扫描线Li,该扫描线有P个采样点;
[0060] 右移接收:对第i+1到第M+i阵元所接收的信号进行波束形成,接收焦点也为VRi,VRi=VSi+Δ,得到第i条右移焦点扫描线LRi,该扫描线有P个采样点。
[0061] 步骤1013,当i=N-M+1时,即第TN-M+1次发射时,用所述阵列中的第N-M+1到第N个阵元发射信号,发射焦点为VSN-M+1,由所述阵列中的第N-M到第N个阵元接收信号(接收阵元的数目较发射阵元的数目多一个),每个阵元的采样点数为P,阵列中的阵元在接收信号时需经历左移接收与正常接收两个阶段;
[0062] 左移接收:对第N-M到第N-1阵元所接收的信号进行波束形成,接收焦点也为VLN-M+1,VLN-M+1=VSN-M+1-Δ,得到第N-M条左移焦点扫描线LL(N-M+1),该扫描线有P个采样点;
[0063] 正常接收:对第N-M+1到第N阵元所接收的信号进行波束形成,接收焦点也为VSN-M+1,得到第N-M+1条第一低分辨率扫描线LN-M+1,该扫描线有P个采样点;
[0064] 上述步骤中所涉及的VS1,VS2,……,VSN-M+1均可称为虚拟源。如图3所示,虚拟源所发出信号的作用范围被称为虚拟源开角α,其计算公式为α=2*arctan(LA/2ZV),其中,zv为虚拟源深度,LA为有效孔径M的长度,即一次发射时使用的孔径的长度。
[0065] 若阵列中阵元间距为Δ,则有VRi=VSi+Δ,VLi=VSi-Δ。
[0066] 步骤1014,依次重复上述步骤,直至得到N-M+1条第一低分辨率扫描线(N为阵元总数)。其中,第一低分辨率扫描线L1,L2,……,LN-M+1构成一幅低分辨率图像LRI;所得到的左移焦点聚焦扫描线LL2,LL3,……,LL(N-M+1)用于与右移焦点聚 焦扫描线LR1,LR2,……,LR(N-M)做运动估计处理。
[0067] 步骤102,对左移焦点聚焦扫描线和右移焦点聚焦扫描线间运动进行估计,对第一低分辨率扫描线进行补偿,从而得到第二低分辨率扫描线;
[0068] 具体地,所述步骤102的具体处理方法如下:
[0069] 步骤1021,计算第一条左移焦点聚焦扫描线LL2与第一条右移焦点聚焦扫描线LR1之间的时延,得到第一时延值D1;
[0070] 步骤1022,计算第i条左移焦点聚焦扫描线LL(i+1)与第i条右移焦点聚焦扫描线LR(i)之间的时延,得到第i个时延值Di;
[0071] 步骤1023,第一条第一低分辨率扫描线L1不做补偿,直接作为第一条第二低分辨率扫描线,记作CL1;
[0072] 步骤1024,对于第二条第一低分辨率扫描线L2,利用第一时延值D1进行延时处理,得到补偿后的第二条第二低分辨率扫描线CL2。
[0073] 步骤1025,利用第1至第i时延值之和D1+D2+…+Di对第i+1条第一低分辨率扫描线Li+1进行延时处理,得到补偿后的第i+1条第二低分辨率扫描线CLn+1,直至i的值达到N-M,得到补偿后的第N-M+1条第二低分辨率扫描线CLN-M+1;
[0074] 其中,CL1、CL2……CLN-M+1称为第二低分辨率扫描线。
[0075] 步骤103,将虚拟源开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与虚拟源对应第二低分辨率扫描线上的相应采样点相加,得到高分辨率扫描线,进而得到高分辨率图像。
[0076] 具体地,所述步骤103的具体处理方法如下:
[0077] 步骤1031,对第一条第二低分辨率扫描线进行第二次波束形成,即将VS1开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与第一条第二低分辨率扫描线上相应采样点相加,得到第一条高分辨率扫描线H1;
[0078] 步骤1032,对第二条第二低分辨率扫描线进行第二次波束形成,即将VS2开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与第二条第二低分辨率扫描线上相应采样点相加,得到第二条高分辨率扫描线H2;
[0079] 步骤1033,重复上述处理,直到对第N-M+1条第二低分辨率扫描线进行第二次波束形成,即将VSN-M+1开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与第N-M+1条第二低分辨率扫描线上相应采样点相加,得到第N-M+1条高分辨率扫描线HN-M+1。
[0080] 在上述步骤中,将虚拟源开角内其它第二低分辨率扫描线上的相应采样点乘以加权值后与虚拟源对应第二低分辨率扫描线上相应采样点相加,得到高分辨率扫描 线的具体计算公式为:
[0081]
[0082] 其中, 表示第二低分辨率扫描线 上延时td(k)对应的值,w为加权函数,其值随深度z和扫描线的水平位置xk而改变;K是深度z的函数,K(z)表示虚拟源开角内的低分辨率扫描线数量;
[0083] K(z)的计算公式如下:
[0084]
[0085] 其中,zv为虚拟源深度,α为虚拟源的开角,L(z)为深度z处的声场宽度,Δ为虚拟源间距;
[0086] 延时td(k)的计算公式为:
[0087]
[0088] 其中,zv为虚拟源深度, 为(x,z)点到线 上的虚拟源 的距离;+号表示(x,z)点深度大于zv,-号表示其深度小于zv,c为声速。
[0089] 最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。