如何在R中访问TIF地图的属性表？

Question

我带有以下文件的TIF格式的downloaded a map：a2_LOt_12.tif.vat.cpg; a2_LOt_12_v170321.tfw; a2_LOt_12_v170321.tif; a2_LOt_12_v170321.tif.aux; a2_LOt_12_v170321.tif.ovr a2_LOt_12_v170321.tif.vat.dbf; a2_LOt_12_v170321.tif（XML doc）。在ArcGIS中，文件打开正常，我可以访问属性表：

但是，我很难在R中打开该表。我希望“levels（）”可以工作，但它返回“NULL”：

r <- raster("a2_LOt_12_v170321.tif")
levels(r)

是否有另一种方法来访问R中的属性表？

Answer 1

看起来Arc似乎更喜欢将RAT存储在单独的文件中，而不是存储在预期的tif中。你应该可以做这样的事情

library(raster)
r <- raster("a2_LOt_12_v170321.tif")

library(foreign)
rat <- read.dbf("a2_LOt_12_v170321.tif.vat.dbf")

colnames(rat)[1] <- "ID"
levels(r) <- rat