问题标签 [r-raster]

我被分配了一项我想在 R 中尝试的特定任务。要导入的 shapefile (SpatialPoints df) 包含多个属性，但最重要的是，特定点坐标 (lat/lon) 的商业捕获权重。

我需要一个脚本：

1) 创建一个网格（可以修改大小和单位） 2) 与导入的文件相交，以便按网格正方形汇总样本（平均值、标准差、范围等）。

我可以通过 ArcGIS 做到这一点，但我对修改网格大小和通过 R 具有可重用算法的易用性感兴趣。下面是正在使用的数据的简短示例。

有人知道该怎么做吗？

我想创建一个具有四个波段的图像子集。因此我在R中使用crop函数。

结果，我在 .tif 文件中得到了一个只有一个波段的栅格。我是否必须定义其他选项才能获得具有 4 个波段的子集图像？

我有一个简单的栅格（使用 R-package: raster 创建）。使用函数“rasterToPolygons”我得到包含值“1”的所有栅格单元的多边形：

我不使用“dissolve = T”来避免所有多边形合并成一个大多边形。相反，我希望获得一个新的 SpatialPolygonsDataFrame，其中合并了所有共享边或点的多边形。明确分离的多边形应可识别为单独的多边形。基于新的 SpatialPolygonsDataFrame 我想分析组合多边形的大小如下：

我的问题是双重的：1）如何组合共享边或点的多边形？2）如何将这些信息转换为允许分析组合多边形区域的格式？非常感谢您的反馈。

这个例子应该是可重现的，第一部分（大多边形）有效，第二部分失败：

耶！一切似乎都很好，是吗？现在有了小多边形：

失败并出现错误：

pmax(ceiling(10/dim(r)[1:2]), c(1, 1)) 中的错误：找不到对象“r”

如何将像素值添加到绘图中？我可以通过使用来获取值，click()但我希望它出现在图中。

我惊讶地发现，我从一个大型栅格创建的矩阵据称占用的内存是父级的 35k 倍。此代码演示：

object.size(region_utm)肯定是低估了，因为 12,128 字节不足以包含 53m 个值，即使考虑在内，尤其是 87%（46m）是唯一值（根据length(unique(region_utm))）。不知道如何获得现实的内存估计..

然而，绘制栅格比使用矩阵要快得多。但我一直认为矩阵大致相当于没有空间数据槽的栅格，但我一定错过了这些数据结构之间的一个重要区别。感谢您澄清什么可以解释这种记忆差异。我经常使用矩阵工作流程，需要了解它们的局限性。

编辑：预期报告请求str()：

大家好，

我有一个包含许多数据列的 data.frame，但是，对于每一行，我只对这些列的子集感兴趣。我想使用另一个特定值来标识我感兴趣的列。然后我将取感兴趣的列和前 5 个列的平均值。

我的 data.frame 包括点位置、收集月份和从 1996-2012 年每个点的一组月度栅格中提取的值。对于每个点，我都对收集日期之前的六个月平均值感兴趣，例如。如果我在 200106 (06/2001) 记录了一个变量，我想要 200101-200106 的栅格平均值。

收集日期值的编码与对应于同一月份提取的值的列名称相同。

给定我的收集日期，有没有办法识别我感兴趣的列？

我的 data.frame 看起来像：

给定 dt (date) 我可以得到相应列的六个月平均值吗？

我不知道如何继续，我想象某种数据转换，但不知道从哪里开始。任何帮助将不胜感激。非常感谢！

干杯，亚当

我有一个从多边形 shapefile 转换的空间线（根据“imagebrick”中的特征手动数字化 - 这意味着“折线”和“imagebrick”在空间上按照我的意愿重叠）

还有一个光栅砖

我试图提取线（空间线）触及的光栅砖中的像素。但我也想提取与这些提取像素相关的 xy 坐标，并将它们一起存储在 SpatialPointsDataFrame 中。我试过了：

然后，我试图使用 RasterToPoints 函数来提取 xy 坐标。这需要一个栅格层，其中包含我想要 xy 的所有像素。我尝试使用“裁剪”来获取该光栅“线”，但它不起作用。

似乎“裁剪”功能不会裁剪栅格的“线”，而是裁剪整个“多边形”。有没有其他方法可以获得我提取的像素的 xy 坐标？

ps 我还试图提取每个提取像素的周围像素。我该怎么做？我注意到提取函数中有一个函数参数。我应该在提取函数中编写一个小函数还是有其他方法可以做到这一点？

pps 我也在考虑使用 Rasterize 函数来栅格化“SpatialLine”，然后使用 RasterToPoints 拉出 xy 坐标。但我无法控制“光栅化空间线”的每个像素等于原始图像砖中“线”所触及的像素。

感谢大家花时间查看我的问题，并提前感谢您的帮助。

我无法使用 R 中的 raster 包对 shapefile 进行光栅化。我遵循了上一篇文章中的代码和提示（“Layer wont rasterize using raster package in R，2013 年 6 月 30 日）。但是，我仍然遇到同样的错误”试图从没有槽的基本类（“NULL”）的对象中获取槽“坐标”。这是与上一篇文章不同的问题吗？非常感谢！

我有一个大RasterBrick对象，其中包含表示每个时间间隔（不规则）的生物量的时间序列栅格层。我需要的是两个连续时间段（difference in biomass = current biomass layer - last weeks biomass layer）之间生物量差异的时间序列。我的想法是使用循环或其中一个应用函数来遍历，RasterBrick并且为每个RasterLayer应用减法函数与时间序列中较早的函数。是有序的RasterBrick，因此实际的时间戳不一定很重要。我试图找到例子，但没有取得进展。任何指针将不胜感激。

我在下面提供了一个关于我的情况的快速示例：

随机设置 10 个栅格图层，放入 RasterStack

现在我需要一个RasterBrick包含 9 层的结果

生物量差异 = 当前生物量层 - 上周生物量层。