接收区域的属性。
·f1 接收区域左下角的x坐标，单位米 ·f2 接收区域左下角的y坐标，单位米 ·f3 接收区域右上角的x坐标，单位米 ·f4 接收区域右上角的y坐标，单位米
·f5 接收区域在x轴的步长，即Δx，单位米 ·f6 接收区域在y轴的步长，即Δy，单位米
#snapshot: i1 f1 f2 f3 f4 f5 f6 f7 file1 c1 或
#snapshot: i1 f1 f2 f3 f4 f5 f6 i2 file1 c1 获取模型一个区域，某一时刻的电磁场的快照。
·i1 第几道的扫描，值介于1和#analysis:命令中的总步数。 ·f1 区域左下角的x坐标，单位米 ·f2 区域左下角的y坐标，单位米 ·f3 区域右上角的x坐标，单位米 ·f4 区域右上角的y坐标，单位米
·f5 区域在x轴方向上的采样间隔，单位米 ·f6 区域在y轴方向上的采样间隔，单位米 ·f7 快照的时刻，单位秒 ·i2 快照的迭代次数
·file1 快照存储的文件名（*.out）
·c1 类型，取值围a（文本文件）或b（二进制文件）
#tx_steps: f1 f2 发射天线的移动步长。
·f1 发射天线沿x轴的移动步长，单位米 ·f2 发射天线沿y轴的移动步长，单位米
#rx_steps: f1 f2 接收天线的移动步长。
·f1 接收天线沿x轴的移动步长，单位米 ·f2 接收天线沿y轴的移动步长，单位米
------------------------------------------------------------------------------------------
#line_source: f1 f2 str1 str2 ·f1 发射天线的振幅 ·f2 发射天线的频率 ·str1 信号波形，有五种： ·cont_sine 连续正弦波 ·sine 脉冲正弦波 ·gaussian 高斯波 ·ricker 雷克子波
·user 用户自定义的激励函数
·str2 发射天线的名称
#excitation_file: str1
选择用户自定义的激励文件。自定义激励文件是文本文件。
------------------------------------------------------------------------------------------
五、MATLAB绘图
笔者使用MATLAB7.0，以该版本为例进行说明，其他版本可参考使用。为方便说明，笔者假设读者已对MATLAB具有最基本的了解和使用能力，同时主要以命令来完成操作。读者若对MATLAB一无所知，请另行学习。
1.设定当前目录和搜索路径。
------------------------------------------------------------------------------------------
cd D:\\GPR %设定当前目录 addpath D:\\GPR %设定本次搜索路径
------------------------------------------------------------------------------------------
2.模型文件绘图
------------------------------------------------------------------------------------------
geo = 'fname1.geo';
[mesh,header,media] = gprmax2g(geo);
[M,N] = size(mesh); %获取模型大小 imagesc((1:N)*header.dx,(1:M)*header.dy,mesh) %绘图
xlabel('Distance/m'); %设置x轴标签 ylabel('Depth/m'); %设置y轴标签
------------------------------------------------------------------------------------------
3.输出文件绘图
------------------------------------------------------------------------------------------
out = 'fname1.out';
[Header,Fields] = gprmax(out);
NN = 1:Header.NSteps; %移动次数 Position = Header.dx*Header.tx+(NN-1)*Header.dx*Header.TxStepX); %天线每次所在位置
Data(:,:) = Fields.ez(:,1,:); %转换数组格式
imagesc(Position,Fields.t*1e9,Data) %绘图
colorbar; %显示色彩条 xlabel('Distance/m'); %设置x轴标签
ylabel('Time/ns'); %设置y轴标签
------------------------------------------------------------------------------------------
绘出的图像请及时保存。 例子：
#medium: 6.0 0.0 0.0 0.01 1.0 0.0 concrete #medium: 20.0 0.0 0.0 0.1 1.0 0.0 wet_sand
--------------------------------------------------------------------- #domain: 2.5 0.5
#dx_dy: 0.0025 0.0025 #time_window: 12e-9
--------------------------------------------------------------------- #box: 0.0 0.0 2.5 0.45 wet_sand
----- Get a cylinder of free space and then put a slab of concrete to cut it in half ------------------------------------- #cylinder: 1.05 0.3 0.2 free_space #box: 0.0 0.3 2.5 0.45 concrete #cylinder: 0.25 0.375 0.0125 pec #cylinder: 0.45 0.375 0.0125 pec #cylinder: 0.65 0.375 0.0125 pec #cylinder: 0.85 0.375 0.0125 pec #cylinder: 1.05 0.375 0.0125 pec #cylinder: 1.25 0.375 0.0125 pec #cylinder: 1.45 0.375 0.0125 pec #cylinder: 1.65 0.375 0.0125 pec #cylinder: 1.85 0.375 0.0125 pec #cylinder: 2.05 0.375 0.0125 pec #cylinder: 2.25 0.375 0.0125 pec
#triangle: 0.25 0.3 0.65 0.3 0.45 0.1 free_space #triangle: 1.45 0.3 1.85 0.3 1.65 0.1 free_space #box: 2.05 0.1 2.25 0.3 free_space
#cylinder: 2.15125 0.10125 0.1 wet_sand
------------------------------------------------------------- #line_source: 1.0 900e6 ricker MyLineSource
------------------------------------------------------------- #analysis: 115 bre2.out b
#tx: 0.0875 0.4525 MyLineSource 0.0 12e-9 #rx: 0.1125 0.4525 #tx_steps: 0.02 0.0 #rx_steps: 0.02 0.0
#end_analysis:
------------------------------------------------------------ #geometry_file: bre2.geo #title: BRE Model 2 #messages: y
GPRMAX2D知识点汇总
