http://www.k-wave.org/forum/topic/simulate-ultrasound-image-from-ct Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 01:36:01 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/simulate-ultrasound-image-from-ct#post-6783 Sun, 17 Mar 2019 22:37:35 +0000 Bradley Treeby 6783@http://www.k-wave.org/forum/ <p>Hi Kunal,</p> <p>The volume data is the CT image in Hounsfield units. Although note an image is just a matrix of numbers. Unfortunately, I don't have any idea about CT data for the feet.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/simulate-ultrasound-image-from-ct#post-6758 Wed, 13 Feb 2019 18:27:34 +0000 khoslaygr26 6758@http://www.k-wave.org/forum/ <p>Hi<br /> In the above code what is exactly in the Volume data.<br /> where can one find this CT data for the feet (specifically the calcenous or heel bone).<br /> Are they images or simply matrices of HU values?<br /> I want to use the data to simulate ultrasound within the heel.<br /> Thanks<br /> Kunal </p> http://www.k-wave.org/forum/topic/simulate-ultrasound-image-from-ct#post-6456 Sun, 13 May 2018 19:59:47 +0000 Bradley Treeby 6456@http://www.k-wave.org/forum/ <p>Hi dhchen,</p> <p>That looks like the correct approach. You will likely also need to add some small scale scatterers (e.g., using random noise) to the background values so you can see contrast on the generated ultrasound image.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/simulate-ultrasound-image-from-ct#post-6441 Sat, 28 Apr 2018 09:07:22 +0000 dhchen 6441@http://www.k-wave.org/forum/ <p>Hi,<br /> I have a CT volume and I want to simulate Ultrasound image from it. I have used hounsfield2density to get the speed and density, but I have no idea how to adjust the parameters in example_us_bmode_linear_transducer.m.</p> <p>% load volumeData from zraw<br /> Nx = size(volumeData,1); % [grid points]<br /> Ny = size(volumeData,2); % [grid points]<br /> Nz = size(volumeData,3); % [grid points]</p> <p>...</p> <p>for scan_line_index = 1:number_scan_lines<br /> ...<br /> [ct_density, ct_sound_speed] = hounsfield2density(volumeData);<br /> medium.sound_speed = ct_sound_speed(:, medium_position:medium_position + Ny - 1, :);<br /> medium.density = ct_density(:, medium_position:medium_position + Ny - 1, :);<br /> % medium.sound_speed = sound_speed_map(:, medium_position:medium_position + Ny - 1, :);<br /> % medium.density = density_map(:, medium_position:medium_position + Ny - 1, :);<br /> ...</p> <p>Is it right? Can anyone give me some suggestions? Thanks. </p>