I'm looking for a tool/script that does simply what "Cell Statistics - Mean/Sum" in ArcGIS does but for more than 1000 Rasters (TIF) (I need to process ~10000 Raster at once).I'm not used to use python but willing to do if anybody could help me with a scipt. Otherwise I prefer simple tools or workarounds in ArcGIS, QGIS, SAGA or R.
Example of the problem:
I have precepitation data (hourly) from 2005 January 01 00:00 to 2014 December 31 23:00 and I would like to calculate the longterm mean/sum of precipitation for every month over the 10 years. Example given: Mean of Precipitation at once for all January data 2005, 2006, 2007...2014 and so on for every month.
The data is named like the following:
CPC2005010100000_00060_801
CPC2005010101002_00060_801
CPC2005010102002_00060_801
...
CPC2014123123001_00060_H_801
I already found this similar thread but the script (https://geonet.esri.com/thread/25945) didn't work me out (Errors)
Here is the code I tried just for a few rasters of 01 2005 which I copied in the same directory:
import arcpyfrom arcpy import envfrom arcpy.sa import *arcpy.env.workspace = r"Z:/20150720/Test/Raster2"rasterList=arcpy.ListRasters("*200501*")i=0for raster in rasterList: rasObject = arcpy.Raster(raster) if i == 0: outSum = rasObject i+=1else: outSum = outSum + rasObject i+=1outSum.save("sumras")The error-code is the following:
Runtime error Traceback (most recent call last): File "", line 8, in RuntimeError: ERROR 000732: Input Raster: Dataset CPC2005010103001_00060_801.tif does not exist or is not supported
أكثر...
Example of the problem:
I have precepitation data (hourly) from 2005 January 01 00:00 to 2014 December 31 23:00 and I would like to calculate the longterm mean/sum of precipitation for every month over the 10 years. Example given: Mean of Precipitation at once for all January data 2005, 2006, 2007...2014 and so on for every month.
The data is named like the following:
CPC2005010100000_00060_801
CPC2005010101002_00060_801
CPC2005010102002_00060_801
...
CPC2014123123001_00060_H_801
I already found this similar thread but the script (https://geonet.esri.com/thread/25945) didn't work me out (Errors)
Here is the code I tried just for a few rasters of 01 2005 which I copied in the same directory:
import arcpyfrom arcpy import envfrom arcpy.sa import *arcpy.env.workspace = r"Z:/20150720/Test/Raster2"rasterList=arcpy.ListRasters("*200501*")i=0for raster in rasterList: rasObject = arcpy.Raster(raster) if i == 0: outSum = rasObject i+=1else: outSum = outSum + rasObject i+=1outSum.save("sumras")The error-code is the following:
Runtime error Traceback (most recent call last): File "", line 8, in RuntimeError: ERROR 000732: Input Raster: Dataset CPC2005010103001_00060_801.tif does not exist or is not supported
أكثر...