ultrasonic.m

The function runs a k-wave acoustic and thermal simulation in CPP interactive mode and water+skull setup. This to avoid extra walltime once it is performed into the ultrasonic_loop.m. As inputs it takes delta_y and delta_z to add to the geometric focus y and z coordinates, while x coordinate is computed as transducer_x_coordinate = geometric_focus_x_coordinate + focal length. It returns the max peak pressure map, Isppa map, and thermal map.

DOWNLOAD ultrasonic.m PARAMETERS Simulation Medium Parameters

ultrasonic_loop.m

The loop allows to use the function ultrasonic.m in multiple simulations, saving the outputs in different matrices.

The code contains also a mask, because sometimes the output values range is too wide to be visualized well on a plot. Then, it is possible to mask the mean and standard deviation maps from a certain value, in order to have a figure of the values of interest and/or regions of interest.

In the terminal on the directory where the file is located, the ultrasonic_loop.m runs with the use of the command ‘matlab_sub ultrasonic_loop.m’. The time and memory limits can be added to the command or chosen by default.

The output data is saved in the same folder where the ultrasonic_loop.m is located. The output data is ultrasonic_output.mat and it can be loaded on the workspace in Matlab.

DOWNLOAD ultrasonic_loop.m

References

  1. B. E. Treeby and B. T. Cox, “k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave-fields,” J. Biomed. Opt., vol. 15, no. 2, p. 021314, 2010.

  2. B. E. Treeby, J. Jaros, D. Rohrbach, and B. T. Cox, “Modelling elastic wave propagation using the k-Wave MATLAB toolbox,” IEEE International Ultrasonics Symposium, pp. 146-149, 2014.