“““
python運(yùn)行g(shù)prmax
讀取.in文件
運(yùn)行api函數(shù)模擬
“““

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“?#項(xiàng)目目錄
filename?=?os.path.join（dmax‘concrete_Ascan_2D.in‘）
#正演??n：仿真次數(shù)（A掃描次數(shù)）->B掃描
api（filename?n=1?geometry_only=True）??#geometry_only：僅幾何圖形
#merge_files（r“.\GprmaxCode\concrete_Ascan_2D“?removefiles=True）

#?獲取回波數(shù)據(jù)
#?A?B掃描時(shí)out文件名不一樣
filename?=?os.path.join（r“.\GprmaxCode\concrete_Ascan_2D.out“）
rxnumber?=?1
rxcomponent?=?‘Ez‘
outputdata?dt?=?get_output_data（filename?rxnumber?rxcomponent）

#?保存回波數(shù)據(jù)
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（）

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????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