I am trying to write one code where I use set null function for different images and also use a simple math function for scaling the images. But when I run this code after three processes it shows Traceback (most recent call last): File "F:\DB_test_data\python_script\FINAL_DB.py", line 192, in arcpy.CopyRaster_management(toa_b2,OutRaster) File "C:\Program Files (x86)\ArcGIS\Desktop10.6\ArcPy\arcpy\management.py", line 14587, in CopyRaster raise e RuntimeError: ERROR 999998: Unexpected Error. here is my code outdir = "F:\\DB_test_data\\TEST_RAY\\TEST1\\c\\" 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\" ) 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" ) ############################################################################# # # 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 # 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 ) I can't understand why it shows this error. please help me, anyone, to resolve this matter. And have any way to use a "foor loop" for set null function of different images? Xander Bakker Luke Webb, Joshua Bixby Joe Borgione Dan Patterson Test_data.rar

Hello Xander Bakker sir, sorry for the late reply. I am confused to write the code using your mentioned function. Can you write one example code with my code?

Hi BIJOY GAYEN , There is no option, without conversion (Raster To Numpy and then Numpy To Raster) to write the code for raster data analysis directly use cos, sin, subtract, plus, minus, division function, etc. I am totally confused to write this code. Give any solution regarding this issue Have a look at these functions that are available in Spatial Analyst: * http://desktop.arcgis.com/en/arcmap/latest/tools/spatial-analyst-toolbox/cos.htm * http://desktop.arcgis.com/en/arcmap/latest/tools/spatial-analyst-toolbox/sin.htm * http://desktop.arcgis.com/en/arcmap/latest/tools/spatial-analyst-toolbox/minus.htm * http://desktop.arcgis.com/en/arcmap/latest/tools/spatial-analyst-toolbox/plus.htm * http://desktop.arcgis.com/en/arcmap/latest/tools/spatial-analyst-toolbox/divide.htm * http://desktop.arcgis.com/en/arcmap/latest/tools/spatial-analyst-toolbox/power.htm

Hello Xander Bakker sir, I rechecked my code as per your suggestion, But when I try to save the output SctAgl and toa_b1 after some processes this code gave me some error in output (the output gave me blank result with unwanted value range in two images) SctAgl: Toa_B1: Second, when I run this code it shows some unwanted messages, I cannot understand why these messages occur F : \DB_test_data\python_script\FINAL_DB . py : 298 : RuntimeWarning : overflow encountered in multiply cosSenZ = np . cos ( senzen_na * np . pi / 180 ) F : \DB_test_data\python_script\FINAL_DB . py : 298 : RuntimeWarning : invalid value encountered in cos cosSenZ = np . cos ( senzen_na * np . pi / 180 ) F : \DB_test_data\python_script\FINAL_DB . py : 300 : RuntimeWarning : overflow encountered in multiply cosSolZ = np . cos ( cosSolZ_na * np . pi / 180 ) F : \DB_test_data\python_script\FINAL_DB . py : 300 : RuntimeWarning : invalid value encountered in cos cosSolZ = np . cos ( cosSolZ_na * np . pi / 180 ) F : \DB_test_data\python_script\FINAL_DB . py : 302 : RuntimeWarning : overflow encountered in multiply sinSenZ = np . sin ( sinSenZ_na * np . pi / 180 ) F : \DB_test_data\python_script\FINAL_DB . py : 302 : RuntimeWarning : invalid value encountered in sin sinSenZ = np . sin ( sinSenZ_na * np . pi / 180 ) F : \DB_test_data\python_script\FINAL_DB . py : 304 : RuntimeWarning : overflow encountered in multiply sinSolZ = np . sin ( sinSolZ_na * np . pi / 180 ) F : \DB_test_data\python_script\FINAL_DB . py : 304 : RuntimeWarning : invalid value encountered in sin sinSolZ = np . sin ( sinSolZ_na * np . pi / 180 ) Third things, when I try to run this code with row data (whole India) after some processing it shows memory issue sinSenZ_na = arcpy . RasterToNumPyArray ( senzen ) File "C:\Program Files (x86)\ArcGIS\Desktop10.6\ArcPy\arcpy\__init__.py" , line 2275 , in RasterToNumPyArray return _RasterToNumPyArray ( * args , ** kwargs ) MemoryError There is no option, without conversion (Raster To Numpy and then Numpy To Raster) to write the code for raster data analysis directly use cos, sin, subtract, plus, minus, division function, etc. I am totally confused to write this code. Give any solution regarding this issue. Dan Patterson Joe BorgioneJoshua BixbyLuke Webb My_full_code: #------------------------------------------------------------------------------- # Name: module1 # Purpose: # # Author: HONEY # # Created: 18-06-2019 # Copyright: (c) HONEY 2019 # 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 glob import itertools import math import numpy as np arcpy . CheckOutExtension ( 'Spatial' ) 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 i in range ( len ( CLDMASK ) ) : AA = CLDMASK [ i ] . split ( "\\" ) [ 5 ] [ 0 : 114 ] print ( "processing_a:" + AA ) BB = AOD [ i ] . split ( "\\" ) [ 9 ] [ 0 : 120 ] print ( "processing_b:" + BB ) CC = TOA_B1 [ i ] . split ( "\\" ) [ 5 ] [ 0 : 88 ] print ( "processing_c:" + CC ) DD = TOA_B2 [ i ] . split ( "\\" ) [ 5 ] [ 0 : 88 ] print ( "processing_d:" + DD ) EE = Height [ i ] . split ( "\\" ) [ 5 ] [ 0 : 71 ] print ( "processing_e:" + EE ) FF = SenZen [ i ] . split ( "\\" ) [ 5 ] [ 0 : 83 ] print ( "processing_f:" + FF ) GG = SenAzm [ i ] . split ( "\\" ) [ 5 ] [ 0 : 85 ] print ( "processing_g:" + GG ) HH = SolZen [ i ] . split ( "\\" ) [ 5 ] [ 0 : 81 ] print ( "processing_h:" + HH ) II = SolAzm [ i ] . split ( "\\" ) [ 5 ] [ 0 : 83 ] print ( "processing_i:" + II ) JJ = Sur_B1 [ i ] . split ( "\\" ) [ 5 ] [ 0 : 121 ] print ( "processing_j:" + JJ ) KK = Sur_B2 [ i ] . split ( "\\" ) [ 5 ] [ 0 : 121 ] print ( "processing_k:" + KK ) LL = Sur_B3 [ i ] . split ( "\\" ) [ 5 ] [ 0 : 121 ] print ( "processing_l:" + LL ) ####################################################################### #Cloud Mask scalling Scale_factor1 = float ( 1.0 ) add_offset1 = float ( 0.0 ) whereClause = "VALUE # Execute SetNull outSetNull = SetNull ( CLDMASK [ i ] , CLDMASK [ i ] , whereClause ) cldmask = ( outSetNull - add_offset1 ) * Scale_factor1 #AOD scalling Scale_factor2 = float ( 0.0010000000474974513 ) add_offset2 = float ( 0.0 ) whereClause = "VALUE = -9999" # Execute SetNull outSetNull2 = SetNull ( AOD [ i ] , AOD [ i ] , whereClause ) aod = ( outSetNull2 - add_offset2 ) * Scale_factor2 aod_1 = Con ( ( aod >= 0.0 ) & ( aod 0.1 ) , 1 ) #TOA_B1 Scalling Scale_factor3 = float ( 0.000053811826 ) add_offset3 = float ( - 0.0 ) whereClause = "VALUE = 65535" # Execute SetNull outSetNull3 = SetNull ( TOA_B1 [ i ] , TOA_B1 [ i ] , whereClause ) toa_b1 = ( outSetNull3 - add_offset3 ) * Scale_factor3 #TOA_B2 Scalling Scale_factor4 = float ( 0.00003255546 ) add_offset4 = float ( - 0.0 ) whereClause = "VALUE = 65535" # Execute SetNull outSetNull4 = SetNull ( TOA_B2 [ i ] , TOA_B2 [ i ] , whereClause ) toa_b2 = ( outSetNull4 - add_offset4 ) * Scale_factor4 #Height Scale_factor5 = float ( 1.0 ) add_offset5 = float ( - 0.0 ) whereClause = "VALUE = -32767" # Execute SetNull outSetNull5 = SetNull ( Height [ i ] , Height [ i ] , whereClause ) height = ( outSetNull5 - add_offset5 ) * Scale_factor5 #Sensor Zenith angle Scale_factor6 = float ( 0.01 ) add_offset6 = float ( 0.0 ) whereClause = "VALUE = -32767" # Execute SetNull outSetNull6 = SetNull ( SenZen [ i ] , SenZen [ i ] , whereClause ) senzen = ( outSetNull6 - add_offset6 ) * Scale_factor6 #Sensor azimuth angle Scale_factor7 = float ( 0.01 ) add_offset7 = float ( 0.0 ) whereClause = "VALUE = -32767" # Execute SetNull outSetNull7 = SetNull ( SenAzm [ i ] , SenAzm [ i ] , whereClause ) senazm = ( outSetNull7 - add_offset7 ) * Scale_factor7 #solar Zenith angle Scale_factor8 = float ( 0.01 ) add_offset8 = float ( 0.0 ) whereClause = "VALUE = -32767" # Execute SetNull outSetNull8 = SetNull ( SolZen [ i ] , SolZen [ i ] , whereClause ) solzen = ( outSetNull8 - add_offset8 ) * Scale_factor8 #solar azimuth angle Scale_factor9 = float ( 0.01 ) add_offset9 = float ( 0.0 ) whereClause = "VALUE = -32767" # Execute SetNull outSetNull9 = SetNull ( SolAzm [ i ] , SolAzm [ i ] , whereClause ) solazm = ( outSetNull9 - add_offset9 ) * Scale_factor9 # OutRaster = os.path.join(outdir,'MOD04_L22.{0}.tif'.format("solazm")) # solazm.save(OutRaster) #(MODO9)surface reflectance band 1 Scale_factor_10 = float ( 0.0001 ) add_offset_10 = float ( 0.0 ) whereClause = "VALUE = -28672" # Execute SetNull outSetNull10 = SetNull ( Sur_B1 [ i ] , Sur_B1 [ i ] , whereClause ) sur_b1 = ( outSetNull10 - add_offset_10 ) * Scale_factor_10 #(MODO9)surface reflectance band 2 Scale_factor_11 = float ( 0.0001 ) add_offset_11 = float ( 0.0 ) whereClause = "VALUE = -28672" # Execute SetNull outSetNull11 = SetNull ( Sur_B2 [ i ] , Sur_B2 [ i ] , whereClause ) sur_b2 = ( outSetNull11 - add_offset_11 ) * Scale_factor_11 #(MODO9)surface reflectance band 3 Scale_factor_12 = float ( 0.0001 ) add_offset_12 = float ( 0.0 ) whereClause = "VALUE = -28672" # Execute SetNull outSetNull12 = SetNull ( Sur_B3 [ i ] , Sur_B3 [ i ] , whereClause ) sur_b3 = ( outSetNull12 - add_offset_12 ) * Scale_factor_12 # Process: Extract by Cloud_Mask tempEnvironment0 = arcpy . env . cellSize arcpy . env . cellSize = "MAXOF" toa_bb1 = arcpy . sa . ExtractByMask ( toa_b1 , cldmask ) toa_bb2 = arcpy . sa . ExtractByMask ( toa_b2 , cldmask ) heightt = arcpy . sa . ExtractByMask ( height , cldmask ) senzenn = arcpy . sa . ExtractByMask ( senzen , cldmask ) senazmm = arcpy . sa . ExtractByMask ( senazm , cldmask ) solzenn = arcpy . sa . ExtractByMask ( solzen , cldmask ) solazmm = arcpy . sa . ExtractByMask ( solazm , cldmask ) sur_b11 = arcpy . sa . ExtractByMask ( sur_b1 , cldmask ) sur_b22 = arcpy . sa . ExtractByMask ( sur_b2 , cldmask ) sur_b33 = arcpy . sa . ExtractByMask ( 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 . sa . ExtractByMask ( toa_bb1 , aod_1 ) toa_b2 = arcpy . sa . ExtractByMask ( toa_bb2 , aod_1 ) height = arcpy . sa . ExtractByMask ( heightt , aod_1 ) senzen = arcpy . sa . ExtractByMask ( senzenn , aod_1 ) senazm = arcpy . sa . ExtractByMask ( senazmm , aod_1 ) solzen = arcpy . sa . ExtractByMask ( solzenn , aod_1 ) solazm = arcpy . sa . ExtractByMask ( solazmm , aod_1 ) sur_b1 = arcpy . sa . ExtractByMask ( sur_b11 , aod_1 ) sur_b2 = arcpy . sa . ExtractByMask ( sur_b22 , aod_1 ) sur_b3 = arcpy . sa . ExtractByMask ( sur_b33 , aod_1 ) arcpy . env . cellSize = tempEnvironment0 ############################################################################# #Calculate the scattering angle from Resample MOD03 products senzen_na = arcpy . RasterToNumPyArray ( senzen ) cosSenZ = np . cos ( senzen_na * np . pi / 180 ) cosSolZ_na = arcpy . RasterToNumPyArray ( solzen ) cosSolZ = np . cos ( cosSolZ_na * np . pi / 180 ) sinSenZ_na = arcpy . RasterToNumPyArray ( senzen ) sinSenZ = np . sin ( sinSenZ_na * np . pi / 180 ) sinSolZ_na = arcpy . RasterToNumPyArray ( solzen ) sinSolZ = np . sin ( sinSolZ_na * np . pi / 180 ) senazm_na = arcpy . RasterToNumPyArray ( senazm ) solazm_na = arcpy . RasterToNumPyArray ( solazm ) RelAzm = np . abs ( ( senazm_na ) - ( solazm_na ) ) index = np . where ( RelAzm > 180.0 ) RelAzm [ index ] = 360.0 - RelAzm [ index ] index = np . where ( RelAzm 180.0 ) RelAzm [ index ] = 180.0 - RelAzm [ index ] cosRelA = np . cos ( np . radians ( RelAzm ) ) SctAgl = np . arccos ( ( - cosSolZ * cosSenZ ) + ( ( sinSolZ * sinSenZ ) * cosRelA ) ) dscc = arcpy . Describe ( senazm ) srr = dscc . SpatialReference extt = dscc . Extent lll = arcpy . Point ( extt . XMin , extt . YMin ) noDataValuee = dscc . noDataValue myArraySumm = SctAgl . sum ( 1 ) myArraySumm . shape = ( SctAgl . shape [ 0 ] , 1 ) myArrayPercc = ( SctAgl * 1.0 ) / myArraySumm newRasterr = arcpy . NumPyArrayToRaster ( myArrayPercc , lll , dscc . meanCellWidth , dscc . meanCellHeight , noDataValuee ) arcpy . DefineProjection_management ( newRasterr , srr ) OutRaster = os . path . join ( outdir , 'MOD022222.{0}.tif' . format ( "SctAgl" ) ) newRasterr . save ( OutRaster ) # # RCR for TOA_BAND_1 and TOA_BAND_2 # TOA_BAND_1 dsc = arcpy . Describe ( toa_b1 ) sr = dsc . SpatialReference ext = dsc . Extent ll = arcpy . Point ( ext . XMin , ext . YMin ) noDataValue = dsc . noDataValue toa_b1_na = arcpy . RasterToNumPyArray ( toa_b1 ) B1_ROD = ( 0.00864 + ( 0.0000065 ) ) * ( 0.67 ) ** ( - ( 3.916 + ( 0.074 * 0.67 ) + ( 0.05 / 0.67 ) ) ) Pr1 = ( 3.0 / 16.0 * np . pi ) * ( 1 + ( np . cos ( SctAgl ) * np . cos ( SctAgl ) ) ) wr1 = 1.0 RayRef_B1 = ( wr1 * B1_ROD * Pr1 ) / ( 4.0 * ( cosSolZ * cosSenZ ) ) myArray = toa_b1_na - RayRef_B1 myArraySum = myArray . sum ( 1 ) myArraySum . shape = ( myArray . shape [ 0 ] , 1 ) myArrayPerc = ( myArray * 1.0 ) / myArraySum newRaster = arcpy . NumPyArrayToRaster ( myArrayPerc , ll , dsc . meanCellWidth , dsc . meanCellHeight , noDataValue ) arcpy . DefineProjection_management ( newRaster , sr ) OutRaster = os . path . join ( outdir , 'MOD022222.{0}.tif' . format ( "TOA_B1" ) ) # toa_b2.save(OutRaster) newRaster . save ( OutRaster ) print arcpy . GetMessages ( ) print 'finished run: %s\n\n' % ( datetime . datetime . now ( ) - start )

nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

RuntimeError: ERROR 999998: Unexpected Error in Set null function

BIJOYGAYEN

I am trying to write one code where I use set null function for different images and also use a simple math function for scaling the images. But when I run this code after three processes it shows

Traceback (most recent call last):
  File "F:\DB_test_data\python_script\FINAL_DB.py", line 192, in <module>
    arcpy.CopyRaster_management(toa_b2,OutRaster)
  File "C:\Program Files (x86)\ArcGIS\Desktop10.6\ArcPy\arcpy\management.py", line 14587, in CopyRaster
    raise e
RuntimeError: ERROR 999998: Unexpected Error.‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

here is my code

outdir<SPAN class="operator token">=</SPAN><SPAN class="string token">"F:\\DB_test_data\\TEST_RAY\\TEST1\\c\\"</SPAN>
<SPAN class="keyword token">for</SPAN> a<SPAN class="punctuation token">,</SPAN>b<SPAN class="punctuation token">,</SPAN>c<SPAN class="punctuation token">,</SPAN>d<SPAN class="punctuation token">,</SPAN>e<SPAN class="punctuation token">,</SPAN>f<SPAN class="punctuation token">,</SPAN>g<SPAN class="punctuation token">,</SPAN>h<SPAN class="punctuation token">,</SPAN>i<SPAN class="punctuation token">,</SPAN>j<SPAN class="punctuation token">,</SPAN>k<SPAN class="punctuation token">,</SPAN>l <SPAN class="keyword token">in</SPAN> zip <SPAN class="punctuation token">(</SPAN>CLDMASK<SPAN class="punctuation token">,</SPAN>AOD<SPAN class="punctuation token">,</SPAN>TOA_B1<SPAN class="punctuation token">,</SPAN>TOA_B2<SPAN class="punctuation token">,</SPAN>Height<SPAN class="punctuation token">,</SPAN>SenZen<SPAN class="punctuation token">,</SPAN>SenAzm<SPAN class="punctuation token">,</SPAN>SolZen<SPAN class="punctuation token">,</SPAN>SolAzm<SPAN class="punctuation token">,</SPAN>Sur_B1<SPAN class="punctuation token">,</SPAN>Sur_B2<SPAN class="punctuation token">,</SPAN>Sur_B3<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_a:"</SPAN><SPAN class="operator token">+</SPAN> a<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_b:"</SPAN><SPAN class="operator token">+</SPAN> b<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_c:"</SPAN><SPAN class="operator token">+</SPAN> c<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_d:"</SPAN><SPAN class="operator token">+</SPAN> d<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_e:"</SPAN><SPAN class="operator token">+</SPAN> e<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_f:"</SPAN><SPAN class="operator token">+</SPAN> f<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_g:"</SPAN><SPAN class="operator token">+</SPAN> g<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_h:"</SPAN><SPAN class="operator token">+</SPAN> h<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_i:"</SPAN><SPAN class="operator token">+</SPAN> i<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_j:"</SPAN><SPAN class="operator token">+</SPAN> j<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_k:"</SPAN><SPAN class="operator token">+</SPAN> k<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_l:"</SPAN><SPAN class="operator token">+</SPAN> l<SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#######################################################################</SPAN>
    <SPAN class="comment token">#Cloud Mask scalling</SPAN>
    Scale_factor1 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1.0</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset1 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull1 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>a<SPAN class="punctuation token">,</SPAN>a<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat1"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" <= 0"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras1<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull1<SPAN class="punctuation token">)</SPAN>
    cldmask<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras1<SPAN class="operator token">-</SPAN>add_offset1<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor1
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"cldmask"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>cldmask<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat1"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#AOD scalling and set lessthan 0.1 AOD </SPAN>
    Scale_factor2 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0010000000474974513</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset2 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull2 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>b<SPAN class="punctuation token">,</SPAN>b<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat2"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -9999"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras2<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull2<SPAN class="punctuation token">)</SPAN>
    aod<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras2<SPAN class="operator token">-</SPAN>add_offset2<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor2
    aod_1 <SPAN class="operator token">=</SPAN> Con<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">(</SPAN>aod <SPAN class="operator token">>=</SPAN> <SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">&</SPAN> <SPAN class="punctuation token">(</SPAN>aod <SPAN class="operator token"><=</SPAN> <SPAN class="number token">0.1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"aod_1"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>aod_1<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat2"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#TOA_B1 scalling</SPAN>
    Scale_factor3 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.000053811826</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset3 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="operator token">-</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull3 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>c<SPAN class="punctuation token">,</SPAN>c<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat3"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = 65535"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras3<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull3<SPAN class="punctuation token">)</SPAN>
    toa_b1<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras3<SPAN class="operator token">-</SPAN><SPAN class="punctuation token">(</SPAN>add_offset3<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor3
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"toa_b1"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>toa_b1<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat3"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#TOA_B2 Scalling</SPAN>
    Scale_factor4 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.00003255546</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset4 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="operator token">-</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull4 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>d<SPAN class="punctuation token">,</SPAN>d<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat4"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = 65535"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras4<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull4<SPAN class="punctuation token">)</SPAN>
    toa_b2<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras4<SPAN class="operator token">-</SPAN><SPAN class="punctuation token">(</SPAN>add_offset4<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor4
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"toa_b2"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>toa_b2<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat4"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#Height</SPAN>
    setnull <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>e<SPAN class="punctuation token">,</SPAN>e<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat5"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -32767"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull<SPAN class="punctuation token">)</SPAN>
    height<SPAN class="operator token">=</SPAN>ras
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"height"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>height<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat5"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#Sensor Zenith angle</SPAN>
    Scale_factor6 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.01</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset6 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull6 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>f<SPAN class="punctuation token">,</SPAN>f<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat6"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -32767"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras6<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull6<SPAN class="punctuation token">)</SPAN>
    senzen<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras6<SPAN class="operator token">-</SPAN>add_offset6<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor6
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"senzen"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>senzen<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat6"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#Sensor azimuth angle</SPAN>
    setnull7 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>g<SPAN class="punctuation token">,</SPAN>g<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat7"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -32767"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras7<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull7<SPAN class="punctuation token">)</SPAN>
    senazm<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras7<SPAN class="operator token">-</SPAN>add_offset6<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor6
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"senazm"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>senazm<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat7"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#solar zenithal angle</SPAN>
    setnull8 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>h<SPAN class="punctuation token">,</SPAN>h<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat8"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -32767"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras8<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull8<SPAN class="punctuation token">)</SPAN>
    solzen<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras8<SPAN class="operator token">-</SPAN>add_offset6<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor6
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"solzen"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>solzen<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat8"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#solar azimuth angle</SPAN>
    setnull9 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>i<SPAN class="punctuation token">,</SPAN>i<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat9"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -32767"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras9<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull9<SPAN class="punctuation token">)</SPAN>
    solazm<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras9<SPAN class="operator token">-</SPAN>add_offset6<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor6
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"solazm"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>solazm<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat9"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#(MODO9)surface reflectance band 1</SPAN>
    Scale_factor_10 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0001</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset_10 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull_10 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>j<SPAN class="punctuation token">,</SPAN>j<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat10"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -28672"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras_10<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull_10<SPAN class="punctuation token">)</SPAN>
    sur_b1<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras_10<SPAN class="operator token">-</SPAN>add_offset_10<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor_10
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"sur_b1"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>sur_b1<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat10"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#(MODO9)surface reflectance band 2</SPAN>
    setnull_11 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>k<SPAN class="punctuation token">,</SPAN>k<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat11"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -28672"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras_11<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull_11<SPAN class="punctuation token">)</SPAN>
    sur_b2<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras_11<SPAN class="operator token">-</SPAN>add_offset_10<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor_10
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"sur_b2"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>sur_b2<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat11"</SPAN><SPAN class="punctuation token">)</SPAN>
   

    <SPAN class="comment token">#(MODO9)surface reflectance band 3</SPAN>
    setnull_12 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>l<SPAN class="punctuation token">,</SPAN>l<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat12"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -28672"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras_12<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull_12<SPAN class="punctuation token">)</SPAN>
    sur_b3<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras_12<SPAN class="operator token">-</SPAN> add_offset_10<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor_10
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"sur_b3"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>sur_b3<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat12"</SPAN><SPAN class="punctuation token">)</SPAN>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍
<SPAN class="comment token">#############################################################################</SPAN>
    <SPAN class="comment token"># # Process: Extract by Cloud_Mask</SPAN>
    <SPAN class="comment token"># tempEnvironment0 = arcpy.env.cellSize</SPAN>
    <SPAN class="comment token"># arcpy.env.cellSize = "MAXOF"</SPAN>
    <SPAN class="comment token"># toa_bb1=arcpy.gp.ExtractByMask_sa(toa_b1, cldmask)</SPAN>
    <SPAN class="comment token"># toa_b22=arcpy.gp.ExtractByMask_sa(toa_b2, cldmask)</SPAN>
    <SPAN class="comment token"># heightt=arcpy.gp.ExtractByMask_sa(height, cldmask)</SPAN>
    <SPAN class="comment token"># senzenn=arcpy.gp.ExtractByMask_sa(senzen, cldmask)</SPAN>
    <SPAN class="comment token"># senazmm=arcpy.gp.ExtractByMask_sa(senazm, cldmask)</SPAN>
    <SPAN class="comment token"># solzenn=arcpy.gp.ExtractByMask_sa(solzen, cldmask)</SPAN>
    <SPAN class="comment token"># solazmm=arcpy.gp.ExtractByMask_sa(solazm, cldmask)</SPAN>
    <SPAN class="comment token"># sur_b11=arcpy.gp.ExtractByMask_sa(sur_b1, cldmask)</SPAN>
    <SPAN class="comment token"># sur_b22=arcpy.gp.ExtractByMask_sa(sur_b2, cldmask)</SPAN>
    <SPAN class="comment token"># sur_b33=arcpy.gp.ExtractByMask_sa(sur_b3, cldmask)</SPAN>
    <SPAN class="comment token"># arcpy.env.cellSize = tempEnvironment0</SPAN>
    <SPAN class="comment token"># # Process: Extract by Mask using AOD less than 0.1 value</SPAN>
    <SPAN class="comment token"># tempEnvironment0 = arcpy.env.cellSize</SPAN>
    <SPAN class="comment token"># arcpy.env.cellSize = cellsize</SPAN>
    <SPAN class="comment token"># toa_b1=arcpy.gp.ExtractByMask_sa(toa_bb1, aod_1)</SPAN>
    <SPAN class="comment token"># toa_b2=arcpy.gp.ExtractByMask_sa(toa_bb2, aod_1)</SPAN>
    <SPAN class="comment token"># height=arcpy.gp.ExtractByMask_sa(heightt, aod_1)</SPAN>
    <SPAN class="comment token"># senzen=arcpy.gp.ExtractByMask_sa(senzenn, aod_1)</SPAN>
    <SPAN class="comment token"># senazm=arcpy.gp.ExtractByMask_sa(senazmm, aod_1)</SPAN>
    <SPAN class="comment token"># solzen=arcpy.gp.ExtractByMask_sa(solzenn, aod_1)</SPAN>
    <SPAN class="comment token"># solazm=arcpy.gp.ExtractByMask_sa(solazmm, aod_1)</SPAN>
    <SPAN class="comment token"># sur_b1=arcpy.gp.ExtractByMask_sa(sur_b11, aod_1)</SPAN>
    <SPAN class="comment token"># sur_b2=arcpy.gp.ExtractByMask_sa(sur_b22, aod_1)</SPAN>
    <SPAN class="comment token"># sur_b3=arcpy.gp.ExtractByMask_sa(sur_b33, aod_1)</SPAN>
    <SPAN class="comment token"># arcpy.env.cellSize = tempEnvironment0</SPAN>



    <SPAN class="comment token">#############################################################################</SPAN>
    <SPAN class="comment token">#Calculate the scattering angle from Resample MOD03 products</SPAN>
      <SPAN class="comment token">#Css = Cos(S2 * !pi / 180.D) IDL</SPAN>
      <SPAN class="comment token">#Css = np.cos(S2 * np.pi / 180) Python</SPAN>
    <SPAN class="comment token"># CosSenZ = np.Cos(senzen * np.pi / 180)</SPAN>
    <SPAN class="comment token"># CosSolZ = np.Cos(solzen * np.pi / 180)</SPAN>
    <SPAN class="comment token"># SinSenZ = np.Sin(senzen * np.pi / 180)</SPAN>
    <SPAN class="comment token"># SinSolZ = np.Sin(solzen * np.pi / 180)</SPAN>

    <SPAN class="comment token"># RelAzm = np.abs((senazm ) - (solazm))</SPAN>
    <SPAN class="comment token"># index = np.where(RelAzm > 180.0)</SPAN>
    <SPAN class="comment token"># RelAzm[index] = 360.0- RelAzm[index]</SPAN>
    <SPAN class="comment token"># index = np.where(RelAzm <= 180.0)</SPAN>
    <SPAN class="comment token"># RelAzm[index] = 180.0- RelAzm[index]</SPAN>
    <SPAN class="comment token"># CosRe = np.Cos(np.radians(RelAzm))</SPAN>

    <SPAN class="comment token"># SctAgl =  np.acos((-CosSolZ * CosSenZ) + ((SinSolZ * SinSenZ) * CosRelA))</SPAN>

    <SPAN class="comment token"># #Ray_correction for TOA_BAND_1 and TOA_BAND_2</SPAN>
    <SPAN class="comment token"># #TOA_BAND_1</SPAN>
    <SPAN class="comment token"># 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</SPAN>
    <SPAN class="comment token"># Pr1 =  (3./16.*np.pi) * (1 + (np.cos(SctAgl) * np.cos(SctAgl)))</SPAN>
    <SPAN class="comment token"># wr1 = 1.0</SPAN>
    <SPAN class="comment token"># RayRef_B1 = (wr1 * B1_ROD * Pr1 ) / (4.0 * (CosSolZ * CosSenZ))</SPAN>
    <SPAN class="comment token"># RCR_B1 = toa_b1 - RayRef_B1</SPAN>
    <SPAN class="comment token"># #TOA_BAND_2</SPAN>
    <SPAN class="comment token"># 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</SPAN>
    <SPAN class="comment token"># Pr2 =  (3./16.*np.pi) * (1 + (np.cos(SctAgl) * np.cos(SctAgl)))</SPAN>
    <SPAN class="comment token"># wr2 = 1.0</SPAN>
    <SPAN class="comment token"># RayRef_B2 = (wr2 * B2_ROD * Pr2 ) / (4.0 * (CosSolZ * CosSenZ))</SPAN>
    <SPAN class="comment token"># RCR_B2 = toa_b2 - RayRef_B2</SPAN>

    <SPAN class="comment token">#Calculating NDVI float raster</SPAN>
    <SPAN class="comment token"># ndvi_raster = Divide(Float(Raster(toa_b2) - Raster(toa_b1)), Float(Raster(toa_b2) + Raster(toa_b1)))</SPAN>

    <SPAN class="comment token"># ndvi = (Float(toa_b2) - toa_b1) / (Float(toa_b2) + toa_b1)</SPAN>


<SPAN class="keyword token">print</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetMessages<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">print</SPAN> <SPAN class="string token">'finished run: %s\n\n'</SPAN> <SPAN class="operator token">%</SPAN> <SPAN class="punctuation token">(</SPAN>datetime<SPAN class="punctuation token">.</SPAN>datetime<SPAN class="punctuation token">.</SPAN>now<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">-</SPAN> start<SPAN class="punctuation token">)</SPAN>
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

I can't understand why it shows this error.

please help me, anyone, to resolve this matter.

And have any way to use a "foor loop" for set null function of different images?

Luke Webb‌,

Find more posts tagged with

Accepted answers

All comments

JoeBorgione

My example uses the in_memory workspace only because I'm writing summary tables there and they are relatively small. My initial approach was to write everything to disc in a file geodatabse; for me that was more overhead than I needed. Processing rasters data is a different animal altogether. Using in_memory is great for speed, but if it gets consumed, things have a tendency to go south....

XanderBakker

Hi BIJOY GAYEN ,

I get the impression that you did not write this script yourself, right? You are trying to adjust it to work with your data.

First of all, checking the result of a glob to check is it is None will not work since it returns a list which could be empty. You could check using for instance if len(CLDMASK) == 0: to see if the list is empty.

Once I fixed that I noticed that I am missing a TOA_B1 and a TOA_B2 raster. So I can't test the code. Also if the idea is to validate if it works in a loop (multiple rasters in a folder), you will have to provide more rasters as test data.

BIJOYGAYEN

I have already uploaded the Test_data kindly check it. Also I am uploading my full code.

<SPAN class="comment token">#-------------------------------------------------------------------------------</SPAN>
<SPAN class="comment token"># Name:        module1</SPAN>
<SPAN class="comment token"># Purpose:</SPAN>
<SPAN class="comment token">#</SPAN>
<SPAN class="comment token"># Author:      HONEY</SPAN>
<SPAN class="comment token">#</SPAN>
<SPAN class="comment token"># Created:     18-06-2019</SPAN>
<SPAN class="comment token"># Copyright:   (c) HONEY 2019</SPAN>
<SPAN class="comment token"># Licence:     <your licence></SPAN>
<SPAN class="comment token">#-------------------------------------------------------------------------------</SPAN>

<SPAN class="keyword token">import</SPAN> datetime
start <SPAN class="operator token">=</SPAN> datetime<SPAN class="punctuation token">.</SPAN>datetime<SPAN class="punctuation token">.</SPAN>now<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">print</SPAN> <SPAN class="string token">'start run: %s\n'</SPAN> <SPAN class="operator token">%</SPAN> <SPAN class="punctuation token">(</SPAN>start<SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">import</SPAN> arcpy <SPAN class="punctuation token">,</SPAN>os <SPAN class="punctuation token">,</SPAN>sys<SPAN class="punctuation token">,</SPAN>csv<SPAN class="punctuation token">,</SPAN>errno
<SPAN class="keyword token">from</SPAN> arcpy <SPAN class="keyword token">import</SPAN> env
<SPAN class="keyword token">from</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa <SPAN class="keyword token">import</SPAN> <SPAN class="operator token">*</SPAN>
<SPAN class="keyword token">import</SPAN> datetime
<SPAN class="keyword token">import</SPAN> re
<SPAN class="keyword token">import</SPAN> numpy <SPAN class="keyword token">as</SPAN> np
<SPAN class="keyword token">import</SPAN> glob
<SPAN class="keyword token">import</SPAN> itertools
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>

cellsize <SPAN class="operator token">=</SPAN> <SPAN class="string token">"F:\\DB_test_data\\TEST_RAY\\TEST1\\MOD02HKM_A2017001_0530_NDVI_AA.img"</SPAN>
<SPAN class="comment token">#Call cloud mask images from directory</SPAN>
d1<SPAN class="operator token">=</SPAN><SPAN class="string token">"F:\\DB_test_data\\VAR2\\cldmask\\"</SPAN>
CLDMASK <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d1 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.Aerosol_Cldmask_Land_Ocean-Aerosol_Cldmask_Land_Ocean.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
CLDMASK<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> CLDMASK <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the CLDMASK raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The CLDMASK raster files was opened successfully'</SPAN>
<SPAN class="comment token">#Call AOD images from directory</SPAN>
d2 <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"F:\\DB_test_data\\VAR2\\MOD04_l2"</SPAN>
AOD <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d2 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*Corrected_Optical_Depth_Land_2-Corrected_Optical_Depth_Land.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
AOD<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> AOD <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the AOD raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The AOD raster files was opened successfully'</SPAN>

<SPAN class="comment token">#Call MOD02HKM_TOA-B1_TOA-B2 images from directory</SPAN>
d3<SPAN class="operator token">=</SPAN><SPAN class="string token">"F:\\DB_test_data\\VAR2\\IDL\\MOD02HKM\\"</SPAN>
TOA_B1 <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d3 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*RefSB_1-EV_250_Aggr500_RefSB.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
TOA_B1<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> TOA_B1 <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the TOA_B1 raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The TOA_B1 raster files was opened successfully'</SPAN>

TOA_B2 <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d3 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*RefSB_2-EV_250_Aggr500_RefSB.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
TOA_B2<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> TOA_B2 <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the TOA_B2 raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The TOA_B2 raster files was opened successfully'</SPAN>

<SPAN class="comment token">#Call MOD03  (Height,SenZen,SenAzm,SolZen,SolAzm) images from directory</SPAN>
d4<SPAN class="operator token">=</SPAN><SPAN class="string token">"F:\\DB_test_data\\VAR2\\IDL\\MOD03_R\\"</SPAN>
Height <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d4 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.Height-Height.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
Height<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> Height <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the  Height raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The Height raster files was opened successfully'</SPAN>

SenZen <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d4 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.SensorZenith-SensorZenith.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
SenZen<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> SenZen <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the SenZen raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The SenZen raster files was opened successfully'</SPAN>

SenAzm <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d4 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.SensorAzimuth-SensorAzimuth.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
SenAzm<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> SenAzm <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the  SenAzm raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The SenAzm raster files was opened successfully'</SPAN>
SolZen <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d4 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.SolarZenith-SolarZenith.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
SolZen<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> SolZen <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the SolZen raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The SolZen raster files was opened successfully'</SPAN>
SolAzm <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d4 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.SolarAzimuth-SolarAzimuth.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
SolAzm<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> SolAzm <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the SolAzm raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The SolAzm raster files was opened successfully'</SPAN>
<SPAN class="comment token">#Call MOD09 images from directory</SPAN>
d5<SPAN class="operator token">=</SPAN><SPAN class="string token">"F:\\DB_test_data\\VAR2\\IDL\\MOD09\\"</SPAN>
Sur_B1 <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d5 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.500m_Surface_Reflectance_Band_1-500m_Surface_Reflectance_Band_1.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
Sur_B1 <SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> Sur_B1  <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the Sur_B1  raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The Sur_B1 raster files was opened successfully'</SPAN>
Sur_B2 <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d5 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.500m_Surface_Reflectance_Band_2-500m_Surface_Reflectance_Band_2.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
Sur_B2 <SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> Sur_B2  <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the Sur_B2  raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The Sur_B2 raster files was opened successfully'</SPAN>
Sur_B3 <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d5 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.500m_Surface_Reflectance_Band_3-500m_Surface_Reflectance_Band_3.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
Sur_B3 <SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> Sur_B3  <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the Sur_B3  raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The Sur_B3 raster files was opened successfully'</SPAN>
Land_cover<SPAN class="operator token">=</SPAN><SPAN class="string token">"F:\\DB_test_data\\MCD12C1\\LAND_COVER_TYPE_1_Grid_2D.img"</SPAN>
<SPAN class="keyword token">if</SPAN> Land_cover  <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the Land_cover raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The Land_cover raster files was opened successfully'</SPAN>

outdir<SPAN class="operator token">=</SPAN><SPAN class="string token">"F:\\DB_test_data\\TEST_RAY\\TEST1\\c\\d\\"</SPAN>
cellsize <SPAN class="operator token">=</SPAN> <SPAN class="string token">"F:\\DB_test_data\\TEST_RAY\\TEST1\\MOD02HKM_A2017001_0530_NDVI_AA.img"</SPAN>  

<SPAN class="keyword token">for</SPAN> a<SPAN class="punctuation token">,</SPAN>b<SPAN class="punctuation token">,</SPAN>c<SPAN class="punctuation token">,</SPAN>d<SPAN class="punctuation token">,</SPAN>e<SPAN class="punctuation token">,</SPAN>f<SPAN class="punctuation token">,</SPAN>g<SPAN class="punctuation token">,</SPAN>h<SPAN class="punctuation token">,</SPAN>i<SPAN class="punctuation token">,</SPAN>j<SPAN class="punctuation token">,</SPAN>k<SPAN class="punctuation token">,</SPAN>l <SPAN class="keyword token">in</SPAN> zip <SPAN class="punctuation token">(</SPAN>CLDMASK<SPAN class="punctuation token">,</SPAN>AOD<SPAN class="punctuation token">,</SPAN>TOA_B1<SPAN class="punctuation token">,</SPAN>TOA_B2<SPAN class="punctuation token">,</SPAN>Height<SPAN class="punctuation token">,</SPAN>SenZen<SPAN class="punctuation token">,</SPAN>SenAzm<SPAN class="punctuation token">,</SPAN>SolZen<SPAN class="punctuation token">,</SPAN>SolAzm<SPAN class="punctuation token">,</SPAN>Sur_B1<SPAN class="punctuation token">,</SPAN>Sur_B2<SPAN class="punctuation token">,</SPAN>Sur_B3<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_a:"</SPAN><SPAN class="operator token">+</SPAN> a<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_b:"</SPAN><SPAN class="operator token">+</SPAN> b<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_c:"</SPAN><SPAN class="operator token">+</SPAN> c<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_d:"</SPAN><SPAN class="operator token">+</SPAN> d<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_e:"</SPAN><SPAN class="operator token">+</SPAN> e<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_f:"</SPAN><SPAN class="operator token">+</SPAN> f<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_g:"</SPAN><SPAN class="operator token">+</SPAN> g<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_h:"</SPAN><SPAN class="operator token">+</SPAN> h<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_i:"</SPAN><SPAN class="operator token">+</SPAN> i<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_j:"</SPAN><SPAN class="operator token">+</SPAN> j<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_k:"</SPAN><SPAN class="operator token">+</SPAN> k<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_l:"</SPAN><SPAN class="operator token">+</SPAN> l<SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#######################################################################</SPAN>
    <SPAN class="comment token">#Cloud Mask scalling</SPAN>
    Scale_factor1 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1.0</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset1 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull1 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>a<SPAN class="punctuation token">,</SPAN>a<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat1"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" <= 0"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras1<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull1<SPAN class="punctuation token">)</SPAN>
    cldmask<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras1<SPAN class="operator token">-</SPAN>add_offset1<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor1
    OutRaster1 <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"cldmask"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>cldmask<SPAN class="punctuation token">,</SPAN>OutRaster1<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat1"</SPAN><SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN>ras1<SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#AOD scalling and set lessthan 0.1 AOD </SPAN>
    Scale_factor2 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0010000000474974513</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset2 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull2 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>b<SPAN class="punctuation token">,</SPAN>b<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat2"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -9999"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras2<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull2<SPAN class="punctuation token">)</SPAN>
    aod<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras2<SPAN class="operator token">-</SPAN>add_offset2<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor2
    aod_1 <SPAN class="operator token">=</SPAN> Con<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">(</SPAN>aod <SPAN class="operator token">>=</SPAN> <SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">&</SPAN>  <SPAN class="punctuation token">(</SPAN>aod <SPAN class="operator token"><=</SPAN> <SPAN class="number token">0.1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
    OutRaster2 <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"aod_1"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>aod_1<SPAN class="punctuation token">,</SPAN>OutRaster2<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat2"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#TOA_B1 scalling</SPAN>
    Scale_factor3 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.000053811826</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset3 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="operator token">-</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull3 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>c<SPAN class="punctuation token">,</SPAN>c<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat3"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = 65535"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras3<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull3<SPAN class="punctuation token">)</SPAN>
    toa_b1 <SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras3<SPAN class="operator token">-</SPAN><SPAN class="punctuation token">(</SPAN>add_offset3<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor3
    OutRaster3 <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"toa_b1"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>OutRaster3<SPAN class="punctuation token">,</SPAN>toa_b1<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat3"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#TOA_B2 Scalling</SPAN>
    Scale_factor4 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.00003255546</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset4 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="operator token">-</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull4 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>d<SPAN class="punctuation token">,</SPAN>d<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat4"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = 65535"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras4<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull4<SPAN class="punctuation token">)</SPAN>
    toa_b2<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras4<SPAN class="operator token">-</SPAN><SPAN class="punctuation token">(</SPAN>add_offset4<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor4
    OutRaster4 <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"toa_b2"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>toa_b2<SPAN class="punctuation token">,</SPAN>OutRaster4<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat4"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#Height</SPAN>
    setnull <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>e<SPAN class="punctuation token">,</SPAN>e<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat5"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -32767"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull<SPAN class="punctuation token">)</SPAN>
    height<SPAN class="operator token">=</SPAN>ras
    OutRaster5 <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"height"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>height<SPAN class="punctuation token">,</SPAN>OutRaster5<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat5"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#Sensor Zenith angle</SPAN>
    Scale_factor6 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.01</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset6 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull6 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>f<SPAN class="punctuation token">,</SPAN>f<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat6"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -32767"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras6<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull6<SPAN class="punctuation token">)</SPAN>
    senzen<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras6<SPAN class="operator token">-</SPAN>add_offset6<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor6
    OutRaster6 <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"senzen"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>senzen<SPAN class="punctuation token">,</SPAN>OutRaster6<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat6"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#Sensor azimuth angle</SPAN>
    setnull7 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>g<SPAN class="punctuation token">,</SPAN>g<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat7"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -32767"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras7<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull7<SPAN class="punctuation token">)</SPAN>
    senazm<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras7<SPAN class="operator token">-</SPAN>add_offset6<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor6
    OutRaster7 <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"senazm"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>senazm<SPAN class="punctuation token">,</SPAN>OutRaster7<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat7"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#solar zenithal angle</SPAN>
    setnull8 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>h<SPAN class="punctuation token">,</SPAN>h<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat8"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -32767"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras8<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull8<SPAN class="punctuation token">)</SPAN>
    solzen<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras8<SPAN class="operator token">-</SPAN>add_offset6<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor6
    OutRaster8 <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"solzen"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>solzen<SPAN class="punctuation token">,</SPAN>OutRaster8<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat8"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#solar azimuth angle</SPAN>
    setnull9 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>i<SPAN class="punctuation token">,</SPAN>i<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat9"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -32767"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras9<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull9<SPAN class="punctuation token">)</SPAN>
    solazm<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras9<SPAN class="operator token">-</SPAN>add_offset6<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor6
    OutRaster9 <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"solazm"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>solazm<SPAN class="punctuation token">,</SPAN>OutRaster9<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat9"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#(MODO9)surface reflectance band 1</SPAN>
    Scale_factor_10 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0001</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset_10 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull_10 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>j<SPAN class="punctuation token">,</SPAN>j<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat10"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -28672"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras_10<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull_10<SPAN class="punctuation token">)</SPAN>
    sur_b1<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras_10<SPAN class="operator token">-</SPAN>add_offset_10<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor_10
    OutRaster10 <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"sur_b1"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>sur_b1<SPAN class="punctuation token">,</SPAN>OutRaster10<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat10"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#(MODO9)surface reflectance band 2</SPAN>
    setnull_11 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>k<SPAN class="punctuation token">,</SPAN>k<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat11"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -28672"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras_11<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull_11<SPAN class="punctuation token">)</SPAN>
    sur_b2<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras_11<SPAN class="operator token">-</SPAN>add_offset_10<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor_10
    OutRaster11 <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"sur_b2"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>sur_b2<SPAN class="punctuation token">,</SPAN>OutRaster11<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat11"</SPAN><SPAN class="punctuation token">)</SPAN>
   

    <SPAN class="comment token">#(MODO9)surface reflectance band 3</SPAN>
    setnull_12 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>l<SPAN class="punctuation token">,</SPAN>l<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat12"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -28672"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras_12<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull_12<SPAN class="punctuation token">)</SPAN>
    sur_b3<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras_12<SPAN class="operator token">-</SPAN> add_offset_10<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor_10
    OutRaster12 <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"sur_b3"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>sur_b3<SPAN class="punctuation token">,</SPAN>OutRaster12<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat12"</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#############################################################################</SPAN>
    <SPAN class="comment token"># # Process: Extract by Cloud_Mask</SPAN>
    <SPAN class="comment token"># tempEnvironment0 = arcpy.env.cellSize</SPAN>
    <SPAN class="comment token"># arcpy.env.cellSize = "MAXOF"</SPAN>
    <SPAN class="comment token"># toa_bb1=arcpy.gp.ExtractByMask_sa(toa_b1, cldmask)</SPAN>
    <SPAN class="comment token"># toa_b22=arcpy.gp.ExtractByMask_sa(toa_b2, cldmask)</SPAN>
    <SPAN class="comment token"># heightt=arcpy.gp.ExtractByMask_sa(height, cldmask)</SPAN>
    <SPAN class="comment token"># senzenn=arcpy.gp.ExtractByMask_sa(senzen, cldmask)</SPAN>
    <SPAN class="comment token"># senazmm=arcpy.gp.ExtractByMask_sa(senazm, cldmask)</SPAN>
    <SPAN class="comment token"># solzenn=arcpy.gp.ExtractByMask_sa(solzen, cldmask)</SPAN>
    <SPAN class="comment token"># solazmm=arcpy.gp.ExtractByMask_sa(solazm, cldmask)</SPAN>
    <SPAN class="comment token"># sur_b11=arcpy.gp.ExtractByMask_sa(sur_b1, cldmask)</SPAN>
    <SPAN class="comment token"># sur_b22=arcpy.gp.ExtractByMask_sa(sur_b2, cldmask)</SPAN>
    <SPAN class="comment token"># sur_b33=arcpy.gp.ExtractByMask_sa(sur_b3, cldmask)</SPAN>
    <SPAN class="comment token"># arcpy.env.cellSize = tempEnvironment0</SPAN>
    <SPAN class="comment token"># # Process: Extract by Mask using AOD less than 0.1 value</SPAN>
    <SPAN class="comment token"># tempEnvironment0 = arcpy.env.cellSize</SPAN>
    <SPAN class="comment token"># arcpy.env.cellSize = cellsize</SPAN>
    <SPAN class="comment token"># toa_b1=arcpy.gp.ExtractByMask_sa(toa_bb1, aod_1)</SPAN>
    <SPAN class="comment token"># toa_b2=arcpy.gp.ExtractByMask_sa(toa_bb2, aod_1)</SPAN>
    <SPAN class="comment token"># height=arcpy.gp.ExtractByMask_sa(heightt, aod_1)</SPAN>
    <SPAN class="comment token"># senzen=arcpy.gp.ExtractByMask_sa(senzenn, aod_1)</SPAN>
    <SPAN class="comment token"># senazm=arcpy.gp.ExtractByMask_sa(senazmm, aod_1)</SPAN>
    <SPAN class="comment token"># solzen=arcpy.gp.ExtractByMask_sa(solzenn, aod_1)</SPAN>
    <SPAN class="comment token"># solazm=arcpy.gp.ExtractByMask_sa(solazmm, aod_1)</SPAN>
    <SPAN class="comment token"># sur_b1=arcpy.gp.ExtractByMask_sa(sur_b11, aod_1)</SPAN>
    <SPAN class="comment token"># sur_b2=arcpy.gp.ExtractByMask_sa(sur_b22, aod_1)</SPAN>
    <SPAN class="comment token"># sur_b3=arcpy.gp.ExtractByMask_sa(sur_b33, aod_1)</SPAN>
    <SPAN class="comment token"># arcpy.env.cellSize = tempEnvironment0</SPAN>



    <SPAN class="comment token">#############################################################################</SPAN>
    <SPAN class="comment token">#Calculate the scattering angle from Resample MOD03 products</SPAN>
      <SPAN class="comment token">#Css = Cos(S2 * !pi / 180.D) IDL</SPAN>
      <SPAN class="comment token">#Css = np.cos(S2 * np.pi / 180) Python</SPAN>
    <SPAN class="comment token"># CosSenZ = np.Cos(senzen * np.pi / 180)</SPAN>
    <SPAN class="comment token"># CosSolZ = np.Cos(solzen * np.pi / 180)</SPAN>
    <SPAN class="comment token"># SinSenZ = np.Sin(senzen * np.pi / 180)</SPAN>
    <SPAN class="comment token"># SinSolZ = np.Sin(solzen * np.pi / 180)</SPAN>

    <SPAN class="comment token"># RelAzm = np.abs((senazm ) - (solazm))</SPAN>
    <SPAN class="comment token"># index = np.where(RelAzm > 180.0)</SPAN>
    <SPAN class="comment token"># RelAzm[index] = 360.0- RelAzm[index]</SPAN>
    <SPAN class="comment token"># index = np.where(RelAzm <= 180.0)</SPAN>
    <SPAN class="comment token"># RelAzm[index] = 180.0- RelAzm[index]</SPAN>
    <SPAN class="comment token"># CosRe = np.Cos(np.radians(RelAzm))</SPAN>

    <SPAN class="comment token"># SctAgl =  np.acos((-CosSolZ * CosSenZ) + ((SinSolZ * SinSenZ) * CosRelA))</SPAN>

    <SPAN class="comment token"># #Ray_correction for TOA_BAND_1 and TOA_BAND_2</SPAN>
    <SPAN class="comment token"># #TOA_BAND_1</SPAN>
    <SPAN class="comment token"># 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</SPAN>
    <SPAN class="comment token"># Pr1 =  (3./16.*np.pi) * (1 + (np.cos(SctAgl) * np.cos(SctAgl)))</SPAN>
    <SPAN class="comment token"># wr1 = 1.0</SPAN>
    <SPAN class="comment token"># RayRef_B1 = (wr1 * B1_ROD * Pr1 ) / (4.0 * (CosSolZ * CosSenZ))</SPAN>
    <SPAN class="comment token"># RCR_B1 = toa_b1 - RayRef_B1</SPAN>
    <SPAN class="comment token"># #TOA_BAND_2</SPAN>
    <SPAN class="comment token"># 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</SPAN>
    <SPAN class="comment token"># Pr2 =  (3./16.*np.pi) * (1 + (np.cos(SctAgl) * np.cos(SctAgl)))</SPAN>
    <SPAN class="comment token"># wr2 = 1.0</SPAN>
    <SPAN class="comment token"># RayRef_B2 = (wr2 * B2_ROD * Pr2 ) / (4.0 * (CosSolZ * CosSenZ))</SPAN>
    <SPAN class="comment token"># RCR_B2 = toa_b2 - RayRef_B2</SPAN>

    <SPAN class="comment token">#Calculating NDVI float raster</SPAN>
    <SPAN class="comment token"># ndvi_raster = Divide(Float(Raster(toa_b2) - Raster(toa_b1)), Float(Raster(toa_b2) + Raster(toa_b1)))</SPAN>

    <SPAN class="comment token"># ndvi = (Float(toa_b2) - toa_b1) / (Float(toa_b2) + toa_b1)</SPAN>


<SPAN class="keyword token">print</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetMessages<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">print</SPAN> <SPAN class="string token">'finished run: %s\n\n'</SPAN> <SPAN class="operator token">%</SPAN> <SPAN class="punctuation token">(</SPAN>datetime<SPAN class="punctuation token">.</SPAN>datetime<SPAN class="punctuation token">.</SPAN>now<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">-</SPAN> start<SPAN class="punctuation token">)</SPAN>
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Test_data.rar

XanderBakker

Hi BIJOY GAYEN ,

It would help to have access to some test data to see if the problem can be resolved. Without data it will be difficult.

BIJOYGAYEN

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

BIJOYGAYEN

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.

BIJOYGAYEN

Thanks for your suggestion.

JoeBorgione

Sorry, I can't without doing what I suggested. You've got your code in a modular fashion; my personal preference is to take what you've got so far and make them into their own functions or defs(). That way you can call each one in order from the main() function. Makes trouble shooting a bit easier as you can turn them on and off in main() and see where the error occurs. I also like to bracket any geoprocessing functions within Try/Except blocks with print statements that tell you where things went wrong, if they do. Here's an example of one I working on today:

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">mergeClientToDelivery</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
    outWS <SPAN class="operator token">=</SPAN> <SPAN class="string token">'in_memory'</SPAN>
    sourceTable <SPAN class="operator token">=</SPAN> <SPAN class="string token">'{}\\ClientSum2'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>outWS<SPAN class="punctuation token">)</SPAN>
    targetTable <SPAN class="operator token">=</SPAN> <SPAN class="string token">'{}\\Deliveries'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>outWS<SPAN class="punctuation token">)</SPAN>
    
    sourceFields <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="string token">'Zipcode'</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="string token">'ClientCount'</SPAN><SPAN class="punctuation token">]</SPAN>
    targetFields <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="string token">'Zipcode'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'ClientCount'</SPAN><SPAN class="punctuation token">]</SPAN>
    
    <SPAN class="comment token">#create dictionary....</SPAN>
    clientDict <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">{</SPAN>r<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">:</SPAN><SPAN class="punctuation token">(</SPAN>r<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">for</SPAN> r <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>sourceTable<SPAN class="punctuation token">,</SPAN>sourceFields<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">}</SPAN>
    
    
    <SPAN class="comment token">#use the dictionary to update the deliveries table...</SPAN>
    <SPAN class="keyword token">try</SPAN><SPAN class="punctuation token">:</SPAN>
        <SPAN class="keyword token">with</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>UpdateCursor<SPAN class="punctuation token">(</SPAN>targetTable<SPAN class="punctuation token">,</SPAN>targetFields<SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">as</SPAN> updateCursor<SPAN class="punctuation token">:</SPAN>
            <SPAN class="keyword token">for</SPAN> updateRow <SPAN class="keyword token">in</SPAN> updateCursor<SPAN class="punctuation token">:</SPAN>
                keyValue <SPAN class="operator token">=</SPAN> updateRow<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN>
                <SPAN class="keyword token">if</SPAN> keyValue <SPAN class="keyword token">in</SPAN> clientDict<SPAN class="punctuation token">:</SPAN>
                    updateRow<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">=</SPAN> clientDict<SPAN class="punctuation token">.</SPAN>get<SPAN class="punctuation token">(</SPAN>keyValue<SPAN class="punctuation token">,</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
                    updateCursor<SPAN class="punctuation token">.</SPAN>updateRow<SPAN class="punctuation token">(</SPAN>updateRow<SPAN class="punctuation token">)</SPAN>
        <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="string token">'Success: mergeClientToDelivery updated fields in deliveries table in memory'</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">except</SPAN> Exception <SPAN class="keyword token">as</SPAN> err<SPAN class="punctuation token">:</SPAN>
        <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="string token">'Error: mergeClientToDelivery could not update deliveries table in memory'</SPAN><SPAN class="punctuation token">)</SPAN>
        <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>err<SPAN class="punctuation token">)</SPAN>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

BIJOYGAYEN

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

JoeBorgione

Sorry, I just don't have the time to wade through your code; I've hundreds of lines of my own code I'm wading through. My I suggest you study what you've done and execute the code block block until you find what giving you the error.

BIJOYGAYEN

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

<SPAN class="comment token">#-------------------------------------------------------------------------------</SPAN>
<SPAN class="comment token"># Name:        module1</SPAN>
<SPAN class="comment token"># Purpose:</SPAN>
<SPAN class="comment token">#</SPAN>
<SPAN class="comment token"># Author:      HONEY</SPAN>
<SPAN class="comment token">#</SPAN>
<SPAN class="comment token"># Created:     18-06-2019</SPAN>
<SPAN class="comment token"># Copyright:   (c) HONEY 2019</SPAN>
<SPAN class="comment token"># Licence:     <your licence></SPAN>
<SPAN class="comment token">#-------------------------------------------------------------------------------</SPAN>

<SPAN class="keyword token">import</SPAN> datetime
start <SPAN class="operator token">=</SPAN> datetime<SPAN class="punctuation token">.</SPAN>datetime<SPAN class="punctuation token">.</SPAN>now<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">print</SPAN> <SPAN class="string token">'start run: %s\n'</SPAN> <SPAN class="operator token">%</SPAN> <SPAN class="punctuation token">(</SPAN>start<SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">import</SPAN> arcpy <SPAN class="punctuation token">,</SPAN>os <SPAN class="punctuation token">,</SPAN>sys<SPAN class="punctuation token">,</SPAN>csv<SPAN class="punctuation token">,</SPAN>errno
<SPAN class="keyword token">from</SPAN> arcpy <SPAN class="keyword token">import</SPAN> env
<SPAN class="keyword token">from</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa <SPAN class="keyword token">import</SPAN> <SPAN class="operator token">*</SPAN>
<SPAN class="keyword token">import</SPAN> datetime
<SPAN class="keyword token">import</SPAN> re
<SPAN class="keyword token">import</SPAN> numpy <SPAN class="keyword token">as</SPAN> np
<SPAN class="keyword token">import</SPAN> glob
<SPAN class="keyword token">import</SPAN> itertools
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>
outdir<SPAN class="operator token">=</SPAN><SPAN class="string token">"F:\\DB_test_data\\TEST_RAY\\TEST1\\c\\"</SPAN>

cellsize <SPAN class="operator token">=</SPAN> <SPAN class="string token">"F:\\DB_test_data\\TEST_RAY\\TEST1\\MOD02HKM_A2017001_0530_NDVI_AA.img"</SPAN>
<SPAN class="comment token">#Call cloud mask images from directory</SPAN>
d1<SPAN class="operator token">=</SPAN><SPAN class="string token">"F:\\DB_test_data\\VAR2\\cldmask\\"</SPAN>
CLDMASK <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d1 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.Aerosol_Cldmask_Land_Ocean-Aerosol_Cldmask_Land_Ocean.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
CLDMASK<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> CLDMASK <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the CLDMASK raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The CLDMASK raster files was opened successfully'</SPAN>
<SPAN class="comment token">#Call AOD images from directory</SPAN>
d2 <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"F:\\DB_test_data\\VAR2\\MOD04_l2"</SPAN>
AOD <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d2 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*Corrected_Optical_Depth_Land_2-Corrected_Optical_Depth_Land.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
AOD<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> AOD <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the AOD raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The AOD raster files was opened successfully'</SPAN>

<SPAN class="comment token">#Call MOD02HKM_TOA-B1_TOA-B2 images from directory</SPAN>
d3<SPAN class="operator token">=</SPAN><SPAN class="string token">"F:\\DB_test_data\\VAR2\\IDL\\MOD02HKM\\"</SPAN>
TOA_B1 <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d3 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*RefSB_1-EV_250_Aggr500_RefSB.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
TOA_B1<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> TOA_B1 <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the TOA_B1 raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The TOA_B1 raster files was opened successfully'</SPAN>

TOA_B2 <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d3 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*RefSB_2-EV_250_Aggr500_RefSB.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
TOA_B2<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> TOA_B2 <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the TOA_B2 raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The TOA_B2 raster files was opened successfully'</SPAN>

<SPAN class="comment token">#Call MOD03  (Height,SenZen,SenAzm,SolZen,SolAzm) images from directory</SPAN>
d4<SPAN class="operator token">=</SPAN><SPAN class="string token">"F:\\DB_test_data\\VAR2\\IDL\\MOD03_R\\"</SPAN>
Height <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d4 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.Height-Height.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
Height<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> Height <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the  Height raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The Height raster files was opened successfully'</SPAN>

SenZen <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d4 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.SensorZenith-SensorZenith.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
SenZen<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> SenZen <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the SenZen raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The SenZen raster files was opened successfully'</SPAN>

SenAzm <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d4 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.SensorAzimuth-SensorAzimuth.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
SenAzm<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> SenAzm <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the  SenAzm raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The SenAzm raster files was opened successfully'</SPAN>
SolZen <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d4 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.SolarZenith-SolarZenith.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
SolZen<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> SolZen <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the SolZen raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The SolZen raster files was opened successfully'</SPAN>
SolAzm <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d4 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.SolarAzimuth-SolarAzimuth.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
SolAzm<SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> SolAzm <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the SolAzm raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The SolAzm raster files was opened successfully'</SPAN>
<SPAN class="comment token">#Call MOD09 images from directory</SPAN>
d5<SPAN class="operator token">=</SPAN><SPAN class="string token">"F:\\DB_test_data\\VAR2\\IDL\\MOD09\\"</SPAN>
Sur_B1 <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d5 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.500m_Surface_Reflectance_Band_1-500m_Surface_Reflectance_Band_1.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
Sur_B1 <SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> Sur_B1  <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the Sur_B1  raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The Sur_B1 raster files was opened successfully'</SPAN>
Sur_B2 <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d5 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.500m_Surface_Reflectance_Band_2-500m_Surface_Reflectance_Band_2.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
Sur_B2 <SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> Sur_B2  <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the Sur_B2  raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The Sur_B2 raster files was opened successfully'</SPAN>
Sur_B3 <SPAN class="operator token">=</SPAN> glob<SPAN class="punctuation token">.</SPAN>glob<SPAN class="punctuation token">(</SPAN>d5 <SPAN class="operator token">+</SPAN> os<SPAN class="punctuation token">.</SPAN>sep <SPAN class="operator token">+</SPAN> <SPAN class="string token">"*.500m_Surface_Reflectance_Band_3-500m_Surface_Reflectance_Band_3.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
Sur_B3 <SPAN class="punctuation token">.</SPAN>sort<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">if</SPAN> Sur_B3  <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the Sur_B3  raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The Sur_B3 raster files was opened successfully'</SPAN>
Land_cover<SPAN class="operator token">=</SPAN><SPAN class="string token">"F:\\DB_test_data\\MCD12C1\\LAND_COVER_TYPE_1_Grid_2D.img"</SPAN>
<SPAN class="keyword token">if</SPAN> Land_cover  <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'Could not open the Land_cover raster files'</SPAN>
    sys<SPAN class="punctuation token">.</SPAN>exit<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="string token">'The Land_cover raster files was opened successfully'</SPAN>

outdir<SPAN class="operator token">=</SPAN><SPAN class="string token">"F:\\DB_test_data\\TEST_RAY\\TEST1\\c\\"</SPAN>
cellsize <SPAN class="operator token">=</SPAN> <SPAN class="string token">"F:\\DB_test_data\\TEST_RAY\\TEST1\\MOD02HKM_A2017001_0530_NDVI_AA.img"</SPAN>  

<SPAN class="keyword token">for</SPAN> a<SPAN class="punctuation token">,</SPAN>b<SPAN class="punctuation token">,</SPAN>c<SPAN class="punctuation token">,</SPAN>d<SPAN class="punctuation token">,</SPAN>e<SPAN class="punctuation token">,</SPAN>f<SPAN class="punctuation token">,</SPAN>g<SPAN class="punctuation token">,</SPAN>h<SPAN class="punctuation token">,</SPAN>i<SPAN class="punctuation token">,</SPAN>j<SPAN class="punctuation token">,</SPAN>k<SPAN class="punctuation token">,</SPAN>l <SPAN class="keyword token">in</SPAN> zip <SPAN class="punctuation token">(</SPAN>CLDMASK<SPAN class="punctuation token">,</SPAN>AOD<SPAN class="punctuation token">,</SPAN>TOA_B1<SPAN class="punctuation token">,</SPAN>TOA_B2<SPAN class="punctuation token">,</SPAN>Height<SPAN class="punctuation token">,</SPAN>SenZen<SPAN class="punctuation token">,</SPAN>SenAzm<SPAN class="punctuation token">,</SPAN>SolZen<SPAN class="punctuation token">,</SPAN>SolAzm<SPAN class="punctuation token">,</SPAN>Sur_B1<SPAN class="punctuation token">,</SPAN>Sur_B2<SPAN class="punctuation token">,</SPAN>Sur_B3<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_a:"</SPAN><SPAN class="operator token">+</SPAN> a<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_b:"</SPAN><SPAN class="operator token">+</SPAN> b<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_c:"</SPAN><SPAN class="operator token">+</SPAN> c<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_d:"</SPAN><SPAN class="operator token">+</SPAN> d<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_e:"</SPAN><SPAN class="operator token">+</SPAN> e<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_f:"</SPAN><SPAN class="operator token">+</SPAN> f<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_g:"</SPAN><SPAN class="operator token">+</SPAN> g<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_h:"</SPAN><SPAN class="operator token">+</SPAN> h<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_i:"</SPAN><SPAN class="operator token">+</SPAN> i<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_j:"</SPAN><SPAN class="operator token">+</SPAN> j<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_k:"</SPAN><SPAN class="operator token">+</SPAN> k<SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">print</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"processing_l:"</SPAN><SPAN class="operator token">+</SPAN> l<SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#######################################################################</SPAN>
    <SPAN class="comment token">#Cloud Mask scalling</SPAN>
    Scale_factor1 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1.0</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset1 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull1 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>a<SPAN class="punctuation token">,</SPAN>a<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat1"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" <= 0"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras1<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull1<SPAN class="punctuation token">)</SPAN>
    cldmask<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras1<SPAN class="operator token">-</SPAN>add_offset1<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor1
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"cldmask"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>cldmask<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat1"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#AOD scalling and set lessthan 0.1 AOD </SPAN>
    Scale_factor2 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0010000000474974513</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset2 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull2 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>b<SPAN class="punctuation token">,</SPAN>b<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat2"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -9999"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras2<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull2<SPAN class="punctuation token">)</SPAN>
    aod<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras2<SPAN class="operator token">-</SPAN>add_offset2<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor2
    aod_1 <SPAN class="operator token">=</SPAN> Con<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">(</SPAN>aod <SPAN class="operator token">>=</SPAN> <SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">&</SPAN> <SPAN class="punctuation token">(</SPAN>aod <SPAN class="operator token"><=</SPAN> <SPAN class="number token">0.1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"aod_1"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>aod_1<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat2"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#TOA_B1 scalling</SPAN>
    Scale_factor3 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.000053811826</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset3 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="operator token">-</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull3 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>c<SPAN class="punctuation token">,</SPAN>c<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat3"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = 65535"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras3<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull3<SPAN class="punctuation token">)</SPAN>
    toa_b1<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras3<SPAN class="operator token">-</SPAN><SPAN class="punctuation token">(</SPAN>add_offset3<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor3
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"toa_b1"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>toa_b1<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat3"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#TOA_B2 Scalling</SPAN>
    Scale_factor4 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.00003255546</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset4 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="operator token">-</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull4 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>d<SPAN class="punctuation token">,</SPAN>d<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat4"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = 65535"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras4<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull4<SPAN class="punctuation token">)</SPAN>
    toa_b2<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras4<SPAN class="operator token">-</SPAN><SPAN class="punctuation token">(</SPAN>add_offset4<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor4
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"toa_b2"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>toa_b2<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat4"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#Height</SPAN>
    setnull <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>e<SPAN class="punctuation token">,</SPAN>e<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat5"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -32767"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull<SPAN class="punctuation token">)</SPAN>
    height<SPAN class="operator token">=</SPAN>ras
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"height"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>height<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat5"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#Sensor Zenith angle</SPAN>
    Scale_factor6 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.01</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset6 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull6 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>f<SPAN class="punctuation token">,</SPAN>f<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat6"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -32767"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras6<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull6<SPAN class="punctuation token">)</SPAN>
    senzen<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras6<SPAN class="operator token">-</SPAN>add_offset6<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor6
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"senzen"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>senzen<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat6"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#Sensor azimuth angle</SPAN>
    setnull7 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>g<SPAN class="punctuation token">,</SPAN>g<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat7"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -32767"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras7<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull7<SPAN class="punctuation token">)</SPAN>
    senazm<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras7<SPAN class="operator token">-</SPAN>add_offset6<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor6
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"senazm"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>senazm<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat7"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#solar zenithal angle</SPAN>
    setnull8 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>h<SPAN class="punctuation token">,</SPAN>h<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat8"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -32767"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras8<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull8<SPAN class="punctuation token">)</SPAN>
    solzen<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras8<SPAN class="operator token">-</SPAN>add_offset6<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor6
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"solzen"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>solzen<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat8"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#solar azimuth angle</SPAN>
    setnull9 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>i<SPAN class="punctuation token">,</SPAN>i<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat9"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -32767"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras9<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull9<SPAN class="punctuation token">)</SPAN>
    solazm<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras9<SPAN class="operator token">-</SPAN>add_offset6<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor6
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"solazm"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>solazm<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat9"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#(MODO9)surface reflectance band 1</SPAN>
    Scale_factor_10 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0001</SPAN><SPAN class="punctuation token">)</SPAN>
    add_offset_10 <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0.0</SPAN><SPAN class="punctuation token">)</SPAN>
    setnull_10 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>j<SPAN class="punctuation token">,</SPAN>j<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat10"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -28672"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras_10<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull_10<SPAN class="punctuation token">)</SPAN>
    sur_b1<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras_10<SPAN class="operator token">-</SPAN>add_offset_10<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor_10
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"sur_b1"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>sur_b1<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat10"</SPAN><SPAN class="punctuation token">)</SPAN>

    <SPAN class="comment token">#(MODO9)surface reflectance band 2</SPAN>
    setnull_11 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>k<SPAN class="punctuation token">,</SPAN>k<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat11"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -28672"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras_11<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull_11<SPAN class="punctuation token">)</SPAN>
    sur_b2<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras_11<SPAN class="operator token">-</SPAN>add_offset_10<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor_10
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"sur_b2"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>sur_b2<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat11"</SPAN><SPAN class="punctuation token">)</SPAN>
   

    <SPAN class="comment token">#(MODO9)surface reflectance band 3</SPAN>
    setnull_12 <SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>gp<SPAN class="punctuation token">.</SPAN>SetNull_sa<SPAN class="punctuation token">(</SPAN>l<SPAN class="punctuation token">,</SPAN>l<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/dat12"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"\"Value\" = -28672"</SPAN><SPAN class="punctuation token">)</SPAN>
    ras_12<SPAN class="operator token">=</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Raster<SPAN class="punctuation token">(</SPAN>setnull_12<SPAN class="punctuation token">)</SPAN>
    sur_b3<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>ras_12<SPAN class="operator token">-</SPAN> add_offset_10<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">*</SPAN>Scale_factor_10
    OutRaster <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outdir<SPAN class="punctuation token">,</SPAN><SPAN class="string token">'MOd0.{0}.img'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"sur_b3"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
    <SPAN class="comment token">#print OutRaster</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>CopyRaster_management<SPAN class="punctuation token">(</SPAN>sur_b3<SPAN class="punctuation token">,</SPAN>OutRaster<SPAN class="punctuation token">)</SPAN>
    arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/dat12"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

JoeBorgione

The error is stating a problem on line 192, but the code you provide only goes to 136: can't help you there.

This stack exchange thread mentions background processing which you asked about in an earlier post: https://gis.stackexchange.com/questions/27190/error-999998-unexpected-error (Dr Google is a good friend when it comes to error codes...)

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.