http://www.k-wave.org/forum/topic/question-about-the-transducercombine_sensor_data-method Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 00:20:41 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/question-about-the-transducercombine_sensor_data-method#post-4770 Thu, 16 Oct 2014 11:33:57 +0000 Bradley Treeby 4770@http://www.k-wave.org/forum/ <p>Hi Shawn,</p> <p>Thanks for your post. After some digging, it seems the problem is that transmit apodization is not currently implemented in the C++ code. If you set <code>transducer.transmit_apodization = &#39;Rectangular&#39;</code>, hopefully everything should match up (please let me know if it doesn't, or if you already have this set). I've added this to the list to implement for the next release. </p> <p>Brad. </p> http://www.k-wave.org/forum/topic/question-about-the-transducercombine_sensor_data-method#post-4769 Thu, 16 Oct 2014 09:51:33 +0000 dsm194 4769@http://www.k-wave.org/forum/ <p>Hello,</p> <p>I don't understand the new "transducer.combine_sensor_data" method. I modified "example_us_bmode_linear_transducer" to additionally generate sensor data using kspaceFirstOrder3DC and to combine the sensor data from kspaceFirstOrder3DC using the above mentioned transducer.combine_sensor_data method. However, the combined sensor data does not appear to match up with the sensor data generated by kspaceFirstOrder3D. For example, when scan_line_index = 60, mean(abs(combined_sensor_data_3DC(:))) = 10504 whereas mean(abs(sensor_data_3D(:))) = 5270.4. </p> <p>Am I missing a step when processing the sensor_data output by kspaceFirstOrder3DC?</p> <p>here are the modifications I made to the example:</p> <pre><code>% preallocate the storage scan_lines_3D = zeros(number_scan_lines, kgrid.Nt); scan_lines_3DC = zeros(number_scan_lines, kgrid.Nt); % set the input settings input_args = {... &#39;PMLInside&#39;, false, &#39;PMLSize&#39;, [pml_x_size, pml_y_size, pml_z_size], ... &#39;DataCast&#39;, DATA_CAST, &#39;DataRecast&#39;, true, &#39;PlotSim&#39;, false}; % run the simulation if set to true, otherwise, load previous results from disk if RUN_SIMULATION % set medium position medium_position = 1; for scan_line_index = 1:number_scan_lines % update the command line status disp(&#39;&#39;); disp([&#39;Computing scan line &#39; num2str(scan_line_index) &#39; of &#39; num2str(number_scan_lines)]); % load the current section of the medium medium.sound_speed = sound_speed_map(:, medium_position:medium_position + Ny - 1, :); medium.density = density_map(:, medium_position:medium_position + Ny - 1, :); % run the simulation sensor_data_3D = kspaceFirstOrder3D(kgrid, medium, transducer, transducer, input_args{:}); sensor_data_3DC = kspaceFirstOrder3DC(kgrid, medium, transducer, transducer, input_args{:}); combined_sensor_data_3DC = transducer.combine_sensor_data(sensor_data_3DC); % extract the scan line from the sensor data scan_lines_3D(scan_line_index, :) = transducer.scan_line(sensor_data_3D); scan_lines_3DC(scan_line_index, :) = transducer.scan_line(combined_sensor_data_3DC); % update medium position medium_position = medium_position + transducer.element_width; end</code></pre> <p>Thanks, Shawn</p> <p>k-wave version 1.1<br /> compiled binary file, SSE3<br /> MATLAB R2013a Student<br /> Ubuntu 13.10 </p>