create_regular_grid

该函数为Gempy中gempy.core.data.Grid类的一个方法，用于创建规则网格。规则网格是以相同大小的正方形或立方体网格排列而成的网格，方便用于建模。

函数参数

create_regular_grid(section=None, spacing=None, resolution=None, center=True)

• section: (optional) 区域的筛选范围，由于模型可能比实际区域大，因此需要指定提取模型的区域。
• spacing: 格网的间距，可以输入标量或长度为3的列表，如果不指定则使用默认值。
• resolution: 格点的分辨率，可以输入标量或长度为3的列表，如果不指定则使用默认值。
• center: 格网中心是否与区域的中心重合，默认为 True。

返回值

(grid, extent)。

• grid: 格网的XYZ坐标；
• extent: XYZ轴上的网格数量。

使用方法

创建并可视化一个二维规则网格

import gempy as gp
from gempy.core.data import Grid

gp.set_data(gp.datasets.get_orebody_2v())
extent = [[0,2000], [0,2000]]
spacing = 50
resolution = [1,1]

grid2D = Grid()
grid, ext = grid2D.create_regular_grid(extent=extent, spacing=spacing, resolution=resolution)
grid2D.centers = grid  

gp.plot_2d(grid2D, show_data=False)

创建并可视化一个三维规则网格

import gempy as gp
from gempy.core.data import Grid

gp.set_data(gp.datasets.get_ica_data())
extent = [[0,1000], [0,1000], [0,900]]

spacing = 50
resolution=[2,2,3]

g3D = Grid()
grid, ext = g3D.create_regular_grid(extent=extent, spacing=spacing, resolution=resolution)
g3D.centers = grid

gp.plot_3d(g3D, show_lith=False)

参考资料

Gempy Documentation: https://www.gempy.org/documentation/core/data.html#gempy.core.data.Grid.create_regular_grid