问题标签 [r-raster]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
631 浏览

r - 光栅包,替代品

我正在尝试重新分类 aRasterLayer以将所有 0 变为 NA。我试过了:

  1. 重新分类:

    错误:无法分配大小为 160.1 Mb 的向量

  2. 潜艇:

    错误:在为函数“as.matrix”选择方法时评估参数“x”时出错

  3. 替换不适用:

    错误:无法分配大小为 160.1 Mb 的向量

我正在使用 R3.0.0。Windows 7 企业版,32 位。我知道内存可能是一个问题,但还没有看到有关如何在 Windows 中更改的任何方向。我也在这段代码中使用过:

我也会在 64 位计算机上尝试这些。

任何人都知道这台机器发生了什么,或者可以推荐一种让这些功能之一工作的方法?
我没有用方法 2(只是 a data.frame)调用任何矩阵,所以也不确定这意味着什么。

0 投票
1 回答
3001 浏览

r - 如何从R中的范围对象中获取坐标?

让我们有一个范围对象:

xmin(e)如何以比仅使用, ymin(e)....更紧凑的方式从中获取范围坐标?

编辑:我实际上需要这个来保存范围以drawExtent供进一步使用来构建extent

将此作为参数复制并粘贴到我的脚本中。完毕。

0 投票
1 回答
2866 浏览

r - 从具有小多边形的栅格中提取()数据 - 舍入的权重太小

使用 R,我正在尝试使用多边形图层从栅格图层中提取数据。多边形比栅格像元小得多

在此处输入图像描述

现在我extract()raster图书馆打电话:

有两个问题——权重是多边形覆盖的单元格区域的比例,权重四舍五入为1/100。在我的例子中,输出中只有左上角的单元格(值 209) - 其他 3 个单元格的权重四舍五入为零,它们被排除在外。但是,左下角的单元格覆盖了多边形的很大一部分,也应该包括在内!

我需要一个适当的加权平均值。这可以用其他方式完成extract()吗?还是有什么其他方式?

PS:请注意:我认为权重extract()设计得不是很好-权重应该是特定单元格覆盖的多边形区域的比例,反之亦然。然后,多边形的加权平均值也将更容易计算(只需将每行中的两个数字相乘并相加),四舍五入到 1/100 不会是一个大问题。

可重现的示例- (下载文件- 简化版,实际数据要大得多):

相关: 对于小多边形和栅格,R 中的提取失败

0 投票
1 回答
3336 浏览

r - 计算网格中的物种出现

我有大约 500,000 个点R 的候鸟在美国各地的出现数据。

我试图在这些点上叠加一个网格,然后计算每个网格中出现的次数。计算完计数后,我想将它们引用到网格单元 ID。

在 R 中,我使用该over()函数来获取范围图中的点,这是一个 shapefile。

上面的工作完全符合我的希望,但没有解决我当前的问题:如何处理 type 的点SpatialPointsDataFrame和栅格的网格。您是否建议对栅格网格进行多边形化,并使用我上面指出的相同方法?或者另一个过程会更有效吗?

0 投票
2 回答
588 浏览

r - 来自矩阵的栅格叠加

我有一个包含 100 个栅格图层的矩阵,我想创建一个新图层,即平均值。我知道如果有两层,我可以简单地使用覆盖功能或者只使用c <- mean (a, b). 但是,我不确定如何处理矩阵。

这是矩阵的示例:

我试过了

我收到错误argument is not numeric or logical: returning NA

0 投票
1 回答
664 浏览

r - 可以为 rasterImage 指定调色板吗?

标题说明了一切,差不多。据我所知,尝试通过使用参数传递一个简单的调色板col,例如, rasterImage(mydata,xleft=1,ybottom=1,xright=2,ytop=2,col=grey(0:4095/4096))失败。更令人困惑的是,默认调色板似乎取决于mydata.

如果我有一个矩阵,rasterImage(as.raster(mymatrix), {etc}则在使用RasterLayer对象时以灰度rasterImage(as.raster(myRaster),{etc}绘制,以拓扑颜色集绘制。

任何关于绘制大图像的替代(快速)工具的建议都非常感激。

0 投票
1 回答
3104 浏览

r - R如何控制彩条/图例的间距

我想绘制一张叶绿素浓度图,但这些值的分散方式使得图例变得不可读(见图)

chl_a

所以我试图控制颜色条/图例中颜色的间距。我想让它均匀分布(同时保持地图本身的不均匀中断)。

下面的示例非常简化,适用于光栅图,但 image.plot 也是如此。

我考虑将数据转换为日志值,但这并没有给我一个令人满意的结果。因此,非常感谢您对传说的任何帮助。

0 投票
2 回答
196 浏览

r - 如何在 r 中的可重现数据中生成 NA

将产生rasterstack在每个像素中都有值的值。

我怎样才能产生一个rasterstack包含NA值的?

0 投票
1 回答
83 浏览

r - 如何搜索非 NA 行、绘图并移至下一行?

我需要在我的 rasterstack 中搜索非 NA 像素,绘制它并移动到下一行。
例如,为了在光栅堆栈中绘制一个像素,我使用了 plot(c(s[2]))。我如何告诉 R 搜索、跳过 s[9] 和 s[10] 但绘制 s[11] 和所有非 NA 像素(我的栅格堆栈很大,我无法手动滚动数据)。

0 投票
1 回答
606 浏览

r - 如何从R中的rasterstack中提取回归系数?

我正在尝试将正弦曲线拟合到数据集并导出系数。由于这是我第一次处理 rasterstack,我不确定如何打印系数。是否可以用栅格层来表示系数?

提前致谢。