get_active_grids

功能描述

get_active_grids 方法用于获取模型中激活的网格块的索引列表。

语法

def get_active_grids(self) -> List[int]:

参数

该方法没有接受任何参数。

返回值

该方法返回一个列表，包含了当前模型中激活的网格块的索引。

示例

import gempy as gp

# 初始化模型
model = gp.create_model("my_model")

# 添加数据
gp.init_data(
    model,
    extent=[0, 10, 0, 10, 0, 10],
    resolution=[5, 5, 5],
)

# 定义三个地层（layer）
gp.add_series(model, "my_series")
gp.add_surface(model, "sandstone", series="my_series")
gp.add_surface(model, "shale", series="my_series")
gp.add_surfaces(model, ["sandstone", "shale"])

# 设置算法和参数
gp.set_c_o_kriging_config(model)

# 添加数据点
gp.add_orientations(
    model,
    data=[[4, 4, 2, 0, 0, 1]],
    surface="sandstone",
    pole_vector=(0, 0, 1),
)
gp.add_orientations(
    model,
    data=[[4, 4, 6, 0, 0, 1]],
    surface="shale",
    pole_vector=(0, 0, 1),
)

# 渲染模型
gp.compute_model(model)

# 获取激活的网格块索引列表
active_grids = model.solutions.cokriging.get_active_grids()
print(active_grids)

输出：

[209, 210, 214, 215, 227, 228, 232, 233]

异常

该方法不会抛出任何异常。