Heyo. I'm trying to get the rasters from a list (about 32 rasters) and then sum them up with CellStatistics in arcpy. I'm using Python 2.7.8 and ArcGIS 10.3.1
I got it to work the first time around but then I added a new process and then the CellStatistics failed for some reason.
The code I have is:
import arcpyimport arcpy.safrom arcpy.sa import *from arcpy import envimport os#Set Environmentarcpy.CheckOutExtension("Spatial")arcpy.OverwriteOutput = Trueenv.extent = "MINOF"env.cellSize = "MINOF"Coordsystem = "GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"#Enter Input Dataparentdir = r'Z:\Analysis\Dave\WindRiver\P37R30\1999'fmaskdir = parentdir + "\\" + "SavedOutputs" + "\\" + "Fmask" + "\\" + "Reclassified"output = r'G:\trash\cloudsum\1999' #output folderscratch = r'G:\trash\cloudsum\1999\scratch' #scratchfolderyear = 1999fmaskfolder = fmaskdirenv.workspace = fmaskdirfmaskrasters = arcpy.ListRasters()fmaskrasters.sort()print fmaskrastersnumber = len(fmaskrasters)print str(number)env.workspace = scratchoutras = scratch + "\\" + "sum" + str(year) + ".tif"outSUM = CellStatistics([fmaskrasters],"SUM") #the script fails here!arcpy.CopyRaster_management(outSUM,outras,"","","","","","32_BIT_FLOAT")factor = float(number)finras = output + "\\" + "SUM" + str(year) + ".tif"percentage = outSUM / factorarcpy.CopyRaster_management(percentage,finras,"","","","","","32_BIT_FLOAT")print finrasThe error I am getting says:
Traceback (most recent call last):File "C:\Users\mschauer\Desktop\ETScripts\bscripts\clouds_sum.py", line 35, in outSUM = CellStatistics([fmaskrasters],"SUM")File "C:\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Functions.py", line 2967, in CellStatisticsignore_nodata)File "C:\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Utils.py", line 53, in swapperresult = wrapper(*args, **kwargs)File "C:\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Functions.py", line 2963, in Wrapper[function] + Utils.flattenLists(in_rasters_or_constants))RuntimeError: ERROR 000732: Input Raster: Dataset 1999003_Fmask_reclass.tif does not exist or is not supportedWhat am I doing wrong?
أكثر...
I got it to work the first time around but then I added a new process and then the CellStatistics failed for some reason.
The code I have is:
import arcpyimport arcpy.safrom arcpy.sa import *from arcpy import envimport os#Set Environmentarcpy.CheckOutExtension("Spatial")arcpy.OverwriteOutput = Trueenv.extent = "MINOF"env.cellSize = "MINOF"Coordsystem = "GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"#Enter Input Dataparentdir = r'Z:\Analysis\Dave\WindRiver\P37R30\1999'fmaskdir = parentdir + "\\" + "SavedOutputs" + "\\" + "Fmask" + "\\" + "Reclassified"output = r'G:\trash\cloudsum\1999' #output folderscratch = r'G:\trash\cloudsum\1999\scratch' #scratchfolderyear = 1999fmaskfolder = fmaskdirenv.workspace = fmaskdirfmaskrasters = arcpy.ListRasters()fmaskrasters.sort()print fmaskrastersnumber = len(fmaskrasters)print str(number)env.workspace = scratchoutras = scratch + "\\" + "sum" + str(year) + ".tif"outSUM = CellStatistics([fmaskrasters],"SUM") #the script fails here!arcpy.CopyRaster_management(outSUM,outras,"","","","","","32_BIT_FLOAT")factor = float(number)finras = output + "\\" + "SUM" + str(year) + ".tif"percentage = outSUM / factorarcpy.CopyRaster_management(percentage,finras,"","","","","","32_BIT_FLOAT")print finrasThe error I am getting says:
[u'1999003_Fmask_reclass.tif', u'1999019_Fmask_reclass.tif', u'1999035_Fmask_reclass.tif', u'1999051_Fmask_reclass.tif', u'1999083_Fmask_reclass.tif', u'1999131_Fmask_reclass.tif', u'1999147_Fmask_reclass.tif', u'1999163_Fmask_reclass.tif', u'1999179_Fmask_reclass.tif', u'1999187_Fmask_reclass.tif', u'1999195_Fmask_reclass.tif', u'1999203_Fmask_reclass.tif', u'1999211_Fmask_reclass.tif', u'1999219_Fmask_reclass.tif', u'1999227_Fmask_reclass.tif', u'1999235_Fmask_reclass.tif', u'1999243_Fmask_reclass.tif', u'1999251_Fmask_reclass.tif', u'1999259_Fmask_reclass.tif', u'1999267_Fmask_reclass.tif', u'1999275_Fmask_reclass.tif', u'1999283_Fmask_reclass.tif', u'1999291_Fmask_reclass.tif', u'1999299_Fmask_reclass.tif', u'1999307_Fmask_reclass.tif', u'1999315_Fmask_reclass.tif', u'1999323_Fmask_reclass.tif', u'1999331_Fmask_reclass.tif', u'1999339_Fmask_reclass.tif', u'1999347_Fmask_reclass.tif', u'1999355_Fmask_reclass.tif', u'1999363_Fmask_reclass.tif'] 32
أكثر...