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