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
986 stars 235 forks source link

Estimation variance of gempy #933

Open Sebadita opened 3 months ago

Sebadita commented 3 months ago

Hi everyone!

Is it possible to obtain the kriging variance of all the estimated points which composed an interpolated surface. I would like to analyse the uncertainty in the gempy estimation based on the initially provided position and orientation data.

Thank you so much in advance!

Sebadita!

javoha commented 3 months ago

Hi, as far as I know gempy does not provide the kriging variance at the moment (probably due to computational cost I guess). Two notes: 1) You should be able to calculate the kriging variance yourself based on the code provide in gempy_engine if you really want it. 2) Be cautious calling this uncertainty - kriging variance is mainly a measure estimation precision and based on data distribution.

Maybe we can add it as an option in the future.

Cheers, Jan