http://www.k-wave.org/forum/topic/sensor-aperture Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 00:43:47 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/sensor-aperture#post-6682 Sun, 09 Dec 2018 21:34:28 +0000 Bradley Treeby 6682@http://www.k-wave.org/forum/ <p>It doesn't work because the two masks don't have the same 'on' value (the ball is 5 and 0, the sensor mask is 1 and 0). Instead, try the following:</p> <p><code>voxelPlot(single(sensor.mask | p0));</code></p> <p>Brad </p> http://www.k-wave.org/forum/topic/sensor-aperture#post-6660 Sat, 24 Nov 2018 15:52:24 +0000 peiliu 6660@http://www.k-wave.org/forum/ <p>hi!<br /> I have used the k-wave for few weeks and have some problems.I want to do a simulation about the aperture of a single element transducer,so I define a transducer and a source use the 'makeBall' in the 3D.The code is shown blow:<br /> % create the computational grid<br /> Nx = 64; % number of grid points in the x direction<br /> Ny = 64; % number of grid points in the y direction<br /> Nz = 64; % number of grid points in the z direction<br /> dx = 100e-3/Nx; % grid point spacing in the x direction [m]<br /> dy = dx; % grid point spacing in the y direction [m]<br /> dz = dx; % grid point spacing in the z direction [m]<br /> kgrid = kWaveGrid(Nx, dx, Ny, dy, Nz, dz);</p> <p>% define the properties of the propagation medium<br /> medium.sound_speed = 1500; % [m/s]</p> <p>% create the time array<br /> kgrid.makeTime(medium.sound_speed);</p> <p>% create a concave sensor<br /> sphere_offset =10; % [grid points]<br /> diameter = Nx/2 + 1; % [grid points]<br /> bowl_pos = [1 + sphere_offset, Ny/2, Nz/2]; % [grid points]<br /> focus_pos = [Nx/2, Ny/2, Nz/2]; % [grid points]<br /> sensor.mask = makeBowl([Nx, Ny, Nz], bowl_pos, inf, diameter, focus_pos);</p> <p>% create initial pressure distribution using makeBall<br /> ball_magnitude = 5; % [Pa]<br /> ball_radius = 5; % [grid points]<br /> p0 = ball_magnitude * makeBall(Nx, Ny, Nz, 30, Ny/2, Nz/2, ball_radius);</p> <p>and I want to show the transducer and the source togather in the codes as follows:</p> <p>voxelPlot(sensor.mask+p0);<br /> view([130, 40]);</p> <p>my question is why the result is wrong.I execute the codes about the transducer and show it ,it's right,and I execute the codes about the source and show it,it's right as well.But when I try show them togather the transducer is can't be seen.I'm very confused,I need your help,thank you in advance! </p>