专利名称：基于分段动态变迹系数的超声成像波束合成方法及其装置的制作方法
技术领域：
本发明涉及超声成像技术领域，具体涉及一种基于分段动态变迹系数的超声成像 波束合成方法及其装置。
背景技术：
在超声成像系统中，通过引入幅度变迹函数来降低合成声束中旁瓣的等级，以抑 制图像中的伪像影响，而通常单一的幅度变迹函数只能作用整个图像的局部区域，使得在 该区域具有较好的图像分辨率，在该区域外，图像的分辨率明显下降。分段动态变迹控制是近年来，针对单一幅度变迹控制的缺点，而提出的改进方式。 现有技术的分段动态变迹控制系数通过实验，将不同变迹函数应用到不同接收焦点范围段 的超声成像技术中，以获得较理想的合成声束，由于使用了不同的变迹函数，而且在合成理 想声束前需要反复实验获得不同变迹函数对应的最佳分段范围，增加了声束合成的运算复 杂度；现有技术的变迹函数的系数中未考虑接收焦点深度对合成声束的影响，对于不同深 度的接收焦点，同一变迹函数获得的一系列变迹系数为一组不变值，无法获得拥有理想旁 瓣和理想主瓣关系的合成声束。
发明内容
有鉴于此，为了解决上述问题，本发明公开了一种基于分段动态变迹系数的超声 成像波束合成方法，使用同一形式变迹函数的系数，在系数中纳入接收焦点深度为声束合 成参考因素，在合成声束前可直接计算出分段动态变迹系数。本发明的目的是这样实现的基于分段动态变迹系数的超声成像波束合成方法， 根据接收焦点深度生成各阵元对应的动态变迹系数，将各阵元的回波数据与该阵元对应的 动态变迹系数相乘，然后对各阵元相乘后的回波数据进行延时叠加，获得该接收焦点深度 的合成波束，所述动态变迹系数通过如下步骤获得1)以阵元列与其中心轴线的交点为起点，沿探测深度方向将阵元列的中心轴线分 段，各分段点为接收焦点，获取接收焦点到阵元列与中心轴线交点的距离，作为各接收焦点 深度；2)计算各接收焦点与阵元列中各阵元的连线和该阵元法线的夹角；3)根据步骤2)获得的夹角构建与夹角大小呈反比关系的分段动态变迹系数。进一步，步骤1)具体为沿探测深度方向将阵元列中心轴均勻分段，接收焦点为等分线与阵元列中心轴线 的交点，得到接收焦点深度为zm = depth*m/M，m = 1，2，3，M ；其中cbpth为探测总深度值，M为均勻分段数，zffl为第m个接收焦点的深度值。进一步，步骤2)具体为21)计算第n个阵元的法线与深度方向第一夹角；
22)计算第m个接收焦点与第n个阵元的连线和深度方向第二夹角；23)第m个接收焦点与第n个阵元的连线和该阵元法线的第三夹角为e = ^ +3 ，其中a 为第一夹角值，3 为第二夹角值，为第三夹角值；24)重复步骤21)至23)遍历所有阵元和所有接收焦点，获取各接收焦点与各阵元 的连线和该阵元法线的第三夹角；进一步，阵元以阵元列中心轴线为对称轴，分布于探头两端并呈轴对称排列；进一步，阵元列所排列的探头为凸阵探头时，步骤21)中第n个阵元的法线与深度方向的第一夹角为 二 r
d ( 7V + 1、
n--
n = 1,2,3,…，N ；
2
“ ‘,其中，d为阵元间距值，R为探头曲率半径值，N为阵元总数；步骤22)中第m个接收焦点与第n个超声探头阵元的连线和深度方向第二夹角 为
nR-sin anPnm = arCtan--TT,-x进一步，步骤3)中利用汉宁窗形式函数计算第n个阵元相对于第m个接收焦点的 分段动态变迹系数为wM = a+bcos O ‘ tan 0 J ；其中，a+b = 1，0 为第m个接收焦点与第n个阵元的连线和该阵元法线的第三 夹角值，w 为分段动态变迹系数值；进一步，a取值为0.5，b取值为0.5。还公开了使用基于分段动态变迹系数的超声成像波束合成方法的装置，包括n个 阵元对应的n个双口 RAM、DSP芯片、动态变迹系数读取控制器、地址译码器、n个乘法器和 聚焦模块；聚焦模块包括n个延时模块和累加器；DSP芯片的数据输出端分别与n个阵元对 应的n个双口 RAM和动态变迹系数读取控制器的数据输入端连接；DSP芯片的地址输出端 与地址译码器的地址输入端连接，地址译码器的地址片选端与n个阵元对应的n个双口 RAM 和动态变迹系数读取控制器的片选端连接；动态变迹系数读取控制器的读取控制端分别与 n个阵元对应的n个双口 RAM的读取被控端连接，n个阵元对应的n个双口 RAM输出的数据 通过n个乘法器和n个延时模块输入累加器，其中n为大于等于1的整数。进一步，装置对超声成像波束合成具体步骤为1)DSP芯片将计算出的各接收焦点对应的各组分段动态变迹系数输出存储到n个 阵元对应的n个双口 RAM中；2)DSP芯片把各组分段动态变迹系数与各接收焦点的对应关系写入到变迹系数读 取控制器中；3)当超声信号发射完毕，n个阵元接收回波数据，动态变迹系数控制器根据回波 数据对应的接收焦点，依次输出各组分段动态变迹系数；4)各组分段动态变迹系数与n个阵元的回波数据通过乘法器进行相乘变迹，生成 变迹数据；5)n个阵元的回波数据经过变迹后，经聚焦模块对n个阵元的变迹数据分别进行
5延时，然后通过累加器进行叠加，实现波束合成。 本发明的有益效果是使用同一形式变迹函数的系数，在系数中纳入接收焦点深 度为声束合成参考因素，在合成声束前直接计算出分段动态变迹系数，提高了声束合成效 率；由于引入接收焦点深度为参考因素，在不同分段所对应的变迹函数系列各不相同，克服 了单一的幅度变迹函数只能作用整个图像的局部区域的缺点，降低合成声束中旁瓣的同时 提高了超声成像的分辨率。
为了使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作进 一步的详细描述图1示出了阵元列以及接收焦点结构关系；图2示出了凸阵探头的分段动态系数值分布示意图；图3示出了动态变迹与波束合成的原理框图。
具体实施例方式以下将对本发明的优选实施例进行详细地描述。本发明公开了基于分段动态变迹系数的超声成像波束合成方法，在接收焦点深度 生成各阵元对应的动态变迹系数，各阵元的回波数据与该阵元对应的动态变迹系数相乘， 对各阵元相乘后的回波数据进行延时叠加，获得该接收焦点深度的合成波束，如图1所示， 所述动态变迹系数通过如下步骤获得1)沿探测深度方向1,将阵元列中心轴Z均勻分段，接收焦点为等分线与阵元列中 心轴线Z的交点，得到接收焦点深度为zm = depth*m/M，m = 1，2，3，...，M;其中depth为探测总深度值，M为均勻分段数，zffl为第m个接收焦点的深度值；2)计算不同接收焦点与阵元列中某阵元的连线和该阵元法线的夹角探头为凸阵 探头时，计算第n个阵元Qn的法线ln与深度方向1第一夹角an:
d ( iV + 1)其中，d为阵元间距值，R为探头曲率半径值，NN为阵元总数；计算第m个接收焦点Fm与第n个阵元Qn的连线1。和深度方向lv的第二夹角3 ■ 为
nisin aPnm =耵汶肪-—r~-——r
第m个接收焦点Fm与第n个阵元Qn的连线1。和该阵元法线ln的第三夹角9 为9 =其中a n为第一夹角值，0 为第二夹角值，0 为第三夹角值；遍历所有阵元和所有接收焦点，获取不同接收焦点与不同阵元的连线和该阵元法线的第三夹角；3)根据步骤2)获得的夹角构建与夹角大小呈反比关系的分段动态变迹系数。利用汉宁窗形式函数计算第n个阵元的单一幅度变迹系数公式w =a + 6cos(-^^) n = 1,2,3, ...，N;
N,其中，N为阵元总数，a+b = 1, wn为单一的幅度变迹系数值；6 替换上式，构建分段动态变迹系数计算公式为wM = a+bcos O ‘ tan 0 J ；其中，Wnm为分段动态变迹系数值；可以取a = 0.5，b = 0.5带入上式有wM = 0. 5+0. 5cos O tan 0 J ；遵循接收焦点与超声探头阵元列中阵元的连线和该阵元法线的第三夹角越大，分 段动态变迹系数应该越小的原则建立的分段动态变迹函数可实现降低合成声束中旁瓣的 等级的同时，不影响合成声束的分辨率。按照这个规律，可以将接收焦点与超声探头阵元列 中阵元的连线和该阵元法线的夹角变形应用到将余弦窗形式、海明窗形式、Blackman窗形 式等单一幅度变迹函数向分段动态变迹函数转换变形中，下列列出了除汉宁窗形式外的， 另外几种动态变迹系数获取式余弦窗海明窗wM = 0. 54+0. 46cos O ‘ tan 0 JBlackman 窗= 0. 42+0. 5cos ( n ‘ tan 0 J +0. 08cos (2 n ‘ tan 0 ■)本发明还公开了使用基于分段动态变迹系数的超声成像波束合成方法的装置，包 括n个阵元对应的n个双口 RAM4、DSP芯片1、动态变迹系数读取控制器3、地址译码器2、n 个乘法器5和聚焦模块6 ；聚焦模块6包括n个延时模块7和累加器8 ；DSP芯片1的数据 输出端分别与n个阵元对应的n个双口 RAM4和动态变迹系数读取控制器3的数据输入端 连接；DSP芯片1的地址输出端与地址译码器2的地址输入端连接，地址译码器2的地址片 选端与n个阵元对应的n个双口 RAM4和动态变迹系数读取控制器3的片选端连接；动态变 迹系数读取控制器3的读取控制端分别与n个阵元对应的n个双口 RAM4的读取被控端连 接，n个阵元对应的n个双口 RAM4输出的数据通过n个乘法器5和n个延时模块7输入累 加器8。装置通过以下步骤完成超声成像波束合成1)DSP芯片1将计算出的各接收焦点对应的各组分段动态变迹系数输出存储到n 个阵元对应的n个双口 RAM4中；2)DSP芯片1把各组分段动态变迹系数与各接收焦点的对应关系写入到变迹系数 读取控制器3中；
3)当超声信号发射完毕，n个阵元接收回波数据，动态变迹系数控制器根据回波 数据对应的接收焦点，依次输出各组分段动态变迹系数；4)各组分段动态变迹系数与n个阵元的回波数据通过乘法器5进行相乘变迹，生 成变迹数据；5)n个阵元的回波数据经过变迹后，经聚焦模块6对n个阵元的变迹数据分别进行 延时，然后通过累加器8进行累加，实现波束合成。图2示出了凸阵探头的分段动态系数值分布。其中探头参数为阵元总数为64， 阵元间距为0. 48mm，探头曲率半径60mm ；超声声速1540m/s，脉冲中心频率3. 5MHz，凸阵探 头在深度20-200mm上分为8段，得到的分段动态变迹系数如图2。图2中横坐标为阵元序 号，纵坐标为变迹系数值，8组曲线由上之下依次对应接收焦点F8F8至&处的变迹系数值。本发明在超声成像技术领域中，通过对探头接收的回波数据进行动态加权变迹控 制，从而有效抑制声束旁瓣，并在整个探测深度上保持良好的分辨率，能够整体提高超声图 像的质量。以上所述仅为本发明的优选并不用于限制本发明，显然，本领域的技术人员可以 对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修 改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和 变型在内。
权利要求
基于分段动态变迹系数的超声成像波束合成方法，其特征在于根据接收焦点深度生成各阵元对应的动态变迹系数，将各阵元的回波数据与该阵元对应的动态变迹系数相乘，然后对各阵元相乘后的回波数据进行延时叠加，获得该接收焦点深度的合成波束，所述动态变迹系数通过如下步骤获得1)以阵元列与其中心轴线的交点为起点，沿探测深度方向将阵元列的中心轴线分段，各分段点为接收焦点，获取接收焦点到阵元列与中心轴线交点的距离，作为各接收焦点深度；2)计算各接收焦点与阵元列中各阵元的连线和该阵元法线的夹角；3)根据步骤2)获得的夹角构建与夹角大小呈反比关系的分段动态变迹系数。
2.如权利要求1所述的基于分段动态变迹系数的超声成像波束合成方法，其特征在 于步骤1)具体为沿探测深度方向将阵元列中心轴均勻分段，接收焦点为等分线与阵元列中心轴线的交 点，得到接收焦点深度为zm = depth 氺 m/M, m = 1,2,3, ".,Μ;其中depth为探测总深度值，M为均勻分段数，Zffl为第m个接收焦点的深度值。
3.如权利要求1所述的基于分段动态变迹系数的超声成像波束合成方法，其特征在 于步骤2)具体为21)计算第η个阵元的法线与深度方向第一夹角；22)计算第m个接收焦点与第η个阵元的连线和深度方向第二夹角；23)第m个接收焦点与第η个阵元的连线和该阵元法线的第三夹角为θΜ = αη+βΜ, 其中α 为第一夹角值，β Μ为第二夹角值，θ Μ为第三夹角值；24)重复步骤21)至23)遍历所有阵元和所有接收焦点，获取各接收焦点与各阵元的连 线和该阵元法线的第三夹角。
4.如权利要求3所述的基于分段动态变迹系数的超声成像波束合成方法，其特征在 于阵元以阵元列中心轴线为对称轴，分布于探头两端并呈轴对称排列。
5.如权利要求4所述的基于分段动态变迹系数的超声成像波束合成方法，其特征在 于阵元列所排列的探头为凸阵探头时，步骤21)中第η个阵元的法线与深度方向的第一夹角为 其中，d为阵元间距值，R为探头曲率半径值，N为阵元总数；步骤22)中第m个接收焦点与第η个超声探头阵元的连线和深度方向第二夹角为
6.如权利要求1至5任一项所述的基于分段动态变迹系数的超声成像波束合成方法， 其特征在于步骤3)中利用汉宁窗形式函数计算第η个阵元相对于第m个接收焦点的分段 动态变迹系数Wnm = a+bcos ( π · tan θ ；其中，a+b = 1，θ ■为第m个接收焦点与第n个阵元的连线和该阵元法线的第三夹角 值，Wnm为分段动态变迹系数值。
7.如权利要求6所述的基于分段动态变迹系数的超声成像波束合成方法，其特征在 于a取值为0. 5，b取值为0.5。
8.使用如权利要求1至7任意一项所述的基于分段动态变迹系数的超声成像波束合 成方法的装置，其特征在于包括η个阵元对应的η个双口 RAM、DSP芯片、动态变迹系数读 取控制器、地址译码器、η个乘法器和聚焦模块；聚焦模块包括η个延时模块和累加器；DSP 芯片的数据输出端分别与η个阵元对应的η个双口 RAM和动态变迹系数读取控制器的数据 输入端连接；DSP芯片的地址输出端与地址译码器的地址输入端连接，地址译码器的地址 片选端与η个阵元对应的η个双口 RAM和动态变迹系数读取控制器的片选端连接；动态变 迹系数读取控制器的读取控制端分别与η个阵元对应的η个双口 RAM的读取被控端连接，η 个阵元对应的η个双口 RAM输出的数据通过η个乘法器和η个延时模块输入累加器，其中 η为大于等于1的整数。
9.如权利要求8所述的装置，其特征在于装置对超声成像波束合成具体步骤为DDSP芯片将计算出的各接收焦点对应的各组分段动态变迹系数输出存储到η个阵元对应的η个双口 RAM中；2)DSP芯片把各组分段动态变迹系数与各接收焦点的对应关系写入到变迹系数读取控 制器中；3)当超声信号发射完毕，η个阵元接收回波数据，动态变迹系数控制器根据回波数据 对应的接收焦点，依次输出各组分段动态变迹系数；4)各组分段动态变迹系数与η个阵元的回波数据通过乘法器进行相乘变迹，生成变迹 数据；5)η个阵元的回波数据经过变迹后，经聚焦模块对η个阵元的变迹数据分别进行延时， 然后通过累加器进行叠加，实现波束合成。
全文摘要
本发明公开了基于分段动态变迹系数的超声成像波束合成方法，沿探测深度方向将阵元列的中心轴线分段生成各接收焦点；根据接收焦点与超声探头阵元列中阵元的连线和阵元法线的夹角构建与夹角大小呈反比的分段动态变迹系数。本发明在系数中纳入接收焦点深度为参考因素，在合成声束前直接生成分段动态变迹系数；不同分段所对应不同变迹函数，克服了单一幅度变迹函数只能作用局部区域的缺点，降低合成声束旁瓣等级，提高超声成像分辨率；还提出使用该方法的装置，通过DSP计算动态变迹系数，初始化存储变迹系数的双口RAM和动态变迹系数读取控制器，实现各阵元回波数据的动态变迹，然后经过聚焦模块，对各个阵元的回波数据进行延时叠加，实现波束合成。
文档编号A61B8/00GK101940479SQ201010228089
公开日2011年1月12日 申请日期2010年7月15日 优先权日2010年7月15日
发明者何为, 张莉, 王平, 许琴, 陈民铀 申请人:重庆大学