, . UTM, Easting Northing , . lat/lon. , :
library(ggplot2)
library(tidyverse)
usa <- map_data("usa", )
shape <- spTransform(shape, CRS("+proj=longlat +datum=WGS84"))
shape_map <- fortify(shape, region="Name")
colnames(shape_map) <- c("long", "lat", "order", "hole", "piece", "region", "group")
prop.test <- proptest.result[which(proptest.result$variable=="Upward N"),]
ggplot() +
geom_map(
data=usa, map=usa, aes(long, lat, map_id=region),
color="#2b2b2b", fill="#00000000"
) +
geom_map(
data=shape_map, map=shape_map,
aes(long, lat, map_id=region)
) +
geom_map(
data=filter(prop.test, season=="DJF"),
map=shape_map, aes(fill=prop.mega, map_id=megaregion)
) +
viridis::scale_fill_viridis(direction=-1) +
coord_map("polyconic") +
ggthemes::theme_map()
: