k-Wave Toolbox |
![]() ![]() |
On this page… |
---|
This example demonstrates the use of k-Wave for the reconstruction of a three-dimensional photoacoustic wave-field recorded over a planar sensor array. The sensor data is simulated and then time-reversed using kspaceFirstOrder3D
. It builds on the 3D FFT Reconstruction For A Planar Sensor and 2D Time Reversal Reconstruction For A Line Sensor examples.
The sensor data is simulated using kspaceFirstOrder3D
in the same way as in the preceeding examples.
The utilised initial pressure distribution and the sensor mask are shown in 3D FFT Reconstruction For A Planar Sensor.
The time-reversal reconstruction is then performed by setting the optional input 'TimeRev'
to true
.
% add time reversal flag to the input options input_args = [input_args, {'TimeRev', true}]; % run the time-reversal reconstruction p0_recon = kspaceFirstOrder3D(sensor_data, kgrid, c, rho, t_array, sensor_mask, input_args{:});
A plot of the reconstructed initial pressure distribution with a positivity condition is shown below.
![]() |
2D Time Reversal Reconstruction For A Circular Sensor | 3D Time Reversal Reconstruction For A Circular Sensor | ![]() |
© 2009 Bradley Treeby and Ben Cox.