R中的栅格合并

时间:2015-02-19 16:06:03

标签: r

我需要一些R语法的帮助来完成(我认为)是一项相当简单的任务 - 希望有人可以提供帮助!

我有一张英国的光栅地图,它被分成邮政编码区域(例如DE,NG,NR等共127个邮政编码)。

我已经安装了包' raster'并已成功将.img绘制在R中。所有工作并使用栅格看起来正确。

我还有一个逗号分隔的CSV文件,其中包含与栅格相同的邮政编码,旁边有另一列,其中包含每个邮政编码的收入。

我想知道是否有人可以帮助我将收入数据合并/绑定到栅格中的正确邮政编码,以便我可以绘制每个邮政编码的收入。

我觉得我应该使用cbind并重新分类来执行此操作,但我无法自行完成

任何帮助将不胜感激。提前谢谢!

这是我到目前为止的代码......还不是火箭科学。

setwd("C:\\Users\\[username]\\Documents\\GIS\\Test Data")
require(raster)
revenue<-read.table("revenue.csv",header=T,row.names=1,sep=",")
postcodes<-raster("C:\\Users\\[username]\\Documents\\GIS\\Test Data\\rasters\\postcodes\\postcodes.img")
trim(postcodes)
plot(postcodes)

1 个答案:

答案 0 :(得分:1)

您应该可以使用'subs'方法执行此操作。您没有向我们展示您的数据(例如头(收入)),但它应该像这样:

library(raster)
setwd("C:\\Users\\[username]\\Documents\\GIS\\Test Data")
postcodes <- raster("rasters\\postcodes\\postcodes.img")
revenue <- read.csv("revenue.csv")

subs(postcodes, revenue, by='code', which='rev')

其中'code'和'rev'是data.frame收入中用于标识邮政编码和收入字段的列名。