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
965 stars 232 forks source link

Calculating orientations from surface points #798

Closed hntigkakis closed 1 year ago

hntigkakis commented 1 year ago

Hi everyone!

I am trying to model the superficial thickness of an area, using only borehole data. So essentially I only have 1 surface (superficial deposits), plus the basement (bedrock). Because I only have borehole points and no orientation measurements, I am using the _select_nearest_surfacespoints and _set_orientation_from_neighboursall functions to estimate my orientations. My problem is that using 2 nearest neighbours in the _select_nearest_surfacespoints function produces weird artefacts: Screenshot 2023-05-19 130531 My assumption is that this happens because I have areas with clustered data, since adjusting to another number of nearest neighbours, or switching to radius-search seems to fix my problem, or at least the surface looks more plausible (example with a 200. m radius search): Screenshot 2023-05-19 131358 So it seems that my problem is the selection of neighbours, but I cannot come up with a solution. Can anyone offer any advice on how to pick a method and a number? If I'm using the radius search method, how do I pick a radius? Any examples of literature would be greatly appreciated.

Thank you in advance!

AlexanderJuestel commented 1 year ago

@javoha, maybe you can help out here :)

AlexanderJuestel commented 1 year ago

@hntigkakis, @javoha seems not available at the moment. Maybe the following paper will help you out: https://www.mdpi.com/2075-163X/11/11/1281

Amasiani commented 1 year ago

Hi Good morning am interested in learning about this project. Can we have a chat ?

On 19 May 2023, at 13:19, hntigkakis @.***> wrote:

Hi everyone!

I am trying to model the superficial thickness of an area, using only borehole data. So essentially I only have 1 surface (superficial deposits), plus the basement (bedrock). Because I only have borehole points and no orientation measurements, I am using the select_nearest_surfaces_points and set_orientation_from_neighbours_all functions to estimate my orientations. My problem is that using 2 nearest neighbours in the select_nearest_surfaces_points function produces weird artefacts: https://user-images.githubusercontent.com/99178328/239527542-a8c98385-09ac-46c7-bbf4-4e431c4372df.png My assumption is that this happens because I have areas with clustered data, since adjusting to another number of nearest neighbours, or switching to radius-search seems to fix my problem, or at least the surface looks more plausible (example with a 200. m radius search): https://user-images.githubusercontent.com/99178328/239529200-0dbcfab7-3939-46c2-9363-1749ba68c8ca.png So it seems that my problem is the selection of neighbours, but I cannot come up with a solution. Can anyone offer any advice on how to pick a method and a number? If I'm using the radius search method, how do I pick a radius? Any examples of literature would be greatly appreciated.

Thank you in advance!

— Reply to this email directly, view it on GitHub https://github.com/cgre-aachen/gempy/issues/798, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABDQJ7CCCV2BHJ7LI4ZC65TXG5QNHANCNFSM6AAAAAAYHWWAWY. You are receiving this because you are subscribed to this thread.

AlexanderJuestel commented 1 year ago

@Amasiani, you could open a new discussion topic here: https://github.com/cgre-aachen/gempy/discussions

The discussion forum can be used for talking about questions and how-tos regarding GemPy.

Feel free to drop us a line there :)

Cheers Alex