rm(list = ls(all.names = TRUE))
library(geojsonio)
##
## Attaching package: 'geojsonio'
## The following object is masked from 'package:base':
##
## pretty
library(mapview)
workdirectory = "~/OneDrive/r-files/"
setwd(paste(workdirectory, "/TCC/", sep = ""))
densidade <- geojson_read("alvaraspoly.json", what = "sp")
densidade$group = NULL; densidade$area = NULL; densidade$ocorrencias = NULL;densidade$id = NULL
length(densidade$box_id)
## [1] 136
mapviewOptions(basemaps = "OpenStreetMap",
na.color = "#ffffff",
layers.control.pos = "topright")
densidade50 = subset(densidade, box_id > quantile(densidade$box_id, 0.50) )
length(densidade50$box_id)
## [1] 68
densidade25 = subset(densidade, box_id > quantile(densidade$box_id, 0.75))
length(densidade25$box_id)
## [1] 34
mapview(densidade, zcol = 'box_id', at=seq(1000, 48000, 100),legend = TRUE, layer.name = 'Alvaras por KM²')
mapview(densidade50, zcol = 'box_id', at=seq(1000, 48000, 500),legend = TRUE, layer.name = 'Alvaras por KM² <br>Áreas com os 50% mais densos')
mapview(densidade25, zcol = 'box_id', at=seq(3000, 48000, 500), legend = TRUE, layer.name = 'Alvaras por KM² <br>Áreas com os 25% mais densos')