使用 raster() 包在 R 中导入 USDA 作物数据层 (GTIFF) 时出现问题

我已经下载了美国农业部作物数据层（此处提供数据，注意大文件：http://www.nass.usda.gov/research/Cropland/SARS1a.htm http://www.nass.usda.gov/research/Cropland/SARS1a.htm）并希望将其导入到R中进行分析。我导入了其他大型 GTIFF 文件（例如 LANDSAT 数据），但在调用该函数时出现以下错误：

r <- raster(file)

Error in .rasterObjectFromFile(x, band = band, objecttype = "RasterLayer",  :
    Cannot create a RasterLayer object from this file.

其中“文件”是 USDA 提供的任何 .tif 文件（我使用的是中西部地区编译）。

我怀疑问题在于 GTiff 是分类数据（土地利用类别）而不是真正的序数/光谱数据，但我不明白为什么它会在图像编辑器（例如预览）中呈现良好，而在 R / GDAL 中表现不佳。

谢谢您的帮助。

编辑：如果这应该出现在 GIS 论坛上，请告诉我。一般来说，空间问题的编程方法似乎就在这里，不是吗？

我能够使用同一个文件重新创建您的错误消息，并且 GDAL 似乎需要比主目录更深的路径名，至少在 Mac 上是这样。也就是说，路径名需要返回到 /Users/myname/cdl_tm_r_nv_2010_utm11.tif 而不是 ~/cdl_tm_r_nv_2010_utm11.tif 作为 USDA 作物文件示例之一。我设法打开文件

library(raster)
r <- raster(file.choose())

然后导航到 GeoTiff。