问题标签 [r-raster]

我有一个rasterstack，我需要使用 shapefile 从每个栅格中提取值。

我知道使用“普通”光栅堆栈很容易，但我有一个RasterStackTS对象。不幸的是，我必须管理一个RasterStackTS对象，因为我已经从以前的计算中获得了 if。

这是一个可重现的例子，取自 Babak Naimi 的网站。

我尝试在 randomForests 的大型 1027206000 单元格栅格上预测森林林分类型（n = 43 个类别，编码地形、地质、水和养分供应等）。

在我用作协变量的许多 DEM 派生参数中，我还有 2 个栅格，其 ID 编号为地质图和土壤图。许多分类映射单元都带有这些 ID。我用数据框训练模型，并通过“合并”将分类映射单元附加到它。

到目前为止一切正常。该模型做了它应该做的事情，我也可以预测数据框中保存的一些测试数据。

但现在我打算制作一些预测地图。但是当使用光栅堆栈或砖运行模型时，只会给出所有 NA 的光栅。我的印象是，我在将因子级别传递给 rasterstack/rasterbrick 时做错了。

她是一些重现问题的代码。

有什么建议么？谢谢！

下面是一个非常有用的例子。我已经拼凑了在 ggplot 中创建的带有插图的地图的代码。

这是我遇到问题的地方。使用 grid.newpage 和 viewport 我的电脑在打印插图时会停止。我究竟做错了什么？您会注意到，我什至努力降低插入插图的国家形状文件的分辨率。

提前致谢。

我有 2 个栅格堆叠在一起：

由于我想稍后为堆栈中的每个特定图层突出显示一些区域，因此我基于栅格值创建了两个点数据集：

现在，我想用 R 中 rasterVis 的levelplot绘制光栅堆栈。我还想用 pts1 覆盖 r1，用 pts2 覆盖 r2。

但是，只要我添加一个点数据集，它就会用于两个图层：

如何在仍使用栅格堆栈的同时使用具有特定图层的特定点数据集？

我想避免制作自己的子图，方法是单独使用其特定点数据集绘制每一层，然后使用 print.trellis。我试过了，但与带有光栅堆栈的 levelplot 相比，结果只是差一些。

关于如何实现这一目标的任何想法？

我目前正在处理我的 R 脚本，似乎 raster 包中的旧 raster() 函数不再起作用。执行 raster("filename") 时出现以下错误 extension(x) 中的错误：4 个参数传递给 .Internal(nchar) 这需要 3 新版本的 raster 包是否发生了变化？或者扩展功能可能已经改变并且光栅功能仍然使用与更新前相同的方式？感谢您的帮助罗比

我正在对 .nc 文件中组织的气候数据进行日常观察。我使用 raster 包的 stack 命令读取它们。每个文件（对应一年）都是一个 RasterStack 元素，具有以下特征：

每一层都是一天的值的栅格。
我想对图层求和以计算每月值。我相信解决方案应该使用 calc 或 stackApply {raster}，但我找不到从第 x 层求和到第 y 层的方法或在求和之前对 RasterStack 进行子集化的方法。

我准备了一个只有 12 层的示例文件（以减小大小）。

我不完全知道如何提出代码，对不起，但它应该是这样的：

谢谢

我尝试将 RasterLayer 16 位转换为 RasterLayer 8 位，但没有成功。有任何想法吗 ？

谢谢！

我想预测randomForest一个巨大的物体RasterStack。
该randomForest对象dataframe使用名为“05_absor_1”、“05_absor_2”的变量进行训练……
有RasterStack189 个层，其名称类似于“X05_absor_1”、“X05_absor_2”……
我试图更改层的名称RasterStack，但是，当它发生时，它们不能以数字开头命名。

是否有可能通过重命名对象的变量名来使这两个数据集相互适应RandomForest？

我正在使用函数 densityplot() 和 grid.arrange()，这是示例：

使用此代码，我的输出是：

没有显示“title 1”，“title 2” .... 与我的数据一起使用的输出是：

equal 不显示密度图的主标题和标签（“Jan” “Feb” “Mar” “Apr” ....）

有可能请指导我。

谢谢你！！

I'm assuming the Raster package has what I need... I'm simply wanting to invert the colors in a Raster image.

The actual scenario is this: I want to invert the raster image returned by a ggmap call:

I want to invert the colors to get a white-on-black version of the Stamen Toner map: