linear.m

超声波扫查中

f0=3e6;
fs=50e6;
c=1540;
width=1/1000;
element_height=5/1000;
kerf=0.1/1000;
focus=[0 0 40]/1000;
N_element=64;
set_field('fs',fs);
emit_aperture=xdc_linear_array(N_element,width,element_height,kerf,1,1,focus);
impulse_response=cos(2*pi*f0*(0:1/fs:2/f0));
impulse_response=impulse_response.*hann(max(size(impulse_response)))';
xdc_impulse(emit_aperture,impulse_response);
excitation=cos(2*pi*f0*(0:1/fs:2/f0));
xdc_excitation(emit_aperture,excitation);
%---------------------
receive_aperture=xdc_linear_array(N_element,width,element_height,kerf,1,1,focus);
xdc_impulse(receive_aperture,impulse_response);

% [phantom_positions,phantom_amplitudes]=cyst_phantom(100);
%------------------------
phantom_positions=[0 0 30;0 0 40;0 0 50]/1000;
phantom_amplitudes=[1;1;1];
no_lines=100;
sector=10/1000;
d_x=sector/no_lines;
z_focus=40/1000;

image_data=zeros(800,no_lines);
x=-sector/2;
for i=1:no_lines;
    xdc_center_focus(emit_aperture,[x 0 0]);
    xdc_focus(emit_aperture,0,[x,0,z_focus]);
    xdc_center_focus(receive_aperture,[x 0 0]);
    xdc_focus(receive_aperture,0,[x,0,z_focus]);
    [v,t1]=calc_scat(emit_aperture,receive_aperture,phantom_positions,phantom_amplitudes);
    image_data(1:max(size(v)),i)=v;
    times(i)=t1;
    x=x+d_x;
end

xdc_free(emit_aperture);
xdc_free(receive_aperture);
min_sample=min(times)*fs;

for i=1:no_lines;
    rf_env=abs(hilbert([zeros(times(i)*fs-min_sample,1);image_data(:,i)]));
    env(1:max(size(rf_env)),i)=rf_env;
end
env=env/max(max(env));
env=log(env+0.1);
env=env-min(min(env));
env=64*env/max(max(env));
figure;
image(env);
colormap(gray(64));


% for i=1:no_lines;
%     rf_env=abs(hilbert([zeros(times(i)*fs-min_sample,1);image_data(:,i)]));
%     env(1:max(size(rf_env)),i)=rf_env;
% end
% D=10;
% env1=env(1:D:max(size(env)),:)/max(max(env));
% env1=env1/max(max(env1));
% env1=log(env1+0.1);
% env1=env1-min(min(env1));
% env1=64*env1/max(max(env1));
% figure;
% image(env1);
% colormap(gray(64));    


% D=10;
% image_data1=image_data(1:D:max(size(image_data)),:);
% 
% for i=1:no_lines;
%     rf_env=abs(hilbert([zeros(times(i)*fs-min_sample,1);image_data1(:,i)]));
%     env(1:max(size(rf_env)),i)=rf_env;
% end
% env1=env1/max(max(env1));
% env1=log(env1+0.1);
% env1=env1-min(min(env1));
% env1=64*env1/max(max(env1));
% figure;
% image(env1);
% colormap(gray(64));