从shapefile中裁剪图

Question

我正在尝试使用欧洲统计局提供的shapefile和数据生成Choroplete地图。 shapefile已下载here：使用this post中的JD Longs代码。

这是重现我在下面发布的数字的最小代码。

library(maptools)
tmpdir <- tempdir()
url <- 'http://ec.europa.eu/eurostat/cache/GISCO/geodatafiles/NUTS_2010_03M_SH.zip'
file <- basename(url)
download.file(url, file)
unzip(file, exdir = tmpdir )
shapeFile <- paste(tmpdir,"/Shape/data/NUTS_RG_03M_2010", sep="")

EU <- readShapeSpatial(shapeFile)
plot(EU)

我的问题是我希望情节区域只关注欧洲，但由于海外地区（法国和西班牙），情节没有正确的焦点。在上面的例子中，有一种简单的方法可以“裁剪”绘图区域吗？

我想摆脱的多边形是“Country_Shape”的一部分，所以过滤它们是没有选择的。 我试图通过在plot命令中定义xlim和ylim参数来实现我的目标，但没有成功。我使用locator（）从图形设备获取坐标，但是插入值并没有提供想要的结果。

Answer 1

shapefile投影文件（Shape / data / NUTS_RG_03M_2010.prj）显示以下内容：

  

GEOGCS [＆＃34; GCS_ETRS_1989＆＃34;，DATUM [＆＃34; D_ETRS_1989＆＃34;，SPHEROID [＆＃34; GRS_1980＆＃34;，6378137.0,298.257222101]，PRIMEM [＆＃34;格林威治＆＃ 34;，0.0]，单位[＆＃34;度＆＃34;，0.0174532925199433]]

.prj文件中包含的信息指定Shapefile中几何数据的地理坐标系。为此坐标系尝试适当的xlim和ylim坐标：

R> plot(EU, ylim=c(30, 70), xlim=c(-10, 44))
R> box()