GemPy 3D点可视化故障排除：Python版本兼容性与正确绘图流程

# 使用Conda创建并激活Python 3.10环境
conda create -n gempy_env python=3.10
conda activate gempy_env

# 在激活的环境中安装GemPy
pip install gempy==2.3.1 # 确保安装指定版本

import gempy as gp
import gempy_viewer as gpv
import numpy as np

# 1. 创建地质模型
# 替换 'your_model_name' 为你的模型名称
geo_model = gp.create_model('my_first_gempy_model')

# 2. 初始化地质数据
# 这里的参数需要根据你的实际数据路径和类型进行调整
# 假设你的数据文件名为 'your_data.csv' 并且包含 x, y, z, surface, series 列
# 如果是内部生成数据，则需要构造相应的DataFrame
# 示例：创建一些模拟数据
np.random.seed(123)
n_points = 50
points = np.zeros((n_points, 5))
points[:, 0] = np.random.uniform(0, 1000, n_points) # x
points[:, 1] = np.random.uniform(0, 1000, n_points) # y
points[:, 2] = np.random.uniform(0, 1000, n_points) # z
points[:, 3] = np.random.choice(['surface1', 'surface2', 'surface3'], n_points) # surface
points[:, 4] = np.random.choice(['Strat_Series'], n_points) # series

# 假设我们将这些点作为地层点
geo_model = gp.init_data(
    geo_model,
    points=points,
    # 其他可能的参数如 surfaces, series, resolution等
    # 例如:
    # surfaces=['surface1', 'surface2', 'surface3'],
    # series=['Strat_Series'],
    # resolution=[50, 50, 50]
)

# 确保数据已正确加载
print(geo_model.surfaces)
print(geo_model.series)

# 3. 映射地质序列到表面
# 替换 'your_fault_series_name' 和 'your_strat_series_name'
# 以及 'your_fault_surface_name', 'your_strat_surface_name' 为你的实际名称
gp.map_stack_to_surfaces(
    geo_model,
    {"Fault_Series": 'your_fault_surface_name',  # 如果有断层
     "Strat_Series": ('surface1', 'surface2', 'surface3')}, # 你的地层表面
)

# 可以查看映射后的地质堆栈结构
gp.plot_stack(geo_model)

# 4. 绘制3D模型
gpv.plot_3d(geo_model)

import gempy as gp
import gempy_viewer as gpv
import numpy as np

# 确保Python版本为3.10，GemPy版本为2.3.1

# 1. 创建地质模型
geo_model = gp.create_model('my_complete_gempy_model')

# 2. 准备模拟数据 (实际应用中请加载你的CSV或DataFrame数据)
# 模拟一些地层点
np.random.seed(42)
n_points = 100
points_data = np.zeros((n_points, 5))
points_data[:, 0] = np.random.uniform(0, 1000, n_points) # x
points_data[:, 1] = np.random.uniform(0, 1000, n_points) # y
points_data[:, 2] = np.random.uniform(0, 1000, n_points) # z
points_data[:, 3] = np.random.choice(['surface1', 'surface2', 'surface3'], n_points) # surface
points_data[:, 4] = 'Strat_Series' # series

# 模拟一些断层点 (如果你的模型包含断层)
n_fault_points = 20
fault_points_data = np.zeros((n_fault_points, 5))
fault_points_data[:, 0] = np.random.uniform(200, 800, n_fault_points)
fault_points_data[:, 1] = np.random.uniform(200, 800, n_fault_points)
fault_points_data[:, 2] = np.random.uniform(100, 900, n_fault_points)
fault_points_data[:, 3] = 'fault_surface'
fault_points_data[:, 4] = 'Fault_Series'

all_points = np.vstack((points_data, fault_points_data))

# 定义表面和序列
surfaces = ['fault_surface', 'surface1', 'surface2', 'surface3']
series = ['Fault_Series', 'Strat_Series']

# 3. 初始化地质数据
geo_model = gp.init_data(
    geo_model,
    points=all_points,
    surfaces=surfaces,
    series=series,
    resolution=[50, 50, 50] # 定义模型网格分辨率
)

# 4. 配置地质堆栈（定义地层和断层关系）
gp.map_stack_to_surfaces(
    geo_model,
    {"Fault_Series": 'fault_surface',
     "Strat_Series": ('surface1', 'surface2', 'surface3')},
)

# 5. 设置地质序列类型（断层、地层等）
geo_model.set_series('Fault_Series', order=1, series_type='Fault')
geo_model.set_series('Strat_Series', order=2, series_type='Stratigraphic')

# 6. 计算地质模型
gp.compute_model(geo_model)

# 7. 绘制3D模型
# 确保你安装了gempy_viewer，它是gp.plot_3d的后端
# 3D视图将弹出，并显示地质点和模型表面
gpv.plot_3d(geo_model, show_points=True, show_surfaces=True, show_scalar=True)

# 也可以绘制2D切片进行验证
gpv.plot_2d(geo_model, direction='y', cell_number=25, show_data=True)