Skip to content

Instantly share code, notes, and snippets.

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

Timothée Giraud rCarto

🛠️
...
255 followers · 69 following
View GitHub Profile
@rCarto
rCarto / iris75.geojson
Last active July 11, 2018 14:08
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@rCarto
rCarto / popcircle.R
Created April 5, 2019 11:48
library(rnaturalearth)
library(sf)
library(wbstats)
library(popcircle)
library(cartography)
# Get countries
ctry <- ne_countries(scale = 50, returnclass = "sf")
ctry <- st_transform(ctry, 54030)
@rCarto
rCarto / tanakapop.R
Created March 12, 2019 14:15
library(raster)
library(sf)
library(viridis)
library(cartography)
library(tanaka)
temp <- tempfile()
data_url <- "http://cidportal.jrc.ec.europa.eu/ftp/jrc-opendata/GHSL/GHS_POP_GPW4_GLOBE_R2015A/GHS_POP_GPW42015_GLOBE_R2015A_54009_1k/V1-0/GHS_POP_GPW42015_GLOBE_R2015A_54009_1k_v1_0.zip"
download.file(data_url, temp)
unzip(temp, exdir = "pop")
@rCarto
rCarto / Penciling.R
Last active May 6, 2019 15:50
https://rgeomatic.hypotheses.org/?p=1333
## Figure 1
#download the dev version of cartography
devtools::install_github("riatelab/cartography")
library(sf)
library(cartography)
mtq <- st_read(system.file("shape/martinique.shp", package="cartography"))
mt1 <- mtq[3,]
mt2 <- getPencilLayer(mt1, size = 1000, lefthanded = T, buffer = 2000)
png(filename = "fig1.png", width = 474, height = 132, res = 100)
par(mfrow=c(1,3), mar = c(0,0,0,0))
@rCarto
rCarto / pencil2019.R
Created July 4, 2019 08:24
library(sf)
library(cartography)
library(png)
# import background image
if (!file.exists("img/background.png")) {
githubURL <- "https://raw.githubusercontent.com/gadenbuie/ggpomological/master/inst/images/pomological_background.png"
download.file(githubURL, "img/background.png")
}
img <- readPNG("img/background.png")
# import Communes of Martinique (sf dataframe within cartography package)
@rCarto
rCarto / disc2019.R
Created July 4, 2019 08:35
library(cartography)
library(sp)
# Load data
data(nuts2006)
# Get a SpatialLinesDataFrame of countries borders
nuts0.contig.spdf <- getBorders(nuts0.spdf)
# Get the GDP per capita
nuts0.df$gdpcap <- nuts0.df$gdppps2008/nuts0.df$pop2008*1000000
@rCarto
rCarto / carto_facet.R
Created November 6, 2019 16:01
# World basemap
download.file(url = "https://raw.githubusercontent.com/riatelab/basemaps/master/World/countries.geojson",
destfile = "country.geojson")
# Graticules layer
download.file(url = "https://raw.githubusercontent.com/riatelab/basemaps/master/World/graticule30.geojson",
destfile = "graticule.geojson")
# Population data base
download.file(url = "https://population.un.org/wup/Download/Files/WUP2018-F12-Cities_Over_300K.xls",
destfile = "citypop.xls")
@rCarto
rCarto / localpeaks.R
Created December 2, 2019 13:15
library(SpatialPosition)
library(raster)
data(hospital)
# Compute Stewart potentials from known points (hospital) on a
# grid defined by its resolution
pot <- stewart(knownpts = hospital, varname = "capacity",
typefct = "exponential", span = 750, beta = 3,
resolution = 100, mask = paris)
# Create a raster of potentials values
ras <- rasterStewart(x = pot)
@rCarto
rCarto / interactiveposition.R
Created February 28, 2020 13:21
library(cartography)
library(sf)
mtq <- st_read(system.file("gpkg/mtq.gpkg",
package="cartography"))
bks <- getBreaks(v = mtq$MED, nclass = 5,
method = "quantile")
cols <- carto.pal(pal1 = "turquoise.pal", n1 = 5)
choroLayer(x = mtq, var = "MED",
breaks = bks, col = cols,
@rCarto
rCarto / mapbox.md
Created April 11, 2020 12:16 
library(cartography)
library(sf)
#> Linking to GEOS 3.7.1, GDAL 2.4.0, PROJ 5.2.0

# create an sf object centered on Tizi Ouzou city
tizi <- data.frame(name = "Tizi Ouzou", longitude = 4.06, latitude = 36.71)
tizi <- st_as_sf(tizi, coords = c("longitude", "latitude"), crs = 4326)
tizi <- st_transform(tizi, 3857)
tizi &lt;- st_buffer(tizi, 5000)