k-Wave

1. asi
   

   Member
   Joined: Oct '14
   Posts: 3

   Hi there,

   can someone explain how the iterative image reconstruction work ?

   Thanks,

   asi.

   Posted 10 years ago #

2. bencox
   

   Developer
   Joined: Mar '10
   Posts: 292

   Hi asi,

   Is this a photoacoustic imaging question? Iterative image reconstruction is much too big a topic to cover here in any general sense. k-Wave doesn't have any built-in iterative photoacoustic reconstruction algorithms, but you could iterate them by using the reconstructed image as a source.p0 in kSpaceFirstOrder3D, saving the output at the positions of the detectors and reconstructing again. If you are detecting over a limited area, ie. not all around the object, then this may improve your image as each time you run the forward model it reimposes the (physically true) condition that the initial particle velocity is zero. There is an example in the Toolbox that does this: Iterative Image Improvement Using Time Reversal Example.

   Ben

   Posted 9 years ago #

3. asi
   

   Member
   Joined: Oct '14
   Posts: 3

   Hi Ben,

   can you please direct me to some references regarding iterative image reconstruction, like papers or an educative book on imaging.

   thanks,

   asi.

   Posted 9 years ago #

4. bencox
   

   Developer
   Joined: Mar '10
   Posts: 292

   asi,

   We're here to help with k-Wave, not anything that you could look up yourself! Try googling it as a start. Please post here only if you have a question related to k-Wave.

   Best wishes,
   Ben

   Posted 9 years ago #

5. wang mengxiao
   

   Member
   Joined: Oct '17
   Posts: 11

   Why i can't simulate Iterative Image Improvement Using Time Reversal in 3D? And can you help me about this?

   Posted 6 years ago #

6. wang mengxiao
   

   Member
   Joined: Oct '17
   Posts: 11

   I modified Iterative Image Improvement Using Time Reversal Example, i changed the kSpaceFirstOrder2D to kSpaceFirstOrder3D, but it shows the error: reference to non-existent field 'p'. Can someone explain why it shows the error. Thank you very much.

   Best,
   Mengxiao

   Posted 6 years ago #

7. bencox
   

   Developer
   Joined: Mar '10
   Posts: 292

   Hi Mengxiao,

   There are two examples of non-iterative 3D time reversal in the toolbox, maybe start with those? Adapt them to be iterative by adding a forward simulation based on the latest estimate of the initial pressure distribution and then time reversing the difference between the measured and modelled time series, then iterating between the time reversal and forward modelling steps. Perhaps use the 2D code as a guide.

   http://www.k-wave.org/documentation/example_pr_3D_tr_planar_sensor.php
   http://www.k-wave.org/documentation/example_pr_3D_tr_spherical_sensor.php

   Best wishes,
   Ben

   Posted 6 years ago #

RSS feed for this topic

Reply

You must log in to post.