Latest Contributions by ShouvikJha

Thank you. Output raster i am trying to save in tif format. I am getting following error massage.  Message File Name Line Position Traceback <module> D:\Arc-GIS-Python-Script\NetCDF_RASTER_TIME_SERIES.py 30 IndexError: string index out of range ‍‍‍‍ ... View more

FC Basson‌, yes exactly ,Good suggestion , i have added the line  ... View more

FC Basson‌. Thank you very much. Code perfectly working  ... View more

jskinner-esristaff‌. I found the reason of error. Actually in between OCTOBER and SEPTEMBER, there were one another file. I deleted it and rerun the program. Now program runs perfectly. Thank you very much for your kind cooperation.   ... View more

Jake Skinner‌. Thank you. I ran the code except SEPTEMBER month all months are calculating correctly. While it execute on month of SEPTEMBER, its showing error massage.  Traceback (most recent call last):   File "<module1>", line 23, in <module> IndexError: list index out of range‍‍‍ ... View more

Hi Jake, Thank you. i have replaced the line but still same issue. there is no output. Actually earlier whatever CSV we converted to point shape for month wise, then we merged  it and save it in different month's folder. Now using above programme i just want to make a loop function that program can take automatically point shape file from different month folder to perform IDW, and  save it as same month folder along with month number and Z filed name  (e.g. For JANUARY: 001_Max_Temper, 001_Min_Temper, 001_Wind, 001_Precipitat, 002_Relative_H, 001_Solar, For FEBRUARY:  002_Max_Temper, 002_Min_Temper, 002_Wind, 002_Precipitat, 002_Relative_H, 002_Solar) Below i have  attached a reference vector. Kindly find it.  ... View more

xander_bakker‌, Thank you. Now code perfectly running  ... View more

Hi Jake, Thank you . I am running the code on Pyscripter. I have updated my code, there is nothing error with my code, but  i am not getting any output raster from IDW, below is my updated code. Whats the possible error  import arcpy, os, calendar from arcpy import env from arcpy.sa import * arcpy.CheckOutExtension("Spatial") topWorkspace = r'D:\NPP_GUJARAT\WEATHER_DATA1\2011' arcpy.env.workspace = topWorkspace # Get dict of months and month-number (i.e. January = 001, March = 003 etc.) months = {calendar.month_name[i].upper(): str(i).zfill(3) for i in range(1, 13)} # Get dict of months and month number (i.e. January = 001, March = 003 etc.) # Step through list of all folders for folderPath in arcpy.ListWorkspaces():        baseName = os.path.basename(folderPath).upper()        if baseName in months: # Test that subfolder is a month name         monthNumber = months[baseName] # Get month-number for use in output filename        arcpy.env.workspace = folderPath         #Loop every month folder ta take shape file from each month (e.g. Merged_001 (JANUARY), Merged_002 (FEBRUARY)....Merged_012 (DECEMBER))        shapeFiles = arcpy.ListFeatureClasses('Merged_{}*.shp') zFields = ["Max_Temper", "Min_Temper", "Precipitat", "Wind", "Relative_H", "Solar"] # IDW variable cellSize = 0.002298707671 power = 2 searchRadius = RadiusVariable(10, 150000) #Mask region of interest mask="D:\Gujarta Shape file\GUJARATSTATE.shp" for shape_num in shapeFiles: Zfield_num = arcpy.ListFeatureClasses(os.path.join(folderPath, shape_num)) # Execute IDW for zField in zFields:   outIDW = Idw(shapeFiles, zField, cellSize, power, searchRadius) # Execute Mask IDWMASk = ExtractByMask(outIDW, mask) if not zField == "Solar":         # Save output, except Solar Radiation raster (e.g for JANUARY:001_Max_Temper, FEBRUARY:002_Max_Temper)         out_name_IDW = os.path.join(topWorkspace, '(months)_.TIF'.format(zField_num))         IDWMASk.save(out_name_IDW) else:         #Only for Solar Radiation raster         PAR = (IDWMASk * 30) * 0.5         #Save output, Only for Solar Radiation raster         out_name_IDW = os.path.join(topWorkspace, '(months)_.TIF'.format(zField))         PAR.save(out_name_IDW) #print done print 'done'‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Xander Bakker‌, In just trying to modify above code, instead of multiplying 0.98, i just want to replace it by constant raster , below is my updated code. but its not working,  i tried both line to replace 0.98 by constant raster.  ras_LUE = getRaster(D:\MOD-REF\NDVI\LUE-CASA.img) ras_LUE = arcpy.raster (D:\MOD-REF\NDVI\LUE-CASA.img)‍‍‍‍‍‍‍‍‍‍‍‍ for ras_num in [str(i).zfill(3) for i in range(1, 13)]:         # get required rasters for current month number         ras_APAR = getRasterFromList(ras_num, lst_ras_APAR, ws_in_apar)         ras_TSCALAR = getRasterFromList(ras_num, lst_ras_TSCALAR, ws_in_tscalar)         ras_WSCALAR = getRasterFromList(ras_num, lst_ras_WSCALAR, ws_in_wscalar)         #ras_LUE         # verify if raster are valid         if all([not ras_APAR is None, not ras_TSCALAR is None, not ras_WSCALAR is None]):                             # 0.98 will replace by constant raster (ras_LUE)             ras_NPP = (ras_APAR * ras_TSCALAR * ras_WSCALAR * 0.98)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Jake Skinner‌, I am trying to solve the above script problem, but i am getting continue error. Error displaying IndentationError unexpected indent (IDW_PAR_Year_WISE.py, line 10) ‍‍ How to resolve the issue. Thanks  ... View more

xander_bakker‌. Excellent. A totally new work i learned from you. Thank you very much   ... View more

Hi Jake,I am trying to modify above code, as keep it IDW function same, But instead of giving path for every month shape file , I just trying to use loop function that program can automatic take shape file from different month folder (Name of point shape file in different month folder : JANUARY: Merged_001, FEBRUARY: Merged_002.....DECEMBER: Merged_012 ) and save as it same folder along with month number and Z field name (e.g for JANUARY : 001_Max_Temper, 001_Min_Temper, 001_Precipitat, FEBRUARY : 002_Max_Temper, 002_Min_Temper, 002_Precipitatetc.) . Below is my working code. But its showing error.  import arcpy, os, calendar arcpy.CheckOutExtension("Spatial") topWorkspace = r'D:\NPP_GUJARAT\WEATHER_DATA\2011' arcpy.env.workspace = topWorkspace # Get dict of months and month-number (i.e. January = 001, March = 003 etc.) months = {calendar.month_name[i].upper(): str(i).zfill(3) for i in range(1, 13)} # Get dict of months and month number (i.e. January = 001, March = 003 etc.) # Step through list of all folders for folderPath in arcpy.ListWorkspaces():     baseName = os.path.basename(folderPath).upper()     if baseName in months: # Test that subfolder is a month name         monthNumber = months[baseName] # Get month-number for use in output filename         arcpy.env.workspace = folderPath         #Loop every month folder ta take shape file from each month (e.g. Merged_001 (JANUARY), Merged_002 (FEBRUARY)....Merged_012 (DECEMBER))         shapeFiles = arcpy.ListFeatureClasses('*Merged_{}.shp') zFields = ["Max_Temper", "Min_Temper", "Precipitat", "Wind", "Relative_H", "Solar"] # IDW variable cellSize = 0.002298707671 power = 2 searchRadius = RadiusVariable(10, 150000) #Mask region of interest mask="D:\Gujarta Shape file\GUJARATSTATE.shp" for shape_num in shapeFiles: Zfield_num = arcpy.ListFeatureClasses(os.path.join(folderPath, shape_num)) # Execute IDW for zField in zFields:   outIDW = Idw(shapeFiles, zField, cellSize, power, searchRadius) # Execute Mask IDWMASk = ExtractByMask(outIDW, mask) if not zField == "Solar":         # Save output, except Solar Radiation raster (e.g for JANUARY:001_Max_Temper, FEBRUARY:002_Max_Temper)         out_name_IDW = os.path.join(topWorkspace, '(months)_.TIF'.format(zField_num))         IDWMASk.save(out_name_IDW) else:         #Only for Solar Radiation raster         PAR = (IDWMASk * 30) * 0.5         #Save output, Only for Solar Radiation raster         out_name_IDW = os.path.join(topWorkspace, '(months)_.TIF'.format(zField))         PAR.save(out_name_IDW) #print done print 'done'‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Thank you. I am trying to do it. I am also not sure Multi-plot facility in Arc GIS has or not. I am taking suggestion from other developer also. Hopefully problem would be be resolve. Which snapshot i have attached, there i have added 3 data frame for reference purpose only. similarly i want to add total 12 raster into one frame and export it to pdf format. ... View more