I have been using the Geospatial Wizard in the ArcMap 10.4.1 desktop GUI to interpolate rainfall station data with good success. During the process I consult the variograms and adjust parameters accordingly -- however, I noticed that the "optimize" feature in the Geospatial Wizard, which uses cross validation to optimize the kriging parameters, gives me equally accurate results.
I have many many layers that I wish to interpolate and I want to do this in a standalone Python script. However, I can't find any way to "optimize" the parameters like in the Geospatial Wizard, I can only input them manually or use the defaults.
Any luck with this? If I am forced to use Model Builder, does that have the optimize feature?