如何在 R 中将 .shp 文件转换为 .csv？

Question

我想将两个 .shp 文件转换成一个数据库，这样我就可以一起绘制地图了。

另外，有没有办法将 .shp 文件转换为 .csv 文件？我希望能够在 .csv 格式下个性化和添加一些对我来说更容易的数据。如果在地图上添加叠加产量数据和降水数据，我有什么想法。

这里是 Morocco 和 Western Sahara 的 shapefile。

绘制两个文件的代码：

# This is code for mapping of CGE_Morocco results

# Loading administrative coordinates for Morocco maps
library(sp)
library(maptools)
library(mapdata)

# Loading shape files
Mor <- readShapeSpatial("F:/Purdue University/RA_Position/PhD_ResearchandDissert/PhD_Draft/Country-CGE/MAR_adm1.shp")
Sah <- readShapeSpatial("F:/Purdue University/RA_Position/PhD_ResearchandDissert/PhD_Draft/Country-CGE/ESH_adm1.shp")

# Ploting the maps (raw)
png("Morocco.png")
Morocco <- readShapePoly("F:/Purdue University/RA_Position/PhD_ResearchandDissert/PhD_Draft/Country-CGE/MAR_adm1.shp")
plot(Morocco)
dev.off()

png("WesternSahara.png")
WesternSahara <- readShapePoly("F:/Purdue University/RA_Position/PhD_ResearchandDissert/PhD_Draft/Country-CGE/ESH_adm1.shp")
plot(WesternSahara)
dev.off()

在研究了@AriBFriedman 和@PaulHiemstra 的建议并随后弄清楚如何合并 .shp 文件后，我设法使用以下代码和数据生成了以下地图（对于 .shp 数据，请参见上面的链接）

代码：

# Merging Mor and Sah .shp files into one .shp file

MoroccoData <- rbind(Mor@data,Sah@data) # First, 'stack' the attribute list rows using rbind() 
MoroccoPolys <- c(Mor@polygons,Sah@polygons) # Next, combine the two polygon lists into a single list using c()

summary(MoroccoData)
summary(MoroccoPolys)

offset <- length(MoroccoPolys) # Next, generate a new polygon ID for the new SpatialPolygonDataFrame object

browser()
for (i in 1: offset)
{
sNew =  as.character(i)
MoroccoPolys[[i]]@ID = sNew
}

ID <- c(as.character(1:length(MoroccoPolys))) # Create an identical ID field and append it to the merged Data component
MoroccoDataWithID <- cbind(ID,MoroccoData)

MoroccoPolysSP <- SpatialPolygons(MoroccoPolys,proj4string=CRS(proj4string(Sah))) #  Promote the merged list to a SpatialPolygons data object

Morocco <- SpatialPolygonsDataFrame(MoroccoPolysSP,data = MoroccoDataWithID,match.ID = FALSE) #  Combine the merged Data and Polygon components into a new SpatialPolygonsDataFrame.

Morocco@data$id <- rownames(Morocco@data)
Morocco.fort <- fortify(Morocco, region='id') 
Morocco.fort <- Morocco.fort[order(Morocco.fort$order), ] 

MoroccoMap <- ggplot(data=Morocco.fort, aes(long, lat, group=group)) + 
geom_polygon(colour='black',fill='white') + 
theme_bw()

结果：

新问题：

1- 如何消除将地图一分为二的边界数据？

2- 如何在 .shp 文件中组合不同的区域？

谢谢大家。

P.S：stackoverflow.com 中的社区非常棒，非常有帮助，尤其是对像 :) 这样的初学者来说，只是想强调一下。

Answer 1

将 shapefile 加载到 Spatial{Lines/Polygons}DataFrame（来自 sp-package 的类）后，您可以使用 fortify 通用函数将它们转换为平面 data.frame 格式。 fortify 泛型的特定函数包含在 ggplot2 包中，因此您需要先加载它。代码示例：

library(ggplot2)
polygon_dataframe = fortify(polygon_spdf)

其中polygon_spdf 是SpatialPolygonsDataFrame。类似的方法适用于SpatialLinesDataFrame's。

我的解决方案与@AriBFriedman 的解决方案之间的区别在于，我的解决方案包括多边形/线的x 和y 坐标，以及与这些多边形/线关联的数据。我真的很喜欢使用 ggplot2 包来可视化我的空间数据。

将数据保存在普通的 data.frame 中后，您可以简单地使用 write.csv 在磁盘上生成 csv 文件。

Answer 2

我想你的意思是你想要每个相关的 data.frame？

如果是这样，可以使用@插槽访问功能进行访问。该插槽称为data：

write.csv( WesternSahara@data, file="/home/wherever/myWesternSahara.csv")

然后当你用read.csv读回它时，你可以尝试分配：

myEdits <- read.csv("/home/wherever/myWesternSahara_modified.csv")
WesternSahara@data <- myEdits

您可能需要对行名等进行一些按摩，以使其接受新的 data.frame 为有效。我可能会尝试将现有的 data.frame 与您在 R 中读取的 csv 合并，而不是破坏性地进行编辑....