3

ウィンドウ（spatstat R）に（三角形、長方形などの）一種のテッセレーションを描画するランダムなセグメントのセットがあります。いくつかのインデックス（面積、形状インデックスなど）を計算するには、それを一連のポリゴン（SpatialPolygons）に変換する必要があります。

これは明らかに簡単ですが、私はそれを行う方法を見つけることができませんでした...

以下は、自己インターセプト セグメントのランダム パターンを生成するCarl Witthoftのコードです。

ranpoly <- function(numsegs=10,plotit=TRUE) {

require(spatstat)
# temp fix: put the first seg into segset. Later make it a constrained random.
segset<-psp(c(0,1,1,0,.25),c(0,0,1,1,0),c(1,1,0,0,1),c(0,1,1,0,.75),owin(c(0,1),c(0,1)) ) #frame the frame
for (jj in 1: numsegs) {
# randomly select a segment to start from, a point on the seg, the slope,and direction
# later... watch for slopes that immediately exit the frame
endx <-sample(c(-0.2,1.2),1)  #force 'x1' outside the frame
# watch that sample() gotcha
if(segset$n<=5) sampset <- c(5,5) else sampset<-5:segset$n
startseg<-sample(sampset,1) #don't select a frame segment
# this is slope of segment to be constructed
slope <- tan(runif(1)*2*pi-pi) # range +/- Inf 
# get length of selected segment
seglen<-lengths.psp(segset)[startseg]
startcut <- runif(1) 
# grab the coords of starting point (similar triangles)
startx<- segset$ends$x0[startseg] + (segset$ends$x1[startseg]-segset$ends$x0[startseg])*startcut #seglen
starty<- segset$ends$y0[startseg] + (segset$ends$y1[startseg]-segset$ends$y0[startseg])*startcut #seglen
# make a psp object with that startpoint and slope; will adjust it after finding intersections
endy <- starty + slope*(endx-startx)
newpsp<-psp(startx,starty,endx,endy,segset$window,check=FALSE)
# don't calc crossing for current element of segset
hits <- crossing.psp(segset[-startseg],newpsp)
segdist <- dist(cbind(c(startx,hits$x),c(starty,hits$y)))
# dig back to get the crosspoint desired -- have to get matrixlike object out of class "dist" object
# And, as.matrix puts a zero in location 1,1 kill that row.
cutx <- hits$x[ which.min( as.matrix(segdist)[-1,1] )]
cuty <- hits$y[which.min(as.matrix(segdist)[-1,1] )]
segset <- superimpose(segset,psp(startx,starty,cutx,cuty,segset$window))

} #end jj loop
if(plotit) plot(segset,col=rainbow(numsegs))
return(invisible(segset))
}

segset=ranpoly()

segset は、SpatialPolygons オブジェクトを作成するために必要な psp オブジェクトです。

4

2 に答える 2

1

spatstatオブジェクトのセットがあると仮定すると、次のようなものを試すことができます（テストされていない）：

require(sp)

# VECTOR OF spatstat OBJECT NAMES
segs <- (seg1,seg2,seg3)

segPolys <- as(segs[1], "SpatialPolygons")
  for( i in 2:length(segs)) {
    y <- as(segs[i], "SpatialPolygons")
      slot(y[[i]], "ID") <- paste(i)
    segPolys <- c(slot(y, "polygons"),segPolys)
   }
于 2013-01-28T17:52:18.383 に答える
1

グーグルで検索するspatstat as spatialPolygonsと、この最初のヒットにたどり着きました。これは、シェープファイルの処理専用のspatstatのビネットです。spクラスをspatstatオブジェクトに変換する方法に多くの時間を費やしています。セクション 3.2.5:クラス SpatialPolygons のオブジェクトとセクション 3.2.6:クラス SpatialPolygonsDataFrame のオブジェクトに最も関心があるかもしれません。

于 2013-01-28T13:05:17.087 に答える