k-Wave Toolbox |

## 3D Time Reversal Reconstruction For A Planar Sensor Example

On this page…
Overview Simulating the sensor data Performing the reconstruction

### Overview

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.

- open the file in the MATLAB Editor
- run the file in MATLAB

Back to Top

### Simulating the sensor data

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.

Back to Top

© 2009 Bradley Treeby and Ben Cox.