0

我正在尝试将多边形的插图添加到地图中。

这个 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()

这是我的代码能走多远,但我仍然无法添加此插图:

4

0 回答 0