DOI: 10.5281/zenodo.7717290
title: Transmission ultrasound data simulated using the k-Wave toolbox as a benchmark for biomedical quantitative ultrasound tomography using a ray approximation to Green's function

description: <strong>Transmission ultrasound data simulated using the k-Wave toolbox as a benchmark for biomedical quantitative ultrasound tomography using a ray approximation to Green's function </strong> The folder ‘’<em>data_ust_kWave_transmission’’ </em>includes the transmission ultrasound data sets used in the project:https://github.com/Ash1362/ray-based-quantitative-ultrasound-tomography ........................................................................................... The ultrasound data were simulated using the k-Wave toolbox (version 1.3.) [5] and using a digital breast phantom [4]. The reported changes in the last version of the k-Wave (version 1.4.) will not affect the simulated signals. The folder ‘’<em>data_ust_kWave_transmission’’ </em> must be added to the path: <em>''…r-Wave/simulation/data/…''</em> For running the Matlab example scripts in the project in the github, the user has two choices: Simulate the k-Wave ultrasound data by setting <em>data_sim=true;</em> in the examples in the project. Upload the already simulated k-Wave ultrasound data according to the description below and load them by setting <em>data_sim=false;</em> in the examples in the project. Please read the description in the example scripts! ………………………………………………………………………………… This folder ‘’<em>data_ust_kWave_transmission’’ </em>includes 2 subfolders, ‘’2D<em>’’ </em> and ‘’3D<em>’’ </em>. ………………………………………………………………………………… The subfolder ‘’2D<em>’’ </em> includes: <strong>data_ust_kWave_transmission/2D/PulsePammoth_1_dx4_cfl1_Nr256_Ne64_Interpoffgrid_Transgeompoint_Absorption1_CodeMatlab/data4_nonsmooth.mat</strong> Two transmission ultrasound data sets were simulated using the k-wave for only water and breast in water according to section <em>‘’6.1. data simulation’’</em> in [1]. 64 emitters and 256 receivers are simulated as off-grid points which are placed on a 2D ring. To simulate the data, each emitter was individually driven by an excitation pulse, and the induced acoustic pressure time series were recorded at all the receivers. The k-Wave simulation was performed on a grid with grid spacing 0.4 mm, and the time spacing was set using a CFL number 0.1. The acoustic absorption and dispersion were accounted for based on the frequency power law. This data set is used for the purpose of image reconstruction, and therefore, the sound speed and absorption coefficients maps are not smoothed, i.e., the original maps are used for simulations. <strong>data_ust_kWave_transmission/2D/PulsePammoth_1_dx4_cfl1_Nr256_Ne64_Interpoffgrid_Transgeompoint_Absorption1_CodeMatlab/<em>data4_smooth_17_1.mat</em></strong> Two transmission ultrasound data sets were simulated using the k-Wave for only water and breast in water as the benchmark for validation of ray approximation to Green’s function in homogeneous and heterogenous media, respectively. The simulation was performed according to section <em>‘’6.2. Numerical validation of the ray approximation to the Green’s function’’</em> in [1]. The pressure field was produced by emitter 1 (of the 64 emitters) and was recorded in time on all 256 receivers. The k-Wave simulation was performed on a grid with grid spacing 0.4 mm, and the time spacing was set using a CFL number 0.1. The acoustic absorption and dispersion were accounted for based on the frequency power law. The sound speed and absorption coefficient maps were smoothed by an averaging window of size 17 grid points. This data set is used as the benchmark for measuring accuracy of ray approximation to Green’s function for computing phase and amplitude of the pressure field on the receivers. <strong>data_ust_kWave_transmission/2D/PulsePammoth_1_dx4_cfl1_Nr256_Ne64_Interpoffgrid_Transgeompoint_Absorption1_CodeMatlab/<em>data4_smooth_17_20.mat</em></strong> This data set is the same as data4_smooth_17_1 except the pressure field is produced by emitter 20. …………………………………………………………………………………………………………………. The subfolder ‘’3D<em>’’ </em> includes: <strong>data_ust_kWave_transmission/3D/PulsePammoth_1_dx5_cfl1_Nr4096_Ne1024_Interpnearest_Transgeompoint_Absorption0_CodeCUDA/data5_nonsmooth_tof_singram.mat</strong> The discrepancy of time-of-flight data for two transmission ultrasound data sets simulated by the k-wave for breast in water and only water according to section 5.2 in [3]. The pressure fields were produced by 1024 emitters separately and were recorded on 4096 receivers. The emitters and receivers were simulated as points which are placed on a 3D hemispherical surface, and are interpolated onto the grid using a neighboring interpolation. The k-Wave simulations were performed on a grid with grid spacing 0.5 mm, and the time spacing was set using a CFL number 0.1. The time-of-flight data were computed and will be used for a refraction-corrected image reconstruction of the sound speed based on the inversion approach proposed in [3]. <strong>References</strong> 1 - A. Javaherian, ❝Hessian-inversion-free ray-born inversion for high-resolution quantitative ultrasound tomography❞, 2022, https://arxiv.org/abs/2211.00316/ . 2 - A. Javaherian and B. Cox, ❝Ray-based inversion accounting for scattering for biomedical ultrasound tomography❞, Inverse Problems vol. 37, no.11, 115003, 2021. https://iopscience.iop.org/article/10.1088/1361-6420/ac28ed/ 3- A. Javaherian, F. Lucka and B. T. Cox, ❝Refraction-corrected ray-based inversion for three-dimensional ultrasound tomography of the breast❞, Inverse Problems, 36 125010. https://iopscience.iop.org/article/10.1088/1361-6420/abc0fc/ 4- Y. Lou, W. Zhou, T. P. Matthews, C. M. Appleton and M. A. Anastasio, ❝Generation of anatomically realistic numerical phantoms for photoacoustic and ultrasonic breast imaging❞, J. Biomed. Opt., vol. 22, no. 4, pp. 041015, 2017. https://anastasio.bioengineering.illinois.edu/downloadable-content/oa-breast-database/ 5 - 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, 021314, 2010. http://www.k-wave.org/
description: Summary: In the linked GitHub project, please add the folder ''data_ust_kWave_transmission'' to the path: ''..../r-Wave/simulation/data/...''
description: {"references": ["A. Javaherian, \u275dHessian-inversion-free ray-born inversion for high-resolution quantitative ultrasound tomography\u275e, 2022, https://arxiv.org/abs/2211.00316/ .", "A. Javaherian and B. Cox, \u275dRay-based inversion accounting for scattering for biomedical ultrasound tomography\u275e, Inverse Problems vol. 37, no.11, 115003, 2021. \u00a0https://iopscience.iop.org/article/10.1088/1361-6420/ac28ed/", "A. Javaherian, F. Lucka and B. T. Cox, \u275dRefraction-corrected ray-based inversion for three-dimensional ultrasound tomography of the breast\u275e, Inverse Problems, 36 125010. \u00a0https://iopscience.iop.org/article/10.1088/1361-6420/abc0fc/", "Y. Lou, W. Zhou, T. P. Matthews, C. M. Appleton and M. A. Anastasio, \u275dGeneration of anatomically realistic numerical phantoms for photoacoustic and ultrasonic breast imaging\u275e, J. Biomed. Opt., vol. 22, no. 4, pp. 041015, 2017. https://anastasio.bioengineering.illinois.edu/downloadable-content/oa-breast-database/", "B. E. Treeby and B. T. Cox, \u275dk-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields\u275e, J. Biomed. Opt. vol. 15, no. 2, 021314, 2010. http://www.k-wave.org/"]}