Latest Contributions by RafaelloBergonse

Hi Xander Bakker, You DID figure it out! It worked perfectly, as you can see by the image below.  Thanks a lot for you help! From now on, I'll use your script whenever I need to calculate percentiles. I wish you a great Christmas!                                                             Best regards,                                                           Rafaello Bergonse ... View more

Hi Xander Bakker,  Thanks a lot for your patience and for the quick response! I tried the latest code you sent, but I get the same results. I am enclosing the shapefile as an annex so you can see what I get. This is the result of the latest code running over the shapefile you sent (so your correct values got overwitten). It may be noteworthy that the result I get for each polygon is invariably a very high value. For example, I get 41 for the first polygon in the shapefile, and 53 for the second one. When I compare these values with the rasters that are generated by the script while running the Extract by Mask tool, I see that there is only 1 pixel in the first raster with a value higher than 41, and only 2 pixels in the second raster with values higher than 53.  There is evidently nothing wrong with the code, so the problem must be related to my own machine or to some configuration issue. I am running ArcGIS desktop 10.7.1. I should also make clear that I cannot program in Python. The only things that I have changed in your code are the inicial statements describing the input raster and feature class, as well as the snap raster and the workspace. I believe the rest is the same independently of the inputs.                                                                                       Best regards,                                                                                   Rafaello Bergonse ... View more

Hi Xander Bakker. Like other users, I followed this conversation with interest and I am using your excellent script to calculate percentiles from a poligon shapefile and a raster. I have tried several iterations of your code, and all of them work without errors, with a big exception: for each polygon, the value obtained is the same, independently of the percentile under calculation. I am aware that other user had the same problem, and I have tried the script you posted in response, but the problem remains. Besides, in my case it is guaranteed that each polygon encompasses many, many raster cells. Would you please take a look at the sample polygon shapefile and raster I am annexing, so that you may find what is wrong? Thanks a lot!                                   Best Regards,                                 Rafaello Bergonse ... View more