CODEWITHSUNDEEP
rr-rasterr-sp
user2507608
user2507608

Reputation: 365

Merging 2 raster images into one plot using the `raster` package

I would like to add an enlarged portion of my map to the original map and have as the final product, one map that shows both the original map, and also the enlarged/zoomed portion. Using the meuse dataset as an example:

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.

I am not sure if there is a command in the raster or rasterVIS packages that would allow one to add the zoomed in part of the raster onto the original map. I have tried the par function but that doesn't work. Any suggestions would be welcomed.

Upvotes: 1

Views: 1451

Answers (1)

Oscar Perpi&#241;&#225;n
Oscar Perpi&#241;&#225;n

Reputation: 4511

This is more or less the same question you asked here. For Raster* objects you have to use the shift function. The result can be combined with the +.trellis function of the latticeExtra package:

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)

shift

Upvotes: 2

Related Questions