r - 使用 raster() 包在 R 中导入美国农业部裁剪数据层 (GTIFF) 时出现问题

Question

我已经下载了 USDA Crop Data Layer(数据可在此处获取，请注意大文件: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.

其中"file"是美国农业部提供的任何 .tif 文件(我使用的是中西部地区汇编)。

我怀疑问题是 GTiff 是分类数据(土地利用类别)而不是真正的序数/光谱数据，但我不明白为什么它会在图像编辑器(例如预览)中渲染得很好并且犹豫不决R/GDAL.

感谢您的帮助。

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

Answer 1

我能够使用同一个文件重新创建您的错误消息，而且 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。