关于超声和光声中的延迟求和算法的详细说明 Detail description of DAS algorithm

背景 background

In LED-base Photoacoustic imaging,DAS algorithm was always used for the image reconstruction. This note is about illustrate the details of DAS algorithm.
LED-based Photoacoustic imaging technology often use a linear probe, generally they have 128 elements in this probe,when the exciting light incident to the samples, if there is an absorb in the sample then it will generate the ultrasound signals and all the 128 elements will detected this signal, but due to the different distance between this absorb to every element, it will have a time-delay in each channel, the small distance will get the signal early and the large distance will get it later.
Therefore if you wanna build the image, you need to know the signal amplitude in every pixel, and each pixel was detected by 128-channels, so you need to sum all the 128-channel’s value together to get the right amplitude of every-pixel.

在LED光声中最基本的图像重建算法就是延迟求和算法DAS，通常我们都会使用128陈元的线性探头。当激励光源照射进成像样品之后，如果样品中存在一个吸收体，那么这个吸收体产生的光声信号会被这128个阵元都检测到，但由于每个阵元到这个吸收体的距离不一样，因此导致每个通道监测到的信号有不同的时间延迟，离吸收体近的阵元将会早一点的接收到超声信号，离得远点的就会晚一点接收到这个信号，但是这个信号就是这个吸收体的，只是在接收的时间上有一定的差异，因此为了正确得到这个吸收体产生的光声信号的幅值应该将这128个阵元监测的信号加起来。

存在问题（个人思考）

There is a condition on this algorithm，because if you just based on the distance difference to determine where is the absorb it will make some confuse problems, we all know that the incident light is scattered as a ball in the tissue, consequently，the generated signal can not existed on a plane, it should positioned as a ball, but in generally we assume all the signal was on a plane, then use the DAS algorithm to build the image.

有一个问题就是入射光进入组织以后会散射向各个方向，因此我们接收到的信号不应该只是来自一个平面上的，而应该是一个立体范围内的，但是为了方便计算成像，我们假设用线阵探头进行光声成像的时候接收到的光声信号来自与一个平面即我们通常所说的B-scan,B-plane.而这种线阵探头之所以成像速度快就是其可以一次就完成一个B-scan,不像Photoacoustic microscopy 需要先扫描A-line 然后组成B-plane.
假设所有的吸收体都来自这个B-plane 以后我们就可以根据这个B-plane每个像素距离每个阵元的距离来判断信号的延迟时间，然后根据每个通道接收到信号的时间差异提取出每个通道处信号的具体位置进行叠加，即可以重建出图像。

算法说明及代码示例

DAS延迟和算法，不难从字面理解，计算延迟+数据求和。

（1）延迟时间：通过计算从目标点P到传感器接收单元的距离，再根据传输介质的速度，以及超声传感器阵列的采集频率，求光声信号从目标点到信号接收所需要的时间，由此定位出所收集的数据点在传感器接收数据集中的位置，找到目标点在某个传感器单元所对应的时间的接收信号。

（2）数据求和：在确定了传感器阵列中每个单元所对应的接收时间及对应的信号强度，将数据全部提前进行求和，即为该目标点的信号强度。

根据这个方法，可以求出图像中每个像素所对应的所有阵列信