如何在R中将csv转换为shp

Question

在过去的几天里,我一直在尝试将csv转换为shapefile.我知道我可以在QGIS或Arc中轻松完成,但是想将此过程添加到我现有的R代码中.

所以我可以在csv中读到没有任何问题

MyData <- read.csv(file="c:/TheDataIWantToReadIn.csv", header=TRUE, sep=",")

我在Packages Shapefile帮助指南中找到了以下代码.但是,我似乎无法找到一种方法来处理我的代码.我的行都是一个点,因此我想创建的shapefile将是所有点.我没有Id列,但是我在两个单独的列中有x和y数据.

dd <- data.frame(Id=c(1,2),X=c(3,5),Y=c(9,6))
ddTable <- data.frame(Id=c(1,2),Name=c("Item1","Item2"))
ddShapefile <- convert.to.shapefile(dd, ddTable, "Id", 1)
write.shapefile(ddShapefile, "c:/test", arcgis=T)

任何帮助将不胜感激.

Answer 1

我建议使用rgdal而不是shapefiles.要使用rgdal,您必须从http://cran.revolutionanalytics.com/web/packages/rgdal/查看系统要求.

以下代码可以帮助您找到正确的方向:

install.packages(c("rgdal", "sp"))
library(rgdal)
library(sp)
MyData <- read.csv(file="c:/TheDataIWantToReadIn.csv", header=TRUE, sep=",")

以下代码片段来自使用ggplot2包的R中的映射.

class(MyData) # data.frame
coordinates(MyData)<-~longitude+latitude # whatever the equivalent is in your 
# data.frame
class(MyData) # [1] "SpatialPointsDataFrame"
          # attr(,"package")
          # [1] "sp"

以下代码片段如何编写带有投影的shapefile - 问题已解决

writeOGR(crest.sp, "c:/test", "layer name", driver = "ESRI Shapefile")