def process_temporal_data(lat, lon, root):
times = [ datetime.fromtimestamp(int(t)) for t in nc.getvar(root, 'time') ]
indexes = range(len(times))
for i in indexes:
show("\rTemporal data: preprocessing image %d / %d " % (i, len(indexes)-1))
dt = times[i]
# Calculate some geometry parameters
# Parameters that need the datetime: gamma, tst_hour, slots, linketurbidity
gamma = geo.getdailyangle(geo.getjulianday(dt),geo.gettotaldays(dt))
tst_hour = geo.gettsthour(geo.getdecimalhour(dt), GREENWICH_LON, lon, geo.gettimeequation(gamma))
declination = geo.getdeclination(gamma)
slots = geo.getslots(dt,IMAGE_PER_HOUR)
omega = geo.gethourlyangle(tst_hour, lat/abs(lat))
solarangle = geo.getzenithangle(declination,lat,omega)
solarelevation = geo.getelevation(solarangle)
excentricity = geo.getexcentricity(gamma)
save_temporal_data(root,i,gamma,tst_hour,declination,solarangle,solarelevation,excentricity,slots)
say("Projecting Linke's turbidity index... ")
linke.cut_projected(root)
say("Calculating the satellital zenith angle... ")
satellitalzenithangle = geo.getsatellitalzenithangle(lat, lon, SAT_LON)
dem.cut_projected(root)
v_satellitalzenithangle = nc.getvar(root,'satellitalzenithangle', 'f4', ('northing','easting',),4)
v_satellitalzenithangle[:] = satellitalzenithangle
nc.sync(root)
def process_temporal_data(lat, lon, root):
times = [ datetime.utcfromtimestamp(int(t)) for t in nc.getvar(root, 'data_time')[:] ]
indexes = range(len(times))
gamma = nc.clonevar(root,'data_time', 'gamma')
nc.sync(root)
tst_hour = nc.clonevar(root,'data', 'tst_hour')
declination = nc.clonevar(root, 'gamma', 'declination')
solarangle = nc.clonevar(root, 'data', 'solarangle')
solarelevation = nc.clonevar(root, 'solarangle', 'solarelevation')
excentricity = nc.clonevar(root, 'gamma', 'excentricity')
slots = nc.getvar(root,'slots', 'u1', ('timing',))
nc.sync(root)
for i in indexes:
show("\rTemporal data: preprocessing image %d / %d " % (i, len(indexes)-1))
dt = times[i]
# Calculate some geometry parameters
# Parameters that need the datetime: gamma, tst_hour, slots, linketurbidity
gamma[i] = geo.getdailyangle(geo.getjulianday(dt),geo.gettotaldays(dt))
tst_hour[i,:] = geo.gettsthour(geo.getdecimalhour(dt), GREENWICH_LON, lon, geo.gettimeequation(gamma[i]))
declination[i] = geo.getdeclination(gamma[i])
slots[i] = geo.getslots(dt,IMAGE_PER_HOUR)
omega = geo.gethourlyangle(tst_hour[i], lat/abs(lat))
solarangle[i] = geo.getzenithangle(declination[i],lat,omega)
solarelevation[i] = geo.getelevation(solarangle[i])
excentricity[i] = geo.getexcentricity(gamma[i])
nc.sync(root)
say("Projecting Linke's turbidity index... ")
linke.cut_projected(root)
say("Calculating the satellital zenith angle... ")
satellitalzenithangle = geo.getsatellitalzenithangle(lat, lon, SAT_LON)
dem.cut_projected(root)
v_satellitalzenithangle = nc.clonevar(root,'lat', 'satellitalzenithangle')
v_satellitalzenithangle[:] = satellitalzenithangle
nc.sync(root)
v_satellitalzenithangle = None
def process_atmospheric_irradiance(lat, lon, data, root):
i0met = geti0met()
dc = nc.getvar(root,'dc')[:]
satellitalzenithangle = nc.getvar(root,'satellitalzenithangle')[:]
excentricity = nc.getvar(root,'excentricity')[:]
say("Calculating atmospheric irradiance... ")
atmosphericradiance = geo.getatmosphericradiance(1367.0, i0met,dc, satellitalzenithangle)
atmosphericalbedo = geo.getalbedo(atmosphericradiance, i0met, excentricity, satellitalzenithangle)
satellitalelevation = geo.getelevation(satellitalzenithangle)
v_atmosphericalbedo = nc.getvar(root, 'atmosphericalbedo', 'f4', ('timing','northing','easting',),4)
v_atmosphericalbedo[:] = atmosphericalbedo
v_satellitalelevation = nc.getvar(root, 'satellitalelevation', 'f4', ('northing','easting',),4)
v_satellitalelevation[:] = satellitalelevation
nc.sync(root)
def process_atmospheric_irradiance(root):
i0met = geti0met()
dc = nc.getvar(root,'dc')[:]
satellitalzenithangle = nc.getvar(root,'satellitalzenithangle')[:]
excentricity = nc.getvar(root,'excentricity')[:]
say("Calculating atmospheric irradiance... ")
atmosphericradiance = geo.getatmosphericradiance(1367.0, i0met,dc, satellitalzenithangle)
atmosphericalbedo = geo.getalbedo(atmosphericradiance, i0met, excentricity, satellitalzenithangle)
satellitalelevation = geo.getelevation(satellitalzenithangle)
v_atmosphericalbedo = nc.clonevar(root, 'data', 'atmosphericalbedo')
v_atmosphericalbedo[:] = atmosphericalbedo
v_satellitalelevation = nc.clonevar(root, 'lon', 'satellitalelevation')
v_satellitalelevation[:] = satellitalelevation
nc.sync(root)
v_atmosphericalbedo, v_satellitalelevation = None, None
