Hi,
maybe I've overlooked something obvious but I'm trying to generate _one_ plot of a world map (actually projected) with an inset of a specific zoomed area based on SpatialDataFrames.
I followed this example which makes usage of the library "maps":
http://wiki.cbr.washington.edu/qerm/sites/qerm/images/7/78/MakingAnInsetMapShorter.r
Thus I wrote this (as a minimal and na?ve example):
library(maptools)
data(wrld_simpl)
# plot world map
plot(wrld_simpl, col = "khaki", bg = "azure2")
# define inset map location and size
par(plt = c(0.57, 0.87, 0.4, 0.7), new = TRUE)
# init inset map view port
plot.window(xlim = c(130, 180), ylim = c(40, 70))
# draw background
rect(130, 40, 180, 70, col = "azure2")
# plot zoomed inset map
plot(wrld_simpl, xlim = c(130, 180), ylim = c(40, 70),
col = "khaki", bg = "azure2", add = TRUE)
The problem is that the inset map won't be cropped, i.e. you see both maps overlapped.
What I'm doing wrong? Or is this with SpatialDataFrames not possible?
I'd be appreciated for any hint!
Kind regards,
--Hans
PS = sessionInfo():
R version 2.14.2 (2012-02-29)
Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] maptools_0.8-14 lattice_0.20-0 sp_0.9-96 foreign_0.8-49
loaded via a namespace (and not attached):
[1] grid_2.14.2 tools_2.14.2