在 R 中导出栅格图层答案

【问题标题】：exporting raster layer in R在 R 中导出栅格图层
【发布时间】：2018-03-28 08:57:17
【问题描述】：

我正在处理bnspatialexample，我想看看导入的数据ConwyData中的实际内容：

source("http://bioconductor.org/biocLite.R")
biocLite("RBGL")
install.packages("bnspatial")
install.packages("raster")
library(raster)
library(bnspatial)
data(ConwyData)

> str(ConwyData)
List of 6
 $ ConwyStatus  :Formal class 'RasterLayer' [package "raster"] with 12 slots
  .. ..@ file    :Formal class '.RasterFile' [package "raster"] with 13 slots
  .. .. .. ..@ name        : chr "E:\\masanda\\Documents\\R\\win-library\\3.4\\bnspatial\\extdata\\ConwyStatus.tif"

因此，我想将 Rater Layer 文件导出为可以导入 ArcGis 或 QGIS 的格式。我尝试了以下方法，但都没有工作：

writeRaster(ConwyData, "raster_bnspatial", format = "GTiff") （函数（类，fdef，mtable）中的错误：无法找到用于签名“列表”、“字符”的函数“writeRaster”的继承方法
ConwyData

1:https://cran.r-project.org/web/packages/bnspatial/vignettes/bnspatial.html

【问题讨论】：

标签： r geospatial r-raster sp geotiff

【解决方案1】：

看起来ConwyData 是一个对象列表，而不是单个栅格。这也是为什么您不能将其作为光栅写入磁盘的原因。

查看example documentation（或?ConwyData），列表包含以下元素：

LandUseChange 谷物类的对象。贝叶斯网络，用于演示目的。

ConwyLU RasterLayer 类的对象。 Conwy 集水区（英国威尔士）当前土地利用地图的简化版本。它包括三个类别：耕地（栅格值 3）、森林（2）、其他（1）。

ConwySlope RasterLayer 类的对象。从分辨率为 50 米的数字高程模型导出的坡度栅格，单位为度。

ConwyStatus 类 RasterLayer 的对象。土地所有权类型（虚拟数据），分为三个可能的类别：公共（栅格值 4）、私有（3）、受保护（1）。

证据矩阵。从流域中的每个位置（即像元）提取的可用空间数据（见上文）的集合，其中后者由栅格对象 ConwyLU 表示。空间数据中的每个值都通过 dataDiscretize 或 bulkDiscretize 函数进行离散化，然后从贝叶斯网络（LandUseChange）分配到相应的状态。

LUclasses 具有输入空间数据分类的列表（其对应的状态和值）。该列表根据 bnspatial 函数要求进行格式化，并由函数 importClasses 和 setClasses 返回。

如您所见，其中只有一小部分是实际的栅格。这些您可以轻松地写入磁盘：

writeRaster(ConwyData$ConwyLU,'/tmp/landuse.tif')

但如果你只是想看看，你可以跳过写入文件，直接用 R 可视化它，例如

base 绘图：

plot(ConwyData$ConwyLU)

或者更通用的rasterVis：

library(rasterVis)

levelplot(ConwyData$ConwySlope,margin=FALSE)

【讨论】：