gprMax3.0建模时如何自定义目标的几何形状

---

简介

gprMax中包含一些矩形、三角形、球体、圆柱体等一些规则形状的构造命令，但是，当我们想要构造一些无规则形状对象或者gprMax没有定义的对象时，该如何做呢？gprMax 提供了允许用户将自定义的几何形状插入模型。例如，gprMax没有直接构造圆锥体对象的命令，下面是自定义了一个圆锥体目标模型。

需要的命令语句为：

#geometry_objects_read: f1 f2 f3 file1 file2

f1 f2 f3 是在整个区域中应放置几何数组的左下角位置的坐标 (x,y,z) 。
file1 是包含自定义几何数组 HDF5 文件的路径和文件名。
file2是包含#material命令的文本文件的路径和文件名。

---

具体方法步骤

1、生成自定义几何数组HDF5文件（file1文件）

我们要把自定义的目标形状数组保存为HDF5文件，数组必须以Numpy数组的形式存储在以data为根目录的位置，且数据类型为np.int16，对于数组中的0元素，它将与file 2文件中的第一个material的材质属性相关联，对于数组中的-1元素，则保留模型中原有的材质属性。生成HDF5文件的Python代码如下：

import h5py
import numpy as np

def cone(r=1):
    a = np.full([64,64,64],-1)
    for c in range(64):
        for i in range(64):
            for j in range(64):
                if (i - 32) ** 2 + (j - 32) ** 2 < r:
                    a[i][j][c] = 0
        r += 1
    return a
#write HDF5 file
f = h5py.File('my_object.h5','w')
dx_dy_dz = (0.005, 0.005, 0.005)
data = cone()
data = np.array(data,dtype=np.int16)
f.attrs['dx_dy_dz'] = dx_dy_dz
f['/data'] = data

2、生成材料属性文件（file2文件）

材料属性文件是一个文本文件，直接打开记事本写入材料命令就行，如：

#material: 3 0.1 1 0 sand

3、在.in文件中加入#geometry_objects_read命令

#title:  A cone in concrete
#domain: 0.60 0.60 0.60
#dx_dy_dz: 0.005 0.005 0.005
#time_window: 5e-9

#material: 6 0.1 1 0 concrete

#waveform: ricker 1 1.5e9 my_ricker
#hertzian_dipole: z 0.32 0.56 0 my_ricker
#rx:  0.34 0.56 0

#box: 0 0 0 0.60 0.56 0.60 concrete

#geometry_objects_read: 0.10 0.10 0 F:\learn_gprmax3.0\define_my_object\my_object.h5 F:\learn_gprmax3.0\define_my_object\my_material.txt

#box: 0 0.56 0 0.60 0.60 0.60 free_space

#geometry_view: 0 0 0 0.60 0.60 0.60 0.005 0.005 0.005 concrete n

总结

理论上来说，通过该方法可以建立任意形状的目标模型，其关键在于HDF5文件中数组的编写，也可以通过读取图像，将图像矩阵转换为特定数值的几何数组来生成HDF5文件。