http://www.k-wave.org/forum/topic/loading-an-external-image-for-reconstruction Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 09:02:42 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/loading-an-external-image-for-reconstruction#post-6914 Tue, 25 Jun 2019 08:48:38 +0000 elahehamini 6914@http://www.k-wave.org/forum/ <p>I found out the problem and it was the name of my m-file which was c and the code tried to execute it as a function and that's why I got an error on line <code>p = c.^2 .* (rhox + rhoy)</code> . Thank you so much for your time </p> http://www.k-wave.org/forum/topic/loading-an-external-image-for-reconstruction#post-6913 Tue, 25 Jun 2019 08:43:27 +0000 Bradley Treeby 6913@http://www.k-wave.org/forum/ <p>If I set</p> <p><code>p0_image = rand(128, 128);</code></p> <p>the simulation runs without errors. Are the values in F real numbers in single/double precision? </p> http://www.k-wave.org/forum/topic/loading-an-external-image-for-reconstruction#post-6912 Tue, 25 Jun 2019 08:10:36 +0000 elahehamini 6912@http://www.k-wave.org/forum/ <p>yes, my image is the same size as the grid and I use it to define the source. the Image is 128*128 pixels. This is the code: </p> <p> <code>p0_magnitude = 1;</code><br /> <code>load d1.mat</code><br /> <code>F =d1;</code><br /> <code>p0_image =p0_magnitude* F;</code></p> <p> % assign the grid size and create the computational grid<br /> <code>PML_size = 20; % size of the PML in grid points</code><br /> <code>Nx = 168 - 2 * PML_size; % number of grid points in the x direction</code><br /> <code>Ny = 168 - 2 * PML_size; % number of grid points in the y direction</code><br /> <code>x = 2.56e-3; % total grid size [m]</code><br /> <code>y = 2.56e-3; % total grid size [m]</code><br /> <code>dx = x / Nx; % grid point spacing in the x direction [m]</code><br /> <code>dy = y / Ny; % grid point spacing in the y direction [m]</code><br /> <code>kgrid = kWaveGrid(Nx, dx, Ny, dy);</code></p> <p> % assign to the source structure<br /> <code>source.p0 = p0_image;</code></p> <p> % define the properties of the propagation medium<br /> <code>medium.sound_speed = 1500; % [m/s]</code></p> <p> % define a centered Cartesian circular sensor<br /> <code>sensor_radius = 1.26e-3; % [m]</code><br /> <code>sensor_angle = pi; % [rad]</code><br /> <code>sensor_pos = [0, 0]; % [m]</code><br /> <code>num_sensor_points = 70;</code><br /> <code>cart_sensor_mask = makeCartCircle(sensor_radius, num_sensor_points, sensor_pos, sensor_angle);</code></p> <p> % assign to sensor structure<br /> <code>sensor.mask = cart_sensor_mask;</code></p> <p> % create the time array<br /> <code>kgrid.makeTime(medium.sound_speed);</code></p> <p> % set the input options<br /> <code>input_args = {&#39;Smooth&#39;, true, &#39;PMLInside&#39;, false, &#39;PlotPML&#39;, false, &#39;PlotSim&#39;, true, &#39;DataCast&#39;, &#39;single&#39;};</code></p> <p> % run the simulation<br /> <code>sensor_data = kspaceFirstOrder2D(kgrid, medium, source, sensor, input_args{:});</code> </p> http://www.k-wave.org/forum/topic/loading-an-external-image-for-reconstruction#post-6906 Tue, 25 Jun 2019 05:50:11 +0000 Bradley Treeby 6906@http://www.k-wave.org/forum/ <p>Hi elahehamini,</p> <p>It's hard to say without seeing more details. Is your image the same size as the grid (Nx, Ny)? Are you using it to define the source or the medium properties?</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/loading-an-external-image-for-reconstruction#post-6905 Mon, 24 Jun 2019 16:55:31 +0000 elahehamini 6905@http://www.k-wave.org/forum/ <p>I have simulated an image using Monte Carlo simulation and saved it as a 2D array and when I load it as an external image to have it reconstructed I get an error in the acoustical forward simulation using kspaceFirstOrder2D which refers it to the line p = c.^2 .* (rhox + rhoy) in the kspaceFirstOrder2D code. what am I doing wrong? </p>