10

我正在尝试在 R 中创建一个 shapefile，稍后将其导入到 Fusion Table 或其他一些 GIS 应用程序中。

首先，我导入了一个包含加拿大所有人口普查区的空白 shapefile。我已根据 CT 的唯一 ID 将其他数据（以表格格式）附加到 shapefile，并映射了我的结果。目前，我只需要温哥华的那些，我想导出一个仅包含温哥华 CT 以及我新附加的属性数据的 shapefile。

这是我的代码（由于隐私原因省略了一些部分）：

shape <- readShapePoly('C:/TEST/blank_ct.shp') #Load blank shapefile
shape@data = data.frame(shape@data, data2[match(shape@data$CTUID, data2$CTUID),]) #data2 is my created attributes that I'm attaching to blank file
shape1 <-shape[shape$CMAUID == 933,] #selecting the Vancouver CTs

我见过其他使用这个的例子：writePolyShape来创建 shapefile。我试过了，它在一定程度上奏效了。它创建了 .shp、.dbf 和 .shx 文件。我缺少 .prj 文件，我不确定如何创建它。有没有更好的方法来创建 shapefile？

对此问题的任何帮助将不胜感激。

4

2 回答 2

8

使用rgdal和writeOGR。rgdal将保留投影信息

就像是

library(rdgal)

shape <- readOGR(dsn = 'C:/TEST', layer = 'blank_ct')
# do your processing
shape@data = data.frame(shape@data, data2[match(shape@data$CTUID, data2$CTUID),]) #data2 is my      created attributes that I'm attaching to blank file
shape1 <-shape[shape$CMAUID == 933,]
writeOGR(shape1, dsn = 'C:/TEST', layer ='newstuff', driver = 'ESRI Shapefile')

请注意，这dsn是包含文件的.shp文件夹，并且是不带扩展layer名的 shapefile 的名称。.shp它将读取 ( readOGR) 并写入 ( writeOGR) 所有组件文件 ( .dbf,.shp等.prj)

于 2013-04-12T03:16:41.077 回答
4

问题解决了！再次感谢帮助过的人！

这是我最终做的事情：

正如 Mnel 所写，这一行将创建 shapefile。

writeOGR(shape1, dsn = 'C:/TEST', layer ='newstuff', driver = 'ESRI Shapefile')

但是，当我运行此行时，它返回此错误：

Can't convert columns of class: AsIs; column names: ct2,mprop,mlot,mliv

这是因为我的属性数据不是数字，而是字符。幸运的是，我的属性数据都是数字，所以我运行了 transform() 来解决这个问题。

shape2 <-shape1
shape2@data <- transform(shape1@data, ct2 = as.numeric(ct2),
mprop = as.numeric(mprop),
mlot = as.numeric(mlot),
mliv = as.numeric(mliv))

我再次尝试了 writeOGR() 命令，但仍然没有得到我正在寻找的 .prj 文件。问题是我在导入文件时没有指定 shapefile 的坐标系。由于我已经知道坐标系是什么，我所要做的就是在导入时定义它。

readShapePoly('C:/TEST/blank_ct.shp',proj4string=CRS("+proj=longlat +datum=WGS84")

之后，我重新运行了我想对 shapefile 和 writeOGR 行进行导出的所有操作。就是这样！

于 2013-04-12T14:54:04.517 回答