该函数用于获取Gempy模型的附加数据,包括节点属性、中心和地质事件。这些数据可以用于可视化、后处理和分析。
gempy.get_additional_data(model_name: str) -> Tuple[pd.DataFrame, pd.DataFrame, pd.DataFrame]
model_name
:字符串类型,表示要获取数据的Gempy模型名称。
该函数返回一个元组,包括以下三个数据帧(DataFrame):
Node_attrs
:节点属性数据,包括坐标、方向、标签等;Center
:模型中心数据,包括模型中心坐标、旋转矩阵、缩放因子等;Events
:地质事件数据,包括地层事件、插入事件等。以下示例演示如何使用get_additional_data
函数获取Gempy模型的附加数据:
import gempy as gp
geo_model = gp.create_model('my_model')
# 填充模型
gp.init_data(geo_model, [0, 100, 0, 100, 0, 100], [2, 2, 2],
path_i="Series.csv", path_o="faults.csv")
gp.map_stack_to_surfaces(geo_model, {"Fault_1": 'fault', "Strat_Series": (
'rock3', 'rock2', 'rock1')})
gp.set_is_fault(geo_model, 'Fault_1')
gp.set_topography(geo_model, [0, 100, 0, 100, 100, 100])
gp.set_section_grid(geo_model, [0, 100, 0, 100], [20, 20])
# 计算模型
gp.compute_model(geo_model)
# 获取附加数据
node_attrs, center_data, event_data = gp.get_additional_data('my_model')
print(node_attrs.head())
print(center_data)
print(event_data)