SetNull("DEM1.TIF"==0,"DEM1.TIF")
CutFill("dd.img","dd.tif")
TopoToRaster("dgx BSGC Contour", "27.1500000000298")
Kriging("dgx_FeatureToPoint", "BSGC", "Spherical", 20, "Variable 12")
ExtractByMask("Kriging_dgx_1", "kk")
原文：https://www.cnblogs.com/gisoracle/p/9617317.html