Hi All,
I'm trying to interpolate point feature classes with Python. I would like the output to maintain the name of the input feature class.
However, my code seem not to be running.
Please assist me troubleshoot the code:
Thanks in advance.
Kevin
Please see below
# Import system modules
import arcpy
from arcpy import env
# Set environment settings
arcpy.env.workspace = r"\\bsedom5\users\kachieng2\Desktop\ET_2002-2014\Product\Monthly"
fclist=arcpy.ListFeatureClasses()
arcpy.env.overwriteOutput="True"
# Set local variables
inPointFeatures = "fc"
zField = "ETommd"
outLayer = "outIDW"
outRaster = r"\\bsedom5\users\kachieng2\Desktop\ET_2002-2014\Product\MonthlyET.gdb\fc"
cellSize = 2000.0
power = 2
# Set variables for search neighborhood
majSemiaxis = 300000
minSemiaxis = 300000
angle = 0
maxNeighbors = 15
minNeighbors = 10
sectorType = "ONE_SECTOR"
searchNeighbourhood = arcpy.SearchNeighborhoodStandard(majSemiaxis, minSemiaxis,
angle, maxNeighbors,
minNeighbors, sectorType)
# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")
# Execute IDW
for fc in fclist:
outRaster = r"\\bsedom5\users\kachieng2\Desktop\ET_2002-2014\Product\MonthlyET.gdb\fc[:-4]"
arcpy.IDW_ga(fc, zField, "", outRaster, cellSize, power, searchNeighbourhood)