Я просто пытаюсь сделать простую карту области исследования, которая также содержит врезку штата, в котором я работаю (Северная Каролина). Я хотел бы преобразовать карту-врезку в объект grob, чтобы нанести его на карту основной области исследования, а затем использовать ggsave, чтобы сохранить карту в виде изображения, pdf и т. д. Я использую шейп-файлы для моей фактической карты, но я буду показать вам, что я пытаюсь использовать map_data:
library(ggplot2)
library(ggmap)
library(maps)
library(mapdata)
library(gridExtra)
library(grid)
# get the NC data:
states <- map_data("state")
nc_df <- subset(states, region == "north carolina")
# study area map:
nc_base <- ggplot() +
geom_polygon(data = nc_df, aes(x = long, y = lat, group = group), fill="grey", color="black") +
coord_fixed(xlim=c(-80, -77.5), ylim=c(33.5, 34.9), ratio = 1.3) +
theme_bw()
nc_base
# inset map:
insetmap<-ggplot() +
geom_polygon(data = nc_df, aes(x = long, y = lat, group = group), fill="grey", color="black") + # get the state border back on top
coord_fixed(ratio = 1.3) +
annotate(geom = "rect", ymax = 34.9, ymin = 33.5, xmax = -77.5, xmin = -80, colour = "red", fill = NA) +
ylab("") +
xlab("") +
theme_nothing()
insetmap
insetmap.grob <- ggplotGrob(insetmap)
final_map <- nc_base + annotation_custom(insetmap.grob, xmin=-79.5, xmax=-79, ymin=33.75, ymax=34)
final_map
Когда я запускаю сценарий для создания окончательной карты, создается только карта изучаемой области. Мне интересно, неправильно ли я использую ggplotGrob или это что-то еще? Возможно, я где-то читал, что функция annotation_custom не работает, если вы не используете функцию coord_cartesian в ggplot2 (и здесь я использую coord_fixed). Если это так, могу ли я увеличить масштаб аналогичным образом с помощью этой функции, или есть другая функция coord_ для увеличения карты моей изучаемой области?
Спасибо, Джей