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

Question

我想将两个.shp文件转换为一个允许我一起绘制地图的数据库。

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

以下是Morocco和Western Sahara的形状文件。

绘制两个文件的代码：

# 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数据，参见上面的链接）< / p>

代码：

# 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} DataFrames（sp-package中的类）后，可以使用fortify泛型函数将它们转换为flat data.frame格式。 fortify通用的具体功能包含在ggplot2包中，因此您需要先加载它。代码示例：

library(ggplot2)
polygon_dataframe = fortify(polygon_spdf)

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

我的解决方案与@AriBFriedman的解决方案之间的区别在于，除了与这些polgons /行关联的数据之外，我还包括多边形/线的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合并，而不是破坏性地编辑....