R中igraph网络的输出shapefile

Question

您好，我在 R 中有一个使用 igraph 库的网络

Vertices: 616 
Edges: 6270 
Directed: TRUE 
No graph attributes.
Vertex attributes: name, Lat, Lon.
Edge attributes: V3.

如何使用顶点中的纬度、经度信息为顶点和边生成两个 shapefile？

Answer 1

您可以使用sp 和maptools 包来执行此操作。 maptools 中有方便的函数 writePointsShape() 和 writeLinesShape() 将写入 ESRI shapefile 格式。

在执行此操作之前，需要从图形顶点中提取纬度/经度信息，并将其放入顶点的SpatialPoints 对象和边的SpatialLinesDataFrame 对象中。

此代码为以下示例生成一个非常简单的igraph 对象：

library(igraph)

## Produce a ring graph with 4 vertices
x <- graph.ring(4)

## Add lat/lon information to vertices
V(x)$lat <- c(50, 50, 51, 51)
V(x)$lon <- c(40, 41, 41, 40)

现在，为顶点创建 SpatialPoints 对象

library(sp)
library(maptools)

## Create SpatialPoints object containing coordinates
xV <- SpatialPoints(cbind(V(x)$lon, V(x)$lat))

## Write vertices to a shapefile
writePointsShape(xV, fn="vertices")

最后，为边缘创建SpatialLinesDataFrame 对象。这有点乱，但我还没有找到一种快速的方法来生成给定坐标的 SpatialLines 对象。

## Create SpatialLinesDataFrame object describing edges
edges <- get.edgelist(x)+1
edges <- cbind(edgeNum=1:nrow(edges), v1=edges[,1], v2=edges[,2])
xE <- apply(edges, 1, function(i) Lines(Line(cbind(c(V(x)$lon[i["v1"]], V(x)$lon[i["v2"]]), c(V(x)$lat[i["v1"]], V(x)$lat[i["v2"]]))), ID=as.character(i["edgeNum"])))
xE <- SpatialLinesDataFrame(SpatialLines(xE), data=data.frame(edgeNum=1:nrow(edges)))

## Write edges to a shapefile
writeLinesShape(xE, fn="edges")