helpful when using mapshot/saveWidget + webshot to ensure map is zoomed to desired level
### zoom to after rendering (relevant for mapshot)
library(mapview)
map = mapview(franconia)@map
tim-salabim
tim-salabim
/ franconia_elevation_code
Last active
December 19, 2019 16:49
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| library(raster) | |
| library(sf) | |
| library(sp) | |
| library(elevatr) | |
| library(mapview) | |
| library(magick) | |
| library(smoothr) | |
| library(RColorBrewer) | |
| # download from http://www.geodaten.bayern.de/opendata/Verwaltungsgebiete_shp_epsg31468.zip |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| ghmd2pdf = function(input, output = NULL, pandoc_opts = NULL) { | |
| # input = "01_analysis/0110_code/01101_r_scripts/030_yield_by_oven.md" | |
| if (is.null(output)) { | |
| output = paste0(tools::file_path_sans_ext(input), ".pdf") | |
| } | |
| lns = readLines(input) | |
| idx = grep("<img src=", lns) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| viewNodes = function(x, cex = 6) { | |
| crds = sf::st_coordinates(x)[, c("X", "Y")] | |
| n = nrow(crds) | |
| pts = sf::st_as_sf( | |
| data.frame( | |
| id = 1:n, | |
| x = crds[, "X"], | |
| y = crds[, "Y"] | |
| ), | |
| coords = c("x", "y"), |
tim-salabim
/ buildings_ch.R
Created
October 15, 2018 16:57
Visualising 1.8 mio buldings in switzerland using leaflet.glify
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| library(sf) | |
| library(leaflet) | |
| library(leaflet.glify) | |
| library(mapview) | |
| bu = read_sf("/media/timpanse/d8346522-ef28-4d63-9bf3-19fec6e13aab/bu_lenovo/software/testing/mapview/switzerland/buildings.shp") | |
| st = seq(1, nrow(bu), 5e5) | |
| nd = st + 5e5 - 1 | |
| nd[length(nd)] = nrow(bu) |
tim-salabim
/ on_render.md
Last active
February 17, 2019 19:47
various snippets for modifying map state after creation
tim-salabim
/ glify fiddle
Created
September 20, 2018 18:02
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| https://jsfiddle.net/7yf6upg5/4/ |
This file has been truncated, but you can view the full file.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| var cols = | |
| [{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.3294},{"r":0.2667,"g":0.0039,"b":0.32 |
tim-salabim
/ leaflet.pixi_20mio
Created
August 28, 2018 21:17
render millions of points with leaflet.pixi
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| # devtools::install_github("tim-salabim/leaflet.pixi") | |
| library(mapview) | |
| library(leaflet) | |
| library(leaflet.pixi) | |
| library(sf) | |
| n = 1e6 | |
| df1 = data.frame(x = runif(n, -100, 30), |
tim-salabim
/ nearest_neighbor_graph
Last active
May 21, 2020 12:55
R version of https://beta.observablehq.com/@mbostock/nearest-neighbor-graph
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| library(sf) | |
| library(data.table) | |
| library(mapview) | |
| n = 1000 | |
| x = runif(n, 0, 2) | |
| y = runif(n, 0, 1) | |
| pts = st_as_sf(data.frame(id_pt = 1:n, x = x, y = y), coords = c("x", "y")) |
tim-salabim
/ merge_sf_buffers
Last active
June 13, 2018 06:15
Solution for https://aghaynes.wordpress.com/2018/06/11/merging-spatial-buffers-in-r/
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| library(osmdata) | |
| library(mapview) | |
| library(sf) | |
| q0 <- opq(bbox = "Bern, Switzerland", timeout = 60) | |
| q1 <- add_osm_feature(q0, key = 'building', value = "hospital") | |
| x <- osmdata_sf(q1) | |
| pols_ctr = st_centroid(x$osm_polygons) | |
| buff = st_buffer(pols_ctr, dist = 0.0015) |