http://www.k-wave.org/forum/topic/save-snapshots-of-the-simulation-every-x-timesteps Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 23:40: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/save-snapshots-of-the-simulation-every-x-timesteps#post-8578 Wed, 13 Jul 2022 18:45:35 +0000 bencox 8578@http://www.k-wave.org/forum/ <p>Thanks guys, noted!</p> <p>Best wishes<br /> Ben </p> http://www.k-wave.org/forum/topic/save-snapshots-of-the-simulation-every-x-timesteps#post-8564 Sat, 09 Jul 2022 23:10:17 +0000 JBAjo-Franklin 8564@http://www.k-wave.org/forum/ <p>Up-voting this! I'm trying to do the same thing and a space-filling sensor array is not what I'm looking for. I'm hoping to capture a quantitative snapshot of the wavefield at one point in time rather than a movie or a sensor trace. Also - great code! </p> http://www.k-wave.org/forum/topic/save-snapshots-of-the-simulation-every-x-timesteps#post-8536 Wed, 18 May 2022 19:35:10 +0000 3dfernando 8536@http://www.k-wave.org/forum/ <p>First, thanks for this code. It's really nice to have a toolbox like k-Wave!</p> <p>I was wondering if there's a way of having kWave save a full snapshot of the acoustic field at predefined times (say, you run for t=linspace(0,1e-3,1000); and you want snapshots for tSnapshot=linspace(0,1e-3,10)). These snapshots would be full pictures of the pressure field, as if there were sensors everywhere.</p> <p>An example with the field filled up with sensors would be something like shown below; but as the code is below it does not skip any time step. This is not great for memory purposes, if the simulation is very long. Right now I'm making do with this, but I wonder if k-wave has something better!</p> <p>%================<br /> Nx = 128; % number of grid points in the x (row) direction<br /> Ny = 128; % number of grid points in the y (column) direction<br /> dx = 1e-3; % grid point spacing in the x direction [m]<br /> dy = 1e-3; % grid point spacing in the y direction [m]</p> <p>t=linspace(0,1e-3,1000);</p> <p>kgrid = kWaveGrid(Nx, dx, Ny, dy);<br /> kgrid.t_array=t;</p> <p>% define the properties of the propagation medium<br /> medium.sound_speed = 343; % [m/s]<br /> medium.alpha_coeff = 0.75; % [dB/(MHz^y cm)]<br /> medium.alpha_power = 1.5;<br /> medium.density = 1.2 * ones(Nx, Ny);</p> <p>% define a single source point<br /> source.p_mask = zeros(Nx, Ny);<br /> source.p_mask(end - Nx/4, Ny/2) = 1;</p> <p>% define a time varying sinusoidal source<br /> source_freq = 10e3; % [Hz]<br /> source_mag = 10; % [Pa]<br /> source.p = source_mag * sin(2 * pi * source_freq * t);</p> <p>sensor.mask=ones(Nx, Ny);</p> <p>sensor_data = kspaceFirstOrder2D(kgrid, medium, source, sensor);<br /> sensor_data=reshape(sensor_data,size(sensor.mask,1),size(sensor.mask,2),length(t));</p> <p>%Plots the simulation<br /> figure;<br /> for i=1:length(t)<br /> imagesc(sensor_data(:,:,i), [-1, 1]);<br /> colormap(getColorMap);<br /> drawnow;pause(0.01)<br /> end</p> <p>%======= </p>