Latest Contributions by BIJOYGAYEN

For miss TOA_B1 and a TOA_B2 raster files, we have to write these two lines: TOA_B1 = glob.glob (d3 + os.sep + "* EV_500_RefSB_1-EV_500_RefSB.tif") TOA_B2 = glob.glob (d3 + os.sep + "* EV_500_RefSB_2-EV_500_RefSB.tif") I wrote that, there is no problem with it. When I check my "test data" (Subset Data,indiviual file size 5-10kb), the code works well but it does not work on row data sets (whole India,indiviual file size 54mb) then it gives me this error: ph.I have attached the data processing screen short, please check it.I've added some files in "test_data" for your test. Please check this problem. There have no option to reduce and modification of this code? Xander Bakker‌ Test_Data processing:working well Row_Data processing:Error is coming Code: #------------------------------------------------------------------------------- # Name: module1 # Purpose: # # Author: HONEY # # Created: 18-06-2019 # Copyright: (c) HONEY 2019 # Licence: <your licence> #------------------------------------------------------------------------------- import datetime start = datetime.datetime.now() print 'start run: %s\n' % (start) import arcpy ,os ,sys,csv,errno from arcpy import env from arcpy.sa import * import datetime import re import numpy as np import glob import itertools arcpy.env.overwriteOutput = True cellsize = "F:\\DB_test_data\\VAR2\\TEST_DATA\\Clip\\MOD02HKM_A2017001_0530_NDVI_AA.img" #Call cloud mask images from directory d1="F:\\DB_test_data\\VAR2\\TEST_DATA\\Clip" CLDMASK = glob.glob(d1 + os.sep + "*.Aerosol_Cldmask_Land_Ocean-Aerosol_Cldmask_Land_Ocean.tif") CLDMASK.sort() if len(CLDMASK) == 0: print 'Could not open the CLDMASK raster files' sys.exit(1) else: print 'The CLDMASK raster files was opened successfully' #Call AOD images from directory d2 = r"F:\\DB_test_data\\VAR2\\TEST_DATA\\Clip" AOD = glob.glob(d2 + os.sep + "*Corrected_Optical_Depth_Land_2-Corrected_Optical_Depth_Land.tif") AOD.sort() if len(AOD) == 0: print 'Could not open the AOD raster files' sys.exit(1) else: print 'The AOD raster files was opened successfully' #Call MOD02HKM_TOA-B1_TOA-B2 images from directory d3="F:\\DB_test_data\\VAR2\\TEST_DATA\\Clip" TOA_B1 = glob.glob(d3 + os.sep + "*EV_500_RefSB_1-EV_500_RefSB.tif") TOA_B1.sort() if len(TOA_B1) == 0: print 'Could not open the TOA_B1 raster files' sys.exit(1) else: print 'The TOA_B1 raster files was opened successfully' TOA_B2 = glob.glob(d3 + os.sep + "*EV_500_RefSB_2-EV_500_RefSB.tif") TOA_B2.sort() if len(TOA_B2) == 0: print 'Could not open the TOA_B2 raster files' sys.exit(1) else: print 'The TOA_B2 raster files was opened successfully' #Call MOD03 (Height,SenZen,SenAzm,SolZen,SolAzm) images from directory d4="F:\\DB_test_data\\VAR2\\TEST_DATA\\Clip" Height = glob.glob(d4 + os.sep + "*.Height-Height.tif") Height.sort() if len(Height) == 0: print 'Could not open the Height raster files' sys.exit(1) else: print 'The Height raster files was opened successfully' SenZen = glob.glob(d4 + os.sep + "*.SensorZenith-SensorZenith.tif") SenZen.sort() if len(SenZen) == 0: print 'Could not open the SenZen raster files' sys.exit(1) else: print 'The SenZen raster files was opened successfully' SenAzm = glob.glob(d4 + os.sep + "*.SensorAzimuth-SensorAzimuth.tif") SenAzm.sort() if len(SenAzm) == 0: print 'Could not open the SenAzm raster files' sys.exit(1) else: print 'The SenAzm raster files was opened successfully' SolZen = glob.glob(d4 + os.sep + "*.SolarZenith-SolarZenith.tif") SolZen.sort() if len(SolZen) == 0: print 'Could not open the SolZen raster files' sys.exit(1) else: print 'The SolZen raster files was opened successfully' SolAzm = glob.glob(d4 + os.sep + "*.SolarAzimuth-SolarAzimuth.tif") SolAzm.sort() if len(SolAzm) == 0: print 'Could not open the SolAzm raster files' sys.exit(1) else: print 'The SolAzm raster files was opened successfully' #Call MOD09 images from directory d5="F:\\DB_test_data\\VAR2\\TEST_DATA\\Clip" Sur_B1 = glob.glob(d5 + os.sep + "*.500m_Surface_Reflectance_Band_1-500m_Surface_Reflectance_Band_1.tif") Sur_B1.sort() if len(Sur_B1) == 0: print 'Could not open the Sur_B1 raster files' sys.exit(1) else: print 'The Sur_B1 raster files was opened successfully' Sur_B2 = glob.glob(d5 + os.sep + "*.500m_Surface_Reflectance_Band_2-500m_Surface_Reflectance_Band_2.tif") Sur_B2.sort() if len(Sur_B2) == 0: print 'Could not open the Sur_B2 raster files' sys.exit(1) else: print 'The Sur_B2 raster files was opened successfully' Sur_B3 = glob.glob(d5 + os.sep + "*.500m_Surface_Reflectance_Band_3-500m_Surface_Reflectance_Band_3.tif") Sur_B3.sort() if len(Sur_B3) == 0: print 'Could not open the Sur_B3 raster files' sys.exit(1) else: print 'The Sur_B3 raster files was opened successfully' Land_cover="F:\\DB_test_data\\VAR2\\TEST_DATA\\Clip\\LAND_COVER_TYPE_1_Grid_2D.img" if len(Land_cover) == 0: print 'Could not open the Land_cover raster files' sys.exit(1) else: print 'The Land_cover raster files was opened successfully' outdir="F:\\DB_test_data\\TEST_RAY\\TEST1\\c\\d\\" cellsize = "F:\\DB_test_data\\VAR2\\TEST_DATA\\Clip\\MOD02HKM_A2017001_0530_NDVI_AA.img" for a,b,c,d,e,f,g,h,i,j,k,l, in zip (CLDMASK,AOD,TOA_B1,TOA_B2,Height,SenZen,SenAzm,SolZen,SolAzm,Sur_B1,Sur_B2,Sur_B3): AA=a.split("\\")[5][0:114] print ("processing_a:"+ AA) BB=b.split("\\")[9][0:120] print ("processing_b:"+ BB) CC=c.split("\\")[5][0:88] print ("processing_c:"+ CC) DD=d.split("\\")[5][0:88] print ("processing_d:"+ DD) EE=e.split("\\")[5][0:71] print ("processing_e:"+ EE) FF=f.split("\\")[5][0:83] print ("processing_f:"+ FF) GG=g.split("\\")[5][0:85] print ("processing_g:"+ GG) HH=h.split("\\")[5][0:81] print ("processing_h:"+ HH) II=i.split("\\")[5][0:83] print ("processing_i:"+ II) JJ=j.split("\\")[5][0:121] print ("processing_j:"+ JJ) KK=k.split("\\")[5][0:121] print ("processing_k:"+ KK) LL=l.split("\\")[5][0:121] print ("processing_l:"+ LL) ####################################################################### # #Cloud Mask scalling Scale_factor1 = float(1.0) add_offset1 = float(0.0) setnull1 =arcpy.gp.SetNull_sa(a,a, "in_memory/dat1", "\"Value\" <= 0") ras1=arcpy.Raster(setnull1) cldmask=(ras1-add_offset1)*Scale_factor1 OutRaster1 = os.path.join(outdir,'{0}.img'.format(AA)) arcpy.CopyRaster_management(cldmask,OutRaster1) arcpy.Delete_management("in_memory/dat1") arcpy.Delete_management(ras1) #AOD scalling and set lessthan 0.1 AOD Scale_factor2 = float(0.0010000000474974513) add_offset2 = float(0.0) setnull2 =arcpy.gp.SetNull_sa(b,b, "in_memory/dat2", "\"Value\" = -9999") ras2=arcpy.Raster(setnull2) aod=(ras2-add_offset2)*Scale_factor2 aod_1 = Con((aod >= 0.0) & (aod <= 0.1),1) OutRaster2 = os.path.join(outdir,'{0}.img'.format(BB)) arcpy.CopyRaster_management(aod_1,OutRaster2) arcpy.Delete_management("in_memory/dat2") #TOA_B1 scalling Scale_factor3 = float(0.000053811826) add_offset3 = float(-0.0) setnull3 =arcpy.gp.SetNull_sa(c,c, "in_memory/dat3", "\"Value\" = 65535") ras3=arcpy.Raster(setnull3) toa_b1 =(ras3-(add_offset3))*Scale_factor3 OutRaster3 = os.path.join(outdir,'{0}.img'.format(CC)) arcpy.CopyRaster_management(toa_b1,OutRaster3) arcpy.Delete_management("in_memory/dat3") #TOA_B2 Scalling Scale_factor4 = float(0.00003255546) add_offset4 = float(-0.0) setnull4 =arcpy.gp.SetNull_sa(d,d, "in_memory/dat4", "\"Value\" = 65535") ras4=arcpy.Raster(setnull4) toa_b2=(ras4-(add_offset4))*Scale_factor4 OutRaster4 = os.path.join(outdir,'{0}.img'.format(DD)) arcpy.CopyRaster_management(toa_b2,OutRaster4) arcpy.Delete_management("in_memory/dat4") #Height setnull =arcpy.gp.SetNull_sa(e,e, "in_memory/dat5", "\"Value\" = -32767") ras=arcpy.Raster(setnull) height=ras OutRaster5 = os.path.join(outdir,'{0}.img'.format(EE)) arcpy.CopyRaster_management(height,OutRaster5) arcpy.Delete_management("in_memory/dat5") #Sensor Zenith angle Scale_factor6 = float(0.01) add_offset6 = float(0.0) setnull6 =arcpy.gp.SetNull_sa(f,f, "in_memory/dat6", "\"Value\" = -32767") ras6=arcpy.Raster(setnull6) senzen=(ras6-add_offset6)*Scale_factor6 OutRaster6 = os.path.join(outdir,'{0}.img'.format(FF)) arcpy.CopyRaster_management(senzen,OutRaster6) arcpy.Delete_management("in_memory/dat6") #Sensor azimuth angle setnull7 =arcpy.gp.SetNull_sa(g,g, "in_memory/dat7", "\"Value\" = -32767") ras7=arcpy.Raster(setnull7) senazm=(ras7-add_offset6)*Scale_factor6 OutRaster7 = os.path.join(outdir,'{0}.img'.format(GG)) arcpy.CopyRaster_management(senazm,OutRaster7) arcpy.Delete_management("in_memory/dat7") #solar zenithal angle setnull8 =arcpy.gp.SetNull_sa(h,h, "in_memory/dat8", "\"Value\" = -32767") ras8=arcpy.Raster(setnull8) solzen=(ras8-add_offset6)*Scale_factor6 OutRaster8 = os.path.join(outdir,'{0}.img'.format(HH)) arcpy.CopyRaster_management(solzen,OutRaster8) arcpy.Delete_management("in_memory/dat8") #solar azimuth angle setnull9 =arcpy.gp.SetNull_sa(i,i, "in_memory/dat9", "\"Value\" = -32767") ras9=arcpy.Raster(setnull9) solazm=(ras9-add_offset6)*Scale_factor6 OutRaster9 = os.path.join(outdir,'{0}.img'.format(II)) arcpy.CopyRaster_management(solazm,OutRaster9) arcpy.Delete_management("in_memory/dat9") #(MODO9)surface reflectance band 1 Scale_factor_10 = float(0.0001) add_offset_10 = float(0.0) setnull_10 =arcpy.gp.SetNull_sa(j,j, "in_memory/dat10", "\"Value\" = -28672") ras_10=arcpy.Raster(setnull_10) sur_b1=(ras_10-add_offset_10)*Scale_factor_10 OutRaster10 = os.path.join(outdir,'{0}.img'.format(JJ)) arcpy.CopyRaster_management(sur_b1,OutRaster10) arcpy.Delete_management("in_memory/dat10") #(MODO9)surface reflectance band 2 setnull_11 =arcpy.gp.SetNull_sa(k,k, "in_memory/dat11", "\"Value\" = -28672") ras_11=arcpy.Raster(setnull_11) sur_b2=(ras_11-add_offset_10)*Scale_factor_10 OutRaster11 = os.path.join(outdir,'{0}.img'.format(KK)) arcpy.CopyRaster_management(sur_b2,OutRaster11) arcpy.Delete_management("in_memory/dat11") #(MODO9)surface reflectance band 3 setnull_12 =arcpy.gp.SetNull_sa(l,l, "in_memory/dat12", "\"Value\" = -28672") ras_12=arcpy.Raster(setnull_12) sur_b3=(ras_12- add_offset_10)*Scale_factor_10 OutRaster12 = os.path.join(outdir,'{0}.img'.format(LL)) arcpy.CopyRaster_management(sur_b3,OutRaster12) arcpy.Delete_management("in_memory/dat12") print arcpy.GetMessages() print 'finished run: %s\n\n' % (datetime.datetime.now() - start)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

I have already uploaded the Test_data kindly check it. Also I am uploading my full code. #------------------------------------------------------------------------------- # Name: module1 # Purpose: # # Author: HONEY # # Created: 18-06-2019 # Copyright: (c) HONEY 2019 # Licence: <your licence> #------------------------------------------------------------------------------- import datetime start = datetime.datetime.now() print 'start run: %s\n' % (start) import arcpy ,os ,sys,csv,errno from arcpy import env from arcpy.sa import * import datetime import re import numpy as np import glob import itertools arcpy.env.overwriteOutput = True cellsize = "F:\\DB_test_data\\TEST_RAY\\TEST1\\MOD02HKM_A2017001_0530_NDVI_AA.img" #Call cloud mask images from directory d1="F:\\DB_test_data\\VAR2\\cldmask\\" CLDMASK = glob.glob(d1 + os.sep + "*.Aerosol_Cldmask_Land_Ocean-Aerosol_Cldmask_Land_Ocean.tif") CLDMASK.sort() if CLDMASK is None: print 'Could not open the CLDMASK raster files' sys.exit(1) else: print 'The CLDMASK raster files was opened successfully' #Call AOD images from directory d2 = r"F:\\DB_test_data\\VAR2\\MOD04_l2" AOD = glob.glob(d2 + os.sep + "*Corrected_Optical_Depth_Land_2-Corrected_Optical_Depth_Land.tif") AOD.sort() if AOD is None: print 'Could not open the AOD raster files' sys.exit(1) else: print 'The AOD raster files was opened successfully' #Call MOD02HKM_TOA-B1_TOA-B2 images from directory d3="F:\\DB_test_data\\VAR2\\IDL\\MOD02HKM\\" TOA_B1 = glob.glob(d3 + os.sep + "*RefSB_1-EV_250_Aggr500_RefSB.tif") TOA_B1.sort() if TOA_B1 is None: print 'Could not open the TOA_B1 raster files' sys.exit(1) else: print 'The TOA_B1 raster files was opened successfully' TOA_B2 = glob.glob(d3 + os.sep + "*RefSB_2-EV_250_Aggr500_RefSB.tif") TOA_B2.sort() if TOA_B2 is None: print 'Could not open the TOA_B2 raster files' sys.exit(1) else: print 'The TOA_B2 raster files was opened successfully' #Call MOD03 (Height,SenZen,SenAzm,SolZen,SolAzm) images from directory d4="F:\\DB_test_data\\VAR2\\IDL\\MOD03_R\\" Height = glob.glob(d4 + os.sep + "*.Height-Height.tif") Height.sort() if Height is None: print 'Could not open the Height raster files' sys.exit(1) else: print 'The Height raster files was opened successfully' SenZen = glob.glob(d4 + os.sep + "*.SensorZenith-SensorZenith.tif") SenZen.sort() if SenZen is None: print 'Could not open the SenZen raster files' sys.exit(1) else: print 'The SenZen raster files was opened successfully' SenAzm = glob.glob(d4 + os.sep + "*.SensorAzimuth-SensorAzimuth.tif") SenAzm.sort() if SenAzm is None: print 'Could not open the SenAzm raster files' sys.exit(1) else: print 'The SenAzm raster files was opened successfully' SolZen = glob.glob(d4 + os.sep + "*.SolarZenith-SolarZenith.tif") SolZen.sort() if SolZen is None: print 'Could not open the SolZen raster files' sys.exit(1) else: print 'The SolZen raster files was opened successfully' SolAzm = glob.glob(d4 + os.sep + "*.SolarAzimuth-SolarAzimuth.tif") SolAzm.sort() if SolAzm is None: print 'Could not open the SolAzm raster files' sys.exit(1) else: print 'The SolAzm raster files was opened successfully' #Call MOD09 images from directory d5="F:\\DB_test_data\\VAR2\\IDL\\MOD09\\" Sur_B1 = glob.glob(d5 + os.sep + "*.500m_Surface_Reflectance_Band_1-500m_Surface_Reflectance_Band_1.tif") Sur_B1 .sort() if Sur_B1 is None: print 'Could not open the Sur_B1 raster files' sys.exit(1) else: print 'The Sur_B1 raster files was opened successfully' Sur_B2 = glob.glob(d5 + os.sep + "*.500m_Surface_Reflectance_Band_2-500m_Surface_Reflectance_Band_2.tif") Sur_B2 .sort() if Sur_B2 is None: print 'Could not open the Sur_B2 raster files' sys.exit(1) else: print 'The Sur_B2 raster files was opened successfully' Sur_B3 = glob.glob(d5 + os.sep + "*.500m_Surface_Reflectance_Band_3-500m_Surface_Reflectance_Band_3.tif") Sur_B3 .sort() if Sur_B3 is None: print 'Could not open the Sur_B3 raster files' sys.exit(1) else: print 'The Sur_B3 raster files was opened successfully' Land_cover="F:\\DB_test_data\\MCD12C1\\LAND_COVER_TYPE_1_Grid_2D.img" if Land_cover is None: print 'Could not open the Land_cover raster files' sys.exit(1) else: print 'The Land_cover raster files was opened successfully' outdir="F:\\DB_test_data\\TEST_RAY\\TEST1\\c\\d\\" cellsize = "F:\\DB_test_data\\TEST_RAY\\TEST1\\MOD02HKM_A2017001_0530_NDVI_AA.img" for a,b,c,d,e,f,g,h,i,j,k,l in zip (CLDMASK,AOD,TOA_B1,TOA_B2,Height,SenZen,SenAzm,SolZen,SolAzm,Sur_B1,Sur_B2,Sur_B3): print ("processing_a:"+ a) print ("processing_b:"+ b) print ("processing_c:"+ c) print ("processing_d:"+ d) print ("processing_e:"+ e) print ("processing_f:"+ f) print ("processing_g:"+ g) print ("processing_h:"+ h) print ("processing_i:"+ i) print ("processing_j:"+ j) print ("processing_k:"+ k) print ("processing_l:"+ l) ####################################################################### #Cloud Mask scalling Scale_factor1 = float(1.0) add_offset1 = float(0.0) setnull1 =arcpy.gp.SetNull_sa(a,a, "in_memory/dat1", "\"Value\" <= 0") ras1=arcpy.Raster(setnull1) cldmask=(ras1-add_offset1)*Scale_factor1 OutRaster1 = os.path.join(outdir,'MOd0.{0}.img'.format("cldmask")) arcpy.CopyRaster_management(cldmask,OutRaster1) arcpy.Delete_management("in_memory/dat1") arcpy.Delete_management(ras1) #AOD scalling and set lessthan 0.1 AOD Scale_factor2 = float(0.0010000000474974513) add_offset2 = float(0.0) setnull2 =arcpy.gp.SetNull_sa(b,b, "in_memory/dat2", "\"Value\" = -9999") ras2=arcpy.Raster(setnull2) aod=(ras2-add_offset2)*Scale_factor2 aod_1 = Con((aod >= 0.0) & (aod <= 0.1),1) OutRaster2 = os.path.join(outdir,'MOd0.{0}.img'.format("aod_1")) arcpy.CopyRaster_management(aod_1,OutRaster2) arcpy.Delete_management("in_memory/dat2") #TOA_B1 scalling Scale_factor3 = float(0.000053811826) add_offset3 = float(-0.0) setnull3 =arcpy.gp.SetNull_sa(c,c, "in_memory/dat3", "\"Value\" = 65535") ras3=arcpy.Raster(setnull3) toa_b1 =(ras3-(add_offset3))*Scale_factor3 OutRaster3 = os.path.join(outdir,'MOd0.{0}.img'.format("toa_b1")) arcpy.CopyRaster_management(OutRaster3,toa_b1) arcpy.Delete_management("in_memory/dat3") #TOA_B2 Scalling Scale_factor4 = float(0.00003255546) add_offset4 = float(-0.0) setnull4 =arcpy.gp.SetNull_sa(d,d, "in_memory/dat4", "\"Value\" = 65535") ras4=arcpy.Raster(setnull4) toa_b2=(ras4-(add_offset4))*Scale_factor4 OutRaster4 = os.path.join(outdir,'MOd0.{0}.img'.format("toa_b2")) arcpy.CopyRaster_management(toa_b2,OutRaster4) arcpy.Delete_management("in_memory/dat4") #Height setnull =arcpy.gp.SetNull_sa(e,e, "in_memory/dat5", "\"Value\" = -32767") ras=arcpy.Raster(setnull) height=ras OutRaster5 = os.path.join(outdir,'MOd0.{0}.img'.format("height")) arcpy.CopyRaster_management(height,OutRaster5) arcpy.Delete_management("in_memory/dat5") #Sensor Zenith angle Scale_factor6 = float(0.01) add_offset6 = float(0.0) setnull6 =arcpy.gp.SetNull_sa(f,f, "in_memory/dat6", "\"Value\" = -32767") ras6=arcpy.Raster(setnull6) senzen=(ras6-add_offset6)*Scale_factor6 OutRaster6 = os.path.join(outdir,'MOd0.{0}.img'.format("senzen")) arcpy.CopyRaster_management(senzen,OutRaster6) arcpy.Delete_management("in_memory/dat6") #Sensor azimuth angle setnull7 =arcpy.gp.SetNull_sa(g,g, "in_memory/dat7", "\"Value\" = -32767") ras7=arcpy.Raster(setnull7) senazm=(ras7-add_offset6)*Scale_factor6 OutRaster7 = os.path.join(outdir,'MOd0.{0}.img'.format("senazm")) arcpy.CopyRaster_management(senazm,OutRaster7) arcpy.Delete_management("in_memory/dat7") #solar zenithal angle setnull8 =arcpy.gp.SetNull_sa(h,h, "in_memory/dat8", "\"Value\" = -32767") ras8=arcpy.Raster(setnull8) solzen=(ras8-add_offset6)*Scale_factor6 OutRaster8 = os.path.join(outdir,'MOd0.{0}.img'.format("solzen")) arcpy.CopyRaster_management(solzen,OutRaster8) arcpy.Delete_management("in_memory/dat8") #solar azimuth angle setnull9 =arcpy.gp.SetNull_sa(i,i, "in_memory/dat9", "\"Value\" = -32767") ras9=arcpy.Raster(setnull9) solazm=(ras9-add_offset6)*Scale_factor6 OutRaster9 = os.path.join(outdir,'MOd0.{0}.img'.format("solazm")) arcpy.CopyRaster_management(solazm,OutRaster9) arcpy.Delete_management("in_memory/dat9") #(MODO9)surface reflectance band 1 Scale_factor_10 = float(0.0001) add_offset_10 = float(0.0) setnull_10 =arcpy.gp.SetNull_sa(j,j, "in_memory/dat10", "\"Value\" = -28672") ras_10=arcpy.Raster(setnull_10) sur_b1=(ras_10-add_offset_10)*Scale_factor_10 OutRaster10 = os.path.join(outdir,'MOd0.{0}.img'.format("sur_b1")) arcpy.CopyRaster_management(sur_b1,OutRaster10) arcpy.Delete_management("in_memory/dat10") #(MODO9)surface reflectance band 2 setnull_11 =arcpy.gp.SetNull_sa(k,k, "in_memory/dat11", "\"Value\" = -28672") ras_11=arcpy.Raster(setnull_11) sur_b2=(ras_11-add_offset_10)*Scale_factor_10 OutRaster11 = os.path.join(outdir,'MOd0.{0}.img'.format("sur_b2")) arcpy.CopyRaster_management(sur_b2,OutRaster11) arcpy.Delete_management("in_memory/dat11") #(MODO9)surface reflectance band 3 setnull_12 =arcpy.gp.SetNull_sa(l,l, "in_memory/dat12", "\"Value\" = -28672") ras_12=arcpy.Raster(setnull_12) sur_b3=(ras_12- add_offset_10)*Scale_factor_10 OutRaster12 = os.path.join(outdir,'MOd0.{0}.img'.format("sur_b3")) arcpy.CopyRaster_management(sur_b3,OutRaster12) arcpy.Delete_management("in_memory/dat12") ############################################################################# # # Process: Extract by Cloud_Mask # tempEnvironment0 = arcpy.env.cellSize # arcpy.env.cellSize = "MAXOF" # toa_bb1=arcpy.gp.ExtractByMask_sa(toa_b1, cldmask) # toa_b22=arcpy.gp.ExtractByMask_sa(toa_b2, cldmask) # heightt=arcpy.gp.ExtractByMask_sa(height, cldmask) # senzenn=arcpy.gp.ExtractByMask_sa(senzen, cldmask) # senazmm=arcpy.gp.ExtractByMask_sa(senazm, cldmask) # solzenn=arcpy.gp.ExtractByMask_sa(solzen, cldmask) # solazmm=arcpy.gp.ExtractByMask_sa(solazm, cldmask) # sur_b11=arcpy.gp.ExtractByMask_sa(sur_b1, cldmask) # sur_b22=arcpy.gp.ExtractByMask_sa(sur_b2, cldmask) # sur_b33=arcpy.gp.ExtractByMask_sa(sur_b3, cldmask) # arcpy.env.cellSize = tempEnvironment0 # # Process: Extract by Mask using AOD less than 0.1 value # tempEnvironment0 = arcpy.env.cellSize # arcpy.env.cellSize = cellsize # toa_b1=arcpy.gp.ExtractByMask_sa(toa_bb1, aod_1) # toa_b2=arcpy.gp.ExtractByMask_sa(toa_bb2, aod_1) # height=arcpy.gp.ExtractByMask_sa(heightt, aod_1) # senzen=arcpy.gp.ExtractByMask_sa(senzenn, aod_1) # senazm=arcpy.gp.ExtractByMask_sa(senazmm, aod_1) # solzen=arcpy.gp.ExtractByMask_sa(solzenn, aod_1) # solazm=arcpy.gp.ExtractByMask_sa(solazmm, aod_1) # sur_b1=arcpy.gp.ExtractByMask_sa(sur_b11, aod_1) # sur_b2=arcpy.gp.ExtractByMask_sa(sur_b22, aod_1) # sur_b3=arcpy.gp.ExtractByMask_sa(sur_b33, aod_1) # arcpy.env.cellSize = tempEnvironment0 ############################################################################# #Calculate the scattering angle from Resample MOD03 products #Css = Cos(S2 * !pi / 180.D) IDL #Css = np.cos(S2 * np.pi / 180) Python # CosSenZ = np.Cos(senzen * np.pi / 180) # CosSolZ = np.Cos(solzen * np.pi / 180) # SinSenZ = np.Sin(senzen * np.pi / 180) # SinSolZ = np.Sin(solzen * np.pi / 180) # RelAzm = np.abs((senazm ) - (solazm)) # index = np.where(RelAzm > 180.0) # RelAzm[index] = 360.0- RelAzm[index] # index = np.where(RelAzm <= 180.0) # RelAzm[index] = 180.0- RelAzm[index] # CosRe = np.Cos(np.radians(RelAzm)) # SctAgl = np.acos((-CosSolZ * CosSenZ) + ((SinSolZ * SinSenZ) * CosRelA)) # #Ray_correction for TOA_BAND_1 and TOA_BAND_2 # #TOA_BAND_1 # B1_ROD = (0.00864 + (0.0000065)) * (0.67)**(-(3.916 + (0.074 * 0.67)+ (0.05/0.67)))#for RED band, please change 0.67 for other bands wavelength # Pr1 = (3./16.*np.pi) * (1 + (np.cos(SctAgl) * np.cos(SctAgl))) # wr1 = 1.0 # RayRef_B1 = (wr1 * B1_ROD * Pr1 ) / (4.0 * (CosSolZ * CosSenZ)) # RCR_B1 = toa_b1 - RayRef_B1 # #TOA_BAND_2 # B2_ROD = (0.00864 + (0.0000065)) * (0.86)**(-(3.916 + (0.074 * 0.86)+ (0.05/0.86)))#for Near Infrared (NIR) band, please change 0.86 for other bands wavelength # Pr2 = (3./16.*np.pi) * (1 + (np.cos(SctAgl) * np.cos(SctAgl))) # wr2 = 1.0 # RayRef_B2 = (wr2 * B2_ROD * Pr2 ) / (4.0 * (CosSolZ * CosSenZ)) # RCR_B2 = toa_b2 - RayRef_B2 #Calculating NDVI float raster # ndvi_raster = Divide(Float(Raster(toa_b2) - Raster(toa_b1)), Float(Raster(toa_b2) + Raster(toa_b1))) # ndvi = (Float(toa_b2) - toa_b1) / (Float(toa_b2) + toa_b1) print arcpy.GetMessages() print 'finished run: %s\n\n' % (datetime.datetime.now() - start) ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

I just Want to ask one question to Esri python developer, there is no option to solve my issue? ... View more

Sorry, I am not able to solve this problem. If any know the way to solve the problem, then please suggest me regarding this code. ... View more

Thanks for your suggestion. ... View more

Luke Webb‌ Sir, can I get your contact information (email address)? My email address is bijoygayen123@gmail.com ... View more

Thank you, sir, for your suggestion. But sir can explain why this error occurs in this case? ... View more

Thank you, sir,I will manipulate my code as per your suggestion. ... View more

Sorry for there some missing line, again I uploading my full code. Please check this issue... Even I did not find any solution according to my problem #------------------------------------------------------------------------------- # Name: module1 # Purpose: # # Author: HONEY # # Created: 18-06-2019 # Copyright: (c) HONEY 2019 # Licence: <your licence> #------------------------------------------------------------------------------- import datetime start = datetime.datetime.now() print 'start run: %s\n' % (start) import arcpy ,os ,sys,csv,errno from arcpy import env from arcpy.sa import * import datetime import re import numpy as np import glob import itertools arcpy.env.overwriteOutput = True outdir="F:\\DB_test_data\\TEST_RAY\\TEST1\\c\\" cellsize = "F:\\DB_test_data\\TEST_RAY\\TEST1\\MOD02HKM_A2017001_0530_NDVI_AA.img" #Call cloud mask images from directory d1="F:\\DB_test_data\\VAR2\\cldmask\\" CLDMASK = glob.glob(d1 + os.sep + "*.Aerosol_Cldmask_Land_Ocean-Aerosol_Cldmask_Land_Ocean.tif") CLDMASK.sort() if CLDMASK is None: print 'Could not open the CLDMASK raster files' sys.exit(1) else: print 'The CLDMASK raster files was opened successfully' #Call AOD images from directory d2 = r"F:\\DB_test_data\\VAR2\\MOD04_l2" AOD = glob.glob(d2 + os.sep + "*Corrected_Optical_Depth_Land_2-Corrected_Optical_Depth_Land.tif") AOD.sort() if AOD is None: print 'Could not open the AOD raster files' sys.exit(1) else: print 'The AOD raster files was opened successfully' #Call MOD02HKM_TOA-B1_TOA-B2 images from directory d3="F:\\DB_test_data\\VAR2\\IDL\\MOD02HKM\\" TOA_B1 = glob.glob(d3 + os.sep + "*RefSB_1-EV_250_Aggr500_RefSB.tif") TOA_B1.sort() if TOA_B1 is None: print 'Could not open the TOA_B1 raster files' sys.exit(1) else: print 'The TOA_B1 raster files was opened successfully' TOA_B2 = glob.glob(d3 + os.sep + "*RefSB_2-EV_250_Aggr500_RefSB.tif") TOA_B2.sort() if TOA_B2 is None: print 'Could not open the TOA_B2 raster files' sys.exit(1) else: print 'The TOA_B2 raster files was opened successfully' #Call MOD03 (Height,SenZen,SenAzm,SolZen,SolAzm) images from directory d4="F:\\DB_test_data\\VAR2\\IDL\\MOD03_R\\" Height = glob.glob(d4 + os.sep + "*.Height-Height.tif") Height.sort() if Height is None: print 'Could not open the Height raster files' sys.exit(1) else: print 'The Height raster files was opened successfully' SenZen = glob.glob(d4 + os.sep + "*.SensorZenith-SensorZenith.tif") SenZen.sort() if SenZen is None: print 'Could not open the SenZen raster files' sys.exit(1) else: print 'The SenZen raster files was opened successfully' SenAzm = glob.glob(d4 + os.sep + "*.SensorAzimuth-SensorAzimuth.tif") SenAzm.sort() if SenAzm is None: print 'Could not open the SenAzm raster files' sys.exit(1) else: print 'The SenAzm raster files was opened successfully' SolZen = glob.glob(d4 + os.sep + "*.SolarZenith-SolarZenith.tif") SolZen.sort() if SolZen is None: print 'Could not open the SolZen raster files' sys.exit(1) else: print 'The SolZen raster files was opened successfully' SolAzm = glob.glob(d4 + os.sep + "*.SolarAzimuth-SolarAzimuth.tif") SolAzm.sort() if SolAzm is None: print 'Could not open the SolAzm raster files' sys.exit(1) else: print 'The SolAzm raster files was opened successfully' #Call MOD09 images from directory d5="F:\\DB_test_data\\VAR2\\IDL\\MOD09\\" Sur_B1 = glob.glob(d5 + os.sep + "*.500m_Surface_Reflectance_Band_1-500m_Surface_Reflectance_Band_1.tif") Sur_B1 .sort() if Sur_B1 is None: print 'Could not open the Sur_B1 raster files' sys.exit(1) else: print 'The Sur_B1 raster files was opened successfully' Sur_B2 = glob.glob(d5 + os.sep + "*.500m_Surface_Reflectance_Band_2-500m_Surface_Reflectance_Band_2.tif") Sur_B2 .sort() if Sur_B2 is None: print 'Could not open the Sur_B2 raster files' sys.exit(1) else: print 'The Sur_B2 raster files was opened successfully' Sur_B3 = glob.glob(d5 + os.sep + "*.500m_Surface_Reflectance_Band_3-500m_Surface_Reflectance_Band_3.tif") Sur_B3 .sort() if Sur_B3 is None: print 'Could not open the Sur_B3 raster files' sys.exit(1) else: print 'The Sur_B3 raster files was opened successfully' Land_cover="F:\\DB_test_data\\MCD12C1\\LAND_COVER_TYPE_1_Grid_2D.img" if Land_cover is None: print 'Could not open the Land_cover raster files' sys.exit(1) else: print 'The Land_cover raster files was opened successfully' outdir="F:\\DB_test_data\\TEST_RAY\\TEST1\\c\\" cellsize = "F:\\DB_test_data\\TEST_RAY\\TEST1\\MOD02HKM_A2017001_0530_NDVI_AA.img" for a,b,c,d,e,f,g,h,i,j,k,l in zip (CLDMASK,AOD,TOA_B1,TOA_B2,Height,SenZen,SenAzm,SolZen,SolAzm,Sur_B1,Sur_B2,Sur_B3): print ("processing_a:"+ a) print ("processing_b:"+ b) print ("processing_c:"+ c) print ("processing_d:"+ d) print ("processing_e:"+ e) print ("processing_f:"+ f) print ("processing_g:"+ g) print ("processing_h:"+ h) print ("processing_i:"+ i) print ("processing_j:"+ j) print ("processing_k:"+ k) print ("processing_l:"+ l) ####################################################################### #Cloud Mask scalling Scale_factor1 = float(1.0) add_offset1 = float(0.0) setnull1 =arcpy.gp.SetNull_sa(a,a, "in_memory/dat1", "\"Value\" <= 0") ras1=arcpy.Raster(setnull1) cldmask=(ras1-add_offset1)*Scale_factor1 OutRaster = os.path.join(outdir,'MOd0.{0}.img'.format("cldmask")) #print OutRaster arcpy.CopyRaster_management(cldmask,OutRaster) arcpy.Delete_management("in_memory/dat1") #AOD scalling and set lessthan 0.1 AOD Scale_factor2 = float(0.0010000000474974513) add_offset2 = float(0.0) setnull2 =arcpy.gp.SetNull_sa(b,b, "in_memory/dat2", "\"Value\" = -9999") ras2=arcpy.Raster(setnull2) aod=(ras2-add_offset2)*Scale_factor2 aod_1 = Con((aod >= 0.0) & (aod <= 0.1),1) OutRaster = os.path.join(outdir,'MOd0.{0}.img'.format("aod_1")) #print OutRaster arcpy.CopyRaster_management(aod_1,OutRaster) arcpy.Delete_management("in_memory/dat2") #TOA_B1 scalling Scale_factor3 = float(0.000053811826) add_offset3 = float(-0.0) setnull3 =arcpy.gp.SetNull_sa(c,c, "in_memory/dat3", "\"Value\" = 65535") ras3=arcpy.Raster(setnull3) toa_b1=(ras3-(add_offset3))*Scale_factor3 OutRaster = os.path.join(outdir,'MOd0.{0}.img'.format("toa_b1")) #print OutRaster arcpy.CopyRaster_management(toa_b1,OutRaster) arcpy.Delete_management("in_memory/dat3") #TOA_B2 Scalling Scale_factor4 = float(0.00003255546) add_offset4 = float(-0.0) setnull4 =arcpy.gp.SetNull_sa(d,d, "in_memory/dat4", "\"Value\" = 65535") ras4=arcpy.Raster(setnull4) toa_b2=(ras4-(add_offset4))*Scale_factor4 OutRaster = os.path.join(outdir,'MOd0.{0}.img'.format("toa_b2")) #print OutRaster arcpy.CopyRaster_management(toa_b2,OutRaster) arcpy.Delete_management("in_memory/dat4") #Height setnull =arcpy.gp.SetNull_sa(e,e, "in_memory/dat5", "\"Value\" = -32767") ras=arcpy.Raster(setnull) height=ras OutRaster = os.path.join(outdir,'MOd0.{0}.img'.format("height")) #print OutRaster arcpy.CopyRaster_management(height,OutRaster) arcpy.Delete_management("in_memory/dat5") #Sensor Zenith angle Scale_factor6 = float(0.01) add_offset6 = float(0.0) setnull6 =arcpy.gp.SetNull_sa(f,f, "in_memory/dat6", "\"Value\" = -32767") ras6=arcpy.Raster(setnull6) senzen=(ras6-add_offset6)*Scale_factor6 OutRaster = os.path.join(outdir,'MOd0.{0}.img'.format("senzen")) #print OutRaster arcpy.CopyRaster_management(senzen,OutRaster) arcpy.Delete_management("in_memory/dat6") #Sensor azimuth angle setnull7 =arcpy.gp.SetNull_sa(g,g, "in_memory/dat7", "\"Value\" = -32767") ras7=arcpy.Raster(setnull7) senazm=(ras7-add_offset6)*Scale_factor6 OutRaster = os.path.join(outdir,'MOd0.{0}.img'.format("senazm")) #print OutRaster arcpy.CopyRaster_management(senazm,OutRaster) arcpy.Delete_management("in_memory/dat7") #solar zenithal angle setnull8 =arcpy.gp.SetNull_sa(h,h, "in_memory/dat8", "\"Value\" = -32767") ras8=arcpy.Raster(setnull8) solzen=(ras8-add_offset6)*Scale_factor6 OutRaster = os.path.join(outdir,'MOd0.{0}.img'.format("solzen")) #print OutRaster arcpy.CopyRaster_management(solzen,OutRaster) arcpy.Delete_management("in_memory/dat8") #solar azimuth angle setnull9 =arcpy.gp.SetNull_sa(i,i, "in_memory/dat9", "\"Value\" = -32767") ras9=arcpy.Raster(setnull9) solazm=(ras9-add_offset6)*Scale_factor6 OutRaster = os.path.join(outdir,'MOd0.{0}.img'.format("solazm")) #print OutRaster arcpy.CopyRaster_management(solazm,OutRaster) arcpy.Delete_management("in_memory/dat9") #(MODO9)surface reflectance band 1 Scale_factor_10 = float(0.0001) add_offset_10 = float(0.0) setnull_10 =arcpy.gp.SetNull_sa(j,j, "in_memory/dat10", "\"Value\" = -28672") ras_10=arcpy.Raster(setnull_10) sur_b1=(ras_10-add_offset_10)*Scale_factor_10 OutRaster = os.path.join(outdir,'MOd0.{0}.img'.format("sur_b1")) #print OutRaster arcpy.CopyRaster_management(sur_b1,OutRaster) arcpy.Delete_management("in_memory/dat10") #(MODO9)surface reflectance band 2 setnull_11 =arcpy.gp.SetNull_sa(k,k, "in_memory/dat11", "\"Value\" = -28672") ras_11=arcpy.Raster(setnull_11) sur_b2=(ras_11-add_offset_10)*Scale_factor_10 OutRaster = os.path.join(outdir,'MOd0.{0}.img'.format("sur_b2")) #print OutRaster arcpy.CopyRaster_management(sur_b2,OutRaster) arcpy.Delete_management("in_memory/dat11") #(MODO9)surface reflectance band 3 setnull_12 =arcpy.gp.SetNull_sa(l,l, "in_memory/dat12", "\"Value\" = -28672") ras_12=arcpy.Raster(setnull_12) sur_b3=(ras_12- add_offset_10)*Scale_factor_10 OutRaster = os.path.join(outdir,'MOd0.{0}.img'.format("sur_b3")) #print OutRaster arcpy.CopyRaster_management(sur_b3,OutRaster) arcpy.Delete_management("in_memory/dat12")‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Thank you Luke Webb sir, your code working well as per my requirement. And also thank full to other persons who continuously helped me. ... View more

Luke Webb sir, Have there any problem within my code which I wrote earlier? ... View more

 Luke Webb sir, are you checking my problem? I am waiting for your feedback. ... View more

 Luke Webb sir, I just test your function and it works well in  but the arguments  minTolerance3 = 0.5 minType3 = ge #Greater Equal maxTolerance3 = None maxType3 = None‍‍‍‍ elif minCheckType == None: #query with just max being checked result = maxCheckType(row[0], maxValue) elif maxCheckType == None: #Query just checking minimums result = minCheckType(row[0], minValue) ‍‍‍‍‍ it does not work as my requirement (I need row[0] => 0.50 ) but it gives all < 0.5 > value and my code is def geoprocess (file_name,Input_Data,base_path,category,land_use,minValue, minCheckType, maxValue, maxCheckType): outdir = os.path.join(base_path,category,land_use) if not os.path.exists(outdir): try: print('\nCreating new output directory!\n'+outdir+'\n') os.makedirs(outdir) except OSError: print ("Creation of the directory failed" ) else: print ("Successfully created the directory" ) ras_name = os.path.join(outdir,("{0}_{1}{2}_{3}".format(category,file_name,land_use,".txt"))) fields = ['grid_code', 'Band1', 'Band2','Band3','Name',"Name2"] with open(ras_name, 'wb') as txtfile: arcpy.AddMessage("fields Writing For:{0}".format(land_use)) txtfile.write(" {0} {1} {2} {3} {4} {5}\n".format(fields[0],fields[1],fields[2],fields[3],fields[4],fields[5])) with arcpy.da.UpdateCursor(Input_Data, fields) as cursor: for row in cursor: if minCheckType == None and maxCheckType == None: #We have no comparison! result = "No comparison to make" elif maxCheckType == None: #Query just checking minimums result = minCheckType(row[0], minValue) row[5] = land_use print ("processing_A:""{0:.8f} {1:.8f} {2:.8f} {3:.8f} {4} {5}".format(row[0],row[1],row[2], row[3],row[4],row[5])) txtfile.write(" {0:.8f} {1:.8f} {2:.8f} {3:.8f} {4} {5}\n".format(row[0],row[1],row[2], row[3],row[4],row[5])) elif minCheckType == None: #query with just max being checked result = maxCheckType(row[0], maxValue) row[5] = land_use print ("processing_B:""{0:.8f} {1:.8f} {2:.8f} {3:.8f} {4} {5}".format(row[0],row[1],row[2], row[3],row[4],row[5])) txtfile.write(" {0:.8f} {1:.8f} {2:.8f} {3:.8f} {4} {5}\n".format(row[0],row[1],row[2], row[3],row[4],row[5])) else: result1 = minCheckType(row[0], minValue) result2 = maxCheckType(row[0], maxValue) if result1 == True and result2 == True: row[5] = land_use print ("processing_C:""{0:.8f} {1:.8f} {2:.8f} {3:.8f} {4} {5}".format(row[0],row[1],row[2], row[3],row[4],row[5])) txtfile.write(" {0:.8f} {1:.8f} {2:.8f} {3:.8f} {4} {5}\n".format(row[0],row[1],row[2], row[3],row[4],row[5])) else: result1 = False txtfile.close() print ("Successfully save the:{0}".format(land_use) ) print "Code is OK" except: print "Error in code" def main(file_name,Input_Data,base_path,category): minTolerance1= 0.10 minType1 = gt #Greater than maxTolerance1 = 0.20 maxType1 = lt #Less than minTolerance2 = 0.2 minType2 = ge #Greater Equal maxTolerance2 = 0.5 maxType2 = lt #Less than minTolerance3 = 0.5 minType3 = ge #Greater Equal maxTolerance3 = None maxType3 = None geoprocess(file_name,Input_Data,base_path,category,'LAND_10',minTolerance1, minType1,maxTolerance1, maxType1) geoprocess(file_name,Input_Data,base_path,category, 'LAND_20',minTolerance2, minType2, maxTolerance2, maxType2) geoprocess(file_name,Input_Data,base_path,category, 'LAND_30', minTolerance3, minType3, maxTolerance3, maxType3) try: if __name__ == '__main__': t1=time.time() main(file_name,AAA,base_path,"OOO") main(file_name,BBB,base_path,"FFF") main(file_name,CCC,base_path,"ddd") print "serial processing took:",time.time()-t1 print "serial processing is OK" except: print "Error in serial processing" print arcpy.GetMessages() print 'finished run: %s\n\n' % (datetime.datetime.now() - start) ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ please check this issue. ... View more

Thank you, Luke Webb sir. I will try to manipulate my code as per your instruction. I think it will work fine. ... View more