k-Wave

1. vat
   

   Member
   Joined: Oct '21
   Posts: 2

   Hello,

   I am trying to resolve some discrepancy/ambiguity about input signal assignment using kWaveArray objects to define my transducer/array geometry.

   Previously I use kWaveTransducer object passing input in m/s to be assigned to particle velocities as follows:

   % create the input signal using toneBurst
   input_signal = toneBurst(1/kgrid.dt, tone_burst_freq, tone_burst_cycles);

   % scale the source magnitude by the source strength in m/s
   input_signal = source_velocity.*input_signal;

   ...

   transducer.input_signal = input_signal;

   % create the transducer using the defined settings
   transducer = kWaveTransducer(kgrid, transducer);

   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:

   karray = kWaveArray;

   % adding disc elements
   for col=1:N
   for row=1:N
   karray.addDiscElement(elem_pos, diameter, [100,0,0]);
   end
   end
   ...
   % assign binary mask
   source.p_mask = karray.getArrayBinaryMask(kgrid);
   ...
   source_sig = repmat(input_signal,count,1);%for assigning to number elements count

   % get distributed source signals (this automatically returns a weighted
   % source signal for each grid point that forms part of the source elements)
   source.p = karray.getDistributedSourceSignal(kgrid, source_sig);

   Based on documentation http://www.k-wave.org/documentation/example_us_defining_transducer.php
   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?

   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?

   Thank you for the attention.

   Cheers,
   Vatche

   Posted 3 years ago #

2. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,643

   Hi Vatche,

   Is some of the difference related to your use of disc elements in kWaveArray, rather than rectangular ones (which are used in kWaveTransducer)?

   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.

   Brad

   Posted 2 years ago #

3. vat
   

   Member
   Joined: Oct '21
   Posts: 2

   Hi Brad,

   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.

   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:

   'function distributed_source_signal_unweighted = getDistributedSourceSignalUnweighted(obj, kgrid, source_signal)

   % check the array has elements
   obj.checkForElements(dbstack);

   % get the binary mask to count how many points contribute to
   % the source
   mask = obj.getArrayBinaryMask(kgrid);
   mask_ind = find(mask);
   num_source_points = sum(mask(:));

   % number of time points in the signal
   Nt = size(source_signal, 2);

   % initialise the source signal
   distributed_source_signal_unweighted = zeros(num_source_points, Nt);

   % loop through the elements
   for ind = 1:obj.number_elements

   % get the offgrid source weights
   source_weights = obj.getElementGridWeights(kgrid, ind);

   % get indices of the non-zero points
   element_mask_ind = find(source_weights ~= 0);

   % convert these to indices in the distributed source
   local_ind = ismember(mask_ind, element_mask_ind);

   % add to distributed source - HERE I ADD SO THERE IS NO
   % WEIGHTING
   distributed_source_signal_unweighted(local_ind, :) = ...
   distributed_source_signal_unweighted(local_ind, :) ...
   + (1)* source_signal(ind, :);

   end

   end'

   But if you can make a recommendation on my first point, I care about that more urgently.

   Thanks for the attention.

   Posted 2 years ago #

RSS feed for this topic

Reply

You must log in to post.