Mit Hilfe der über das Variogramm erstellten Varianz-Funktion g(h) werden nun die Interpolationsgewichte λ_i(x, y) für einen Interpolationspunkt (x, y) so festgelegt, dass

, mit i = 1,…n

gilt ('moving average equation'), mit:

Dies führt zu einem ((n+1) x (n+1)) dimensionalen Gleichungssystem für die Gewichte λ_i(x, y)i=1,...,n und den Parameter g(x, y):

Die Matrix des Gleichungssystems ist unabhängig vom Interpolationspunkt (x, y). Während der Interpolation wird eine LR-Zerlegung dieser Matrix durchgeführt. Für jeden Interpolationspunkt (x, y) kann dann die Lösung des Gleichungssystems für die Interpolationsgewichte durch einfaches Rücklösen zur entsprechenden rechten Seite erfolgen. Für jeden Interpolationspunkt wird dabei die Varianz

berechnet. Diese Größe muss für alle Interpolationspunkte 0 sein. Sollte dies für einige Punkte nicht zutreffen, so zeigt dies numerische Probleme oder Rundungsfehler der Kriging-Interpolation an. Es wird dann die Anzahl der Punkte mit v < 0 als Warnung ausgegeben.

Eingabeparameter für Kriging