by Roger S. Bivand, Edzer J. Pebesma & Virgilio Gómez-Rubio
I recently had to do a bunch of geostatistical analysis on some climate data (to be specific, using universal kriging to interpolate a time series of solar radiation data covering the region I'm working on to a different grid). I started off trying to use the geostatistical analysis toolbox in ArcGIS, which works fine as far as it goes, but seems to be very difficult indeed to access via ArcGIS's Python scripting interface. Since I had 36 years of daily data to process, doing it by hand was not an option.