从shapefile中裁剪图

Question

我正在尝试使用欧洲统计局提供的shapefile和数据生成Choroplete地图.shapefile已在此处下载:使用此帖子中的 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 [ "GCS_ETRS_1989",DATUM [ "D_ETRS_1989",SPHEROID [ "GRS_1980",6378137.0,298.257222101],PRIMEM [ "格林威治",0.0],单位[ "度",0.0174532925199433]]

.prj文件中包含的信息指定Shapefile中几何数据的地理坐标系.尝试使用适当的坐标xlim和ylim该坐标系:

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