我正在尝试将多边形的插图添加到地图中。
这个pdf(https://journal.r-project.org/archive/2013-1/kahle-wickham.pdf)的第11页显示了我想要获得的内容。不幸的是我是R的新手并且不知道它是如何工作的,即使是这个例子...在我的情况下它显示了这个错误:f(...)中的错误:找不到函数" viewport&#34 ;我不明白。 我的代码
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()
这是我可以用我的代码获得的,我还是不能添加这个插图: