get_surfaces

描述

get_surfaces 是 Gempy 库中的一个函数，用于获取指定模型中的地层面。地层面是在样本中定义的，它们刻画了岩层的分层结构及物性分布。

语法

get_surfaces(model=None, name=None, indices=None, as_array=False)

参数

• model ：一个 GeMpy Model 类型的对象，在该对象中包含了定义地质场的所有信息，包括地层面、采样数据、三维网格等。默认值为 None。

• name ：一个字符串类型的变量，指定要获取的地层面的名称。默认值为 None。如果不指定 name，那么将返回所有的地层面。

• indices ：一个 tuple 类型的变量，用于指定获取哪个采样点或者空间网格点的地层面。默认值为 None。如果 indices 为 None，那么将返回整个模型的地层面。

• as_array ：一个布尔型变量，表示返回数据类型。如果为 True，那么将返回一个二维数组类型的数据；如果为 False，那么将返回一个带有名称、颜色、位置等数据的 Pandas DataFrame 格式数据。

返回值

• 当 as_array 为 False 时，返回一个 Pandas DataFrame 数据格式，包括以下几列：

  • Surface：地层面的名称。

  • X, Y, Z：地层面在三维空间中的位置坐标。

  • values：地层面对应的属性值。

  • colors：地层面对应的颜色值。

• 当 as_array 为 True 时，返回一个二维数组，其中每一行都代表着一个采样点或者空间网格上的点。数组的第一列为地层面名称，后续各列为坐标、属性值以及颜色等数据。

示例

import gempy as gp

# 创建一个模型
geo_model = gp.create_model('My Model')

# 定义地质场模型
gp.init_data(geo_model, [0, 100, 0, 100, 0, 100], [50, 50, 50],
             surfaces=['surface1', 'surface2', 'basement'])

# 计算模型
gp.compute_model(geo_model)

# 获取指定地层面的数据
surfaces_data = gp.get_surfaces(model=geo_model, name='surface1')

# 获取整个模型的数据
all_surfaces_data = gp.get_surfaces(model=geo_model)

# 将数据转换为数组格式
surfaces_array = gp.get_surfaces(model=geo_model, name='surface2', as_array=True)

参考文献

1. Di Giuseppe, E., Ialongo, S., Caumon, G., & Lécuyer, C. (2019). Gempy 2.0: a Python package for quantitative analysis of geological structures from implicit 3D geological modeling. Geoscientific Model Development, 12(3), 1165-1179.