根据R中的属性表值导出栅格

我有一个大型栅格文件，它具有关联的栅格属性表 (RAT)：

class       : RasterLayer 
dimensions  : 32550, 13900, 452445000  (nrow, ncol, ncell)
resolution  : 1, 1  (x, y)
extent      : 0, 13900, 0, 32550  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
data source : H:\My Drive\raster.grd 
names       : mr200_2012 
values      : -318, 151078  (min, max)
attributes  :

ID COUNT  Value_of_interest 
1   4        45
2   7        1
3   94       5
4   5        7
...

栅格显示在从 -318 到 151078 的值上。我希望能够导出以属性值之一作为“值”的栅格。 例如，将“value”更改为“Value_of_interest”，以便栅格显示 1 到 45 的（最小值、最大值）。

这与如何Lookup http://desktop.arcgis.com/en/arcmap/latest/tools/spatial-analyst-toolbox/lookup.htm工具在 ArcMap 中工作。我只是在寻找一种编程方式来做到这一点。

我没有太多使用 RAT，并且不熟悉它的一些功能。我尝试过使用data.table and data.frame方法，但生成的表确实很大，并且大多数时候由于大小而不会生成。我的想法是我可以创建一个data.frame with xy坐标，然后更改z值“Value_of_interest”。

我认为这将是足够的可重复数据：

xy <- matrix(rnorm(200),20,20)
xy<-raster(xy)
xy
rat<-data.frame("ID" = runif(10,0,10),
            "b" = runif(10,0,10),
            "variable of interest" = runif(10,10,20))


levels(xy)[[1]]<-rat

我找到了解决这个问题的方法。使用deratify使用该函数，您可以分离出某些属性并创建新的单波段栅格。

New.raster <- deratify(Raster.with.attributes, "Some_Column_in_raster")

新栅格将是具有指定列值的单波段。