spatstatの頂点からネットワークを構築する際のエラー

Question

spatstatで線形ネットワークとして読み込むための非常に大きな（シェイプファイル）道路ネットワークがあります。だから私は、この本の第17章で説明したように頂点と辺を読んでから、基本的なネットワークを構築しようとしています - 空間点パターンをバデリーらによって

私はエラーを取得するの下に、私はこのコードを使用して私のデータhere

を添付Error: length(x0) == length(x1) is not TRUE。問題を見つけるためにx0とx1が何であるかはわかりません。

library(maptools) 
library(spatstat) 
setwd("~/documents/rwork/traced/a") 

pt <- readShapePoints("collected.shp") #read vertices from a shapefile. 

edgeRecords<-read.delim("edgelist.txt") #read edge connectivity list 

ed<-data.frame(from=edgeRecords$from,to=edgeRecords$to) 

xx<-pt@bbox[1,]#read x bounds of owin 
yy<-pt@bbox[2,]#read y bounds of owin 

v<-ppp(x=pt@coords[,1], y=pt@coords[,2], xx,yy) #read list of vertices 


edg<-as.matrix(ed) # read node pairs as matrix 

built_network<-linnet(v,edges = edg) 

これは、上記のコメントの1つに記載のエラー

Error: length(x0) == length(x1) is not TRUE 

Answer 1

もたらします。 Rインデックスが1から始まる間にGISインデックスが0から始まることに気付きました。

この問題を解決するために、エッジマトリックスに+1を追加しました。 GISソフトウェアからエッジ行列を収集した場合、from_nodeまたはto_nodeのいずれかでノード0への参照があるためです。 Rのエッジ行列がemの場合、+1を追加します。つまり、em+1です。サンプルコードはこのようになります

edgelist <- read.delim("edgelist.txt") 
em <- matrix(c(edgelist$from, edgelist$to), ncol=2) +1 
net <- linnet(n,edges = em) 
plot(net) 

これは私にとって問題を解決しました。誰かを助けることを願っています。または、他の解決策がある場合は、自由に共有してください。