Выходной шейп-файл для сети igraph в R

Вы можете сделать это с помощью пакетов sp и maptools. В maptools есть удобные функции writePointsShape() и writeLinesShape(), которые будут записывать в формат шейп-файла ESRI.

Прежде чем сделать это, необходимо извлечь информацию о широте и долготе из вершин графа и поместить ее в объект 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")