k-Wave User Forum » Topic: Drive the input signals

k-Wave User Forum » Topic: Drive the input signals http://www.k-wave.org/forum/topic/drive-the-input-signals Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 03:07:07 +0000 http://bbpress.org/?v=1.0.2 <![CDATA[Search]]> q http://www.k-wave.org/forum/search.php JOKEjoe on "Drive the input signals" http://www.k-wave.org/forum/topic/drive-the-input-signals#post-4955 Fri, 23 Jan 2015 03:53:21 +0000 JOKEjoe 4955@http://www.k-wave.org/forum/ Hi Brad , thanks for your answer. joe Bradley Treeby on "Drive the input signals" http://www.k-wave.org/forum/topic/drive-the-input-signals#post-4953 Thu, 22 Jan 2015 22:27:05 +0000 Bradley Treeby 4953@http://www.k-wave.org/forum/ Hi JOKEjoe, If you want to use multiple source transducers in the same simulation, you will need to define these individually using <code>source.p_mask</code> and <code>source.p</code> (rather than using the <code>kWaveTransducer</code> class). Essentially, you define the position and shape of the transducers using the source mask, and then assign the time-varying signals used to drive the grid points within the source mask to <code>source.p</code>. There is a lot more information on how to get started in the <a href="http://www.k-wave.org/documentation/k-wave_time_varying_source_problems.php">Time Varying Source Problems</a> examples (these are accessible in the MATLAB help browser), and in the <a href="http://www.k-wave.org/documentation.php">k-Wave manual</a>. Brad. JOKEjoe on "Drive the input signals" http://www.k-wave.org/forum/topic/drive-the-input-signals#post-4944 Thu, 22 Jan 2015 14:00:26 +0000 JOKEjoe 4944@http://www.k-wave.org/forum/ HI, If my understanding to the question is accurate,maybe i face the same question.i am a college student,and i have a homework recendly:i kown the simple use of k-wave 3D,and my goal is use the tool to achieve two point communicate with one receiver,the receiver can receive the two signal at the same time,but the k-wave only support one transducer communicate with one receiver,i try to define two different transducer(different coordinate,transduce different signal),then maybe i should run kspaceFirstOrder3D,now i do not kown how to make the two transducer work in the same kgrid.later,i want to use frequency demodulator to separate the two signal.It puzzle me for a long time. Please, any suggestion is welcome. THX gigi on "Drive the input signals" http://www.k-wave.org/forum/topic/drive-the-input-signals#post-4881 Tue, 09 Dec 2014 21:01:19 +0000 gigi 4881@http://www.k-wave.org/forum/ Hi Brad, thank you, but if I had 32 transmitting elements and one receiving how is the received signal? it is the result of a sort of superposition? It is possible to download the example_us_bmode_linear_transducer_modified.m from here: <a href="https://www.dropbox.com/s/v5bfybo4i3n36f5/example_us_bmode_linear_transducer_modified.m?dl=0" rel="nofollow">https://www.dropbox.com/s/v5bfybo4i3n36f5/example_us_bmode_linear_transducer_modified.m?dl=0</a> I did not change many things.. thank you so much. Gg Bradley Treeby on "Drive the input signals" http://www.k-wave.org/forum/topic/drive-the-input-signals#post-4873 Fri, 05 Dec 2014 18:10:56 +0000 Bradley Treeby 4873@http://www.k-wave.org/forum/ Hi Gg, When an object of the <code>kWaveTransducer</code> class is used as a sensor, the output signals are those received by the individual elements. No beamforming is performed unless you explicitly call the <code>scan_line</code> method. For example, if you run the same simulation with 32 active elements or 1 active element, the signal received by the 1 active element would be the same in both cases. I am not sure about the NaN. If you post your full code, I will try to run it and see what the problem is. Brad. gigi on "Drive the input signals" http://www.k-wave.org/forum/topic/drive-the-input-signals#post-4861 Thu, 27 Nov 2014 18:51:33 +0000 gigi 4861@http://www.k-wave.org/forum/ Hi Brad, Thanks a lot for you answer, but for me is not clear if the output signal registered by one element is the result of the superposition of 32 signals that come from 32 transmitting elements at once or no. Can you explain me this ? To change the transmit and receive configuration I defined two different objects of the kWaveTransducer, I posted the code before your answer, but I have a lot of NAN, I do not understand why. I used a step of one element. Please, any suggestion is welcome. Gg Bradley Treeby on "Drive the input signals" http://www.k-wave.org/forum/topic/drive-the-input-signals#post-4858 Tue, 25 Nov 2014 18:37:48 +0000 Bradley Treeby 4858@http://www.k-wave.org/forum/ Hi Gg, To answer your first question, the simulation is performed a little differently than how you describe. All 32 elements are switched on at once for both transmit and receive, and are used to form a single scan line. There is a detailed explanation of how this works in the html help for the example in MATLAB (you can also access this <a href="http://www.k-wave.org/documentation/example_us_bmode_linear_transducer.php">online</a>). If you want to transmit and receive on different numbers or configurations of elements, then you need to define two different objects of the <code>kWaveTransducer</code> class. Hope that helps, Brad. gigi on "Drive the input signals" http://www.k-wave.org/forum/topic/drive-the-input-signals#post-4855 Sun, 23 Nov 2014 10:39:45 +0000 gigi 4855@http://www.k-wave.org/forum/ Dear All, below there is the code that I wrote, that has some problem because my output (sensor_data) is full of NAN.... <pre><code>for scan_line_index = 1:number_scan_lines % update the command line status disp(''); disp(['Computing scan line ' num2str(scan_line_index) ' of ' 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 for index_active_element=1:transducer.number_elements transducer1.active_elements(index_active_element) = 1; sensor_data = kspaceFirstOrder3D(kgrid, medium, transducer1, transducer, input_args{:}); if save_mat_files file_name = sprintf('examples_us_pre_beamforming_sensor_data_full_scan_aperture_%s_tr_%s', num2str(scan_line_index), num2str(index_active_element)); cmd=['save ',file_name,'.mat sensor_data']; disp(cmd) eval(cmd); end transducer1.active_elements(index_active_element) = 0; end % extract the scan line from the sensor data % scan_lines(scan_line_index, :) = transducer.scan_line(sensor_data); % update medium position medium_position = medium_position + transducer.element_width; end</code></pre> gigi on "Drive the input signals" http://www.k-wave.org/forum/topic/drive-the-input-signals#post-4853 Sat, 22 Nov 2014 19:05:50 +0000 gigi 4853@http://www.k-wave.org/forum/ I apologize, but I wrote transducer1.mask(i) = 1; but was transducer_input.mask thanks a lot. Gg gigi on "Drive the input signals" http://www.k-wave.org/forum/topic/drive-the-input-signals#post-4851 Sat, 22 Nov 2014 15:36:06 +0000 gigi 4851@http://www.k-wave.org/forum/ I apologize but I forget to tell that what are N and M in the NxM matrix: N is the number of active elements and M is the number of samples thanks a lot. Gg gigi on "Drive the input signals" http://www.k-wave.org/forum/topic/drive-the-input-signals#post-4850 Sat, 22 Nov 2014 15:28:40 +0000 gigi 4850@http://www.k-wave.org/forum/ Hi All, Hi Dr. Cox and Dr. Treeby, The example "example_us_b_mode_linear_transducer", % ...illustrates how k-Wave can be used for the simulation of % B-mode ultrasound images (including tissue harmonic imaging) analogous to % those produced by a modern diagnostic ultrasound scanner. It builds on % the Defining An Ultrasound Transducer, Simulating Ultrasound Beam % Patterns, and Using An Ultrasound Transducer As A Sensor examples. just to use the same words of the authors :-) The script in the RUN SIMULATION session is principally constituted by a loop on the number of scan lines (in this case 96), inside it is performed the generation of a phantom and contextually it is applied the scan_lines method; At the input of kspaceFirstOrder3D function we have as source and as output the same transducer object kspaceFirstOrder3D(kgrid, medium, transducer, transducer, input_args{:}); For what concern the output we have a sensor_data matrix of nxm elements Now my question, every one single signal received it is generated by only one active element at a time? I understand that the simulation is performed as below. If we have 32 active elements, the first one element release the input pulse and register the returned pulse, after this, in sequence, the second element release the input pulse and register the returned pulse, the third element release the input pulse and register the returned pulse and so on... until all 32 signals were recorded, it is true? I would like to change the things in way to register with all 32 elements the signal that come from only one element, and choosing from time to time which element or which group of elements will be active in input. At the limit it could have the received signal for each combination of the elements transmitters and receivers (openings full scan), where it will sequentially for each transmitting element all the receiving elements are illuminated. That is, if I for example, an array of 32 elements, for each of them which emits a pulse, it will have 32 that receive up to cover all combinations to have a total of 1024 combinations. I think that it will be possible to do this putting the kspaceFirstOrder3D inside a inner loop (32 or different number_of_active_elements iterations) after defining a new object transducer ie using a transducer_input (that will change inside at the loop) as input and a transducer_output object as output, it will be something like this: for i=1:step:number_of_elements transducer1.mask(i) = 1; sensor_data = kspaceFirstOrder3D(kgrid, medium, transducer_input(i), transducer_output, input_args{:}). end but maybe there is a solution more smart than this mine... How can I do it? I hope that my questions are clear. thanks a lot Gg