Построение растра за шейпфайлом

Как я могу построить объект "растр" за объектом шейп-файла? Оба объекта строятся нормально сами по себе, но точки не строятся поверх растрового объекта:

require(rgdal)
require(maptools)
require(raster)

myproj = "+proj=utm +zone=12 +north +ellps=WGS84 +units=m"
shp = readShapeSpatial(fn.shp, proj4string = CRS(myproj))
ras = raster(fn.tif)

plot(ras)
plot(shp, bg="transparent", add=TRUE)
11
задан Benjamin 14 February 2012 в 17:17
поделиться