我一直在尝试将数据帧中的一列数值数据添加到 SpatialPolygon
数据帧,但每次尝试都会导致后者数据帧被转换为类似于前者的标准数据帧。我需要添加该列,以便我可以创建一个以该列的变量为焦点的等值线图。显然标准数据框不好,因为我正在尝试使用 tmap
创建 map 。
这就是我尝试添加列的方式(其中 shapefilecomb
是空间数据框,wardturnout
是包含相关列的变量):
shapefilecomb <- c(wardturnout)
最佳答案
通过赋值运算符将列添加到 SpatialPolygonsDataFrame 的数据槽中 shapefilecomb$wardturnout <- wardturnout
有效,但这不是完成这项工作的最安全方法。它仅依赖于位置(第一个数据项转到第一个多边形,第二个到第二个等等)。它可能会变得困惑。
最好保留给计算字段——shapefile$valuepercapita <- shapefile$value / shapefile$population
一种作业。
对于来自外部来源的数据,最好按键分配值。函数append_data
来自 tmap
package 做得非常好,不仅在发生错误时给你一条消息,而且在所有数据完美匹配时给你一个确认(我发现在处理大量不完美的数据时这是一个很好的接触)。
outShape <- append_data(srcShape, frmData, key.shp = "KOD_LAU1", key.data = "LAU1")
编辑(截至 9/2019):这个答案似乎仍然很强大……不过世界已经改变了。
tmap::append_data()
已移至tmaptools::append_data()
并且现在已经弃用了sf
已替换sp
作为 R 中空间数据的首选包
在
sf
世界空间数据存储在经过修改的 data.frames
中, 和
通过键分配数据项的最合适的方法是
*_join()
来自 dplyr
的功能- dplyr::left_join()
成为
安全起见,或dplyr::inner_join()
如果两侧过滤实际上是所需的行为。
关于r - 将数据框中的列添加到 SpatialPolygon 数据框中,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/43505467/