2

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

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

这是摩洛哥和西撒哈拉的 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 文件中组合不同的区域？

谢谢大家。

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

4

2 回答 2

8

将 shapefile 加载到 Spatial{Lines/Polygons}DataFrames（sp 包中的类）后，您可以使用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 文件。

于 2013-04-28T18:29:12.470 回答
6

我想你的意思是你想要每个相关的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 合并，而不是破坏性地进行编辑......

于 2013-04-28T18:25:14.167 回答