,
Nicolas Lambert [aut]| Title: | Create Map Insets |
|---|---|
| Description: | Map insets are small zoom-in or zoom-out maps that focus on particular territories. mapinsetr provides a set of functions that helps to create such insets. |
| Authors: | Timothée Giraud [cre, aut] ,
Nicolas Lambert [aut] |
| Maintainer: | Timothée Giraud <[email protected]> |
| License: | GPL-3 |
| Version: | 0.3.0 |
| Built: | 2025-03-20 12:28:53 UTC |
| Source: | https://github.com/riatelab/mapinsetr |
Create a mask based on a bounding box (bbox), a simple
feature collection (sf) extent, a simple feature geometry list column
(sfc) extent or an interactively defined rectancle.
create_mask(bb, prj, interactive = FALSE, add = FALSE)create_mask(bb, prj, interactive = FALSE, add = FALSE)
bb |
either a bounding box, a numeric vector (xmin, ymin, xmax, ymax), an sf or sfc object. |
prj |
a CRS string. |
interactive |
define the mask interactively. |
add |
add the mask to the current plot. |
An sf object is returned.
library(sf) nc <- st_read(system.file("shape/nc.shp", package="sf")) nc <- st_transform(nc, 32119) plot(st_geometry(nc)) bb <- st_bbox(nc[nc$CNTY_ID %in% c('2030', '1989', '1938'),]) mask <- create_mask(bb = bb, add = TRUE) plot(st_geometry(nc)) bb <- nc[nc$CNTY_ID %in% c('2030', '1989', '1938'),] mask <- create_mask(bb = bb, add = TRUE) plot(st_geometry(nc)) bb <- c(589912, 159757, 694332, 257053) mask <- create_mask(bb = bb, prj=32119, add = TRUE)library(sf) nc <- st_read(system.file("shape/nc.shp", package="sf")) nc <- st_transform(nc, 32119) plot(st_geometry(nc)) bb <- st_bbox(nc[nc$CNTY_ID %in% c('2030', '1989', '1938'),]) mask <- create_mask(bb = bb, add = TRUE) plot(st_geometry(nc)) bb <- nc[nc$CNTY_ID %in% c('2030', '1989', '1938'),] mask <- create_mask(bb = bb, add = TRUE) plot(st_geometry(nc)) bb <- c(589912, 159757, 694332, 257053) mask <- create_mask(bb = bb, prj=32119, add = TRUE)
Takes a list of sf polygons or multipolygons objects and output a single sf polygon or multipolygon object.
inset_rbinder(l = list())inset_rbinder(l = list())
l |
a list of sf POLYGON or MULTIPOLYGON objects. |
An sf polygon or multipolygon object is returned.
library(sf) nc <- st_read(system.file("shape/nc.shp", package="sf")) nc <- st_transform(nc, 32119) plot(st_geometry(nc)) mask1 <- st_buffer(st_centroid(st_geometry(nc[nc$CNTY_ID == 2026,])),dist = 30000) mask2 <- st_buffer(st_centroid(st_geometry(nc[nc$CNTY_ID == 2016,])),dist = 30000) plot(st_geometry(mask1), border = "red", lwd = 2, add = TRUE) plot(st_geometry(mask2), border = "red", lwd = 2, add = TRUE) inset1 <- move_and_resize(nc, mask1, xy = c(200000, 5000), k = 2) inset2 <- move_and_resize(nc, mask2, xy = c(200000 + 130000, 5000), k = 2) plot(st_geometry(inset1), add = TRUE) plot(st_geometry(inset2), add = TRUE) nc_insets <- inset_rbinder(l = list(nc, inset1, inset2)) plot(st_geometry(nc_insets)) plot(nc_insets[,9])library(sf) nc <- st_read(system.file("shape/nc.shp", package="sf")) nc <- st_transform(nc, 32119) plot(st_geometry(nc)) mask1 <- st_buffer(st_centroid(st_geometry(nc[nc$CNTY_ID == 2026,])),dist = 30000) mask2 <- st_buffer(st_centroid(st_geometry(nc[nc$CNTY_ID == 2016,])),dist = 30000) plot(st_geometry(mask1), border = "red", lwd = 2, add = TRUE) plot(st_geometry(mask2), border = "red", lwd = 2, add = TRUE) inset1 <- move_and_resize(nc, mask1, xy = c(200000, 5000), k = 2) inset2 <- move_and_resize(nc, mask2, xy = c(200000 + 130000, 5000), k = 2) plot(st_geometry(inset1), add = TRUE) plot(st_geometry(inset2), add = TRUE) nc_insets <- inset_rbinder(l = list(nc, inset1, inset2)) plot(st_geometry(nc_insets)) plot(nc_insets[,9])
do stuff, all layers must use a cartographic projection, no lon/lat.
m_r(x, mask, y, return_k = FALSE)m_r(x, mask, y, return_k = FALSE)
x |
the layer to cut, resize and move, sf |
mask |
the targeted area in x, sf or sfc |
y |
destination, sf or sfc |
return_k |
return the k factor |
a layer
library(sf) library(mapsf) nc <- st_read(system.file("shape/nc.shp", package="sf")) # Ne fonctionne qu'avec 2 fonds projetés nc <- st_transform(nc, 32119) # Créer des boites y <- st_as_sfc(st_bbox( c(xmin = 270000, ymin = 00000, xmax = 550000, ymax = 100000), crs = 32119)) y <- st_make_grid(y, n = c(3,2)) mf_map(nc) mf_map(y, add = TRUE) # plusieurs objets agrandis inset1 <- m_r(x = nc, mask = nc[1,], y = y[1]) mf_map(nc[1, ], col = 2, add = TRUE) mf_map(inset1, col = 2, add = TRUE) # plusieurs objets agrandis inset2 <- m_r(x = nc, mask = nc[2,], y = y[2]) mf_map(nc[2, ], col = 3, add = TRUE) mf_map(inset2, col = 3, add = TRUE) # un seul objet, réduit inset3 <- m_r(x = nc[nc$CNTY_ID == 2000, ], mask = nc[nc$CNTY_ID == 2000, ], y = y[3]) mf_map(nc[nc$CNTY_ID == 2000, ], col = 4, add = TRUE) mf_map(inset3, col = 4, add = TRUE) # plusieurs objets dans une autre proj mtq <- mf_get_mtq() inset4 <- m_r(x = mtq, mask = mtq, y = y[4]) mf_map(inset4, col = 5, add = TRUE) # bouger des points pts <- st_as_sf(st_sample(x = nc[1, ], 10)) inset5 <- m_r(x = pts, mask = nc[1,], y = y[1]) mf_map(pts, cex = .2, add=TRUE) mf_map(inset5, cex = .2, add = TRUE) # MULTIPOINT pts <- (st_sample(x = nc[3, ], 30)) pts <- st_as_sf(st_combine(x = pts)) inset6 <- m_r(x = pts, mask = nc[3,], y = y[5]) mf_map(pts, cex = .5, add = TRUE) mf_map(inset6, cex = .5, add = TRUE) # MULTILINESTRING line <- st_cast(nc, 'MULTILINESTRING') inset7 <- m_r(x = line, mask = line[4,], y = y[6]) mf_map(line[4, ], col = 3, add = TRUE, lwd = 3) mf_map(inset7, add = TRUE, col = 3, lwd = 3)library(sf) library(mapsf) nc <- st_read(system.file("shape/nc.shp", package="sf")) # Ne fonctionne qu'avec 2 fonds projetés nc <- st_transform(nc, 32119) # Créer des boites y <- st_as_sfc(st_bbox( c(xmin = 270000, ymin = 00000, xmax = 550000, ymax = 100000), crs = 32119)) y <- st_make_grid(y, n = c(3,2)) mf_map(nc) mf_map(y, add = TRUE) # plusieurs objets agrandis inset1 <- m_r(x = nc, mask = nc[1,], y = y[1]) mf_map(nc[1, ], col = 2, add = TRUE) mf_map(inset1, col = 2, add = TRUE) # plusieurs objets agrandis inset2 <- m_r(x = nc, mask = nc[2,], y = y[2]) mf_map(nc[2, ], col = 3, add = TRUE) mf_map(inset2, col = 3, add = TRUE) # un seul objet, réduit inset3 <- m_r(x = nc[nc$CNTY_ID == 2000, ], mask = nc[nc$CNTY_ID == 2000, ], y = y[3]) mf_map(nc[nc$CNTY_ID == 2000, ], col = 4, add = TRUE) mf_map(inset3, col = 4, add = TRUE) # plusieurs objets dans une autre proj mtq <- mf_get_mtq() inset4 <- m_r(x = mtq, mask = mtq, y = y[4]) mf_map(inset4, col = 5, add = TRUE) # bouger des points pts <- st_as_sf(st_sample(x = nc[1, ], 10)) inset5 <- m_r(x = pts, mask = nc[1,], y = y[1]) mf_map(pts, cex = .2, add=TRUE) mf_map(inset5, cex = .2, add = TRUE) # MULTIPOINT pts <- (st_sample(x = nc[3, ], 30)) pts <- st_as_sf(st_combine(x = pts)) inset6 <- m_r(x = pts, mask = nc[3,], y = y[5]) mf_map(pts, cex = .5, add = TRUE) mf_map(inset6, cex = .5, add = TRUE) # MULTILINESTRING line <- st_cast(nc, 'MULTILINESTRING') inset7 <- m_r(x = line, mask = line[4,], y = y[6]) mf_map(line[4, ], col = 3, add = TRUE, lwd = 3) mf_map(inset7, add = TRUE, col = 3, lwd = 3)
Map insets are small zoom-in or zoom-out maps that focus on particular
territories. mapinsetr provides a set of functions that helps to
create such insets.
Maintainer: Timothée Giraud [email protected] (ORCID)
Authors:
Nicolas Lambert [email protected]
Useful links:
Move and resize a simple feature collection of polygons, multipolygons or points.
move_and_resize(x, mask = NULL, xy, prj, k = 1)move_and_resize(x, mask = NULL, xy, prj, k = 1)
x |
an sf POINT, POLYON or MULTIPOLYGON object to resize and move. |
mask |
an sf or sfc POLYGON or MULTIPOLYGON object used to select the area to move an resize. |
xy |
coordinates used to move the inset, bottomleft corner of the inset. |
prj |
CRS string of the output projection of the inset. |
k |
factor used to resize. |
An sf object is returned.
library(sf) nc <- st_read(system.file("shape/nc.shp", package="sf")) nc <- st_transform(nc, 32119) plot(st_geometry(nc)) bb <- st_bbox(nc[nc$CNTY_ID == '2030',]) mask <- create_mask(bb = bb, add = TRUE) inset <- move_and_resize(nc, mask, xy = c(190000, 1000), k = 3) plot(st_geometry(inset), add = TRUE) plot(st_geometry(nc)) mask <- st_buffer(st_centroid(st_geometry(nc[nc$CNTY_ID == 2026,])), dist = 30000) plot(st_geometry(mask), border = "red", lwd = 2, add = TRUE) inset <- move_and_resize(nc, mask, xy = c(270000, 5000), k = 2.5) plot(st_geometry(inset), add = TRUE)library(sf) nc <- st_read(system.file("shape/nc.shp", package="sf")) nc <- st_transform(nc, 32119) plot(st_geometry(nc)) bb <- st_bbox(nc[nc$CNTY_ID == '2030',]) mask <- create_mask(bb = bb, add = TRUE) inset <- move_and_resize(nc, mask, xy = c(190000, 1000), k = 3) plot(st_geometry(inset), add = TRUE) plot(st_geometry(nc)) mask <- st_buffer(st_centroid(st_geometry(nc[nc$CNTY_ID == 2026,])), dist = 30000) plot(st_geometry(mask), border = "red", lwd = 2, add = TRUE) inset <- move_and_resize(nc, mask, xy = c(270000, 5000), k = 2.5) plot(st_geometry(inset), add = TRUE)