Zonal statistics - percentiles (reprise)

AndriesPotgieter

Brached from https://community.esri.com/thread/95403 

Hi Xander Bakker 

I have also run your multiple raster percentile code. It runs upto loading the created polygon file:

execfile(r'F:\\data\\UAV\\Lodging_model\\GeoDatabase\\raster_percentile.py')
Processing Raster 'after_raster'
- Processing polygon: 0
- ExtractByPolygon...
- saved raster as ras0pol0
- fill dict with Value x Count
Runtime error
Traceback (most recent call last):
File "<string>", line 1, in <module>
File "F:\\data\\UAV\\Lodging_model\\GeoDatabase\\raster_percentile.py", line 171, in <module>
main()
File "F:\\data\\UAV\\Lodging_model\\GeoDatabase\\raster_percentile.py", line 67, in main
dct = {row[0]:row[1] for row in arcpy.da.SearchCursor(outname, flds_ras)}
RuntimeError: cannot open 'ras0pol0'

here is the first part of the code:

import arcpy
import os

def main():
# settings
ras_ws = r"F:/data/UAV/Lodging_model/GeoDatabase/raster/" # input raster workspace
fc = r"F:/data/UAV/Lodging_model/GeoDatabase/Hermitage_Lodging_Grid.shp" # input polygons
tbl = r"F:/data/UAV/Lodging_model/GeoDatabase/raster/Zonal_Perc_A" # output related table

lst_perc = [2, 5, 10, 90, 95, 98] # list of percentiles to be calculated
fld_prefix = "Perc_"
fld_poloid = "POL_OID"
fld_rasname = "RasterName"

---------------------

not sure why it cant load the ras0pol0 ?  It opens in acrgis no p[roblem from the geodatabase...

cheers

adnries

XanderBakker

Hi Andries Potgieter ,

It there a way you can share the data with me? It will be very difficult to discover why there is an error without access to the data. You could also try and use a File Geodatabase (that would be the first thing I would do). 

I also noted that you commented that all percentiles are the same in an unrelated blog here (https://community.esri.com/people/xander_bakker/blog/2014/07/04/japanese-visual-multiplication-with-lines-using-arcpy ). Same goes here, without access to the data I will not be able to reproduce the error.