R / GIS / AIR POLLUTION

Titulo en inglés, texto en espanol desde un teclado en portugues. Lo que importa en el codigo. Me pase medio dia resolviendo algo que deveria ser muy simple pero al final no lo es: Convertir un SpatialPolygonsDataFrame en RasterStack . Despues de ver varias respuestas en StackOverFlow , me di cuenta que deberia crear un loop. La linea 4 rasteriza el SpatialPolygonsdataFrame: g_NOx_NMHC La Linea 5 es el raster utilizado en la linea 4 que incluye en numero de columnas y filas Las Lineas 2, 3 y 6 son lapply para que crear un raster iterando para cada columna del SpatialGridDatFrame. La ṕrimera columna es "id" por tanto se parte desde la segunda linea. Finalmente la Linea 1 es do.call ("stack" para convertir la lista de rastaer en un RasterStack de una vez. Demora aproximadamente 5 min. y aqui el codigo:

Un simple script para generar plotear las salidas WRF

Convertí las salidas .DAT de Calpuff en un vector de puntos con Arcgis. La fuente es   esta . Luego cargué el shapefile en R, luego lo convertí a raster, y rasterizé con varias configuraciones de filas y columnas, con promedios.

Buscando Modelos Digitales de Elevación ( Digital Elevation Model; DEM) encontré el  The ASTER Global Digital Elevation Mode l (ASTER GDEM), con resolución de celdas de 30 mts!!! File Name(Example of N35~36,E135~136) ASTGTM_ N35E135_dem.tif ASTGTM_ N35E135_num.tif Output format GeoTIFF, signed 16 bits Geographic coordinates Geographic latitude and longitude Tile Size 3601-by-3601 pixels (1-by-1 degree) Posting interval 1 arc-second DN values 1m/DN referenced to the WGS84/EGM96 geoid-9999 for void pixels, and 0 for sea water body Stacking number or reference DEM SRTM3:-1,-2 CDED:-6 Alaska DEM:-11 etc. Coverage North 83 degrees to south 83 degrees, 22,600 tiles Algunas características: Resolución celda 30 mts, con precisión de 7 a 14 mts,  Para bajar la información hay que clickear en "search", en http://www.gdem.aster.ersdac.or.jp/search.jsp específicamente. Lo mas fácil es contar con shapefiles en coordenadas geográficas del área a descargar. As...