library(raster)
library(sf)
library(ggplot2)
library(viridis)
library(rgdal)
canada_raw <- readOGR(dsn = "Canada_Map_1/Canada.shp", layer = "Canada",
use_iconv=TRUE, encoding="CP1250")
can_adm <- sf::st_as_sf(canada_raw)
can_simp <- rmapshaper::ms_simplify(can_adm)
str(can_simp)
can_proj <- "+proj=lcc +lat_1=49 +lat_2=77 +lat_0=40 +lon_0=-91.86666666666666 +x_0=0 +y_0=0"
can_simp <- sf::st_transform(can_simp, can_proj)
can_simp$Popu<-1:13
ggplot() +
geom_sf(data = can_simp, aes(fill = Popu)) +
scale_fill_continuous(low="red",high="blue")+
theme_minimal()
image.png
ggplot() +
geom_sf(data = can_simp, aes(fill = Popu)) +
scale_fill_continuous(low="red",high="blue")
image.png
我原来用的R版本是3.5.1,
sf
包安装的时候一直遇到报错。
后来换成3.6.1就直接安装好了