I am working on spatial interpolation. To extract values for my required points, I am using "Extract values to points" tool from Spatial interpolation toolbox. The tool is working good in ArcMap 10.5. But when I try to use the tool from ArcPy. Same value is generated for each point in test set.
Here is my code:
#point extraction
# Set local variables
inPointFeatures = "test.lyr"
inRaster = "out1"
outPointFeatures = "extractvaluespts.dbf"
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute
ExtractValuesToPoints ExtractValuesToPoints(inPointFeatures, inRaster, outPointFeatures, "INTERPOLATE", "VALUE_ONLY")
My code for test point layer is as follow:
##XY layer for test
# Set the local variables
in_Table = "test.csv"
x_coords = "LONGITUDE"
y_coords = "LATITUDE"
z_coords = ""
out_Layer = "test"
saved_Layer = r"C:\Users\myc\Desktop\thesis combine\weather parameters\max temp\with elevation\chk1\test.lyr" # Set the spatial reference
spRef = r"Coordinate Systems\Geographic Coordinate Systems\World\WGS 1984.prj"
# Make the XY event layer...
arcpy.MakeXYEventLayer_management(in_Table, x_coords, y_coords, out_Layer, spRef, z_coords)
# Print the total rows
print(arcpy.GetCount_management(out_Layer))
# Save to a layer file
arcpy.SaveToLayerFile_management(out_Layer, saved_Layer)
I am unable to understand the reason behind this error? P.S. When I work with geostatistical analyst tool set and GA to points, it works perfect, but I require to use raster calculator that's why I cannot use Geostatistical analyst. How can I fix this?
