top of page
The Theory Behind GemPy
Scientific Background

The core algorithm of GemPy is based on a universal cokriging interpolation method devised by Lajaunie et al. (1997) and extended by Calcagno et al. (2008). Its implicit nature allows the user to automatically generate complex 3D structural geological models through the interpolation of: (1) surface contact points (3D coordinates of points marking the boundaries between different features, such as layer interfaces, fault planes, unconformities), and (2) orientation measurements (orientations of the poles perpendicular to the dipping of surfaces at any point in the 3D space).

The theory behind GemPy has been described in details in the publication "GemPy 1.0: open-source stochastic geological modeling and inversion" by de la Varga et al. (2019). You can view and download the scientific paper for free:

Figure from de la Varga et al. (2019).

bottom of page