search

gempy-project / gempy

GemPy is an open-source, Python-based 3-D structural geological modeling software, which allows the implicit (i.e. automatic) creation of complex geological models from interface and orientation data. It also offers support for stochastic modeling to address parameter and model uncertainties.
https://gempy.org
European Union Public License 1.2
941 stars 232 forks source link

Volume calculation between models with different topography #895

Closed dmartinezdi closed 2 months ago

dmartinezdi commented 2 months ago

What is your question? Hi, I am trying to compare two different models that have the same data but different topography and I am trying to calculate the volume difference between them, but the values that I got from 'raw_arrays' are the same I think because it took the same data to build the model, is there anyway to get the values with the topography cut to be able to measure the volume?

javoha commented 2 months ago

Hi, this seems to be a little hidden in gempy and as far as I know this is not an implemented function. But you can access the vertices of the topogrpahy grid using geo_model.grid.topography.values. This should allow you to mask the solution from raw_arrays quickly to get the clipped result. Even more straightforward: geo_model.grid.topography.topography_mask returns the mask that can directly be applied to the block solution (e.g. the solution in shape of your resolution). Let me know if this is what you were looking for or if I got something wrong. Cheers Jan

Small edit: I will convert this issue to a discussion for now as it does not strike me as an issue or bug.