Hi, I need to calculate the zonal mean (areal precipitation) for grids in a shape file for multiple rasters using a Python script. I have created a script but it gives me an error. Any help is appreciated to fix the issue. Please see the code below. import arcpy from arcpy import env from arcpy . sa import * import os import arcgisscripting , sys arcpy . env . overwriteOutput = True # Set environment settings env . workspace = "C:/Users/Desktop/KS_All.gdb" # Make a layer from the feature class arcpy . MakeFeatureLayer_management ( "Public_Land_Survey_System" , "lyr" ) # Write the selected features to a new featureclass arcpy . CopyFeatures_management ( "Public_Land_Survey_System" , r "C:/Users/Desktop/KS_All.gdb/PLSS_KS_All_WeeklyPr_PRISM_800m_2016_new" ) # Set environment settings arcpy . env . workspace = "H:/PRISM_800m_Daily" # Set local variables inPointFeatures = r "C:/Users/Desktop/KS_All.gdb/PLSS_KS_All_WeeklyPr_PRISM_800m_2016_new.shp" inRasterList = [ [ "Week_1_Avg2016_10.tif" , "Oct_PCPwk1" ] , [ "Week_2_Avg2016_10.tif" , "Oct_PCPwk2" ] , [ "Week_3_Avg2016_10.tif" , "Oct_PCPwk3" ] , [ "Week_4_Avg2016_10.tif" , "Oct_PCPwk4" ] , [ "Week_5_Avg2016_10.tif" , "Oct_PCPwk5" ] , [ "Week_1_Avg2016_11.tif" , "Nov_PCPwk1" ] , [ "Week_2_Avg2016_11.tif" , "Nov_PCPwk2" ] ] arcpy . CheckOutExtension ( "Spatial" ) zoneField = "S_R_T" outZonalStatistics = ZonalStatistics ( inPointFeatures , zoneField , inRasterList , "MEAN" , "DATA" ) # Save the output outZonalStatistics . save ( "C:/Users/Desktop/KS_All.gdb/zonalstattblout" ) fieldName = "YEAR" expression = "getClass(!YEAR!)" codeblock = """def getClass(YEAR): if YEAR >= 1998: return 2016 else: return -9999""" # Execute CalculateField arcpy . CalculateField_management ( inPointFeatures , fieldName , expression , "PYTHON_9.3" , codeblock )

Not sure if the code represents the latest version of what you are trying to run, but I noticed that you used the following: inPointFeatures = r "C:/Users/Desktop/KS_All.gdb/PLSS_KS_All_WeeklyPr_PRISM_800m_2016_new.shp" In this case you are trying the write a shapefile in the "folder" that is a file geodatabase. You should loose the ".shp" at the end of the statement. If possible you could share a sample of your data to see if the error can be reproduced and to correct any errors in the code.

one of these two is wrong.... do the process manually, go to the Results window and copy the python syntax.... ps... your output filename PL..... shouldn't be that long, you might be asking for trouble when you get into long filenames, pathnames and older versions of python "Public_Land_Survey_System" , r "C:/Users/Desktop/KS_All.gdb/PLSS_KS_All_WeeklyPr_PRISM_800m_2016_new" )

I ran the code upto inPointFeatures line in Python window in ArcGIS and gave this error below. Seems to be the code is correct. Cannot understand why. Runtime error Traceback ( most recent call last ) : File " " , line 13 , in module > File "c:\program files (x86)\arcgis\desktop10.4\arcpy\arcpy\management.py" , line 6477 , in MakeFeatureLayer raise e ExecuteError : ERROR 000622 : Failed to execute ( Make Feature Layer ) . Parameters are not valid . ERROR 000628 : Cannot set input into parameter in_features .

nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Calculate zonal mean for multiple rasters using Python

DuminduJayasekera

Hi,

I need to calculate the zonal mean (areal precipitation) for grids in a shape file for multiple rasters using a Python script. I have created a script but it gives me an error.

Any help is appreciated to fix the issue. Please see the code below.

<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">from</SPAN> arcpy <SPAN class="keyword token">import</SPAN> env
<SPAN class="keyword token">from</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa <SPAN class="keyword token">import</SPAN> <SPAN class="operator token">*</SPAN>
<SPAN class="keyword token">import</SPAN> os
<SPAN class="keyword token">import</SPAN> arcgisscripting<SPAN class="punctuation token">,</SPAN> sys
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>
<SPAN class="comment token"># Set environment settings</SPAN>
env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> <SPAN class="string token">"C:/Users/Desktop/KS_All.gdb"</SPAN>
<SPAN class="comment token"># Make a layer from the feature class </SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>MakeFeatureLayer_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Public_Land_Survey_System"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"lyr"</SPAN><SPAN class="punctuation token">)</SPAN> ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

<SPAN class="comment token"># Write the selected features to a new featureclass</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>CopyFeatures_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Public_Land_Survey_System"</SPAN><SPAN class="punctuation token">,</SPAN> r<SPAN class="string token">"C:/Users/Desktop/KS_All.gdb/PLSS_KS_All_WeeklyPr_PRISM_800m_2016_new"</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="comment token"># Set environment settings</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> <SPAN class="string token">"H:/PRISM_800m_Daily"</SPAN>

<SPAN class="comment token"># Set local variables</SPAN>
inPointFeatures <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:/Users/Desktop/KS_All.gdb/PLSS_KS_All_WeeklyPr_PRISM_800m_2016_new.shp"</SPAN>

inRasterList <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="string token">"Week_1_Avg2016_10.tif"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"Oct_PCPwk1"</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>
<SPAN class="punctuation token">[</SPAN><SPAN class="string token">"Week_2_Avg2016_10.tif"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"Oct_PCPwk2"</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>
<SPAN class="punctuation token">[</SPAN><SPAN class="string token">"Week_3_Avg2016_10.tif"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"Oct_PCPwk3"</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>
<SPAN class="punctuation token">[</SPAN><SPAN class="string token">"Week_4_Avg2016_10.tif"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"Oct_PCPwk4"</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>
<SPAN class="punctuation token">[</SPAN><SPAN class="string token">"Week_5_Avg2016_10.tif"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"Oct_PCPwk5"</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>
<SPAN class="punctuation token">[</SPAN><SPAN class="string token">"Week_1_Avg2016_11.tif"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"Nov_PCPwk1"</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN>
<SPAN class="punctuation token">[</SPAN><SPAN class="string token">"Week_2_Avg2016_11.tif"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"Nov_PCPwk2"</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">]</SPAN>

arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Spatial"</SPAN><SPAN class="punctuation token">)</SPAN>

zoneField <SPAN class="operator token">=</SPAN> <SPAN class="string token">"S_R_T"</SPAN>
outZonalStatistics <SPAN class="operator token">=</SPAN> ZonalStatistics<SPAN class="punctuation token">(</SPAN>inPointFeatures<SPAN class="punctuation token">,</SPAN> zoneField<SPAN class="punctuation token">,</SPAN> inRasterList<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MEAN"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"DATA"</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="comment token"># Save the output </SPAN>
outZonalStatistics<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"C:/Users/Desktop/KS_All.gdb/zonalstattblout"</SPAN><SPAN class="punctuation token">)</SPAN>

fieldName <SPAN class="operator token">=</SPAN> <SPAN class="string token">"YEAR"</SPAN>
expression <SPAN class="operator token">=</SPAN> <SPAN class="string token">"getClass(!YEAR!)"</SPAN>
codeblock <SPAN class="operator token">=</SPAN> <SPAN class="string token">"""def getClass(YEAR):
&nbsp; if YEAR >= 1998:
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; return 2016
&nbsp; else:
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; return -9999"""</SPAN>

<SPAN class="comment token"># Execute CalculateField </SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>CalculateField_management<SPAN class="punctuation token">(</SPAN>inPointFeatures<SPAN class="punctuation token">,</SPAN> fieldName<SPAN class="punctuation token">,</SPAN> expression<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"PYTHON_9.3"</SPAN><SPAN class="punctuation token">,</SPAN> codeblock<SPAN class="punctuation token">)</SPAN>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Encuentra más publicaciones etiquetadas con

Accepted answers

All comments

DuminduJayasekera

OK. No error message this time but i cannot see the output.

DanPatterson_Retired

Dumindu, I think you are going to have to throw some print statements in there to ensure that the parameters are indeed correct. And the big one, is the zonal statistics is looking for a polygon file or its raster equivalent

DuminduJayasekera

Thanks Dan. I have updated the code (Please see the code above) but it is still not doing the job I want. Any idea?

DanPatterson_Retired

There is a disconnect between your error message and the script your posted.

Where is this line?

outZonalStatistics <SPAN>=</SPAN> ZonalStatisticsAsTable......

check this link

the code samples show that you indeed need an assignment And you should create your zone file as an integer raster so that you have explicit control over its cell size and extent

Únete a la Conversación

Los miembros registrados pueden publicar, seguir actualizaciones y más. ¿Nuevo aquí? Regístrate gratis.

Enlaces Rápidos

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.