计算正则网格上的Marching Cubes解决方案。
mesh
: gempy.core.grid.TensorGrid
的实例property_name
: 包含地质属性值的网格中的属性名称。threshold
: float, 用于地质属性值的 Marching Cubes 截取参数。默认值为0.5。mode
: 髙级别筛选算法的模式, 参见 skimage.measure.marching_cubes_lewiner
。默认值为'L'
。在指定的网格上执行 Marching Cubes 算法以计算地质体的等值面。
sol = Solution(grid)
sol.compute_model()
vertices, faces = sol.compute_marching_cubes_regular_grid(mesh=grid, property_name='lith', threshold=0.5, mode='L')
Lewiner, Thomas, et al. "Marching cubes: A high resolution 3D surface construction algorithm." ACM Siggraph Computer Graphics. Vol. 26. No. 3. ACM, 1992.
Lewiner, Thomas, Helio Lopes, and Antonio Wilson Vieira. "A survey of surface reconstruction from points." Computer Graphics Forum. Vol. 30. No. 1. Oxford, UK: Blackwell Publishing Ltd, 2011. 1-28.