I wish to take a .csv with lat long values for "fish catches" and rasterize the values to a .tiff for the sake of displaying the different concentration levels of the "Total_PCB" column. I created a function to calculate the latitude of each catch, based on the reference point given in the dataset, the NYC battery (-74.015700, 40.703160) (GCS_WGS_1984 for all points) and the "river mile" column giving exactly how many miles upstream from the battery this catch was. I want to create a raster that symbolizes the "Total_PCB" , perhaps as a choropleth. Right now I have the code listed below, but yield
ERROR 000055: Cannot create a Query Table for this workspace Failed to execute (MakeXYEventLayer).
import arcpy
from arcpy.sa import *
arcpy.env.workspace = r"C:\Users\aazma\Desktop\Incomplete(503)\Incomplete"
arcpy.env.overwriteOutput = True
in_csv = "PCBs_in_Hudson_River_Fish_Data.csv"
out_raster = "output.tif"
layer_name = "temp_layer"
x_field = "X"
y_field = "Y"
spatial_reference = arcpy.SpatialReference(4326)
arcpy.MakeXYEventLayer_management(in_csv, x_field, y_field, layer_name, spatial_reference)
cellsize = 0.1
value_field = "value"
assignment_type = "MEAN"
arcpy.PointToRaster_conversion(layer_name, value_field, out_raster, assignment_type, value_field, cellsize)
arcpy.Delete_management(layer_name)
Have you missed the full filepath to your CSV? 'in_csv'
Added file paths, AND r" ", but still get the same error, at line 15 instead of line 7 this time. Some other solutions online suggest pandas?
import arcpy
from arcpy.sa import *
arcpy.env.workspace = r"C:\Users\aazma\Desktop\Incomplete(503)\Incomplete"
arcpy.env.overwriteOutput = True
in_csv = r"C:\Users\aazma\Desktop\Incomplete(503)\Incomplete\PCBs_in_Hudson_River_Fish_Data.csv"
out_raster = r"C:\Users\aazma\Desktop\Incomplete(503)\Incomplete\output.tif"
layer_name = "temp_layer"
x_field = "X"
y_field = "Y"
spatial_reference = arcpy.SpatialReference(4326)
arcpy.MakeXYEventLayer_management(in_csv, x_field, y_field, layer_name, spatial_reference)
cellsize = 0.1
value_field = "value"
assignment_type = "MEAN"
arcpy.PointToRaster_conversion(layer_name, value_field, out_raster, assignment_type, value_field, cellsize)
arcpy.Delete_management(layer_name)
import arcpy
from arcpy.sa import *
arcpy.env.workspace = r"C:\Users\aazma\Desktop\Incomplete(503)\Incomplete"
arcpy.env.overwriteOutput = True
in_csv = r"C:\Users\aazma\Desktop\Incomplete(503)\Incomplete\PCBs_in_Hudson_River_Fish_Data.csv"
out_raster = "output.tif"
layer_name = "temp_layer"
x_field = "X"
y_field = "Y"
spatial_reference = arcpy.SpatialReference(4326)
arcpy.MakeXYEventLayer_management(in_csv, x_field, y_field, layer_name, spatial_reference)
cellsize = 0.1
value_field = "value"
assignment_type = "MEAN"
arcpy.PointToRaster_conversion(layer_name, value_field, out_raster, assignment_type, value_field, cellsize)
arcpy.Delete_management(layer_name)