BITeam | Our Planet, Our Home
Albacete, Castilla-La Mancha
Team Updates
Next Steps:
We will carry on about how to collect more information by several geographical worldwide zones. Trying to do an average of pollution measurement around a specific zone until complete all the world.
P
Pablo Palazon SanchezProject Resources:
P
Pablo Palazon SanchezJSON Pollution's city format:
"manchester": { "center": { "lat": 53.467, "lng": -2.233 }, "ch4": 3.779105907159125e+19, "ch4_per": 95.51937264091019, "city": "Manchester", "co": 5.537136305022727e-10, "co2": 397.56538211116316, "co2_per": 63.18718816857975, "co_per": 3.005560852106043, "dust": 1.6923117073019966e-05, "dust_per": 2.3738678003459066, "lw_flux": -45.75823974609375, "lw_flux_per": 77.03036410685885, "no2": 4823121035001856.0, "no2_per": 24.110171253212517, "sw_flux": 86.77875518798828, "sw_flux_per": 31.76687690180772, "wikipedia": "Manchester"}P
Pablo Palazon Sanchez
Panoply CH4 Plot over WorldMap
P
Pablo Palazon Sanchez| from netCDF4 import Dataset | |
| import numpy as np | |
| import json | |
| import sys,getopt | |
| opts, args = getopt.getopt(sys.argv[1:],':ir',['nc=', 'cities=', 'variable=', 'attr=', 'output=']) | |
| for opt,arg in opts: | |
| if opt=='--nc': | |
| filename=arg | |
| if opt=='--variable': | |
| variable=arg | |
| if opt=='--cities': | |
| cities=arg | |
| if opt=='--attr': | |
| attr=arg | |
| if opt=='--output': | |
| output=arg | |
| with open(cities, 'r') as infile: | |
| cities = json.loads(infile.read()) | |
| dataset = Dataset(filename, 'r') | |
| limit = { | |
| 'min': None, | |
| 'max': None | |
| } | |
| def atributos(group): | |
| for name in group.ncattrs(): | |
| print "Global attr", name, "=", getattr(group,name) | |
| def location(city, limit): | |
| indlat=0 | |
| indlon=0 | |
| city[attr] = None | |
| lat1 = dataset.variables['lat'] | |
| lon1 = dataset.variables['lon'] | |
| center = city['center'] | |
| for longitud in lon1: | |
| if longitud < center['lng']: | |
| indlon+=1 | |
| for latitud in lat1: | |
| if latitud < center['lat']: | |
| indlat+=1 | |
| print variable + ' ' + city['city'] | |
| print limit | |
| if indlon < lon1[:].size and indlat < lat1[:].size : | |
| city[attr]=dataset.variables[variable][indlat][indlon].tolist() | |
| if(city[attr]!=None) : | |
| if(limit['min']==None or city[attr] < limit['min']) : | |
| limit['min'] = city[attr] | |
| if(limit['max']==None or city[attr] > limit['max']) : | |
| limit['max'] = city[attr] | |
| def percentlocation(city, limit): | |
| value = city[attr] | |
| if(value != None): | |
| print 'Percent values: '+variable + ' ' + city['city'] | |
| if (limit['max'] - limit['min']) > 0 : | |
| city[attr+'_per'] = abs(((city[attr] - limit['min']) * 100) / (limit['max'] - limit['min'])) | |
| else : | |
| city[attr+'_per'] = 100 | |
| for city in cities.values(): | |
| location(city, limit) | |
| for city in cities.values(): | |
| percentlocation(city, limit) | |
| with open(output, 'w') as outfile: | |
| json.dump(cities, outfile, sort_keys=True, indent=2) | |
| dataset.close() |
P
Pablo Palazon Sanchez
Giovani NASA Data Resource
P
Pablo Palazon Sanchez
CO Pollution's city over map - Relative to others cities
P
Pablo Palazon Sanchez
Panoply netCDF, HDF and GRIB Data Viewer
P
Pablo Palazon SanchezThe code is in this repository https://github.com/alejandrorueda/spaceapp
R
Rubén Donate Serrano
SpaceApps is a NASA incubator innovation program.