问题标签 [r-raster]

我想使用分类 Landsat 图像和参考数据集之间的混淆矩阵进行准确性评估。地图具有相同的分辨率和范围。我想逐像素评估协议。在许多研究中，我发现他们使用移动内核窗口（理想的 3x3）来“处理”Landsat 像素配准错误。但是我找不到任何方法来使用这个移动窗口在 R 中进行混淆矩阵评估，通常它用于值插值。

您对如何在分类准确性评估中实施移动窗口有任何想法吗？还是我误解了这种方法？

非常感谢，

示例：

给所有可能相关的人：

这是源代码：

我正在尝试将输出写入空的 TIFF 栅格，但我不断收到以下错误消息：

我想知道这是否与滥用 RGDAL 或 RASTER 包中的函数有关。

有人可以帮助我吗？

预先感谢您的慷慨。

干杯，广告

我在一个文件夹中有大约 40 个 .tiff 格式的空间栅格。我正在尝试从 R 中的每个栅格生成直方图，并将每个直方图作为 jpeg 保存在单独的文件夹中。我编写了代码来遍历每个光栅，创建一个直方图并使用“jpeg”包保存它。

我希望每个生成的 jpeg 都具有不同的名称，最好是原始栅格名称的子集。例如，如果栅格的原始名称为“MODIS101_265”，则 jpeg 的名称应为 265。这里，265 是年份中的儒略日期。我假设这可能涉及在 C 中使用像 %d 这样的格式说明符，但我不确定这在 R 中是如何工作的。

当我运行上面的代码时，我只得到一个直方图。似乎代码在原始栅格上正确循环，但将所有生成的直方图保存到单个 jpeg。

任何建议都会有所帮助！谢谢！

我是 R 和空间分析的初学者。我有一个栅格，其中栖息地多边形的单元格等于 1。

从这个栅格中，我如何为每个多边形构建一个栅格，其中位于给定多边形周围 1 公里缓冲区中的所有多边形单元都等于 1，所有其他矩阵单元都等于 0。

非常感谢您的帮助。

我已经将一个光栅文件导入到 R 中——它是分类的，有两个值，0 和 1。我需要将它存储为一个因子。

我尝试了以下方法：

我不能完全得到我想要的，这是一个在调用另一个函数（spdedp 包中的 joincount.multi 函数）时被识别为一个因素的层。任何提示将非常感谢！

我目前正在学习有关栅格图层的知识，我可以使用 inset_raster() 在 ggmap 上放置栅格图层。下面是一个示例代码：

请忽略地图上的点是否合理。我知道如何在使用plot(). 我也想知道如何改变ggmap中点的颜色。现在颜色总是白色、粉色、黄色和绿色（如果我希望颜色是蓝色黄色和红色怎么办）。是否有可以像 ggplot 中一样指定颜色的函数？

提前致谢

Does any one knows that ,how to remove background gray color when using gplot method in Rastervis? I already try to remove the theme color but it did not work.

我将 10 多个光栅文件（单波段文件，不是很大）从我的工作目录之外的文件夹中引入 R。单独调用文件有效，但它变得乏味。我见过人们使用类似这样的代码：

但是，它给我带来了各种错误。包括，“ .rasterObjectFromFile(x, band = band, objecttype = "RasterLayer", : 无法从此文件创建 RasterLayer 对象。（文件不存在） "

请注意，我尝试了几种不同的格式（.grd、.img），并且在运行命令之前通常还需要（rgdal）。提前感谢您的想法！

如何将栅格对象提取或子集化为一行数据？我有一个全球温度数据栅格，并想计算每个经度/列的每年平均月值。

所以基本上我想为剩余的 72 个网格单元中的每一个提供 12 个月的平均值。

在对一堆栅格进行 PCA 之后（类似于 this & in the 2014 Raster Package documentation），我想查看我的特征值、特征向量和加载...

对返回碎石图、变化比例、累积比例的 princomp 的典型调用 - summary()、print()、hist() - 似乎没有从我的 RasterBrick 输出中提取信息。这是示例代码：

summary()返回似乎是栅格图层的汇总统计数据，而不是分析中的值；print()显示栅格的最小值和最大值等。

感谢您的想法，尤其是关于如何找到与 PCA 相关的特征值和特征向量的任何清晰说明。