I am generating sensor data through the simulation 2D Time Reversal Reconstruction For A Circular Sensor and then I enter them into a network to improve the synogram, however when reconstructing the image, the circle of the sensors is seen, I want to know if you can reconstruct only the central section of the image so that the sensors are not shown, I have been trying to make a kgrid to generate the data as follows:
PML_size = 20; % size of the PML in grid points.
Nx = 256 - 2 * PML_size; % number of grid points in the x direction
Ny = 256 - 2 * PML_size; % number of grid points in the y direction
x = 50.1e-3; % total grid size [m]
y = 50.1e-3; % total grid size [m] dx = 2.3194e-3; % total grid size [m]
dx = x/Nx; % grid point spacing in the x direction [m] dy = 2.3194e-3; % total grid size [m]
dy = y/Ny; % grid point spacing in the y direction [m] kgrid = kWaveGrid [m] kgrid = kWaveGrid [m
kgrid = kWaveGrid(Nx, dx, Ny, dy);
I have subsequently generated a reconstruction kgrid that is smaller with:
Nx = 300; % number of grid points in the x direction.
Ny = 300; % number of grid points in the y direction
xx=20.1e-3;
yy=20.1e-3;
dx = xx/Nx; % grid point spacing in the x direction [m]; dy = 6.7e-5; % number of grid points in the y direction
dy = yy/Ny; % grid point spacing in the y direction [m]
kgrid_recon = kWaveGrid(Nxx, dxx, Nyy, dyy);
but when I do it I get the following error
Error using cart2grid (line 118)
Cartesian points must lie within the grid defined by kgrid.
Error in kspaceFirstOrder_inputChecking (line 616)
[sensor.mask, order_index, reorder_index] = cart2grid(kgrid, sensor.mask,
flags.axisymmetric);
Error in kspaceFirstOrder2D (line 537)
kspaceFirstOrder_inputChecking;
Error in Reconstruccion_Sensor_circular_ejemplo (line 87)
p0_recon = kspaceFirstOrder2D(kgrid_recon, medium, source, sensor,input_args{:});
i.e. my general grid is 50.1mm x50.1mm and the reconstruction grid I need to be 20.1mmx20.1mm.