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 cut_projected_terrain(filename): from libs.dem import dem root = nc.open(filename)[0] lat = nc.getvar(root, 'lat') lon = nc.getvar(root, 'lon') data = nc.getvar(root, 'data') time = nc.getvar(root, 'data_time') root_cut = nc.open('wterrain.' + filename)[0] nc.getdim(root_cut, 'timing', data.shape[0]) nc.getdim(root_cut, 'northing', data.shape[1]) nc.getdim(root_cut, 'easting', data.shape[2]) lat_cut = nc.getvar(root_cut, 'lat', 'f4', ('northing','easting',),4) lon_cut = nc.getvar(root_cut, 'lon', 'f4', ('northing','easting',),4) data_cut = nc.getvar(root_cut, 'data', 'f4', ('timing','northing','easting',),4) time_cut = nc.getvar(root_cut, 'data_time', 'f4', ('timing',),4) lat_cut[:] = lat[:] lon_cut[:] = lon[:] data_cut[:] = data[:] time_cut[:] = time[:] dem.cut_projected(root_cut) nc.close(root) nc.close(root_cut)