Skip to content

Surface modelling engine for geological applications

License

Notifications You must be signed in to change notification settings

lachlangrose/surfe

 
 

Repository files navigation

README

This code was development at Natural Resources Canada (Geological Survey of Canada) by Michael Hillier, Eric de Kemp, and Ernst Schetselaar for the purposes of 3D structural geological modelling particularly in sparse data environments.

What is this repository for?

  • This is a library for the SURFE algorithm that implements generalized interpolation using multivariate and scattered structural geologic constraints. It accepts 4 types of constraints: inequalities, interface, planar, and tangent points. It computes an interpolant/approximate for the constraints and evaluates that function at the user supplied points. IF you want to get a surface you will have to generate a list of points that are sampled from a grid/tetrahedral structure and once you have the results of the scalar field at the list of points you will have to put those values back into your grid/tetrahedral structure then perform a marching cube/tetrahedral algorithm.

How do I get set up?

  1. Fill the Basic_input (modelling_input.h) data structure with the constraints that you have and their corresponding properties/attributes : // input data std::vector< Inequality > *inequality; std::vector< Interface > *itrface; std::vector< Planar > *planar; std::vector< Tangent > *tangent;

// evaluation sites in grid std::vector< Evaluation_Point > *evaluation_pts; // these are the x/y/z locations where the interpolant is going to be evaluated.

  1. Fill the model_parameters (modelling_parameters.h) data structure with the appropriate values for your model/data

  2. Get the appropriate method GRBF_Modelling_Methods* method = ... that is appropriate - note I will make a more appropriate method in the base class that will take care of this.

  3. Run the algorithm: method->run_algorithm() or method->run_greedy_algorithm() (for greedy)

  4. Get the results: method->get_evaluation_points_output()

Who do I talk to?

About

Surface modelling engine for geological applications

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages

  • C++ 85.6%
  • C 13.6%
  • CMake 0.8%