1、软件
1.spyder
2.Gprmax
2、代码及解释
1.首先在spyder中编辑好python部分代码，代码如下：主要是确定含水量、半径、深度三个部分，随机生成需要的范围内的数据；
import scipy.interpolate
import numpy as np
from random import random as rd
m=np.array([12,9.3,6.2,5.5,2.8,0.2])
e_stat = np.array([12.84,11.19,9.14,8.63,6.75,4.814])
e_inf = np.array([7.42,7.2,5.93,6.023,5.503,4.507])
t = np.array([0.611,0.73,0.8,1,2.28,0.82])*10**-9
sigma = np.array([20.6,23,6.7,5.15,2.03,0.606])*10**-3
es = scipy.interpolate.interp1d(m, e_stat)
ei = scipy.interpolate.interp1d(m, e_inf)
t0 = scipy.interpolate.interp1d(m, t)
si = scipy.interpolate.interp1d(m, sigma)
mv=0.2+11.8*rd()
print("#material:{} {} 1 0 concrete n".format(ei(mv), si(mv)))
print("#add_dispersion_debye: 1 {} {} concrete n".format(es(mv)-ei(mv), t0(mv)))
print("#box: 0 0 0 0.5 0.3 0.3 concrete n")
r=0.005 + 0.0245*rd()
z=(0.3-r)*rd()
print("#cylinder: 0.25 0 {} 0.25 0.3 {} {} pec n".format(z,z,r))
2.生成结果：每次都随机生成含水量、半径、深度，含水量有Debye模型来确定；
3.整合到gprmax中可批量生成，代码如下：
#domain: 0.5 0.3 0.4
#dx_dy_dz: 0.001 0.001 0.001
#time_window: 3000
#python:
import scipy.interpolate
import numpy as np
from random import random as rd
m=np.array([12,9.3,6.2,5.5,2.8,0.2])
e_stat = np.array([12.84,11.19,9.14,8.63,6.75,4.814])
e_inf = np.array([7.42,7.2,5.93,6.023,5.503,4.507])
t = np.array([0.611,0.73,0.8,1,2.28,0.82])*10**-9
sigma = np.array([20.6,23,6.7,5.15,2.03,0.606])*10**-3
es = scipy.interpolate.interp1d(m, e_stat)
ei = scipy.interpolate.interp1d(m, e_inf)
t0 = scipy.interpolate.interp1d(m, t)
si = scipy.interpolate.interp1d(m, sigma)
mv=0.2+11.8*rd()
print("#material: {} {} 1 0 concrete n".format(ei(mv), si(mv)))
print("#add_dispersion_debye: 1 {} {} concrete n".format(es(mv)-ei(mv), t0(mv)))
print("#box: 0 0 0 0.5 0.3 0.3 concrete n")
r=0.005 + 0.0245*rd()
z=(0.3-r)*rd()
print("#cylinder: 0.25 0 {} 0.25 0.3 {} {} pec n".format(z,z,r))
filename="concrete_labels"
f=open(filename,'a')
f.write("{}: {} {} {}".format(current_model_run, r,z,mv))
#定义天线
from user_libs.antennas import GSSI
GSSI.antenna_like_GSSI_1500(0.25,0.15,0.3,resolution=0.001)
#end_python:
#geometry_view: 0 0 0 0.5 0.3 0.4 0.002 0.002 0.002 geo n
4.输入指令在gprmax中执行改输入文件：
- python -m gprMax ×××（文件名）-n 2000 ,这个两千是指模型数量
最后结果会生成你所定义的模型数量的out文件；
