linefoki.blogg.se

1. #OPTIMAL WAVES NX SETTINGS HOW TO#
2. #OPTIMAL WAVES NX SETTINGS UPDATE#
3. #OPTIMAL WAVES NX SETTINGS SOFTWARE#

Vmcboundary.lightsource(boundary_with_lightsource) = uneisen_parameter = 0.02*ones(Nx, Ny) % Create a light source % Set a light source with a width of 2 mm and cosinic directional profile % in -x directionīoundary_with_lightsource = findBoundaries(vmcmesh, 'direction'.

Vmcmedium.refractive_index(:,Ny/2:end) = 1.4 % refractive index % Define the Gruneisen parameter describing photoacoustic efficiency Vmcmedium.refractive_index = 1.0*ones(Nx, Ny) Vmcmedium.scattering_anisotropy = 0.9 % scattering anisotropy parameter Vmcmedium.absorption_coefficient(disc_indices) = 0.01 Vmcmedium.absorption_coefficient = 0.001*ones(Nx, Ny) vmcmedium.scattering_coefficient = 0.01*ones(Nx, Ny)

#OPTIMAL WAVES NX SETTINGS HOW TO#

See the example 'Working with pixel and voxel data' on how to achieve similar assignments using one dimensional indexing. If one of the arrays defining the medium is given as multidimensional array to ValoMC, the code will assume that the mesh was created using 'createGridMesh' and the output fluence will also given as a two dimensional array in id_fluence. meshgrid) and ValoMC uses matrix(Y,X) Therefore x and y should be swapped when moving between ValoMC arrays and k-Wave arrays vmcmesh = createGridMesh(kgrid.y_vec*1e3, kgrid.x_vec*1e3) % Define optical coefficientsįor users accustomed to k-Wave, the optical coefficients can be set in similar fashion as in k-Wave, i.e. Note that K-Wave uses matrices in the format matrix(X,Y), whereas MATLAB t(c.f.

createGridMesh can be used to create a straightforward mapping between the triangles and pixels. ValoMC uses triangles and tetrahedrons as the basis shape, whereas in k-Wave pixels and voxels are used.

#OPTIMAL WAVES NX SETTINGS UPDATE#

Update codes to use the syntax kgrid = kWaveGrid(.). WARNING: makeGrid will be deprecated in a future version of k-Wave. Medium.sound_speed(disc_indices) = 1800 % Medium.sound_speed = 1500*ones(Nx, Ny) % % Create two internal structures using makeDiskĭiscs = makeDisc(Nx, Ny, 55, 55, 5) + makeDisc(Nx, Ny, 75, 85, 10) Ny = 150 % number of grid points in the y (column) directionĭx = 0.1e-3 % grid point spacing in the x direction ĭy = 0.1e-3 % grid point spacing in the y direction Nx = 150 % number of grid points in the x (row) direction Care must be taken at the initialization ValoMC, to make a matching computational simulation area for (see ValoMC initialization) the photon propagation simulation. The initialisation is done as normal in k-Wave.

#OPTIMAL WAVES NX SETTINGS SOFTWARE#

The software is designed for time domain acoustic and ultrasound simulations in complex and tissue-realistic media. K-Wave is an open source acoustics toolbox for MATLAB and C++ developed by Bradley Treeby and Ben Cox (University College London) and Jiri Jaros (Brno University of Technology). ) and ValoMC works in millimetre-scale (e.g. The example also shows how the computation grid of k-Wave and mesh of ValoMC can be made compatible. The light propagation is simulated using ValoMC and the propagation and detection of pressure wavefield is simulated using k-Wave, see. This example demonstrates simulation of a pressure field generated through the absorption of an externally introduced light pulse.

Simulating the photoacoustic effect using k-Wave: kwavetest.m