http://www.k-wave.org/forum/topic/kwavearray-distributed_source_signal-and-alternative-source-assignment-methods Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 01:09:09 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/kwavearray-distributed_source_signal-and-alternative-source-assignment-methods#post-8397 Tue, 21 Dec 2021 04:27:43 +0000 vat 8397@http://www.k-wave.org/forum/ <p>Hi Brad,</p> <p>For the example I indicate, is there a recommendation for routine or approach to most simply assign all the grid points on the disc element a particular velocity of given time varying amplitude in m/s? This amplitude is desired of the velocity component normal to the disc area, and all other velocity components would be 0 m/s. </p> <p>On your first question, there is a geometry difference but I do not understand why I still would not see the specified velocity on a monitoring node located at the array adding a modified function below in which changing the weights to 1 in the karray class:</p> <p>'function distributed_source_signal_unweighted = getDistributedSourceSignalUnweighted(obj, kgrid, source_signal)</p> <p> % check the array has elements<br /> obj.checkForElements(dbstack);</p> <p> % get the binary mask to count how many points contribute to<br /> % the source<br /> mask = obj.getArrayBinaryMask(kgrid);<br /> mask_ind = find(mask);<br /> num_source_points = sum(mask(:));</p> <p> % number of time points in the signal<br /> Nt = size(source_signal, 2);</p> <p> % initialise the source signal<br /> distributed_source_signal_unweighted = zeros(num_source_points, Nt);</p> <p> % loop through the elements<br /> for ind = 1:obj.number_elements</p> <p> % get the offgrid source weights<br /> source_weights = obj.getElementGridWeights(kgrid, ind);</p> <p> % get indices of the non-zero points<br /> element_mask_ind = find(source_weights ~= 0);</p> <p> % convert these to indices in the distributed source<br /> local_ind = ismember(mask_ind, element_mask_ind);</p> <p> % add to distributed source - HERE I ADD SO THERE IS NO<br /> % WEIGHTING<br /> distributed_source_signal_unweighted(local_ind, :) = ...<br /> distributed_source_signal_unweighted(local_ind, :) ...<br /> + (1)* source_signal(ind, :);</p> <p> end</p> <p> end' </p> <p>But if you can make a recommendation on my first point, I care about that more urgently.</p> <p>Thanks for the attention. </p> http://www.k-wave.org/forum/topic/kwavearray-distributed_source_signal-and-alternative-source-assignment-methods#post-8375 Thu, 16 Dec 2021 11:12:46 +0000 Bradley Treeby 8375@http://www.k-wave.org/forum/ <p>Hi Vatche,</p> <p>Is some of the difference related to your use of disc elements in <code>kWaveArray</code>, rather than rectangular ones (which are used in <code>kWaveTransducer</code>)?</p> <p>One point to note, if you're using velocity sources and the physical array elements are not aligned with the grid, you'll have to consider how to assign the vector components of the velocity so the resulting field is still normal to the element surface.</p> <p>Brad </p> http://www.k-wave.org/forum/topic/kwavearray-distributed_source_signal-and-alternative-source-assignment-methods#post-8327 Wed, 13 Oct 2021 19:20:37 +0000 vat 8327@http://www.k-wave.org/forum/ <p>Hello,</p> <p>I am trying to resolve some discrepancy/ambiguity about input signal assignment using kWaveArray objects to define my transducer/array geometry.</p> <p>Previously I use kWaveTransducer object passing input in m/s to be assigned to particle velocities as follows:</p> <p>% create the input signal using toneBurst<br /> input_signal = toneBurst(1/kgrid.dt, tone_burst_freq, tone_burst_cycles);</p> <p>% scale the source magnitude by the source strength in m/s<br /> input_signal = source_velocity.*input_signal;</p> <p>...</p> <p>transducer.input_signal = input_signal;</p> <p>% create the transducer using the defined settings<br /> transducer = kWaveTransducer(kgrid, transducer);</p> <p>Now a similar geometry defined kWaveArray has a different pressure output and I assume this relates partially to my use but perhaps to implementation? Here is essentially what I tried:</p> <p>karray = kWaveArray;</p> <p>% adding disc elements<br /> for col=1:N<br /> for row=1:N<br /> karray.addDiscElement(elem_pos, diameter, [100,0,0]);<br /> end<br /> end<br /> ...<br /> % assign binary mask<br /> source.p_mask = karray.getArrayBinaryMask(kgrid);<br /> ...<br /> source_sig = repmat(input_signal,count,1);%for assigning to number elements count </p> <p>% get distributed source signals (this automatically returns a weighted<br /> % source signal for each grid point that forms part of the source elements)<br /> source.p = karray.getDistributedSourceSignal(kgrid, source_sig);</p> <p>Based on documentation <a href="http://www.k-wave.org/documentation/example_us_defining_transducer.php" rel="nofollow">http://www.k-wave.org/documentation/example_us_defining_transducer.php</a><br /> I assume the particle velocities/grid points are assigned the input_signal in m/s on my attempt to use kWaveTransducer but here the weighting scheme reduces the values at each of the grid points?</p> <p>Is there a recommended approach that assures grid points for all elements of the kWaveArray are assigned specific particle velocity time traces I would like to pass to them in source_sig or some other structure? </p> <p>Thank you for the attention.</p> <p>Cheers,<br /> Vatche </p>