Latest Contributions by ShouvikJha

Hi All,  I have one NetCDF file, that contain monthly value for the different year (1910 - 2014). I want to extract few years data (e.g 2000 - 2014), date range reading from a csv file  and save it as a raster map as the  name format of year-month-date.  I am working on below script but I am getting error  Below i have attached the csv file which i am using for date range purpose.  Any help would be greatly appreciated! Thanks  The image showing the format of Year-month-date Error massage  Traceback (most recent call last):   File "<module4>", line 18, in <module> ValueError: 'Date' is not in list‍‍‍‍‍‍ Here what i have  so far  import os, sys import arcpy # Check out any necessary licenses arcpy.CheckOutExtension("spatial") arcpy.env.overwriteOutput = True # Script arguments netCDF = "C:\\Users\\Downloads\\SPIE" rainfall = "C:\\Users\\Downloads\\SPIE\\New folder" arcpy.env.workspace = netCDF # Read Date from csv file eveDate = open ("C:\\Users\\Downloads\\SPIE\\Name1.csv") headerLine = eveDate.readline() valueList = headerLine.split(",") dateValueIndex = valueList.index("Date") eventList = [] for line in eveDate.readlines():     segmenLine = line.split(",")     variable = "spei"     x_dimension = "lon"     y_dimension = "lat"     band_dimension = "" for year in range(2010):     nc_name = 'spei.%d.nc' % (year,)     valueSelectionMethod = "BY_VALUE"     outFile = "Pre"     # extract dimensionValues from csv file     arcpy.MakeNetCDFRasterLayer_md("spei.nc", variable, x_dimension, y_dimension, outFile, band_dimension, segmenLine[dateValueIndex], valueSelectionMethod)     print "layer done"     #copy and save as raster tif file     arcpy.CopyRaster_management(outFile, rainfall + segmenLine[dateValueIndex] + ".tif" , "", "", "", "NONE", "NONE", "")     print "raster done"‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Hi,  I am trying to execute an equation into python program considering monthly folder, but I am getting error, Below i have attached reference of the equation also, Where P and R  for rasterlist1 and rasterlist2 respectively,  Error message is. Traceback (most recent call last):   File "<module3>", line 23, in <module> TypeError: can't multiply sequence by non-int of type 'list'‍‍‍‍‍‍‍‍‍‍‍‍ Script i am working,  import arcpy, os, calendar import numpy from arcpy.sa import * arcpy.env.overwriteOutput = True arcpy.CheckOutExtension("Spatial") topWorkspace = r'D:\NPP_GUJARAT_MONSOON' ws_out = r'D:\NPP_GUJARAT_MONSOON\EVA3' arcpy.env.workspace = topWorkspace 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         # list rasters (Erlier:*Max_Temper.tif, We changed it to below)         rasterList1 = arcpy.ListRasters('Precipitat*.tif')         rasterList2 = arcpy.ListRasters('Solar*.tif')         print rasterList1, rasterList2         for ras_name in rasterList1:             ras_mean = arcpy.Raster(os.path.join(folderPath, ras_name))         for ras_name in rasterList2:             ras_mean = arcpy.Raster(os.path.join(folderPath, ras_name))         outCellStatistics = ((rasterList1 * rasterList2) * ((Square(rasterList1)) + (Square(rasterList2)) + rasterList1 * rasterList2)) /rasterList1 + rasterList2 * ((Square(rasterList1)) + (Square(rasterList2)))         #Save the output         outRasterName = os.path.join(ws_out, "EO_{}.tif".format(monthNumber))         outCellStatistics.save(outRasterName)         print outRasterName #print done print 'done' ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍                                     ... View more