を使用してポリゴン間の小さな「断片」を排除する方法はありR
ますか?望ましいソリューションはSpatialPolygonsDataFrame
、ポリゴン間の共有境界が一致する新しいものを作成します。R
ArcMapやQGISではなく、を使用するソリューションに特に興味があります。
そもそもなぜこれらのギャップが存在するのかについての説明も聞きたいです。
これは、私が使用している空間データの再現可能な例です。
library(rgdal)
library(sp)
library(tigris)
library(magrittr)
library(leaflet)
library(gplots)
# This project will use WGS 84 projected coordinate system
crs_proj <- CRS("+init=epsg:4326")
# These are the FIPS codes of the specific block groups in my study area
sel <- c("530330079005", "530330079001", "530330079004",
"530330085002", "530330085003", "530330086003",
"530330087003", "530330085001", "530330090001",
"530330091001", "530330091002", "530330092001",
"530330092002", "530330086001", "530330090002",
"530330086002", "530330079003", "530330079002",
"530330087002", "530330087001")
# Create polygons
polygons <- tigris::block_groups(state = "WA",county = "King") %>%
.[.@data$GEOID %in% sel,] %>%
spTransform(CRSobj = crs_proj)
# Map the result
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = polygons,
stroke = F,
fillColor = col2hex("red"), fillOpacity = 1)
上のスクリーンショットからわかるように、国勢調査ブロックグループのポリゴン間に小さなギャップがあります。これらのギャップのこれらの位置は、ズームレベルに応じて異なる方法でレンダリングされますが、常にいくつかのギャップが表示されます。
誰かR
がプログラムでポリゴンを調整してこれらのギャップをなくすための関数(または関数の組み合わせ)を推奨できますか?