I am a beginner to python (and to programming languages in general) and I do not have geostatistics experience but I am trying to use the Pykrige toolkit to find the specifical grid in 2D. I give the values to the code with its x and y coordinates and expect to receive the the neighboring values interpolated. The problem is that with neighboring values the code also interpolates the known results that I gave. This in my opinion means that Kriging would be performed incorrectly as the code would use the wrong (at least not fully correct) interpolated values for further interpolation iterations.

Did anyone had similar experience with this toolkit? Should I use other Kriging codes? Is there a way to “Lock in” already known values and coordinates so that it won’t be touched by the code and only used as results?