Add plot of stack of 3D sections for geological interpretations about gemgis HOT 6 CLOSED

@AlexanderJuestel I've implemented something like this using pyvista on the lithblock. I think there is a larger idea of making a generic pyvista compatible grid class that could take "any" pyvista object, using pyvista to manage the resulting solution as well

from gemgis.

Hey @AndrewAnnex,

Would you mind sharing that to have a look if it is practical for my needs?

Best regards
Alex

from gemgis.

hey @AlexanderJuestel, I'm pretty busy for the next two weeks but I think I could post some psuedo code, below I post a written pseudo method that should be easy to follow. My current bad approach involves building a huge regular grid (code was posted into a wip branch some time ago, I think the branch is still in my fork which is public) and then slicing the resulting grid object using pyvista, but this approach is actually very inefficient as you end up throwing away a lot of the data that was computed by the model.

My (better) suggestion would be to use pyvista to create the sliced "view" that you see in the figure you posted, and then use the .points numpy array of the object to create a new custom grid object in gempy. You would then use the solution object to get the scalar field values and lithology (aka formation) values to assign pyvista attributes to the pyvista mesh you made to start with (https://docs.pyvista.org/getting-started/what-is-a-mesh.html#what-are-attributes).

I think with a little effort, a new grid object for gempy could be made that would allow for "any" pyvista mesh object to be used directly within gempy for evaluating models, once that becomes stable, all the other grid objects (cross sections, topography, regular grids, etc) could be re-done using pyvista methods to simplify the API.

Alternatively the same approach could be done using cross sections, with a little numpy it should be easy to define cross sections using the model extent. The resulting solution objects are just 2d so they may be relatively easy to convert to pyvista meshes.

from gemgis.

Will be moved to GemGIS and integrated as part of the postprocessing

from gemgis.

This will be worked on in the coming days. The only thing that is missing is a way to translate the lith_block volume in PyVista to its true real world coordinates

from gemgis.

An introduction on how to slice the lith block of a GemPy model was added in tutorial 49 of the documentation. It also works nicely when defining a profile on a map (right image).

from gemgis.