[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] 最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。