B模式超声成像仿真(MATLAB k-Wave仿真)

1. 仿真方法

将超声换能器中的阵元总数设置为32，对96条扫描线进行B模式成像仿真。

2. 成像模拟

定义超声换能器，依次模拟每个扫描线。对于每个扫描线，更新计算域内medium.sound_speed和medium.density的值。然后返回并存储每个换能器阵元（在本例中为32）处的压力。

for scan_line_index = 1: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 = kspaceFirstOrder3D(kgrid, medium, transducer, transducer, input_args{:});
 	
	% 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

3. 数据处理

在超声成像模拟之后，需要多个步骤来将原始回波信号处理成超声图像。

4. 成像结果

谐波图像提高了分辨率和对比度。

相关代码获取方式：请关注微信公众号：超声前沿，回复关键词：kwave