超声B mode聚焦成像仿真（ultrasound B mode fcous imaging simulation）

上篇文章超声平面波复合成像仿真（ultrasound plane_wave compound imaging simulation）使用Field II仿真了平面波成像，相比平面波成，聚焦成像像在焦区处分辨率、信噪比更强、通过聚焦穿透力也有提升，平面波发射一次获取整个深度信息，发射聚焦需要对每一条线上进行聚焦，然后做波束合成获取一条线上的信息，需要发射N次，帧率比平面波慢。

在仿真之前先介绍一下聚焦成像，通过对发射和接收的信号进行延时叠加求和（delay and sum ,DAS）可以达到聚焦成像的目的，聚焦成像中的波束合成包含两种模式，发射聚焦及接收聚焦。

理想上希望对每一个深度上的点进行发射聚焦即发射动态聚焦，由于声波传播需要时间，一副图像的形成包含N*M个点，需要上万次发射，无法实时成像，而实时成像又是我们超声成像（相比CT、MRI成像）最重要的优势，因此发射聚焦往往只在感兴趣区域进行1~3次，超过一个发射焦点就叫多焦点成像。对于有些厂家宣传的十几个以上的焦点只是为了满足招标要求去做的，临床上很少甚至不会使用。

根据互易定理，接收聚焦与发射聚焦的原理是一样的，都是通过对信号进行延时达到聚焦目的。对每一个深度上的点进行聚焦称为接收动态聚焦。如果只对一个感兴趣深度做接收聚焦就叫做接收定点聚焦，目前商业机器上使用的大多是接收动态聚焦。

下图是一副大家经常看到的波束合成的示意图

从图上可以看出发射延时控制各阵元发出的波束同时到达聚焦点(x,z)(x,z)(x,z)，对接收到的回波信号通过计算接收延时进行接收聚焦，获取波束合成信号。

本文仿真的是定点发射聚焦，发射聚焦通过发射延时控制，当然也可以通过Field II自带函数去计算，感兴趣的话可以看下超声成像发射声场仿真（Ultrasound Emit Field Simulation），这篇文章介绍了通过控制发射延时形成不同成像模式的发射声场；对于动态接收聚焦可以通过Field II自带函数去计算，为了更深刻了解波束合成过程，这里使用仿真代码去实现。在此处的仿真中仅仿真正常发射图像，对于偏转仿真仅提供结果以验证结果是否正确，因为偏转图像需要另一种技术空间复合（spatial compound imaging），这里不做过多介绍。

与平面波成像仿真一样，本文仅对聚焦波束合成进行仿真，不涉及额外处理，如发射变迹、接收变迹等。仿真工具使用 Filed II，探头包括线阵、凸阵、相控阵。

Field II仿真包含四个部分，

参数配置
AD数据获取
波束合成
图像显示

代码与上篇文章超声平面波复合成像仿真（ultrasound plane_wave compound imaging simulation）流程一样，只是发射延时、接收延时计算方式不同，下面对线阵的延时计算做说明，相应代码略过。

一线阵聚焦成像（fcous imaging of linear array）

上图为线阵非偏转发射延时计算，波束合成中的延时由上图中的两部分构成发射延时TxdelayTxdelayTxdelay和接收延时RvdelayRvdelayRvdelay，在仿真过程中需要减去数据采样起始时间stimestimestime，在实际工程中stimestimestime可以代表波束合成参数准备时间、声束在透镜中传播时间、校正发射脉冲产生时间等，具体根据工程实现的方式去计算。对于延时时间最终转化为：
delay=Txdelay+Rxdelay−stimedelay =Txdelay+Rxdelay-stime delay=Txdelay+Rxdelay−stime
在Field II中stimestimestime由下面的函数返回

[scat, stime] = calc_scat_multi (Th, Rh, points, amplitudes);

Txdelay=(sqrt((xn−xline)2+(zn−zline2))/c=zn/cTxdelay = (sqrt((x_n-x_line)^2+(z_n-z_line^2))/c = z_n/c Txdelay=(sqrt((xn−xline)2+(zn−zline2))/c=zn/c

Rxdelay=(sqrt((xn−x)2+(zn−z)2))/cRxdelay = (sqrt((x_n-x)^2+(z_n-z)^2))/c Rxdelay=(sqrt((xn−x)2+(zn−z)2))/c

上述即延时计算公式，ccc为声速，(xn,yn)(x_n,y_n)(xn,yn)为接收线位置，(x,z)(x,z)(x,z)为阵元坐标位置。

Field II仿真代码：

参数设置
AD数据获取

这里对于每一条线的发射孔径固定为64，通道并没有全部使用，下图是每一条线对应的发射孔径，注意左右两侧

下图是位于最左侧的线 64通道发射。128通道接收和64通道发射，32通道接收时获取到的的AD数据示意图，从图像可以看出通道打开的位置是正确的。
波束合成
图像显示

结果为：

这里附上偏转结果:

图像上边分别为-15, 0, 15偏转角度波束合成后图像，下侧为坐标系转换后结果，最后经过空间复合：

从仿真效果上看靶点位置正确，3 cm 聚焦处比其他位置分辨率及信噪比更好，偏转图像经过空间复合后伪像减少。

二凸阵聚焦成像（fcous imaging of curve array）

仿真流程与线阵一样，不同地方在于发射延时计算方式与线阵不一样，波束合成中延时的距离计算需要由极坐标系转化为笛卡尔坐标系，坐标系转化后应该如下图：

由于凸阵探头发射频率低，成像深度比线阵深，聚焦点位置放在 6 cm，仿真的靶点位置重新定义

按照线阵的流程配置，最终结果为：

偏转图像空间复合后结果：

三相控阵聚焦成像（fcous imaging of phase array）

相控阵探头发射线永远位于中心线处，聚焦延时修改为

xdc_center_focus(Th,[0,0,0]);
xdc_focus(Th,0,[fcous_x,0,fcous_z]);%  每条发射线聚焦处位置

按照线阵的流程配置，最终结果为：

相控阵一般是对心脏部位进行成像，帧率要求高，在聚焦成像中不会对相控阵进行空间复合。

总结

通过Filed II仿真了线阵、凸阵、相控阵的聚焦波束合成，有助于对聚焦成像有个深入的了解。对于超声机器中的偏转成像、线阵梯形成像、凸阵扩展成像都可以通过本文的方法进行相应扩展仿真。

上文仅针对接收动态聚焦做了仿真，对于定点接受聚焦只需要计算一次接受延时就可以确定一条线上所有点的位置。

目前商用机器中的聚焦波束合成方法不断发展，聚焦成像不仅能在焦区获取良好的分辨率，在焦区外也可以获得优异的图像，如philips的nSIGHT、迈瑞的Dynamic Pixel Focusing