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)
