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