k-Wave

1. HustYjq
   

   Member
   Joined: Jul '20
   Posts: 20

   Hi,
   I want to simulate Photoacoustic Waveforms in 3D.In the examples,the initial pressure distribution is defined as a ball (filled sphere) using makeBall.However,I have a 3DDisk initial pressure.I don’t know whether its shape can be changed.So can anyone tell me how to define it?

   Any help will be appreciated.
   Thanks!

   Posted 4 years ago #

2. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,643

   Hi HustYjq,

   The initial pressure distribution is just a matrix, so you can define it to be any shape. One way to make a disc in 3D is to assign a 2D disc to one slice. For example:

   disc2D = makeDisc(50, 50, 25, 25, 20);
   disc3D = zeros(50, 50, 50);
   disc3D(:, :, 25) = disc2D;
   voxelPlot(disc3D);

   Hope that helps,

   Brad

   Posted 4 years ago #

3. HustYjq
   

   Member
   Joined: Jul '20
   Posts: 20

   Hi Bradley Treeby，
   Thanks for your help.I still have a question.If i have a matrix of initial pressure ,how do i assign it to source.p?In the examples,the pressure in each position is the same.I have not found the intput matrix of initial pressure.

   Thanks

   Posted 4 years ago #

4. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,643

   If you're modelling an initial value problem (i.e., photoacoustics), you can assign the matrix to source.p0.

   Posted 4 years ago #

RSS feed for this topic

Reply

You must log in to post.