ggplot2を使用して複数のサブポリゴンと穴を持つ空間ポリゴンを描画する

Question

ライブラリspにggplot2の穴があるSpatialPolygonsを描画します。
http://stackoverflow.com/questions/12047643/geom-polygon-with-multiple-hole/12051278#12051278
実際 broom::tidyを使用してSpatialPolygonsを変換する際に、（ggplot2::fortifyを置き換える）、穴のポリゴンが時計回りに保存されているとして描画する：StackOverflowの上の他の人の質問に おかげで、私は時計回りに書かれたポリゴンを扱う一方で、これが許可されていることを知っています穴。
ggplot2では、穴があるポリゴンを強制的に描画して、fillを使用して一度描画し、もう一度colourを使用して描画します。そうしないと、ポリゴンと交差する線が表示されることがあります。 複数のサブポリゴンを扱う場合、穴があるものもありますが、これはより難しく、broom::tidyで定義されている点のフィーチャの順番でポリゴンを塗りつぶすことはできません（下の画像を参照）。
あなたはこの充填問題の問題を解決するための解決策をお持ちですか？ここで

は再現例です。

library(sp) 
library(ggplot2) 

# Create two polygons: second would be a hole inside the first 
xy = cbind(
    x = c(13.4, 13.4, 13.6, 13.6, 13.4), 
    y = c(48.9, 49, 49, 48.9, 48.9) 
    ) 
hole.xy <- cbind(
    x = c(13.5, 13.5, 13.45, 13.45, 13.5), 
    y = c(48.98, 48.92, 48.92, 48.98, 48.98) 
) 

# Transform as SpatialPolygons with holes 
xy.sp <- SpatialPolygons(list(
    Polygons(list(Polygon(xy), 
       Polygon(hole.xy, hole = TRUE)), "1"), 
    Polygons(list(Polygon(xy + 0.2), 
       Polygon(xy + 0.35), 
       Polygon(hole.xy + 0.2, hole = TRUE)), "2") 
)) 

# Transform SpatialObject to be used by ggplot2 
xy.sp.l <- broom::tidy(xy.sp) 

ggplot(xy.sp.l) + 
    geom_polygon(aes(x = long, y = lat, group = id, fill = id)) 

Answer 1

はsfパッケージに "オーバー行く" ための良い機会でした。 sfオブジェクトでの作業ggplotで実際にははるかに簡単です、geom_sfジオメトリのおかげ：行を追加

library("sf") 
library("rgeos") 
sf_poly <- as(xy.sp, "sf") 
sf::st_crs(sf_poly) <- 4326 
sf_poly$id <- c(1,2) 
ggplot(sf_poly) + 
    geom_sf(aes(fill = as.factor(id))) 

Answer 2

は、問題の原因を示しています。青色の "ポリゴン"が描画され、下から上の穴が描画されます。 （非常に、エレガントに残念ではありません）

このコードは、パスが3位に進む前に戻って最初の作品の出発点に戻ります。

library(dplyr) 
    extra <- xy.sp.l %>% 
     filter(piece != 1) %>% 
     group_by(id, group) %>% 
     summarise(last_pt = max(order)) 


for (n in 1:nrow(extra)) { 
    id_ex <- as.character(extra[n,"id"]) 
    x <- subset(xy.sp.l, id == id_ex & piece == 1 & order == 1) 
    x$order <- as.numeric(extra[n,"last_pt"]) + 0.5 
    xy.sp.l <- rbind(xy.sp.l,x) 
} 

xy.sp.l <- xy.sp.l[order(xy.sp.l$id, xy.sp.l$order),] 

Answer 3

この投稿は良い質問です、すでに偉大な答えを受けました。私は人々がsfオブジェクトを扱う方法を学ぶべきだとも信じています。それはRの空間データ型の次世代であるためです。しかし、この場合、ggspatialのgeom_spatialパッケージをSpatialPolygonsをプロットするオプションにすることができます。

library(sp) 
library(ggplot2) 
library(ggspatial) 

ggplot() + 
    geom_spatial(xy.sp, aes(fill = id)) 
# Ignoring argument 'mapping' in geom_spatial.SpatialPolygons 
# Autodetect projection: assuming lat/lon (epsg 4326)