Наложение растровых графиков с точками, линиями и полигонами должно работать нормально, как показано в следующем примере.
Мое лучшее предположение состоит в том, что объекты Spatial*
, которые вы пытаетесь нанести поверх растра, попадают за пределы отображаемой области. Вы проверили, что оба объекта raster
и Spatial*
находятся в одной и той же CRS и (при условии, что они есть) перекрываются ли ограничивающие рамки? (т.е. попробуйте bbox(shp)
и bbox(ras)
и сравните результаты).
library(rgdal)
library(raster)
# Create a raster
ras <- raster(ncols=36, nrows=18)
ras[] <- runif(ncell(ras))
# Create a SpatialPoints object
shpPts <- spsample(Spatial(bbox=bbox(ras)), 20, type="random")
# Create a SpatialPolygons object
p1 <- rbind(c(-10,0), c(140,60), c(160,0), c(140,-55), c(-10,0))
shpPolys <- SpatialPolygons( list(Polygons(list(Polygon(p1)), 1)))
# Plot them, one layer after another
plot(ras)
plot(shpPts, pch=16, col="red", add=TRUE)
plot(shpPolys, col="yellow", add=TRUE)
person
Josh O'Brien
schedule
14.02.2012