r - R Levelplot 在 RasterLayer 投影上覆盖两个 SpatialPolygonsDataFrame

Question

也在gis.stackexchange上发布

我在 R 的水平图中绘制 RasterLayer 和 SpatialPolygonsDataFrame 时遇到了一些问题。投影有问题，但我不明白是什么。

在这里您可以找到我 可重现的示例数据和代码：

library(sp)
library(rgdal)
library(rgeos)
library(raster)
library(rasterVis)
library(maptools)

setwd("C:/...path_to/test2")
data<-read.csv("test.csv", header=TRUE)

#creating the raster from a data.frame and giving projection
raster<-rasterFromXYZ(data, crs="+proj=longlat")
raster_proj<-projectRaster(raster, crs="+proj=longlat +datum=WGS84 
+ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs")

#loading shapes
nuts<-readShapePoly("NUTS_RG_WGS84.shp", proj4string=CRS("+proj=longlat"))
countries<-readShapePoly("countries.shp", proj4string=CRS("+proj=longlat"))

#subsetting the raster
raster_clip<- crop(raster_proj, countries, snap="near")

#plot
p.strip <- list(cex=1.5, lines=1, fontface='bold') 
x.scale <- list(cex=1.5)
y.scale <- list(cex=1.5)
label <- list(labels=list(width=1, cex=1.5), height=0.95)
levelplot(raster_clip, par.settings = RdBuTheme, margin=FALSE,
      at=seq(min(na.omit(values(raster_clip))), max(na.omit(values(raster_clip))), length.out=15),scales=list(x=x.scale, y=y.scale), 
      par.strip.text=p.strip, colorkey=label, xlab=list(label="Longitude", cex=2), ylab=list(label="Latitude", cex=2))+ 
layer(sp.polygons(nuts))+layer(sp.polygons(countries))

错误信息是：

Error: Attempted to create layer with no geom.

但是栅格和多边形在同一个投影中。我还尝试使用以下方法重新投影多边形：

nuts <- spTransform(nuts, CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs"))
countries <- spTransform(countries, CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs"))

任何想法？

Answer 1

如果您不想或不能分离 ggplot2 因为您需要它来做其他事情，您也可以layer使用以下方法从 latticeExtra 包中专门调用该函数：

levelplot(myRaster) + 
  latticeExtra::layer(sp.polygon(myPolygon))

Answer 2

detach("package:ggplot2", unload=TRUE)

Answer 3

我通过从我的库中分离 ggplot2 解决了这个问题。