问题标签 [terra]

我正在尝试复制我在raster使用 new(ish) terraR 包时所做的情节，但alphain 中的论点terra::plot似乎与在raster.

我正在尝试使用 alpha 层来指示基础（彩色）层中的不确定性。下面是一些表明我的问题的简化代码。使用该raster包，alpha 层会根据每个单独像素的值更改透明度。使用时terra似乎不起作用。

我想我正在使用最新版本的terra（1.3.4），并且我在 Mac 上（版本 15.5.11）。

对此的任何帮助将不胜感激。

我想使用包中的aggregate函数terra R以分位数方法聚合栅格作为聚合函数。下面，我使用quantile函数 fromR base使用本地包目录中的栅格计算第 50 个百分位数（即中位数）。我选择了第 50 个百分位数与中位数进行比较，但我的目标确实是计算其他分位数......

我的电脑花了大约。6秒做20次。

时差 6.052727 秒

当我aggregate使用中值内置函数运行相同的运行时，大约需要。执行相同的 20 次迭代的时间减少了 40 倍！

时间差 0.1456101 秒

由于我想计算第 50 位以外的其他百分位数，有人可以提供一些建议来加快aggregate使用自定义函数的速度吗？

我有多个多边形和一个栅格来组成我的地图。我想“放大”每个多边形，以便更好地可视化每个多边形内的栅格信息。我想我最终可以将facet_grid和的某种组合ggdraw用于最终地图，但现在我需要帮助创建地图的“放大”部分。

例子：

我知道raster和terra包都有一个zoom功能，但我无法将它与 ggplot 一起使用。欢迎任何指导，谢谢！

我想在使用r 包的freq函数输出中添加层名称来代替层号。terra这是一个可重现的示例

从输出层可以看出，代表层数。现在如何freq使用 terra 包在函数的输出中添加图层名称？

假设我们有一个数组modeldata（数据来自陆地模型），其维度为：

第一个维度包括网格的单元格，第二个维度是时间序列1500:2017
。第一个维度的不寻常长度是由于单独存在陆地单元格以节省空间。

在raster包中，我通过以下方式处理它：

请不要以这种方式评判我，
我对使用 terra 包的可比较（高性能）方法更感兴趣。
另一种很好的方法是使用SpatRaster遮罩层而不是坐标。

谢谢 ：-）

更新- 这个问题与 相关terra 1.4-1，现在已经过时terra 1.5-12（或更早？）。

我正在寻找terra相当于raster::gridDistance(..., origin = my_origin, omit = my_omit).

terra::gridDistance我从terra v0.2-8 这里发现了一个旧网页，但我可以收集terra::distance到的是当前的替换raster::gridDistance（替换功能列表terra's）。

但是，我不知道omit = my_omit在terra::distance. 从文档页面看来，任何非 NA 都被视为origin，但没有提及 ，也没有omit将 更改origin为特定值的选项（与raster::gridDistance下面的示例不同）。

这是来自的示例raster::gridDistance：

任何人都可以使用复制这个例子terra::distance吗？

我可以使用 做我需要做的事情raster，但我仍在学习如何正确处理空间栅格数据，并且正在努力学习新terra包。

我有一个全球 fAPAR 数据的栅格块b ，它是 35 年的日常数据（即大数据）。

目的

有没有办法可以循环遍历这块砖并一次在 365 层（天）上运行一个函数？最终目标是运行一个函数来识别每个网格单元的值超过该网格单元每年最大值的 20% 的阈值的层（日期）。

所以代码需要：

1. 计算出一年内每个网格单元的最大值
2. 返回每个网格单元超过该最大值的 20% 的图层（日期）

代码

以下代码查找每年的最大值，但由于其大小（每日、全球、0.5 度数据）而不适用于我的数据。关键可能是一次使用 365 层，而不是整个砖块（如果可能的话）。

我的目标与this和this非常相似。考虑到大小问题，我还想知道是否需要将我的数据划分为单独的年份并在多个文件上运行某些东西。

我在一个文件夹中有一组多波段栅格，我想使用单个 shapefile 进行遮罩。我希望屏蔽的输出（光栅）进入一个单独的文件夹。不久前，我在这里获得了帮助，以for循环完成此操作。它工作得很好，但是当我放大它（增加 shapefile 中的多边形数量，并添加更多栅格）时，它就不再工作了。具体来说，它有时会正确掩盖栅格，有时则不会。我无法辨别任何模式，因为我已经多次运行此代码，并且每次没有被屏蔽的栅格集都是不同的。

到目前为止的代码：

需要明确的是，我知道上面的代码运行正确。出于某种原因，它不能与我较长的数据集一起运行，我想知道是否有人可以阐明这种for循环的任何潜在缺陷。

我正在尝试使用R'sterra::nearest()函数，但我无法理解输出格式以及如何实际找出最近的点！这是一个简单的例子：

好的，所以to_x和to_y列distance对我来说很有意义并且看起来正确。但to_id对我来说没有任何意义。我会假设该列会有索引dest，但它显然没有，因为dest[1, ]不是。(0, 0)(5, 0)

我的最终目标是找出最接近每个点attr的点的值。如何找到最接近的点以便查找它的值？destsrcdestattr

（这是使用 terra 版本1.3-4）

terra使用地理坐标（经度/纬度）时，如何计算面积？

1 - 在球体或椭球体上？

2 - 使用在包装中完成的单元格的恒定宽度（中值宽度）raster，还是使用更准确的算法？

我的问题涉及expanse和cellSize功能。

我试图在terra手册和帮助页面中找到答案，但没有成功。谢谢您的帮助。