k-Wave

1. Dustin Qi
   

   Member
   Joined: Aug '16
   Posts: 3

   Now I'm modelling the complete problem of photoacoustic tomography, which means that I need to simulate the light propagation(using TOAST software) before I use k-wave to model acoustic problem. But k-Wave and TOAST use different numerical techniques (regular mesh in k-Wave and irregular mesh in TOAST).So I have to interpolate the function from the nodes of the TOAST basis onto the points in the k-Wave basis in Matlab.
   Is this way is reliable and accurate? And how to do it? I mean how to interpolate functions from one basis to another? Does anybody have done that before or are there some information or papers about it? I really have no ideas about how to do it.
   Any answer will be helpful to me:) Thanks~

   Posted 6 years ago #

2. Felix Lucka
   

   Member
   Joined: Sep '14
   Posts: 1

   Hi!

   First about how to do it:
   - Have a look at the toast Matlab documentation "TOAST > Definitions > Basis definitions": This is one build-in way to map from unstructured to regular meshes within toast. The following paper might be helpful: http://onlinelibrary.wiley.com/doi/10.1002/nme.5271/full
   - Another way is to use Matlab's build in functions: Use scatteredInterpolant to interpolate data from an unstructured mesh onto a structured one.

   The second thing is about whether it is reliable and accurate. If you smooth the initial pressure in k-Wave, I don't see any potential instability. Of course you can think more carefully about the coupling in mathematical terms, i.e., which functions spaces are used and whether the interpolation fits to this function space setting (but I would expect that it won't make much difference in practical terms). About the accuracy: Of course having a fine mesh in toast (especially around the source) will lead to a higher numerical accuracy, but I would estimate that for experimental data, the more problematic part is the model error that you get from using the diffusion approximation for modelling light propagation.

   I hope that helps!

   Posted 6 years ago #

RSS feed for this topic

Reply

You must log in to post.