10
如何繪製shapefile對象後面的「柵格」對象?對自己的兩個情節不錯,但點不超過光柵情節:在shapefile背後繪製柵格
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)
良好的通話,我的光柵被錯誤地定位。感謝您確認這可以正常工作。 – Benjamin 2012-02-14 20:15:18