kidwellj/mapping_environmental_action
1
0
Fork 0
mirror of https://github.com/kidwellj/mapping_environmental_action.git synced 2024-11-01 07:52:21 +00:00
Code Issues Packages Projects Releases Wiki Activity

transitioning to sf, working on plots

Browse source
This commit is contained in:
Jeremy Kidwell 2018-11-26 15:18:27 +00:00
parent cd6793b4f9
commit dac4a9c3a4
1 changed files with 10 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

15
mapping_draft.Rmd
View file

@ -261,15 +261,20 @@ Perhaps the first important question to ask of these groups is, where are they?
# Row 1 plot using polygons from admin_lev1 and row 2 plot using ploygons from admin_lev2 # Row 1 plot using polygons from admin_lev1 and row 2 plot using ploygons from admin_lev2
# 3. Need to clip choropleth polygons to buildings shapefile # 3. Need to clip choropleth polygons to buildings shapefile
savePlot <- function(myPlot) {
pdf(file="figures/admin_choropleth_ecs.pdf")
print(myPlot)
dev.off()
}
myplot <- ggplot() + geom_sf(data = admin_lev1_sf) + geom_point(data=as.data.frame(ecs), aes(x=X, y=Y))
# + geom_point(data=as.data.frame(ecs), aes(x=X, y=Y))
savePlot(myplot)
ggplot() +
geom_sf(data = admin_lev1_sf)
# geom_point(data=as.data.frame(ecs), aes(x=X, y=Y))
pdf(file="figures/admin_choropleth_ecs.pdf", width=4, height=4)
var01 <- admin_lev1$ecs_count var01 <- admin_lev1$ecs_count
bins <- unique(quantile(var01, seq(0,1,length.out=30))) bins <- unique(quantile(var01, seq(0,1,length.out=30)))
admin_lev1$binId01 <- findInterval(var01, bins) admin_lev1$binId01 <- findInterval(var01, bins)