http://www.k-wave.org/forum/topic/understanding-beamforming Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 00:48:03 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/understanding-beamforming#post-4368 Fri, 14 Mar 2014 13:04:41 +0000 Bradley Treeby 4368@http://www.k-wave.org/forum/ <p>Hi Margot,</p> <p>If you are using the C++ code, the situation is slightly different. In this case, <code>sensor_data</code> will contain one time trace for every grid point that forms part of the active elements of your transducer (rather than one time trace per element). If your transducer elements contain more than one grid point (which they normally will), then further processing is required to put the C++ output data in the correct format before using the <code>scan_line</code> method. An additional method called <code>combine_sensor_data</code> has been added to the <code>kWaveTransducer</code> class ready for the next release. I have also copied a beta version of it below. You should add this to the "general class methods". After running the C++ simulation, you can then call <code>sensor_data = transducer.combine_sensor_data(sensor_data)</code>.</p> <p>Let me know if you have any problems,</p> <p>Brad.</p> <pre><code>% function to average the individual time series recorded at each % grid point of each transducer element (as returned by the C++ % code) and return a single time series per active transducer % element (as returned by the MATLAB code) function sensor_data_sum = combine_sensor_data(obj, sensor_data) % check the data is the correct size if size(sensor_data, 1) ~= (obj.number_active_elements * obj.element_width * obj.element_length) error(&#39;The number of time series in the input sensor_data must match the number of grid points in the active tranducer elements.&#39;); end % get index of which element each time series belongs to ind = obj.indexed_active_elements_mask(obj.indexed_active_elements_mask &gt; 0); % create empty output sensor_data_sum = zeros(obj.number_active_elements, size(sensor_data, 2)); % check if an elevation focus is set if isinf(obj.elevation_focus_distance) % loop over time series and sum for ii = 1:length(ind) sensor_data_sum(ind(ii), :) = sensor_data_sum(ind(ii), :) + sensor_data(ii, :); end else % get the elevation delay for each grid point of each % active transducer element (this is given in units of grid % points) dm = obj.delay_mask(2); dm = dm(obj.active_elements_mask ~= 0); % loop over time series, shift and sum for ii = 1:length(ind) sensor_data_sum(ind(ii), 1:end - dm(ii)) = sensor_data_sum(ind(ii), 1:end - dm(ii)) + sensor_data(ii, 1 + dm(ii):end); end end % divide by number of time series in each element sensor_data_sum = sensor_data_sum .* (1/ (obj.element_width * obj.element_length)); end</code></pre> http://www.k-wave.org/forum/topic/understanding-beamforming#post-4362 Fri, 07 Mar 2014 16:24:24 +0000 mgt_p 4362@http://www.k-wave.org/forum/ <p>I can't edit my previous post, but in fact, I am interested in understanding lines 693 to 701 in kWaveTransducer.m . It is just before the summation to get an A-line. </p> <p>To me, it really doesn't look like the whole data gets through apodization and delays.<br /> Is it because only the first line is linked to the element number? Or am I totally mistaken on sensor_data?</p> <p>Thanks ! </p> http://www.k-wave.org/forum/topic/understanding-beamforming#post-4361 Fri, 07 Mar 2014 14:55:51 +0000 mgt_p 4361@http://www.k-wave.org/forum/ <p>Hello Bradley,</p> <p>Thanks for your answer.</p> <p>Yes, I am using the C++ code (so much faster !). </p> http://www.k-wave.org/forum/topic/understanding-beamforming#post-4343 Tue, 25 Feb 2014 14:32:42 +0000 Bradley Treeby 4343@http://www.k-wave.org/forum/ <p>Hi Margot,</p> <p>The transducer object that is created by <code>makeTransducer</code> has two fields related to the number of elements. The first is <code>number_elements</code> which sets the physical number of elements in the transducer. The second is <code>active_elements</code> which defines which of the elements are actually used to generate or detect sound. When using a transducer as a sensor, the number of time series in <code>sensor_data</code> should be <code>active_elements</code>. Thus, the processing only applies to these time series.</p> <p>Note, if you're using the C++ code, this doesn't return a time trace per transducer element, but rather, a time trace per grid point. In this case, some additional processing is required. Let me know if this is the case and I'll give you the additional code required.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/understanding-beamforming#post-4332 Thu, 20 Feb 2014 14:21:22 +0000 mgt_p 4332@http://www.k-wave.org/forum/ <p>Hi,</p> <p>I don't quite understand how delays are applied to sensor_data in the kWaveTransducer.m.</p> <p>I understand the calculation, but it seems to me that only the first number_active_elements lines of it is modified. So a wide majority of values in sensor_data isn't modified...</p> <p>I suspect I don't understand exactly the meaning of sensor_data, but could somebody please take a minute to explain this to me?</p> <p>Thanks a lot,</p> <p>Margot </p>