create_section_grid

该方法用于创建一个“section grid”，即在地质体中提取出一个平面，从而创建一个网格。具体而言，该方法通过在指定区域内定义一个“section plane”，并且沿着该平面创建一个网格，从而获得一个三维视角的地质体剖面。

输入参数

• resolution：float类型，表示网格分辨率。
• extent：list类型，包含两个元素，表示section plane的起点和终点坐标。
• bounding_box：list类型，包含六个元素，表示该section plane所处的局部区域的范围。这些元素为xmin、ymin、zmin、xmax、ymax、zmax。

输出结果

• SectionGrid：该对象包括提取出来的平面地质体，包括它的模拟属性和所在位置。可以读取或更新该属性来交互地重新构建图形。

示例

import gempy as gp

# Define the resolution
resolution = 50

# Define the extent of the section
extent = [0, 200]

# Define the bounding box of the section
bounding_box = [0, 0, 0, 200, 200, 1000]

# Instantiate a GeoModel
geo_model = gp.create_model("Section Example")

# Add a lithology to the model
gp.init_data(
    geo_model,
    extent=bounding_box,
    resolution=[resolution, resolution, resolution],
    default_values=1,
)

# Create the section grid
section_grid = geo_model.grid.create_section_grid(resolution=resolution, extent=extent, bounding_box=bounding_box)

# Interact with the SectionGrid object
section_grid.update_grid_properties(
    property="lithology",
    values=geo_model.solutions.lithology[:, section_grid.extent[0]:section_grid.extent[1], :],
    resolution=res,
)