|
33 | 33 | mesh = ss.visualization.to_pyvista_mesh(obj) |
34 | 34 | ss.visualization.pv_plot([mesh]) |
35 | 35 |
|
36 | | -# |
37 | | -# # %% md |
38 | | -# # Four Layers |
39 | | -# # ----------- |
40 | | -# |
41 | | -# # Pull gempy model |
42 | | -# model_files_2 = pooch.retrieve( |
43 | | -# url="https://github.com/cgre-aachen/gempy_data/raw/master/data/gempy_models/Kim.zip", |
44 | | -# known_hash="f530a88351ed0e38673c6937161c59a2f69df92202e14c1e5d5729ed5d72a323", |
45 | | -# processor=pooch.Unzip() |
46 | | -# ) |
47 | | -# |
48 | | -# # %% md |
49 | | -# # Triangular meshes |
50 | | -# # +++++++++++++++++ |
51 | | -# |
52 | | -# # %% |
53 | | -# fname, = [i for i in model_files_2 if "meshes.nc" in i] |
54 | | -# gempy_tri_mesh_unstruct = ss.reader.read_netcdf.read_unstruct(fname) |
55 | | -# tri_surf = ss.TriSurf(gempy_tri_mesh_unstruct) |
56 | | -# mesh = ss.visualization.to_pyvista_mesh(tri_surf) |
57 | | -# ss.visualization.pv_plot([mesh]) |
58 | | -# |
59 | | -# # %% md |
60 | | -# # Regular grid |
61 | | -# # ++++++++++++ |
62 | | -# |
63 | | -# fname, = [i for i in model_files_2 if "regular_grid.nc" in i] |
64 | | -# gempy_struct = ss.reader.read_netcdf.read_struct(fname) |
65 | | -# regular_grid = ss.StructuredGrid(gempy_struct) |
66 | | -# |
67 | | -# # %% |
68 | | -# pyvista_mesh = ss.visualization.to_pyvista_grid(regular_grid, |
69 | | -# data_set_name='property_matrix', |
70 | | -# attribute_slice={'Properties': 'id'} |
71 | | -# ) |
72 | | -# |
73 | | -# ss.visualization.pv_plot([pyvista_mesh]) |
74 | | -# |
75 | | -# |
76 | | -# |
77 | | -# |
78 | | -# # %% |
79 | | -# |
80 | | -# pyvista_mesh = ss.visualization.to_pyvista_grid(regular_grid, |
81 | | -# data_set_name='block_matrix', |
82 | | -# attribute_slice={'Properties': 'id', |
83 | | -# 'Features': 'Default series'}) |
84 | | -# |
85 | | -# ss.visualization.pv_plot([pyvista_mesh]) |
| 36 | + |
| 37 | +# %% md |
| 38 | +# Four Layers |
| 39 | +# ----------- |
| 40 | + |
| 41 | +# Pull gempy model |
| 42 | +model_files_2 = pooch.retrieve( |
| 43 | + url="https://github.com/cgre-aachen/gempy_data/raw/master/data/gempy_models/Kim.zip", |
| 44 | + known_hash="f530a88351ed0e38673c6937161c59a2f69df92202e14c1e5d5729ed5d72a323", |
| 45 | + processor=pooch.Unzip() |
| 46 | +) |
| 47 | + |
| 48 | +# %% md |
| 49 | +# Triangular meshes |
| 50 | +# +++++++++++++++++ |
| 51 | + |
| 52 | +# %% |
| 53 | +fname, = [i for i in model_files_2 if "meshes.nc" in i] |
| 54 | +gempy_tri_mesh_unstruct = ss.reader.read_netcdf.read_unstruct(fname) |
| 55 | +tri_surf = ss.TriSurf(gempy_tri_mesh_unstruct) |
| 56 | +mesh = ss.visualization.to_pyvista_mesh(tri_surf) |
| 57 | +ss.visualization.pv_plot([mesh]) |
| 58 | + |
| 59 | +# %% md |
| 60 | +# Regular grid |
| 61 | +# ++++++++++++ |
| 62 | + |
| 63 | +fname, = [i for i in model_files_2 if "regular_grid.nc" in i] |
| 64 | +gempy_struct = ss.reader.read_netcdf.read_struct(fname) |
| 65 | +regular_grid = ss.StructuredGrid(gempy_struct) |
| 66 | + |
| 67 | +# %% |
| 68 | +pyvista_mesh = ss.visualization.to_pyvista_grid(regular_grid, |
| 69 | + data_set_name='property_matrix', |
| 70 | + attribute_slice={'Properties': 'id'} |
| 71 | + ) |
| 72 | + |
| 73 | +ss.visualization.pv_plot([pyvista_mesh]) |
| 74 | + |
| 75 | + |
| 76 | + |
| 77 | + |
| 78 | +# %% |
| 79 | + |
| 80 | +pyvista_mesh = ss.visualization.to_pyvista_grid(regular_grid, |
| 81 | + data_set_name='block_matrix', |
| 82 | + attribute_slice={'Properties': 'id', |
| 83 | + 'Features': 'Default series'}) |
| 84 | + |
| 85 | +ss.visualization.pv_plot([pyvista_mesh]) |
0 commit comments