520489f3fa80b0538a60b4f96ea53b56.zip

520489f3fa80b0538a60b4f96ea53b56.zip

“““
python运行gprmax
读取.in文件
运行api函数模拟
“““

import os
import numpy as np
import matplotlib.pyplot as plt
from gprMax.gprMax import api
from tools.outputfiles_merge import get_output_data merge_files



#文件路径+文件名
dmax=r“.\GprmaxCode“ #项目目录
filename = os.path.join（dmax‘concrete_Ascan_2D.in‘）
#正演  n：仿真次数（A扫描次数）->B扫描
api（filename n=1 geometry_only=True）  #geometry_only：仅几何图形
#merge_files（r“.\GprmaxCode\concrete_Ascan_2D“ removefiles=True）

# 获取回波数据
# A B扫描时out文件名不一样
filename = os.path.join（r“.\GprmaxCode\concrete_Ascan_2D.out“）
rxnumber = 1
rxcomponent = ‘Ez‘
outputdata dt = get_output_data（filename rxnumber rxcomponent）

# 保存回波数据
np.savetxt（‘concrete_Ascan_2D.txt‘outputdatadelimiter=‘ ‘）

## B扫描绘图
#from tools.plot_Bscan import mpl_plot
#plt = mpl_plot（filenameoutputdata dt*1e9 rxnumber rxcomponent）
#plt.ylabel（‘Time [ns]‘）
#plt.show（）

## A扫描绘图
from tools.plot_Ascan import mpl_plot
from gprMax.receivers import Rx
outputs = Rx.defaultoutputs
outputs = [‘Ez‘]
print（outputs）
plt = mpl_plot（filename outputs）
plt.show（）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2020-07-03 10:38  .idea\
     目录           0  2020-05-30 19:37  .idea\inspectionProfiles\
     文件         228  2020-05-30 19:37  .idea\inspectionProfiles\profiles_settings.xml
     文件         197  2019-12-07 09:11  .idea\misc.xml
     文件         267  2019-12-07 09:11  .idea\modules.xml
     文件        5660  2020-07-03 10:38  .idea\workspace.xml
     文件       15759  2020-07-03 10:38  concrete_Ascan_2D.txt
     目录           0  2020-07-03 10:38  GprmaxCode\
     文件         398  2020-05-20 11:57  GprmaxCode\concrete_Ascan_2D.in
     文件       26576  2020-05-15 20:07  GprmaxCode\concrete_Ascan_2D.out
     文件      121050  2020-07-03 10:38  GprmaxCode\concrete_model.vti
     文件      121050  2020-05-30 20:01  GprmaxCode\concrete_model1.vti
     文件      121050  2020-05-30 20:01  GprmaxCode\concrete_model2.vti
     文件        1234  2020-07-03 10:38  model_building.py
     文件        2740  2019-12-07 19:59  wavedata处理.py