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