我需要一些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)
答案 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收入中用于标识邮政编码和收入字段的列名。