k-Wave

1. HustYjq
   

   Member
   Joined: Jul '20
   Posts: 20

   Hi!
   I am doing a simulation about proton-induced acoustic emission.I don't know how to set the initial pressure distribution in 2D grid.I am confused about the Example:Photoacoustic Waveforms In 2D.It shows:
   % create initial pressure distribution
   source.p0 = makeDisc(Nx, Nx, Nx/2, Nx/2, source_radius);
   But in makeDisc, it shows Discdisc = makeDisc(Nx, Ny, cx, cy, radius, plot_disc),i don't know what Nx/2 means.
   And i am also confused about Loading the image map,can i import initial pressure data by this way?I am new in k-wave.

   Any help will be appreciated!

   Posted 3 years ago #

2. RuiLiu
   

   Member
   Joined: Jul '20
   Posts: 6

   Hi, I am also new in K wave, I am not really understand your question, but if we want we could discuss in wechat: liurui194114

   Posted 3 years ago #

3. HustYjq
   

   Member
   Joined: Jul '20
   Posts: 20

   Hi，Liu
   I don't know how to set the initial pressure distribution in 2D (not a time varying pressure source distribution).On the same 2d disc, the pressures at different positions are different.Maybe i should use source.p0?

   Posted 3 years ago #

4. prabodhphyee
   

   Member
   Joined: Aug '20
   Posts: 5

   Hi,

   I would suggest running the example: 'example_pr_2D_TR_circular_sensor' [it loads the initial pressure distribution from an image and scales it to the grid. Afterwrds, it generates 2D PA measurements and then uses time reversal to reconstruct the initial pressure source].

   Indeed source.p0 is the initial pressure distribution. You can load it from an image (as demonstrated in the previous example), use the functions such as makedisc to make your pressure source. You can also write your own script to create (various shapes) the pressure source.

   For the makedisc function : look into the function file.
   disc = makeDisc(Nx, Ny, cx, cy, radius, plot_disc)
   You will find:-
   % INPUTS:
   % Nx, Ny - size of the 2D grid [grid points]
   % cx, cy - centre of the disc [grid points]
   % radius - disc radius [grid points]

   Therefore the command "source.p0 = makeDisc(Nx, Nx, Nx/2, Nx/2, source_radius);" simply means that your pressure source is a disc centered at (Nx/2, Nx/2) [center of the grid] and radius = source_radius.

   Hope this helps.

   Regards
   Prabodh

   Posted 3 years ago #

5. HustYjq
   

   Member
   Joined: Jul '20
   Posts: 20

   Hi，prabodhphyee

   Thanks for your reply.Can i make a 3DDisk source in k-wave?In the example,it used makeball.

   Posted 3 years ago #

RSS feed for this topic

Reply

You must log in to post.