Я нашел похожие ответы на такие вопросы, как этот, но большинство из них использует пакеты rworldmap, ggmap, ggsubplot или geom_subplot2d. См. Например здесь или здесь.
Я хотел бы знать, как я могу нарисовать другие объекты ggplot, такие как диаграмма на карте, которая создается из шейп файла. Тот, который я использую, можно загрузить здесь.
ИЗМЕНИТЬ
Как правильно указал @beetroot, новый файл, который можно загрузить по ссылке выше, значительно изменился. Поэтому изменяются имена шейп файла и т.д.
library(rgdal)
library(ggplot2)
library(rgeos)
library(maptools)
map.det<- readOGR(dsn="<path to your directory>/swissBOUNDARIES3D100216/swissBOUNDARIES3D/V200/SHAPEFILE_LV03", layer="VECTOR200_KANTONSGEBIET")
map.kt <- map.det[[email protected]$KANTONSNUM=="CH01000000"|[email protected]$KANTONSNUM=="CH19000000",]
#get centroids
map.test.centroids <- gCentroid(map.kt, byid=T)
map.test.centroids <- as.data.frame(map.test.centroids)
map.test.centroids$KANTONSNR <- row.names(map.test.centroids)
#create df for ggplot
kt_geom <- fortify(map.kt, region="KANTONSNUM")
#Plot map
map.test <- ggplot(NULL)+
geom_polygon(data=kt_geom, aes(long, lat, group=group), fill="white")+
coord_fixed()+
geom_path(data=kt_geom, color="gray48", mapping=aes(long, lat, group=group), size=0.2)+
geom_point(data=map.test.centroids, aes(x=x, y=y), size=9, alpha=6/10)
mapp
В результате получается такая карта. Все идет нормально.
Однако у меня возникают трудности с объединением двух графиков, таких как карта map.test и, например, эта:
geo_data <- data.frame(who=rep(c(1:2), each=2),
value=as.numeric(sample(1:100, 4, replace=T)),
KANTONSNR=rep(c(1,19), 2))
bar.testplot <- ggplot()+
geom_bar(data=geo_data, aes(factor(id),value,group=who),position='dodge',stat='identity')
Баррикады должны лежать в центре двух полигонов, т.е. где находятся две точки. Я мог бы создавать штрих-коды и рисовать их на карте отдельно, если это облегчает процесс.
