Dear Ben Cox & Bradley Treeby
I am working on harmonic ultrasound imaging. I am using k-Wave for my simulations, but I have two main questions about k-Wave,
The first is about the excitation signals. I want to know that can k-Wave do the simulations with other types of excitation signals (input signals) except the Gaussian pulse?
The second is about the cyst phantom. I does my simulations on wire targets and I have get good images of them. I want to evaluate my simulations for cyst phantom and I want to have a cyst phantom with low intensity inside (black) and homogeneous, and random scatterers outside the cyst phantom with high intensity than inside (bright), but in simulations I can’t see the cyst clearly.
My cyst phantom code is:
background_map_mean = 1;
background_map_std = 0.0008;
background_map = background_map_mean + background_map_std*randn([Nz, Nx]);
% define a random distribution of scatterers for the highly scattering
% region
scattering_map =randn([Nz, Nx]);
scattering_c0 = c0 + 25 + 75*scattering_map;
scattering_c0(scattering_c0 > 1600) = 1600;
scattering_c0(scattering_c0 < 1400) = 1400;
scattering_rho0 = scattering_c0/1.5;
% define properties
sound_speed_map = c0*ones(Nz, Nx).*background_map;
density_map = rho0*ones(Nz,Nx).*background_map;
% define a sphere for a highly scattering region
scattering_region1 =sqrt(xx.^2+(zz-40e-3).^2)<5e-3;
scattering_region1=ones(Nz,Nx)-scattering_region1;
% assign region
sound_speed_map(scattering_region1 == 1) = scattering_c0(scattering_region1 == 1);
density_map(scattering_region1 == 1) = scattering_rho0(scattering_region1 == 1);
% assign to the medium inputs
medium.sound_speed = sound_speed_map;
medium.density = density_map;
figure
imagesc(medium.density)
colormap(gray(128))
Thanks in advanced for your help
Best regards
Sepand