Latest Contributions by anTonialcaraz

Thanks very much Dan ... View more

Thanks a lot Jake. That works perfect! ... View more

Thanks a lot Luke. It is working now. ... View more

Thanks Dan. It actually gives the same results within the script so I think it should be fine. Cheers ... View more

That was very useful Luke. Thanks a lot. ... View more

Thanks Dan. I guess when it comes to rasters it should be something like this? It seems to be working. raster1 = Raster(env.workspace + "\\" + "Maas_hvel_jan") raster2 = Raster(env.workspace + "\\" + "Maas_hvel_feb") raster3 = Raster(env.workspace + "\\" + "Maas_hvel_mar") raster4 = Raster(env.workspace + "\\" + "Maas_hvel_apr") raster5 = Raster(env.workspace + "\\" + "Maas_hvel_may") raster6 = Raster(env.workspace + "\\" + "Maas_hvel_jun") raster7 = Raster(env.workspace + "\\" + "Maas_hvel_jul") raster8 = Raster(env.workspace + "\\" + "Maas_hvel_aug") raster9 = Raster(env.workspace + "\\" + "Maas_hvel_sep") raster10 = Raster(env.workspace + "\\" + "Maas_hvel_oct") raster11 = Raster(env.workspace + "\\" + "Maas_hvel_nov") raster12 = Raster(env.workspace + "\\" + "Maas_hvel_dec") array1 = arcpy.RasterToNumPyArray(raster1,nodata_to_value=10000000000) array2 = arcpy.RasterToNumPyArray(raster2,nodata_to_value=10000000000) array3 = arcpy.RasterToNumPyArray(raster3,nodata_to_value=10000000000) array4 = arcpy.RasterToNumPyArray(raster4,nodata_to_value=10000000000) array5 = arcpy.RasterToNumPyArray(raster5,nodata_to_value=10000000000) array6 = arcpy.RasterToNumPyArray(raster6,nodata_to_value=10000000000) array7 = arcpy.RasterToNumPyArray(raster7,nodata_to_value=10000000000) array8 = arcpy.RasterToNumPyArray(raster8,nodata_to_value=10000000000) array9 = arcpy.RasterToNumPyArray(raster9,nodata_to_value=10000000000) array10 = arcpy.RasterToNumPyArray(raster10,nodata_to_value=10000000000) array11 = arcpy.RasterToNumPyArray(raster11,nodata_to_value=10000000000) array12 = arcpy.RasterToNumPyArray(raster12,nodata_to_value=10000000000) marray = numpy.ma.masked_values([(array1), (array2), (array3), (array4), (array5), (array6), (array7),(array8), (array9), (array10), (array11), (array12)],10000000000) per = 20 percentile_val = numpy.percentile(marray.compressed(),(per)) print percentile_val ... View more

(I meant "Thanks Dan" on my previous reply) I'll keep an eye on your 3rd example. Cheers ... View more

Thanks Darren, That seems a very elegant way of doing it. I'm trying to adapt your script to my set of 12 rasters using the RasterToNumpyArray. Thanks again! ... View more

Thanks Darren, That works nicely. The only thing is that the Max value you get is 24 so that needs to be set back to 12 in the final output. Other than that, that's good. Thanks again Tony ... View more

It is something that I have in excel but it seems not that straight forward to code in Python (or maybe it is...) A function in python that calculates the "last occurrence"; in this case the last occurrence of 0 so it can start counting from there, I guess it would do the job. However, I cannot find in the documentation anything related to that. ... View more

Thanks Neil/Darren, Darren, WHEN running the code it is not relevant here. I'll try to graphically explain again: I have a stack of 12 cells each one representing presence or absence of data (1/0) for 12 months. I need to produce an output raster cell that tells me the maximum number of cumulative months with data. In example 1 above the output is 5 because there is data for sep, oct, nov, dec and jan In example 2 above the output is 3 because there is data for apr, may and jun Cheers ... View more

Thanks Neil. Sorry, I probably didn't explain myself properly. I'm not dealing with years here; just months of any year. I need to calculate cumulative number of months that hold data (in this case value 1). The calculation should loop through the 12 months with no fixed starting point. Hope it is a bit more clear now. Thanks again! ... View more