CPC G06T 11/006 (2013.01) [A61B 5/0095 (2013.01); G06T 2210/52 (2013.01); G06T 2211/424 (2013.01); G06T 2211/428 (2013.01)] | 7 Claims |
1. A method for photoacoustic imaging comprising:
measuring radiofrequency (RF) time samples from transducer elements wherein the transducer elements are arranged in a 3D shape;
reconstructing an estimated initial pressure image from the measured RF time samples and a pre-calculated forward model matrix wherein the pre-calculated forward model matrix takes into account the arrangement of the 3D shape;
wherein reconstructing comprises minimizing the difference between the pre-calculated forward model matrix applied to the estimated initial pressure image and the measured RF time samples by implementing a superiorized modified conjugate gradient least squares algorithm to minimize the total variation norm;
wherein the pre-calculated forward model matrix factors each of the transducer elements into sub-wavelength mathematical elements.
|