prepare_marchinig_cubes_args

描述

prepare_marching_cubes_args是Gempy中gempy.core.solution.Solution类的一个函数。它的主要作用是根据模型数据准备Marching Cubes算法的输入参数。

语法

prepare_marching_cubes_args(self, extent, resolution, force_vector_extent=None, output='vertices_and_transform')

参数

• extent（tuple）：表示模型的尺寸，格式为(x_min, x_max, y_min, y_max, z_min, z_max)。
• resolution（tuple）：表示Marching Cubes算法的分辨率，格式为(nx, ny, nz)。
• force_vector_extent（tuple，可选）：表示力向量场的范围，格式为(x_min, x_max, y_min, y_max, z_min, z_max)。如果不提供，则默认使用模型范围。
• output（str，可选）：表示输出的数据类型。默认值为"vertices_and_transform"，表示返回顶点和坐标变换矩阵。

返回值

根据output参数的值返回不同类型的数据：

• 如果output='vertices_and_transform'，则返回(vertices, transform_matrix)，其中vertices是表示网格顶点坐标的3D数组，transform_matrix是坐标变换矩阵。

• 如果output='full_data'，则返回一个字典，包含以下数据：

  • vertices：表示网格顶点坐标的3D数组。
  • triangles：表示三角形顶点索引的2D数组。
  • voxel_centers：表示体素中心坐标的3D数组。
  • voxel_values：表示每个体素的值的1D数组。
  • voxel_indices：表示每个体素编号的1D数组。
  • input_args：保存了原始输入参数的字典。

示例

import gempy as gp

# 创建模型
geo_model = gp.create_model('Model1')

# 定义模型参数
gp.init_data(geo_model, extent=[0, 1000, 0, 1000, 0, 1000], resolution=[50, 50, 50],
             path_points=[[100, 100, 999], [800, 800, 1]], surface_points=[[400, 500, 450], [600, 500, 650]],
             orientations=[[-90, 0, 0], [90, 0, 0]])

# 创建Solution对象
sol = gp.compute_model(model=geo_model)

# 准备Marching Cubes参数
vertices, transform_matrix = sol.prepare_marching_cubes_args(extent=[0, 1000, 0, 1000, 0, 1000], resolution=[50,50,50])
print(vertices.shape)  # (29842, 3)
print(transform_matrix)

参考文献

• Gempy官方文档