http://www.k-wave.org/forum/topic/delays-for-recorded-sensor_data-in-synthetic-aperture-ultrasound Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 00:22:28 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/delays-for-recorded-sensor_data-in-synthetic-aperture-ultrasound#post-5213 Fri, 21 Aug 2015 12:05:55 +0000 Bradley Treeby 5213@http://www.k-wave.org/forum/ <p>Hi sepand,</p> <p>I would suggest using an additive source unless you specifically want to enforce a boundary condition (for example, to perform time reversal using experimental data). Using an additive source in k-Wave corresponds to the normal way mass and force sources are defined in the wave equation.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/delays-for-recorded-sensor_data-in-synthetic-aperture-ultrasound#post-5189 Sat, 01 Aug 2015 07:07:30 +0000 sepand 5189@http://www.k-wave.org/forum/ <p>Hi Brad</p> <p>Thank you so much for your help. By reducing the grid spacing and the source frequency most of unknown fluctuations has been removed, But for removing all of them I have to change the boundary condition from dirichlet to additive(the default). Now it is a question for me that, is it true to use the additive boundary condition for ultrasound simulations? </p> <p>Thank you in advanced for your help<br /> Sepand </p> http://www.k-wave.org/forum/topic/delays-for-recorded-sensor_data-in-synthetic-aperture-ultrasound#post-5174 Tue, 28 Jul 2015 12:32:44 +0000 Bradley Treeby 5174@http://www.k-wave.org/forum/ <p>Hi Sepand,</p> <p>Your source center frequency is very close to the maximum frequency supported by the grid. I would suggest reducing the grid spacing (dx/dy) or reducing the source frequency. For any simulation, it's always a good idea to run a convergence test by reducing dx and CFL and checking you get the same answer.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/delays-for-recorded-sensor_data-in-synthetic-aperture-ultrasound#post-5167 Thu, 23 Jul 2015 10:26:22 +0000 sepand 5167@http://www.k-wave.org/forum/ <p>Dear Brad<br /> in addition to unknown fluctuations the mentioned echo signals are on an ascending curve around zero axis(y=0). I mean that the echo signals should be zero for other parts of homogeneous medium and they should have amplitudes for phantom position. But my received echoes are not on zero axis(y=0). </p> http://www.k-wave.org/forum/topic/delays-for-recorded-sensor_data-in-synthetic-aperture-ultrasound#post-5166 Thu, 23 Jul 2015 09:02:06 +0000 sepand 5166@http://www.k-wave.org/forum/ <p>Hi Brad<br /> My codes are as following:</p> <p>clear all<br /> close all<br /> clc<br /> data_cast = 'single';<br /> run_simulation = true;<br /> pml_x_size = 20; % [grid points]<br /> pml_y_size = 20; % [grid points]</p> <p>% set total number of grid points not including the PML<br /> sc = 1;<br /> Nx = 256/sc - 2*pml_x_size; % [grid points]<br /> Ny = 256/sc - 2*pml_y_size; % [grid points]</p> <p>% set desired grid size in the x-direction not including the PML<br /> x = 50e-3; % [m]</p> <p>% calculate the spacing between the grid points<br /> dx = x/Nx; % [m]<br /> dy = dx; % [m]</p> <p>kgrid = makeGrid(Nx, dx, Ny, dy);</p> <p>% define the properties of the propagation medium</p> <p>medium.sound_speed= 1540; % [m/s]<br /> medium.density = 1000; % [kg/m^3]<br /> medium.alpha_coeff = 0.75; % [dB/(MHz^y cm)]<br /> medium.alpha_power = 1.5;<br /> medium.BonA = 6;</p> <p>% create time array</p> <p>t_end = (Nx*dx)*2.2/medium.sound_speed; % [s]<br /> [kgrid.t_array, dt]= makeTime(kgrid, medium.sound_speed, [], t_end);</p> <p>point_magnitude =100*randn(Nx,Ny); % [Pa]<br /> x_pos = 3*Nx/4; % [grid points]<br /> y_pos = Ny/2; % [grid points]<br /> radius = 1; % [grid points]<br /> arc_angle = 2*pi; % [radians]</p> <p>Disc = point_magnitude.*makeDisc(Nx, Ny, x_pos, y_pos, radius, arc_angle);<br /> imagesc(double(Disc));<br /> %colormap(gray)<br /> density_map = 1000*ones(Nx, Ny);<br /> medium.density = density_map (:,:);<br /> medium.density = medium.density + Disc;<br /> imagesc(double(medium.density));<br /> colormap(gray)</p> <p>sound_speed_map = 1540*ones(Nx, Ny);<br /> medium.sound_speed = sound_speed_map (:,:);<br /> %medium.sound_speed = medium.sound_speed + Disc;<br /> imagesc(double(medium.sound_speed));<br /> colormap(gray)</p> <p>number_elements=64;<br /> element_width=1;<br /> element_spacing=0;<br /> transducer_width = number_elements*element_width ...<br /> + (number_elements - 1)*element_spacing;</p> <p>sensor.mask = zeros(Nx, Ny);<br /> sensor.mask(1, Ny/2-(number_elements/2-1):Ny/2+number_elements/2) = 1;<br /> %sensor.record_start_index=800;<br /> imagesc(double(sensor.mask));</p> <p>sensor_element_width=1;</p> <p>%%================================================<br /> %--synthetic aperture simulation------------------<br /> %=================================================</p> <p>%matlabpool('open',4); % start of Parallel Computing<br /> j=1;<br /> for i=Ny/2-(number_elements/2-1):Ny/2+number_elements/2;<br /> j<br /> source.p_mask = zeros(Nx, Ny);<br /> source.p_mask(1, i) = 1;<br /> imagesc(source.p_mask);<br /> colormap(gray);</p> <p>tone_burst_cycles = 1;<br /> tone_burst_freq=3e6; </p> <p>% % create the input signal using toneBurst<br /> input_signal = toneBurst(1/dt, tone_burst_freq, tone_burst_cycles);%,'Plot', true);</p> <p>source_strength=1e6;<br /> source.p =(source_strength./(1000*1540)).*input_signal;<br /> source.p_mode='dirichlet';</p> <p>input_args = {...<br /> 'PMLInside', false, 'PMLSize', [pml_x_size, pml_y_size], ...<br /> 'DataCast', data_cast, 'DataRecast', true, 'PlotSim', false};</p> <p>sensor_data{j,1} = kspaceFirstOrder2D(kgrid, medium, source, sensor,input_args{:});</p> <p>j=j+1;<br /> end</p> <p> thank you in advanced for your help<br /> Sepand </p> http://www.k-wave.org/forum/topic/delays-for-recorded-sensor_data-in-synthetic-aperture-ultrasound#post-5159 Sun, 19 Jul 2015 13:06:44 +0000 Bradley Treeby 5159@http://www.k-wave.org/forum/ <p>Hi Sepand,</p> <p>Without seeing your code it's hard to say. Have you tried running a convergence test, i.e., reducing dx/dy and the CFL number and seeing if you get the same answer?</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/delays-for-recorded-sensor_data-in-synthetic-aperture-ultrasound#post-5144 Wed, 15 Jul 2015 09:12:37 +0000 sepand 5144@http://www.k-wave.org/forum/ <p>Hi Brad</p> <p>Thank you very much for your help. As you suggested I have eliminated the beginning of received signal and The echo signal could be seen in remained signal, but the big problem is some unknown oscillations in addition to echos. these oscillations have variable amplitudes and they sounds to have a constant frequency. I can not really understand the origin of them. Because,as you know, for a single scatterer in a homogeneous medium, we should have a echo signal similar to input signal with a specific time delay according to the distance from emitter. but I have gotten a long signal with unknown fluctuations. How can I get rid of these additional oscillations? </p> <p>In my 2D simulations the input signal is a 3MHz tune_burst with one cycle.</p> <p>thank you in advanced</p> <p>Sepand </p> http://www.k-wave.org/forum/topic/delays-for-recorded-sensor_data-in-synthetic-aperture-ultrasound#post-5124 Wed, 08 Jul 2015 11:46:48 +0000 Bradley Treeby 5124@http://www.k-wave.org/forum/ <p>Hi Sepand,</p> <p>In k-Wave, sources are implemented as the injection of force or mass into your medium. This means there is no physical interface at the transducer surface. Similarly, the sensor mask just tells the simulation where to sample the pressure field at each time step (so there is no physical sensor in the field). Consequently, what you are seeing in your simulations is the input signal propagating directly to the sensor positions. As you would expect, this arrives at the closest element first. To avoid this, you need to time-gate out the beginning of the received signal. This should allow you to see the reflections from the scatterer. Indeed, the same is true in a real ultrasound machine, as the elements cannot both transmit and receive at the same time.</p> <p>Hope that helps,</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/delays-for-recorded-sensor_data-in-synthetic-aperture-ultrasound#post-5114 Tue, 30 Jun 2015 14:42:06 +0000 sepand 5114@http://www.k-wave.org/forum/ <p>Hi</p> <p>Dear all , Bencox and Bradley Treeby</p> <p>I am using source and sensor class in 2D for ultrasound synthetic aperture imaging. I have defined 32 grid points as transducer elements with no spacing between them. For each emission, just one of the source points is active but all 32 sensor elements receive the echoes. Source and sensor are at the same place and with same geometry. The scattering medium is a homogeneous medium with a point target at the middle of medium (a single pixel with a little bit more density than background). As I know from ultrasound, this point target should reflect the echo signals to the sensors (receiving elements). But after plotting the sensor_data, the echo signals arranged for each transmission, so that, the nearest element to the emitter has the smallest delay. While, as you know, for synthetic aperture imaging (for a point target at the middle of medium) it should have the smallest delay for middle element for all transmissions regardless of location of emitter. So, I think that, for each transmission, all sensor points just sense the transmitted signal not echo signals from point target. How can I get to the true echo signals from point target (PSF)?<br /> Best regards<br /> Sepand </p>