Dear all,
I am going to do the 2D photoacoustic simulation about a rectangle homogeneous phantom.
But I found the channel data of the simulation is weird, like two lines, so as the beamforming image.
Here is my MATLAB code:
% Transducer Setting
f0 = 11e6; % Center frequency [Hz]
fs = 4*f0; % Sampling frequency [Hz]
soundv = 1540;
lambda = soundv/f0;
BW = 60;
% 2D Simulation Setting
Nx = 128;
dx = 0.12e-3; %[m]
F_number = 1;
depth = F_number*Nx*dx;
dz = lambda/3;
Nz = 2*round(depth/dz);
kgrid = kWaveGrid(Nx,dx,Nz,dz);
% Medium
medium.sound_speed = soundv;
% Source
source.p0 = zeros(Nx,Nz);
hypo = zeros(Nx,Nz);
hypo(:, (Nz/2-50 : Nz/2+50)) = 1e6; % Rectangle phantom
source.p0 = hypo;
% Sensor
sensor.mask = zeros(Nx,Nz);
sensor.mask(:,1) = 1;
senor.frequency_response = [f0,BW];
% 2D Simulation
sensor_data_2D = kspaceFirstOrder2D(kgrid,medium,source,sensor);Is there any stupid mistake I miss?
Thanks.