The goal of this package is to reduce the friction when working with maps & time-series.
It provides two main functionalities: the visualisation of gridded data & choropleths, and the aggregation from gridded data to time-series using geographical borders.
The package is not available on CRAN, you can get the development version from github using devtools:
devtools::install_github("matteodefelice/panas")
Load the panas
package:
library(panas)
The function get_region_from_coordinates
returns the NUTS region where a set of points lie.
> x = data.frame(lat = c(48.5, 49), lon = c(12, 4))
> get_region_from_coordinates(x, shapefile = 'NUTS1')
[1] DE2 FRF
121 Levels: AL0 AT1 AT2 AT3 BE1 BE2 BE3 BG3 BG4 CH0 CY0 ... UKN
This package gives you the possibility to visualise gridded data with the function plot_field_discrete
and choropleths using NUTS classification (this means only Europe).
> data(ncep)
> g = plot_field_discrete(z, lon, lat, latlim = c(-30, 70), breaks = c(5e-5, 1e-4), color_scale = 'PuBu', varname = 'prec.', grid_step = 60)
> print(g + coord_equal())
> my_data = tibble(area = c('IT', 'ITC', 'UKG', 'ITC1'), value = c('a','b','b', 'c'))
> g = get_european_choropleth(my_data)
> print(g + coord_map('lambert', 35, 58, ylim = c(35, 68), xlim = c(-15, 25)))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.