set_centered_grid

set_centered_grid()方法是一个Project类中的成员方法。该方法用于设置Gempy模型的中心点以及栅格大小，以便模拟领域的地质属性。

语法

set_centered_grid(self, extent=None, resolution=None, origin=None, **kwargs)

参数

• extent：array_like，可选参数，表示网格被涵盖的空间范围，这是一个四元组，从左到右分别为最小x坐标、最大x坐标、最小y坐标和最大y坐标。
• resolution：array_like，可选参数，表示网格中的单元格数量，从左到右按顺序为x方向单元格数、y方向单元格数和z方向单元格数。当没有给出时，默认按照默认坐标系的尺寸进行网格划分。
• origin：array_like，可选参数，定义网格中心的坐标，分别是x-axis, y-axis和z-axis的浮点坐标。如果未给出，则默认为地图的中心点。
• **kwargs：可变关键字参数，可以添加其他选项参数。

返回值

该方法无返回值。

用法

import gempy as gp
import numpy as np

# 创建Gempy 3D模型对象
geo_model = gp.create_model("NewModel")

# 加载地图数据
gp.init_data(geo_model, [0, 1000, 0, 1000, 0, 1000], [50, 50, 50])

# 定义模型中不同地层的属性
gp.map_series_to_surfaces(geo_model,
    {"Flat_Series": 'Flat', 'Tilted_Series': ('Tilt_1', 'Tilt_2')})

# 设置模型中心点和栅格大小
geo_model.set_centered_grid(
     [
         0, 1000,   # x轴
         0, 1000,   # y轴
         0, 1000,   # z轴
     ], 
     [50, 50, 50],
     [500, 500, 500]
)

# 打印模型中点、边和面的数量
gp.get_surfaces(geo_model)
print("Surfaces - {}".format(geo_model.surfaces.df.shape[0]))
print("Edges    - {}".format(geo_model.edges.df.shape[0]))
print("Points   - {}".format(geo_model.pts.df.shape[0]))

参考文献

• Gempy文档，https://www.gempy.org/documentation/gempy/core/model/
• Gempy示例，https://www.gempy.org/examples/tutorials/1_Getting_started.html