http://www.k-wave.org/forum/topic/time-reversal-image-reconstruction-with-kwavearray-elements Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 21:25: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/time-reversal-image-reconstruction-with-kwavearray-elements#post-9243 Tue, 27 Jan 2026 11:00:38 +0000 Oskar Sadowski 9243@http://www.k-wave.org/forum/ <p>Hello,<br /> I'm still somewhat new to kWave and am trying to work out Time Reversal image reconstruction for a segmented ring array, created using kWaveArray (i.e. <code>addArcElement</code> and <code>addCustomElement</code>). During testing around, I noticed a unexpected behaviour, and would be grateful for some insights on the matter.</p> <p>I've put together a simple minimal example issued below to explain my problem. The geometry is simply a couple arc shaped elements equally distributed in a circular pattern, all with the same parameters and focus on (0,0). Identical signals (here just <code>toneBurst</code>, the exact signal doesn't matter) are assigned to all elements using <code>sensor.time_reversal_boundary_data</code></p> <p>I would expect each elements to send out the same signal, as they are identical and the geometry is symmetrical. However, the simulation process shows that the signal intensity differs quite heavily between different elements, with (as far as I can tell) the ones closer to the cartesian axes having a higher intensity assigned than the ones off-axis. The number of elements does not affect this behavior</p> <p>While I'm sure there is a proper explanation for this, that I'm just missing, I can not make sense of this behavior currently, and would be grateful for any advise on how to get the proper results. I should note that a identical result can be achieved when constructing the sensor as a source and setting <code>source.p_mode = &#39;dirichlet&#39;</code></p> <p>Minimal example:</p> <p>clear; clearvars;</p> <p>% =========================================================================<br /> % PARAMETERS<br /> % =========================================================================</p> <p>numElements = 7;</p> <p>radius = 1E-2;<br /> sizeArc = 1E-3;<br /> focus = [0,0];</p> <p>gridSize = [2.5E-2, 2.5E-2];<br /> gridN = [512, 512];<br /> gridStep = [gridSize(1)/gridN(1), gridSize(2)/gridN(2)];</p> <p>% =========================================================================<br /> % ELEMENT DEFINITION<br /> % =========================================================================</p> <p>karray = kWaveArray;</p> <p>for i=1:numElements</p> <p> pos = [radius*cosd(i*360/numElements), ...<br /> radius*sind(i*360/numElements)];<br /> karray.addArcElement(pos, radius, sizeArc, focus);</p> <p>end</p> <p>% =========================================================================<br /> % GRID &amp; MEDIUM PROPERTIES<br /> % =========================================================================</p> <p>medium.sound_speed = 1520;<br /> medium.density = 1040;</p> <p>kgrid = kWaveGrid(gridN(1), gridStep(1), gridN(2), gridStep(2));<br /> kgrid.makeTime(medium.sound_speed)</p> <p>% =========================================================================<br /> % ASSIGN SOURCE &amp; SENSOR<br /> % =========================================================================</p> <p>signal = repmat(toneBurst(1/kgrid.dt, 100E4, 24).*15, numElements, 1);</p> <p>source.p0 = 0;</p> <p>sensor.mask = karray.getArrayBinaryMask(kgrid);<br /> sensor.time_reversal_boundary_data = karray.getDistributedSourceSignal(kgrid, signal);</p> <p>% =========================================================================<br /> % RUN SIMULATION<br /> % =========================================================================</p> <p>simulation = kspaceFirstOrder2D(kgrid, medium, source, sensor, ...<br /> 'PMLInside', false, 'DataCast', 'gpuArray-single');</p> <p>figure;<br /> imagesc(kgrid.y_vec, kgrid.x_vec, simulation)<br /> axis image; colormap("gray")</p> <p>hold on;<br /> karray.plotArray(false); </p>