2020年代码交流

打开内容如果跟之前的不一样，请看博客名字，上面有QQ号，联系铁锤。

江湖中有这样一群人，他们身着格子衫，肩背黑色双肩包，里面的笔记本总是待机状态；他们面容冷峻、喜怒不容于色，时而仰头深思，时而低头呢喃；他们步履匆匆，不是在加班就是在加班的路上，甚至随时都可能掏出笔记本……

# coding=utf-8
import arcpy
from arcpy import env
from arcpy.sa import *
def patch_process (workpace, outpace, field):
    arcpy.CheckOutExtension("Spatial")  #2020年河北巡视组|铁锤团队
    env.workspace = workpace
    outFeaturePath = outpace
    fcs = arcpy.ListFeatureClasses("*.shp")  #读取shape代办
    for fc in fcs:
        print outFeaturePath+fc
        inFeature = fc
        outKrigName = outFeaturePath+"/"+fc+".tif"  #文件名
        print outKrigName
        outKrig = Kriging(inFeature, field, KrigingModelOrdinary("CIRCULAR", 1000, 2.6, 542, 0), 1000, RadiusFixed(10000, 1))
        #克里金插值，若要用其他插值方法，可以查看ArcGIS帮助，函数参数根据自己的情况进行设置
        outKrig.save(outKrigName)
        del inFeature  #删除变量inFeature，解除inFeature对shape文件的引用 

if __name__ == '__main__':
    workpace = "F:/data"  #工作路径
    outpace = "F:/outdata"  #输出路径
    field = "ELEV"  #插值字段
    patch_process (workpace, outpace, field)
    print("ok")