本文使用的软件gprMax
MATLAB R 2021a
gprMax是电磁波数值模拟软件,通过模型的编写(.in/.txt文件)建立相应的模型。in文件的编写网上已经有非常详细的教程了。
.in文件编写好后,通过cmd命令行输入相关的命令,可生成.vti(模型)、.out(单道波形模拟数据)、_merged.out(合并单道波形后的数据)
通过cmd命令行,可查看.out文件的图像。
(rx1)
(model_merged.out_-_rx1)
将merged.out文件直接转换成常见的雷达数据格式,这样更加方便使用者对数据进行分析和处理。
转换中需要用到的是MATLAB和gprMax中的文件。
1、在gprMax中找到outputfile_converter.m(转换工具)
plot_Ascan.m(查看A剖面工具)
plot_Bscan.m(查看B剖面工具)
路径在 \gprMax\tools\MATLAB_scripts
2、打开MATLAB
3、点击主页/打开,然后选择outputfile_converter.m 文件
4、运行命令
5、运行后弹窗选择.out文件。一般是合并单道波形的_merged.out文件
6、选择后弹窗设置相关参数。
(第一个选择分量、探地雷达模拟是向下的,一般选择Ez
第二行选择天线中心频率,单位Mhz,填写与模型中设置的频率一致
第三行天线间距,填写模型中设置的间距,激发源与接收天线的距离
第四行步长,填写模型中的设置步长)
7、点击确认,弹窗选择要转换的文件类型
(1是mala的格式
2是劳雷的格式
3是加拿大ekko的格式
4是迪普瑞达的格式)
8、转换成功后,在上面路径中找到转换后的文件。使用Reflexw可正常打开。MATLAB查看A剖面
打开matlab后打开plot_Ascan.m文件,然后点击执行当前节,选择单道波形的out文件即可。
MATLAB查看B剖面
打开matlab后打开plot_Bscan.m文件,然后点击执行当前节,选择合并单道波形后的out文件即可。这里要注意的时,选择out文件后,会在下方命令行窗口提示选择哪个分量。这里打上Ez即可。