シェイプファイルvignetteをspatstatで使用すると、シェイプファイルの各フィーチャが属性としてこのpspリストの要素になるようにpspオブジェクトのリストを作成できました。具体的には、pspオブジェクトのリストをlinnetに変換する方法
> x<- readShapeSpatial("includes/roadshapefile.shp")
> out <- lapply([email protected], function(z) { lapply([email protected], as.psp) })
> dat <- [email protected]
> for(i in seq(nrow(dat)))
+ out[[i]] <- lapply(out[[i]], "marks<-", value=dat[i, , drop=FALSE])
> roadlist <- do.call("c", out)
を次のように私はビネットにセクション3.2.4を使用今私はPSPのオブジェクトのリストから、ネットワーク全体linnet
を作成したいと思います。直接関数がありますか、またはそれぞれのpspオブジェクトをループして、それぞれをlinnetにして、何とかそれを完全なネットワークに結合する必要がありますか? as.linnet
は、そのクラスがプレーン・リストであるため、pspのこのリストを有効な入力として受け取ることはありません。
私はいくつかの小さなステップが欠けていると感じています。誰かが私を助けてくれますか。
私はlinnetを取得した後、linfun
を使用してポイントパターンを作成し、これらの属性(マーク)にアクセスして非均質なポイントを作成したいと考えています。
おかげ
最近Adrian BaddeleyがRoger Bivandに 'linnet'コードを送りました。このコードは、R-forgeで利用できる' maptools'の開発版に入れました。http://r-forge.r-project.org/R/?group_id = 943おそらく、これを使用してより簡単で素早くインポートされたlinnetを得ることができます。どのように属性を保存して、それらを 'linfun'や' linim'に入れることができるのかよく分かりません。私はそれを調べる必要があります。 –
ありがとうございます。私はあなたの参照を読んでいます。あなたはおそらく、私がspatstatを使ってネットワーク上で非同質の強度をどのように生成するか教えてください。シェイプファイルから少なくとも2つの属性を取り込みたい。列が破棄される前に、属性を1つだけ使用できますか? – BKS
ありがとうございます。 SpatialLinesオブジェクトからlinnetを作成する新しい関数が見つかりましたが、あなたが言ったように、私は属性を失った:( – BKS