forked from robsalasco/censo_2017_geojson_chile
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathscript.R
45 lines (29 loc) · 2.26 KB
/
script.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
library(sf)
library(dplyr)
library(rmapshaper)
simplygeojson <- function(x, name, keep){
shp <- ms_simplify(x, keep=keep)
write_sf(shp, dsn = paste0(name, "/","R", sprintf("%02d", as.numeric(unique(x$REGION))),".geojson"), delete_dsn = FALSE, overwrite=TRUE)
}
data1 <- lapply(paste0("R",sprintf("%02d",1:15)), function(x){read_sf(paste0("/Users/robsalasco/Documents/CENSO2017/Cartografia/",x,"/Comuna.shp"))})
lapply(data1, function(x) {simplygeojson(x, "Comunas", 0.2)})
data2 <- lapply(paste0("R",sprintf("%02d",1:15)), function(x){read_sf(paste0("/Users/robsalasco/Documents/CENSO2017/Cartografia/",x,"/Limite_Urbano_Censal.shp"))})
lapply(data2, function(x) {simplygeojson(x, "Limite_Urbano", 0.2)})
data3 <- lapply(paste0("R",sprintf("%02d",1:15)), function(x){read_sf(paste0("/Users/robsalasco/Documents/CENSO2017/Cartografia/",x,"/Provincia.shp"))})
lapply(data3, function(x) {simplygeojson(x, "Provincias", 0.2)})
data4 <- lapply(paste0("R",sprintf("%02d",1:15)), function(x){read_sf(paste0("/Users/robsalasco/Documents/CENSO2017/Cartografia/",x,"/Manzana_Precensal.shp"))})
lapply(data4, function(x) {simplygeojson(x, "Manzanas", 0.2)})
data5 <- lapply(paste0("R",sprintf("%02d",1:15)), function(x){read_sf(paste0("/Users/robsalasco/Documents/CENSO2017/Cartografia/",x,"/Region.shp"))})
lapply(data5, function(x) {simplygeojson(x, "Regiones", 0.2)})
data6 <- lapply(paste0("R",sprintf("%02d",1:15)), function(x){read_sf(paste0("/Users/robsalasco/Documents/CENSO2017/Cartografia/",x,"/Distrito_Censal.shp"))})
lapply(data6, function(x) {simplygeojson(x, "Distritos", 0.2)})
data7 <- lapply(paste0("R",sprintf("%02d",1:15)), function(x){read_sf(paste0("/Users/robsalasco/Documents/CENSO2017/Cartografia/",x,"/Zona_Censal.shp"))})
lapply(data7, function(x) {simplygeojson(x, "Zonas", 0.2)})
chl <- lapply(paste0("R",sprintf("%02d",1:15)), function(x){read_sf(paste0("/Users/robsalasco/Documents/CENSO2017/Cartografia/",x,"/Comuna.shp"))})
chl <- do.call(rbind, chl)
s.chl1 <- ms_simplify(chl, keep=0.001)
write_sf(s.chl1, dsn = "Extras/CHILE1.geojson", delete_dsn = FALSE)
s.chl2 <- ms_simplify(chl, keep=0.01)
write_sf(s.chl2, dsn = "Extras/CHILE2.geojson", delete_dsn = FALSE)
s.chl3 <- ms_simplify(chl, keep=0.1)
write_sf(s.chl3, dsn = "Extras/CHILE3.geojson", delete_dsn = FALSE)