Latest Contributions by 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. 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)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Hi, I have multiple rasters named as follows for each year. But, I need to calculate mean for weekly ONLY for 09, 10 and 11 months based on the following week definition. Week 1: Day 01, 02, 03, 04, 05, 06 Week 2: Day 07, 08, 09, 10, 11, 12 Week 3: Day 13, 14, 15, 16, 17, 18 Week 4: Day 19, 20, 21, 22, 23, 24 Week 5: Day 25, 26, 27, 28, 29, 30, 31 Can anybody suggest me a way to process months 9, 10 and 11 in one go for the following script.?   Thanks in advance, import arcpy import arcpy.sa from arcpy.sa import * from arcpy import env import os arcpy.CheckOutExtension("Spatial") arcpy.OverwriteOutput = True arcpy.env.workspace = r'H:\Test_ONLY' scratch = r'H:\Test_ONLY\scratch' year = 1998 month = 01 fmaskrasters = arcpy.ListRasters() fmaskrasters.sort() if year % 4 > 0 and month == 9:      fmaskrasters1 = fmaskrasters[243:249]      fmaskrasters2 = fmaskrasters[249:255]      fmaskrasters3 = fmaskrasters[255:261]      fmaskrasters4 = fmaskrasters[261:267]      fmaskrasters5 = fmaskrasters[267:273] elif year % 4 == 0 and month == 9:      fmaskrasters1 = fmaskrasters[244:250]      fmaskrasters2 = fmaskrasters[250:256]      fmaskrasters3 = fmaskrasters[256:262]      fmaskrasters4 = fmaskrasters[262:268]      fmaskrasters5 = fmaskrasters[268:274] elif year % 4 > 0 and month == 10:      fmaskrasters1 = fmaskrasters[273:279]      fmaskrasters2 = fmaskrasters[279:285]      fmaskrasters3 = fmaskrasters[285:291]      fmaskrasters4 = fmaskrasters[291:297]      fmaskrasters5 = fmaskrasters[297:304] elif year % 4 == 0 and month == 10:      fmaskrasters1 = fmaskrasters[274:280]      fmaskrasters2 = fmaskrasters[280:286]      fmaskrasters3 = fmaskrasters[286:292]      fmaskrasters4 = fmaskrasters[292:298]      fmaskrasters5 = fmaskrasters[298:305] elif year % 4 > 0 and month == 11:      fmaskrasters1 = fmaskrasters[304:310]      fmaskrasters2 = fmaskrasters[310:316]      fmaskrasters3 = fmaskrasters[316:322]      fmaskrasters4 = fmaskrasters[322:328]      fmaskrasters5 = fmaskrasters[328:334] elif year % 4 == 0 and month == 11:      fmaskrasters1 = fmaskrasters[305:311]      fmaskrasters2 = fmaskrasters[311:317]      fmaskrasters3 = fmaskrasters[317:323]      fmaskrasters4 = fmaskrasters[323:329]      fmaskrasters5 = fmaskrasters[329:335] else:      print("Not configured in the script!") number = len(fmaskrasters) # Week 1 finras1 = scratch + "\\" + "Week_1_Avg" + str(year) + "_" + str(month) + ".tif" calc1 = arcpy.sa.CellStatistics([fmaskrasters1], statistics_type = "MEAN") arcpy.CopyRaster_management(calc1,finras1,"","","","","","32_BIT_FLOAT") # Week 2 finras2 = scratch + "\\" + "Week_2_Avg" + str(year) + "_" + str(month) + ".tif" calc2 = arcpy.sa.CellStatistics([fmaskrasters2], statistics_type = "MEAN") arcpy.CopyRaster_management(calc2,finras2,"","","","","","32_BIT_FLOAT") # Week 3 finras3 = scratch + "\\" + "Week_3_Avg" + str(year) + "_" + str(month) + ".tif" calc3 = arcpy.sa.CellStatistics([fmaskrasters3], statistics_type = "MEAN") arcpy.CopyRaster_management(calc3,finras3,"","","","","","32_BIT_FLOAT") # Week 4 finras4 = scratch + "\\" + "Week_4_Avg" + str(year) + "_" + str(month) + ".tif" calc4 = arcpy.sa.CellStatistics([fmaskrasters4], statistics_type = "MEAN") arcpy.CopyRaster_management(calc4,finras4,"","","","","","32_BIT_FLOAT") # Week 5 finras5 = scratch + "\\" + "Week_5_Avg" + str(year) + "_" + str(month) + ".tif" calc5 = arcpy.sa.CellStatistics([fmaskrasters5], statistics_type = "MEAN") arcpy.CopyRaster_management(calc5,finras5,"","","","","","32_BIT_FLOAT") ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Hi, I have a python script that works fine with ArcMap 10.4.1 version and read some comments on how to modify to suit to run with ArcGIS Pro version using "arcpy.AnalyzeToolsForPro_management" tool. After doing some modifications to the script I ran it in ArcGIS PRO but initially it worked fine and showed the progress up to 45% and it does not progress after 45% and remains in the same position for a long time. I did it several times and it does not complete running the script.  Basically, it converts to process netCDF file to extract a data variable, then project it and save it a .tif file. This process repeats for about 180 netCDF files. The original script works fine in ArcMap 10.4.1 version without any problem.  Any comments or suggestions of what might be happening.  Thanks. ... View more