我正在尝试将多边形的插图添加到地图中。
这个 pdf 的第 11 页(https://journal.r-project.org/archive/2013-1/kahle-wickham.pdf)显示了我想要得到的东西。不幸的是,我是 R 的新手,即使使用这个示例也无法理解它是如何工作的......在我的例子中,它显示了这个错误: f(...) 中的错误:找不到我不理解的函数“viewport” . 我的代码
swd <- read.csv2("C:/MapyR/rozmieszczenieswd/swdwsp.csv",header= TRUE, sep = ";", encoding = "unknow")
swd$Lat <- as.numeric((strtrim(swd$Lat,8)))
swd$Long <- as.numeric((strtrim(swd$Long,8)))
polska <- get_map(
location=c("Olsztyn"),
zoom=12,
maptype="roadmap"
)
olsztyn <- subset(swd, Lat>=53.72 & Lat<= 53.84 & Long>=20.40 & Long<=20.55,select = rodzaj:Long)
mapa <- ggmap(polska, extent = "device", legend ="bottomright")
punkty <- stat_density_2d( aes(x=Long,y=Lat,fill= ..level..,alpha=..level..)
,data=olsztyn
,size=2
,geom = "polygon"
)
nowa <- mapa + inset(
grob = ggplotGrob(ggplot()+punkty + theme_inset()),
xmin = 20.50, xmax = 20.55, ymin = 53.72, ymax = 53.76 #I tryed many versions of x,y
)
sciezka <- file.path("C:","MapyR","proby density",paste("inset",".png"), sep="")
png(file=sciezka, width = 2500 , height=1567)
print(nowa)
dev.off()
这是我的代码能走多远,但我仍然无法添加此插图: