Field II 超声线阵成像系列1——复合平面波成像

       平面波成像与传统超声聚焦成像相比，可通过一次发射可以获取整个感兴趣区域图像信息，具有很高的帧频，可用来做血流成像、ARFI成像等。但由于无聚焦，图像对比度和分辨率降低。为了克服这一限制，可以通过发射多个不同角度的平面波，通过相干叠加得到完整图像，在一定角度范围内平面波复合的数量越多，其成像质量越好，但相应的也会是其帧频下降。

       本文对Filed II对线阵复合平面波成像进行仿真，通过发射延时控制平面波角度偏转，对于发射延时的计算可参考Field II 超声相控阵成像系列1——复合平面波成像(更正)。

 图1

       图1为相控阵平面波偏转发射的示意图，对于无偏转角度，可使各阵元发射延时全为0，对于偏转发射（包括正常发射）延时计算为：

tx_d=pitch*[0:number_of_ele-1]*sin(steer);

     偏转角度表示为负往左偏转，角度为正往右偏转。

     在波束合成中时，采用最常用的延时叠加算法，延时由两部分构成发射延时Tx 和接收延时Rx ，在成像过程中需要减去起始采样时间tstart，在实际工程中tstart可以代表波束合成参数准备时间、声束在透镜中传播时间等，最终成像的延时表示为：

delay_t=Tx+Rx-tstart;

Tx=(z(k)*cos(TXangle) + (x(k)+ halfaper)*sin(TXangle))/c; % TX distance

Rx = (sqrt((xT-x(k)).^2 + z(k).^2))/c; % RX distance

 halfaper = sign(TXangle)*xT(end); 

    其中tstart可有Field II中 calc_scat_multi 函数提供, c为声速，[x(k),z(k)] 为像素点坐标。xT为换能器阵元坐标。

仿真程序：

1.计算平面波数据
 

for i=1:userset.angnum
    angle=ang(i);
    disp(['now angled  ',num2str(angle),' is being made']);
   
    tx(i).delay=calc_delay(trans,c,angle);
    xdc_focus_times(emit,0,tx(i).delay);
    xdc_focus_times(rcv,0,zeros(1,trans.numele));
%% set point target
    point_pos=[
               0 0 15;
               0 0 20;
               0 0 25;
               -10 0 20;
               10 0 20
                ]/1000;

    point_amp=20*ones(1,max(size(point_pos)));
    [v_temp,tstart(i)]=calc_scat_multi(emit,rcv,point_pos,point_amp');
    rf_data(1:max(size(v_temp)),:,i)=v_temp;
    disp(['angle ',num2str(angle),' has been done.'])
end

2 波束形成

for ii=1:userset.angnum
    rf_an=rf(:,:,ii);
    TXangle = ang(ii);             
    halfaper = sign(TXangle)*xT(end);      
    
%% Do Imaging
   
        % Delay
        dTX = z*cos(TXangle) + (x+ halfaper)*sin(TXangle); % TX distance
        dRX = sqrt((xT-x).^2 + z.^2); % RX distance
        tau = (dTX + dRX) / trans.c; % TX+RX travel time
        dasdata(:,ii) = das(rf_an,tau,trans);
    end
      disp(['Angle ',num2str(ang(ii)),' has been done']);
end

3. 11角度复合 

 4. 21角度复合