Read data from gprMAX result in MATLAB

1 H5disp

 see the detail of the FDTD result

h5disp('ZCLW_merged.out')

--------------------------

HDF5 ZCLW_merged.out 
Group '/' 
    Attributes:
        'Title':  'Simulation of zhangchao's paper'
        'gprMax':  '3.1.5'
        'Iterations':  3393
        'dt':  0.000000
        'nrx':  1
    Group '/rxs' 
        Group '/rxs/rx1' 
            Dataset 'Ex' 
                Size:  385x3393
                MaxSize:  385x3393
                Datatype:   H5T_IEEE_F32LE (single)
                ChunkSize:  []
                Filters:  none
                FillValue:  0.000000
            Dataset 'Ey' 
                Size:  385x3393
                MaxSize:  385x3393
                Datatype:   H5T_IEEE_F32LE (single)
                ChunkSize:  []
                Filters:  none
                FillValue:  0.000000
            Dataset 'Ez' 
                Size:  385x3393
                MaxSize:  385x3393
                Datatype:   H5T_IEEE_F32LE (single)
                ChunkSize:  []
                Filters:  none
                FillValue:  0.000000
            Dataset 'Hx' 
                Size:  385x3393
                MaxSize:  385x3393
                Datatype:   H5T_IEEE_F32LE (single)
                ChunkSize:  []
                Filters:  none
                FillValue:  0.000000
            Dataset 'Hy' 
                Size:  385x3393
                MaxSize:  385x3393
                Datatype:   H5T_IEEE_F32LE (single)
                ChunkSize:  []
                Filters:  none
                FillValue:  0.000000
            Dataset 'Hz' 
                Size:  385x3393
                MaxSize:  385x3393
                Datatype:   H5T_IEEE_F32LE (single)
                ChunkSize:  []
                Filters:  none
                FillValue:  0.000000

---------------------

2 h5read

Read the result after  merging all the A-scans.

% %%读取.out文件
% relative path
data_zc=h5read('ZCLW_merged.out','/rxs/rx1/Ez');
FDTD.GPRdata=data_zc';
% % 绘制B-scan图谱
[m,n]=size(FDTD.GPRdata);
tw=4;
FDTD.Fs=10^9*m/tw;
figure(1)
Plotfun(FDTD.GPRdata,FDTD.Fs,0.0075,1,[],0.5)
box off