3

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 つの多角形領域にマージするコードを作成する方法について、誰かアドバイスがあれば、よろしくお願いします。

4

1 に答える 1