1

我想将我的地图的放大部分添加到原始地图中,并作为最终产品,一张显示原始地图和放大/缩放部分的地图。以meuse数据集为例:

library(sp)
library(lattice) 
data(meuse)
coordinates(meuse)=~x+y
gridded(meuse)<-TRUE
rasters.m<-list()
for (i in 1:12){
rasDF <- raster(meuse, layer=i)
rasters.m[[i]]<-rasDF
}
stack.sp<-stack(rasters.m)
 plot(stack.sp) # gives a gridded view of the stacked rasters. But now I would like to zoom in..
 zoom.ent<-zoom(stack.sp,1) #  The zoomed in portion appears as a new window, with the boundaries of the zoomed area highlighted in red on the original map.

我不确定rasterorrasterVIS包中是否有命令可以将栅格的放大部分添加到原始地图上。我已经尝试过该par功能,但这不起作用。任何建议都会受到欢迎。

4

1 回答 1

2

这或多或少与您在此处提出的问题相同。对于Raster*对象,您必须使用该shift功能。结果可以结合包的+.trellis功能latticeExtra

library(raster)
library(rasterVis)

f <- system.file("external/test.grd", package="raster")
r <- raster(f)

rZoom <- crop(r, extent(180000,  181000, 330000, 331500))

displaced <- shift(rZoom, x = -1200, y = 2000)

levelplot(r) + levelplot(displaced)

转移

于 2014-11-10T08:20:07.013 回答