surface_points

surface_points是Gempy中Project类的一个重要属性，表示地质图面的点。在Gempy中，surface_points通常用于定义地层中的不同岩石单元的表面。这些点既可以由用户手动定义，也可以从地形数据中提取出来。

基本属性

surface_points属性是一个二维的numpy数组，存储了所有岩石单元表面的点坐标。数组的每一行都表示一个点的三维坐标，列分别表示x，y，z坐标轴。

添加surface_points

可以使用以下方法在Project对象中添加surface_points属性：

import numpy as np
import gempy

geo_model = gempy.create_model('New_model')
x = np.random.rand(10)
y = np.random.rand(10)
z = np.zeros(10)
points = np.column_stack((x, y, z))
geo_model.add_surface_points(points)

在以上代码中，我们首先创建了一个Project对象，并且随机生成了10个平面点坐标。接着，将平面坐标合并成一个numpy数组，并调用add_surface_points方法将这个数组添加到Project对象中。

使用surface_points

在Gempy中，surface_points的最重要的用途是作为模型构建中其他模块的输入。例如：

• 在数据预处理模块中，surface_points可以用于创建插值网格。
• 在块模型构建模块中，surface_points可以用于定义地层的多个模型单元。
• 在模型解析模块中，surface_points可以用于估算地下岩石单元的物理属性。

结论

总之，surface_points是Gempy中的一个重要属性，用于定义地层中的不同岩石单元的表面。了解surface_points的基本属性和添加方法，可以帮助用户更好地使用Gempy进行地质建模分析。