シェープファイル オブジェクトの背後に「ラスター」オブジェクトをプロットするにはどうすればよいですか? 両方とも単独でうまくプロットされますが、ポイントはラスタ上にプロットされません:
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)
次の例に示すように、点、線、および多角形を使用したラスター プロットのオーバープロットは問題なく機能するはずです。
私の最善の推測は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)
