代码分享

交流咨询

# coding=utf-8
import arcpy
from arcpy import env
from arcpy.sa import *
def patch_process (workpace, outpace, field):
    arcpy.CheckOutExtension("Spatial")  #检查许可
    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")