use Dr.Tzanis matGPR to transfer FDTD data to Reflexw by MATLAB

Target:

  Aim to use Reflexw to analyze the FDTD data.

 

Reason:

I am lazy to understand the detail of the format of GPR system, fortunetely, matGPR( Andreas Tzanis,University of Athens) has this funciton, so I use its function to transfer data to GPR raw data and use reflexw to analyze it, especially for FDTD analyze.

Prepare

   use GPRmax to generate GPR DATA.

matGPR

  (This softare give you chance do analyze GPR data in Matlab. Great!)

 

  use matGPR to transfer your FDTD data to GPR data.

   matGPR read data into IPD struct varaiable. use it to trasfer your data.

 

1 use matGPR to read a GPR data, such as .rd3 format.

 

2  read your FDTD data into workspace as Data1.

 

3 replace the IPD

If you want to download, go: https://download.csdn.net/download/jaz/12149839

 

 

The followings are the detail.

% FDTD
GPRData=Data1; %% set GPR data
[ns,ntr]=size(GPRData);
% ns- number of samples in a trace;
%ntr- number of traces
timewindow= 8.924*10^(-9);% totaol time for a trace
dx=0.02;% sample distance
%r=1892;
fs=1/(timewindow/ns);% sample frequency
Position=(1:270)*0.002;% sample length for each trace;
traceTime=10^9*(0:ns-1)/fs;% time step of trace;


NIPD=IPD;% new data ,to replace IPD
NIPD.ns=ns; % sample number of traces;
NIPD.ntr=ntr; %trace numbers.
NIPD.dx=dx;% sample distance
NIPD.x=Position;%(1-1:270-1)*0.02;% distance of each trace
NIPD.d=GPRData;%Data of traces;
NIPD.Antenna='[User] HF 1.6 GHz shield';
NIPD.dt=10^9*1/fs;% dt in ns
NIPD.tt2w=traceTime;
NIPD.sigpos=1;

%Replace original data
IPDtemp=IPD;% bak
IPD=NIPD;% repalce

4 use matGPR to export data.

 

 

5 Now, read in ReflexW

use this type for GSSI radar data.

 

Now, we realize our target.

 

Analyze in ReflexW