r - 如何在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=",")

以下代码片段来自 Mapping in R using the ggplot2 package 。

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

以下代码片段来自 How to write a shapefile with projection - problem solved

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