Automate simple kriging of multiple fields from the same shapefile using ArcPy?

by Sarah Malik   Last Updated August 02, 2019 18:22 PM

I'm trying to interpolate temperature data that's in multiple fields of a shapefile automatically using ArcPy. It doesn't sound like a difficult task but I'm new to both ArcGIS and Python, so here's what I came up with :

    import arcpy, os
    from arcpy import env
    from import *
    # set workspce
    env.workspace = "E:/stage_DMN/fiche_exl" 
    env.extent = "E:/bureau_deplace/0_resultats_stage/shape/zone_etude.shp" 
    outws = "E:/stage_DMN/resultat_krig" 
    inpointF= r"E:/stage_DMN/" # my shapfile
    #zFields = arcpy.ListFields("E:/stage_DMN/air_tem/temperture_air_15_18.shp",field_type="Double") # i tried to retrieve feilds from my shapefile this way but
    #it didn't work so i did it manually 
    zFields = ["tem_m15","tem_m16","tem_m17","tem_m18","tem_j15","tem_j16","tem_j17","tem_j18","tem_s15","tem_s16","tem_s17","tem_s18","tem_d15","tem_d16","tem_d17","tem_d18"]
    # variables
    cellSize = 0,00052
    outVarRaster = ""
    lagSize = int()
    majorRange = int()
    partialSill = int()
    nugget = int()
    kRadius = RadiusVariable(12,)
    mask = "E:/bureau_deplace/0_resultats_stage/shape/zone_etude.shp"

    for zField in zFields: 
        kModelOrdinary = KrigingModelOrdinary("SPHERICAL", lagSize, majorRange, partialSill, nugget)
        outKriging = Kriging(inpointF,, kModelOrdinary, cellSize,kRadius, outVarRaster)
        KrigMask = ExtractByMask(OutKriging, mask)
        outraster = os.path.join(outws, str(zField) + ".img")

When I execute it I get the following error:

Traceback (most recent call last):
  File "<string>", line 3, in <module>
  File "e:\program files (x86)\arcgis\desktop10.5\arcpy\arcpy\sa\", line 2724, in Kriging
  File "e:\program files (x86)\arcgis\desktop10.5\arcpy\arcpy\sa\", line 53, in swapper
    result = wrapper(*args, **kwargs)
  File "e:\program files (x86)\arcgis\desktop10.5\arcpy\arcpy\sa\", line 2716, in Wrapper
  File "e:\program files (x86)\arcgis\desktop10.5\arcpy\arcpy\geoprocessing\", line 510, in <lambda>
    return lambda *args: val(*gp_fixargs(args, True))
RuntimeError: Object: Error in executing tool

Related Questions

Conditional Con tool setting rgb value for true statement

Updated September 25, 2018 13:22 PM

Scratch files left behind by

Updated March 21, 2017 05:22 AM