r - 从 shapefile 裁剪图

Question

我正在尝试使用 Eurostat 提供的 shapefile 和数据制作 Choroplete map 。 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”的一部分，因此无法过滤掉它们。我试图通过在绘图命令中定义 xlim 和 ylim 参数来实现我的目标，但没有成功。我使用 locator() 从图形设备获取坐标，但插入值并没有提供想要的结果。

Answer 1

shapefile 投影文件 (Shape/data/NUTS_RG_03M_2010.prj) 显示如下:

GEOGCS["GCS_ETRS_1989",DATUM["D_ETRS_1989",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]

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

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