http://www.k-wave.org/forum/topic/new-grid-for-reconstruction Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 13:33:59 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/new-grid-for-reconstruction#post-5854 Mon, 20 Mar 2017 21:07:08 +0000 Bradley Treeby 5854@http://www.k-wave.org/forum/ <p>Hi nfreundl,</p> <p>You can find the complete m-file for these examples from within the MATLAB help browser (or in the examples folder in the toolbox). Only parts of the m-file are included in the html description, so it's likely that you are missing some lines of code.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/new-grid-for-reconstruction#post-5853 Sun, 19 Mar 2017 06:44:02 +0000 nfreundl 5853@http://www.k-wave.org/forum/ <p>No I'm using the example of k-Wave documentation website<br /> <a href="http://www.k-wave.org/documentation/k-wave.php" rel="nofollow">http://www.k-wave.org/documentation/k-wave.php</a> </p> http://www.k-wave.org/forum/topic/new-grid-for-reconstruction#post-5851 Mon, 13 Mar 2017 12:13:04 +0000 Bradley Treeby 5851@http://www.k-wave.org/forum/ <p>Hi nfreundl,</p> <p>Are you using the example m-file that is part of k-Wave (you can open this from the MATLAB help browser by clicking on "open the file")?</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/new-grid-for-reconstruction#post-5842 Sat, 25 Feb 2017 21:09:38 +0000 nfreundl 5842@http://www.k-wave.org/forum/ <p>Hi,<br /> I'm new in K-Wave, I've been trying the initial value example with the 2-D with the sensor disc. The simulation worked well. I want to make the reconstruction of the the initial pressure distribution.<br /> After setting manually the time array, I add those lines</p> <p>sensor.time_reversal_boundary_data=sensor_data;</p> <p>source.p0=0<br /> %run the reverse simulation<br /> p0_recon = kspaceFirstOrder2D(kgrid, medium, source, sensor);</p> <p>So far, it works. They suggested to add noise and change the computational grid:</p> <p>noise_ratio=2.0;</p> <p>for i=1:num_sensor_points<br /> sensor_data(i,:)=addNoise(sensor_data(i,:),noise_ratio);<br /> end</p> <p>It's still okay but when it goes to the grid change, it doesn't work. The interpolation works, but not the reverse simulation : I used<br /> Nx_recon=100;<br /> Ny_recon=100;<br /> dx_recon=dx*Nx/Nx_recon;<br /> dy_recon=dy*Ny/Ny_recon;</p> <p>kgrid_recon=makeGrid(Nx_recon,dx_recon,Ny_recon,dy_recon);</p> <p>%interpolation of the different k-grids<br /> sensor_radius_grid_points = round(sensor_radius/kgrid_recon.dx);<br /> binary_sensor_mask = makeCircle(kgrid_recon.Nx, kgrid_recon.Ny, kgrid_recon.Nx/2, kgrid_recon.Ny/2, sensor_radius_grid_points);<br /> source.p0=0;<br /> sensor.time_reversal_boundary_data = interpCartData(kgrid_recon, sensor_data, sensor.mask, binary_sensor_mask);</p> <p>p0_recon = kspaceFirstOrder2D(kgrid_recon, medium, source, sensor);</p> <p>and I got these errors</p> <p>Error in kspaceFirstOrder_inputChecking (line 628)<br /> sensor.time_reversal_boundary_data(:, new_col_pos) = order_index;</p> <p>Error in kspaceFirstOrder2D (line 496)<br /> kspaceFirstOrder_inputChecking;</p> <p>Error in changeGridKwave (line 80)<br /> p0_recon = kspaceFirstOrder2D(kgrid_recon, medium, source, sensor);</p> <p>So if someone has a solution, thanks :) </p>