add_surface_values()方法用于将地表数据添加到模型中。更具体地说,可以在模型中添加样本点的位置和地表值(例如形貌)。
X: 一个numpy.array,包含每个样本点的位置信息。数组的大小必须是[n_samples, n_dimensions],其中n_dimensions是组成每个位置向量的坐标数。values: 一个numpy.array,包含每个样本点的表面值。数组的大小必须是[n_samples, n_additional_data],其中n_additional_data是给定的额外数据(例如地表形貌)的数量。name: 表示给定表面值的名称。如果不指定,该方法默认使用字符串"surface_value"。此方法返回一个三元组(surface_points, surface_values, surface_name):
surface_points:一个包含新加入的表面点的numpy.array。大小必须是[n_new_surface_points, n_dimensions],其中n_new_surface_points表示添加了多少个新的表面点。surface_values:一个包含新加入的表面值的numpy.array。大小是[n_new_surface_points, n_additional_data]。如果没有添加新的表面点,则返回None。surface_name:一个字符串,表示新加入表面值的名称。import numpy as np
import gempy
# 定义位置和表面值
X = np.array([[0, 0, 0], [1, 1, 1], [2, 2, 2], [3, 3, 3]])
values = np.array([[10], [20], [30], [40]])
# 创建一个ImplictCoKriging模型实例
model = gempy.create_model("example_model")
# 添加surfaces_values到模型中,并指定名称
surface_points, surface_values, surface_name = model.add_surface_values(X, values, name="my_surface")
# 打印结果
print(surface_points)
print(surface_values)
print(surface_name)
输出结果:
[[0 0 0]
[1 1 1]
[2 2 2]
[3 3 3]]
[[10]
[20]
[30]
[40]]
my_surface