用空间线连接 xy 点

Question

我有一个如下所示的 XY 点数据集

x<-c(2,4,6,3,7,9,1)
y<-c(6,4,8,2,9,6,1)
id<-c("a","b","c","d","e","f","g")
dataset<-data.frame(cbind(x,y,id))

我想用空间线连接所有点的所有组合，用它们连接的点的组合命名的线

在输出结果的“属性表”中，空间线的名称可能如下所示：

a_b
a_c
a_d
a_e
a_f
a_g
b_a
b_c
b_d
b_e
b_f
b_g
c_a
etc.

Answer 1

我在此推测您到底想要什么，但我认为您希望可视化从任何点到其他点的连接。如果是这样的话，那么这可能会奏效。 但首先，一些假设：

• 您的x 和y 坐标是起点。因此，id 是 id.origin
• 所有其他点都需要成为“目的地”，然后它们自己的坐标将成为x_destination等等。

应该有更好、更优雅的方式来做到这一点。如果有经验丰富的人可以加入并告诉我任何*ply 方法，我将不胜感激。 免责声明>

复制数据框以涵盖所有可能的组合

dataset<-do.call(rbind, replicate(7, dataset, simplify=FALSE))

现在，创建一个包含所有相同目标点的矩阵，混合：

nm=matrix(ncol=3)
for (i in 1:7){
    nm<-rbind(nm,do.call(rbind,replicate(7,as.matrix(dataset[i,]),simplify=FALSE)))
}
nm<-nm[-1,]

重命名矩阵的列，使其有意义，并将现有数据框与新矩阵绑定

colnames(nm)<-c("x2","y2","id.dest")
newds<-cbind(dataset,as.data.frame(nm))

删除重复的轨迹：

newds<-newds[-which(newds$id.origin==newds$id.dest),]

并使用geom_segment绘制结果：

p<-ggplot(newds,aes(x=x,y=y))+geom_segment(aes(xend=x2,yend=y2))

有一种方法可以命名片段，但从观察情节来看，我不建议这样做。相反，您可以考虑使用geom_text 命名这些点（其他选项可用，请参阅?annotate 之一）。

p<-p + geom_text(aes(x=1.8,y=6.1,label="a"))

这将产生一个像这里这样的情节：

Answer 2

整个解决方案如下所示：

plot(dataset$x,dataset$y)

复制数据框以涵盖所有可能的组合

dataset<-do.call(rbind, replicate(7, dataset, simplify=FALSE))

现在，创建一个包含所有相同目标点的矩阵，混合：

nm=matrix(ncol=3)
for (i in 1:7){
  nm<-rbind(nm,do.call(rbind,replicate(7,as.matrix(dataset[i,]),simplify=FALSE)))
}
nm<-nm[-1,]

重命名矩阵的列，使其有意义，并将现有数据框与新矩阵绑定

colnames(nm)<-c("x2","y2","id.dest")
newds<-cbind(dataset,as.data.frame(nm))

删除重复的轨迹：

newds1<-newds[-which(newds$id==newds$id.dest),]

library(ggplot2)

将目标 x 和 y 从因子转换为数值

newds1$x2<-as.numeric(as.character(newds1$x2)) #converting from factor to numeric
newds1$y2<-as.numeric(as.character(newds1$y2))

绘制目标点。 . .与原点相同

plot(newds1$x, newds1$y) 
plot(newds1$x2, newds1$y2, col="red") 

现在使用此答案中的代码：

Convert Begin and End Coordinates into Spatial Lines in R

存储 Lines 对象的原始列表：

l <- vector("list", nrow(newds1)) #

这个l 现在是一个空向量，其中包含由 newds1 的长度 (nrow) 定义的行数

分割起点和终点坐标以便我可以运行这个脚本：

origins<-data.frame(cbind(newds1$x, newds1$y))
destinations<-data.frame(cbind(newds1$x2, newds1$y2))

library(sp)
for (i in seq_along(l)) {
  l[[i]] <- Lines(list(Line(rbind(origins[i, ], destinations[i,]))), as.character(i))
}

l.spatial<-SpatialLines(l)
plot(l.spatial, add=T)