在Gempy中,create_orientation_from_surface_points
函数是用于根据地质表面上的点数据创建方位数据的函数。该函数可以让用户直接根据二维坐标系中的点数据创建出方位数据,并将该数据应用到模型中。
surface_points
:包含表面点坐标的numpy.ndarray
类型数组,形状为(n,2)
,其中n
表示点的数量。rescale
:如果设置为True
,则将表面坐标重新缩放到[0, 1]
的范围内。如果没有指定,则默认为False
。create_orientation_from_surface_points
函数返回一个包含方位数据的pandas.DataFrame
类型数据框。数据框中包含以下列:
X
:表面点的水平坐标。Y
:表面点的垂直坐标。Z
:模型中对应位置的垂直坐标。azimuth
:对应点的方位角,单位为度,以N-S方向为0度,向东逆时针旋转为正方向。dip
:对应点的倾角,单位为度,以水平面为0度,朝上为正方向。import numpy as np
from gempy.core.data_modules.geometric_data import create_orientation_from_surface_points
# 创建二维坐标系中的表面点
surface_points = np.asarray([[0,0], [0,1], [1,0], [1,1]])
# 创建方位数据
orientation_data = create_orientation_from_surface_points(surface_points)
# 显示输出
print(orientation_data)
输出结果:
X Y Z azimuth dip
0 0.0 0.0 0.161290 90.0 90.0
1 0.0 1.0 0.161290 0.0 90.0
2 1.0 0.0 0.161290 180.0 90.0
3 1.0 1.0 0.161290 270.0 90.0
上述示例中,我们创建了一个四个点的表面坐标系,并使用create_orientation_from_surface_points
函数创建了对应的方位数据。函数返回的数据框中包含了每个点的方位角和倾角。