python脚本运⾏gprMax3.0批量仿真GPR数据
python脚本运⾏gprMax3.0批量仿真GPR B-scan图像
1.引⾔
当归的副作用探地雷达（GPR）结合深度学习通常需要⼤量的训练数据集，对于GPR仿真数据集的获取，我们⼀般通过gprMax⽣成，⽽gprMax3.0仿真数据时需要通过cmd命令提⽰符窗⼝⼈⼯⼀条⼀条地输⼊指令（通过cmd命令⽣成GPR B-scan图像：.），对于批量⽣成GPR数据⾮常不⽅便。因此，有必要写⼀些Python脚本，实现⾃动化批量⽣成GPR数据集。
2.Python脚本
该Python脚本可以⼀次性批量⽣成指定⽂件夹中多个GPR输⼊⽂件（.in⽂件或者.txt⽂件）的输出⽂件（.out⽂件），并且可以把输出⽂件（.out⽂件）保存为数据格式（.txt⽂件）和图像格式（.png⽂件）。运⾏该脚本之前需要安装好gprMax3.0（）
import sys
sys.path.append('D:/my_gprmax/gprMax')#把gprMax安装路径添加⾄系统，使import可以找到gprMax模块
import os
from gprMax.gprMax import api
import numpy as np
from tools.plot_Bscan import get_output_data, mpl_plot
from tools.outputfiles_merge import merge_files
宁正import matplotlib.pyplot as plt
num_scan =120#正演仿真次数（A扫描次数）->B扫描
geo_only =Fal#是否只⽣成模型图
path = os.getcwd()#获得当前脚本所在的⽂件路径
#print(path)
root = path+'/in_data'#指定输⼊⽂件的路径
files = os.listdir(root)#得到路径下的⽂件夹名或者⽂件名，形成列表
for file in files:
dswith('.txt'):#找到.txt⽂件
对话框英文# print(file)
filename = root+'/'+file#得到⽂件名的绝对路径
元代历史fi = filename[0:-4]#去掉⽂件名后的.txt后缀，注意保留了前⾯的路径
api(filename, n=num_scan, geometry_only=geo_only)#geometry_only：仅⼏何图形
merge_files(fi, removefiles=True)
"""B扫描绘图"""清华和北大
filename_b = fi+'_merged.out'
rxnumber =1
rxcomponent ='Ez'
#获取回波数据
outputdata, dt = get_output_data(filename_b, rxnumber, rxcomponent)
#保存回波数据
fi2 = fi.split('/')[-1]#得到⽂件名，去掉了前⾯的路径
out_path = path+'/out_data'
if(ists(out_path)):
pass
平安用英语怎么说el:
os.mkdir(out_path)
np.savetxt(out_path+'/'+fi2+'.txt',outputdata,delimiter=' ')#未去除直达波
#绘图
plt.imshow(outputdata, extent=[0, outputdata.shape[1], outputdata.shape[0],0], interpolation='nearest', aspect='auto', cmap='gray', vmin=-np.amax(np.abs(outputdata)), vmax=np.amax(np.abs(outputdata)))
img_path = path+'/img_data'
if(ists(img_path)):
三衢山pass
el:
os.mkdir(img_path)
plt.savefig(img_path+'/'+fi2+'.png',dpi=300)#保存图⽚
# plt.show()
3.可能出现的报错
运⾏之后，如果出现报错：“no module named terminaltables"，则打开cmd命令提⽰符窗⼝，使⽤“pip install terminaltables ”指令，安装好这个模块。若运⾏之后还出现报错：“no module named gprMax.fields_updates_ext"，则gprMax3.0安装的时候可能出现了错误，我的解决办法是，下载: （提取码：3ss8 ），并将这个⽂件替换gprMax⽂件夹下的gprMax模块。
4.数据展⽰
实验批量⽣成了3个输⼊⽂件的输出⽂件，结果如下：
刮大风
