4

一部だけをプロットしたいシェープファイル (ポリゴン) があるので、ggplot の geom_path を使用して、表示されるパスを xlim と ylim で制限します。シェープファイルはこちら (2.5 MB) から入手できます。

library(rgdal)
library(ggplot2)


coast <- readOGR(dsn=".", layer="coast-rgdal")
coast.df <- fortify(coast, region="STATE_NAME")
p <- ggplot(coast.df, aes(x=long, y=lat, group=group)) +
  geom_path(colour="black", size=1/4) +
  xlim(146, 148) + ylim(-39.25, -37.5) + coord_fixed()
p

次のようになります。

ビックコーストオープン

私がしたいのは、閉じたパスを使用した次のようなものです(gimpで実行):

ビックコースト閉鎖

ggplot コマンドでそれを行うことはできますか? それともRの何か？

参考までに、完全なポリゴンは次のようになります。

ここに画像の説明を入力

4

1 に答える 1

3

パッケージから使用gIntersectionして、rgeosこの領域を抽出できます。

library(rgdal)
library(ggplot2)
library(rgeos)

coast <- readOGR(dsn=".", layer="coast-rgdal")

lim <- cbind(c(146, 148, 148, 146, 146),
             c(-39.25, -39.25, -37.5, -37.5, -39.25))
x <- SpatialPolygons(list(Polygons(list(Polygon(lim)), ID="1")))
proj4string(x) <- proj4string(coast)
res <- gIntersection(coast, x)

coast.df <- fortify(res)
p <- ggplot(coast.df, aes(x=long, y=lat, group=group)) +
     geom_path(colour="black", size=1/4) +
     xlim(146, 148) + ylim(-39.25, -37.5) + coord_fixed()
p

R プロット

于 2013-09-24T09:54:13.527 に答える