3

Eurostatが提供するシェープファイルとデータを使用してコロプレテマップを作成しようとしています。シェープファイルはここからダウンロードされています：この投稿のJDLongsコードを使用しています。

これは、以下に投稿した図を再現するための最小限のコードです。

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（）を使用してグラフィックデバイスから座標を取得しましたが、値をプラグインしても目的の結果が得られませんでした。

最小限の例

4

1 に答える 1

6

シェープファイルプロジェクションファイル（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ファイルに含まれる情報は、シェープファイル内の幾何学的データの地理座標系を指定します。この座標系に適切な座標xlimを試してください。ylim

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

プロット例

于 2012-10-10T14:09:28.587 に答える