Rで空間データとポリゴンを扱うのは初めてです。
Google Earth から抽出した 2 つのポリゴンの 2 つの個別のシェイプ ファイルがあります。したがって、基本的に最初のシェイプ ファイルは場所 (ショッピング モールなど) であり、2 番目のシェイプ ファイルは最初の場所を中心とした半径 3 キロメートルです。両方のシェイプ ファイルを SpatialPolygonsDataFrames として R に読み込みます。次のコードを使用します。
library(maptools)
library(sp)
library(spatstat)
options(digits=10)
# Read polygon a
a <- readShapeSpatial(file.choose())
class(a)
spatstat.options(checkpolygons=FALSE)
r <- slot(a,"polygons")
r <- lapply(r, function(a) { SpatialPolygons(list(a)) })
windows <- lapply(r, as.owin)
Ploy_One <- tess(tiles=windows)
# Read polygon b
b <- readShapeSpatial(file.choose())
class(b)
spatstat.options(checkpolygons=FALSE)
s <- slot(b,"polygons")
s <- lapply(s, function(b) { SpatialPolygons(list(b)) })
windows <- lapply(s, as.owin)
Poly_Two <- tess(tiles=windows)
# Read polygon b
Combined_Region <- intersect.tess(Poly_One, Poly_Two)
plot(Combined_Region)
ただし、2 つのポリゴンを組み合わせたビュー (他のポリゴン内の 1 つのポリゴンのビュー) は得られません。
R で 2 つの多角形領域を 1 つの多角形領域にマージするコードを作成する方法について、誰かアドバイスがあれば、よろしくお願いします。