GGPLOT2     HEALTH     R     EMISSIONS     STATISTICS     GIS   WRF Open Air Air Pollution  

jueves, 23 de febrero de 2017

From SpatialPolygonsDataFrame to RasterStack in R

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: