Skip to content

Instantly share code, notes, and snippets.

View rCarto's full-sized avatar
🛠️
...

Timothée Giraud rCarto

🛠️
...
View GitHub Profile
@rCarto
rCarto / mtq.R
Last active August 21, 2023 19:06
View mtq.R
library(raster)
library(cartography)
library(sf)
library(SpatialPosition)
mtq <- st_read(system.file("shape/martinique.shp", package="cartography"))
# use WGS84 proj
mtq_latlon <- st_transform(mtq, 4326)
# this call throw an error but seems to work...
getData('SRTM', lon=-61, lat=14)
@rCarto
rCarto / zoom.R
Last active March 3, 2023 22:04
View zoom.R
# load libraries
library(sf)
library(cartography)
library(rnaturalearth)
# import world countries and coastlines from Natural Earth
country <- ne_countries(scale = 50, returnclass = "sf",
type="sovereignty")
coastline <- ne_coastline(scale = 50, returnclass = "sf")
View roads.R
library(sf)
library(cartography)
library(osrm)
# set osrm to my own server
options(osrm.server = "http://address.of.my.server/", osrm.profile = "driving")
# destination point (useR2019 conference in Toulouse)
dst <- data.frame(id="dst", x = 1.4344, y = 43.6113)
dst <- st_as_sf(dst, coords = c('x','y'), crs = 4326)
@rCarto
rCarto / cartomix.R
Last active July 13, 2022 02:19
Script to build the cartomix figure
View cartomix.R
library(cartography)
library(sp)
library(sf)
# Load data
data(nuts2006)
# Save image
sizes <- getFigDim(x = nuts0.spdf, width = 700, mar = c(0,0,0,0), res = 100)
png('./img/map8.png', width = sizes[1], height = sizes[2], res = 100)
View fisheye-world.R
library(rnaturalearth)
library(mapsf)
library(fisheye)
ct <- ne_download(scale = 50, returnclass = "sf")
oc <- ne_download(type = "ocean", category = "physical",
scale = 50, returnclass = "sf")
x <- st_transform(ct, "+proj=natearth2")
z <- st_transform(oc, "+proj=natearth2")
View graticules.md
library(mapsf)
#> Le chargement a nécessité le package : sf
#> Linking to GEOS 3.7.1, GDAL 3.1.2, PROJ 7.1.0
mtq <- mf_get_mtq()
mf_theme("nevermind")
mf_map(mtq, graticule = st_crs(4326), 
       axes = TRUE, 
       lon = seq(-60,-62, by=-.2), 
 lat = seq(14, 15, by = .2))
View osrm.R
library(sf)
library(osrm)
library(maptiles)
# build a bbox for Paris
bb <- st_bbox(c(xmin = 643069, ymin = 6857478,
xmax = 661079, ymax = 6867081),
crs = 2154)
# get map tiles
osm <- get_tiles(x = bb, provider = "CartoDB.PositronNoLabels",
crop = TRUE, zoom = 13)
@rCarto
rCarto / bang.R
Last active November 27, 2020 12:22
View bang.R
library(sf)
library(elevatr)
library(raster)
library(tanaka)
# a polygon of bangalore
bangalore <- st_read("export.geojson")
elevation <- get_elev_raster(bangalore, z = 10)
bangalore_elevation <- trim(mask(elevation, bangalore))
# inspect min and max values
@rCarto
rCarto / points.R
Last active November 6, 2020 09:57
View points.R
library(units)
library(sf)
library(cartography)
# library(osmdata)
#
# # define a bounding box
# q0 <- opq(bbox = c(2.2247, 48.8188, 2.4611, 48.9019))
#
# # extract Paris boundaries
View lakes.R
library(rnaturalearth)
library(cartography)
library(sf)
lakes <- ne_download(scale = 10, type = "lakes", category = c("physical"),
destdir = tempdir(), load = TRUE, returnclass = c("sf"))
countries <- ne_download(scale = 10, type = "countries",
category = c("cultural"), destdir = tempdir(),
load = TRUE, returnclass = c("sf"))
countries <- st_transform(countries, "ESRI:54017")