将 R 绘图导出到 shapefile

Question

我对 R 相当陌生，但对 ArcView 不熟悉。我正在绘制一些双模式数据，并希望将绘图转换为 shapefile。具体来说，如果可能的话，我想转换顶点和边，以便在 ArcView 中显示相同的图以及属性。

我已经安装了“shapefiles”包，并且看到了 convert.to.shapefile 命令，但是帮助中没有提到如何将 XY 坐标分配给顶点。

谢谢，

提姆

Answer 1

好的，我在这里做了几个假设，但我阅读了这个问题，因为您希望将空间坐标分配给二分图并将顶点和边导出为点 shapefile 和折线以在 ArcGIS 中使用。

这个解决方案有点笨拙，但会生成坐标限制 xmin、ymin 和 xmax、ymax 为 -0.5 和 +0.5 的 shapefile。由您决定图形布局算法（例如 Kamada-Kawai），并根据 @gsk3 的建议将 shapefile 投影到所需的坐标系中。可以在创建 points.data 和 edge.data 数据框的位置添加顶点和边的附加属性。

library(igraph)
library(shapefiles)

# Create dummy incidence matrix
inc <- matrix(sample(0:1, 15, repl=TRUE), 3, 5)
colnames(inc) <- c(1:5) # Person ID
rownames(inc) <- letters[1:3] # Event 

# Create bipartite graph
g.bipartite <- graph.incidence(inc, mode="in", add.names=TRUE)

# Plot figure to get xy coordinates for vertices
tk <- tkplot(g.bipartite, canvas.width=500, canvas.height=500)
tkcoords <- tkplot.getcoords(1, norm=TRUE) # Get coordinates of nodes centered on 0 with +/-0.5 for max and min values

# Create point shapefile for nodes
n.points <- nrow(tkcoords)
points.attr <- data.frame(Id=1:n.points, X=tkcoords[,1], Y=tkcoords[,2])
points.data <- data.frame(Id=points.attr$Id, Name=paste("Vertex", 1:n.points, sep=""))
points.shp <- convert.to.shapefile(points.attr, points.data, "Id", 1)
write.shapefile(points.shp, "~/Desktop/points", arcgis=TRUE)


# Create polylines for edges in this example from incidence matrix
n.edges <- sum(inc) # number of edges based on incidence matrix
Id <- rep(1:n.edges,each=2) # Generate Id number for edges.
From.nodes <- g.bipartite[[4]]+1 # Get position of "From" vertices in incidence matrix
To.nodes <- g.bipartite[[3]]-max(From.nodes)+1 # Get position of "To" vertices in incidence matrix

# Generate index where position alternates between "From.node" to "To.node"
node.index <- matrix(t(matrix(c(From.nodes, To.nodes), ncol=2))) 

edge.attr <- data.frame(Id, X=tkcoords[node.index, 1], Y=tkcoords[node.index, 2])
edge.data <- data.frame(Id=1:n.edges, Name=paste("Edge", 1:n.edges, sep=""))
edge.shp <- convert.to.shapefile(edge.attr, edge.data, "Id", 3)
write.shapefile(edge.shp, "~/Desktop/edges", arcgis=TRUE)

希望这会有所帮助。

Answer 2

我将根据您的数据是什么样子的疯狂猜测来尝试一下。

基本上，您需要将数据强制转换为 data.frame，其中两列包含 x 和 y 坐标（或纬度/经度等）。

library(sp)
data(meuse.grid)
class(meuse.grid)
coordinates(meuse.grid) <- ~x+y
class(meuse.grid)

一旦你将它作为 SpatialPointsDataFrame，sp 提供了一些不错的功能，包括导出 shapefile：

writePointsShape(meuse.grid,"/home/myfiles/wherever/myshape.shp")

相关帮助文件示例来自：

• 坐标
• SpatialPointsDataFrame
• 读取ShapePoints

至少在几年前，当我最后一次使用 sp 时，投影非常好，而将投影信息写入 shapefile 非常糟糕。所以最好保持坐标不变，并手动告诉 Arc 它是什么投影。或者使用writeOGR 而不是writePointsShape。