Agustin Lobo

Recently Published

Problems witth lmPerm::lmp() for regression
Some problems with lmp described
Simple Markdown syntax for formatting your R scripts into WEB pages using knitr::spin()
Most Basic use of knitr::spin()
Comparisson of Copernicus Time Series of NDVI V1 and V2 in Valderejo Test Site (Val2013_V1V2comp_log.R)
Exercise R Lab 2.1: Pollutants in soils of the Kola Peninsula
Length of Day An example of computing with dates in R
R Lab1
talleR ICTJA RLab1
to be removed
Overlay your raster layer on a background GM layer in R
GM tiles can be downloaded as R objects using either ggmap:get_map() or dismo:gmap(). We prefer dismo:gmap() because, unlike ggmap:get_map(), it returns a raster layer from package raster including complete CRS information
Variables that describe local displacement, such as wind fields, can be very effectively represented with function rasterVis:vectorplot() by Oscar Perpinan. We will use a climatological wind field dowloaded from the NOAA Blended Sea Winds
Plotting vector maps overlaid on top of a raster background in R
Objects returned by dismo:gmap() are raster layers of package *raster* with the correct CRS information ("Google PseudoMercator" epsg:3857) that can be used as background for displaying other geographic layers in R provided their CRS are consistent. We show here the different alternatives.
Introduction to rasterVis
rasterVis is an R package for the visualization of raster layers authored by Oscar Perpinan ( This introductory material that does not substitute the standard documentation provided with the package and the documents provided by Oscar Perpinan:
Finding out the CRS of the object returned by get_map()
The CRS of objects returned by get_map() is not documented. Here we use an empirical approach to confirm that, as it is usual in tools interacting with GM, the map is in so called "Google PseudoMercator" while coordinates are provided in "Longitude, Latitude" (datum WGS84 in both cases) and show how to produce correct plots of overlaid maps