我正在尝试制作一个 SDM（物种分布模型），为此我构建了一个 GAM model（假设物种 A： 在这种情况下，“存在”指的是二元因子结构，其中 1 表示存在，0 表示不存在（背景点）。 此二元因子结构与一组数字栅格预测变量（例如，高程、坡度和各种生物气候变量）相结合。 注意：存在/不存在数据是全局 ...
我有一个包含大量层的堆栈，我想在总和变为 <0 时通过重置跨层运行像素累积和 function。 我编译了以下代码： ，效果很好，但我有一大堆光栅（就 nlyrs 而言）。 我想知道是否有更快的方法来执行此任务。 ...
我有一个 .netCDF 文件（“SSP119.nc”），其中包含 2015 年到 2100 年的土地利用变化信息，包括在内（来源： https://luh.umd.edu/data.shtml ）。 我还有一个包含全球关键生物多样性区域 (KBA) 多边形的 shapefile。 我想知道在某些年 ...
我想要通过瑞士的 map 绘制浮雕的官方来源。 我下载了 https://cms.geo.admin.ch/ogd/topography/DHM25_BM_SHP.zip来使用 dhm25_p.shp-File 现在使用代码 我会得到错误 ...
我想要通过瑞士的 map 绘制浮雕的官方来源。 在 swisstopo 的其中一页上，我找到了以下文件： https://data.geo.admin.ch/ch.swisstopo.pixelkarte-farbe-pk1000.noscale/swiss-map-raster1000_2022_ ...
从 19 个气候栅格中，我尝试仅裁剪多边形中的栅格区域（在法国）。 为此，我尝试生成一个循环，用多边形逐一裁剪每个栅格，如下所示： 它有效，但现在我试图找到一个循环来一次裁剪和遮罩 19 个栅格，并从我会命名的这个裁剪和遮罩中创建新栅格。 有什么建议么？ 我不知道是否清楚，如果您需要更多信息，请告 ...
我的问题遵循之前的问题，该问题已得到部分回答： Reading a grib2 dataset with 4 dimensions and 2 variables with R 我正在尝试读取带有 R 的 GRIB2 文件。此文件是概率气象预报：10 个变量、1 个提前期、17 个经度、23 个纬度 ...
我有一个范围从 0（最小值）到 100（最大值）的栅格图层，我想将其值更改为 1 到 100 的范围。重要的是数据不要改变，基本上我想要相同的栅格但不同范围的值。 我可以在 R 和 QGIS 上工作。 谢谢大家的帮助：:) ...
我需要对很多大栅格（28000 个像元，181 层）进行计算。 我在一小部分（24 个单元，181 层）上尝试了我的代码。 我从这个论坛获得了帮助，以尽可能多地进行优化。 现在我在光栅 package 中使用了砖块，因为我读到砖块被加载到 memory 中并且处理起来更快。 但后来有人认为 ter ...
我有一个测深栅格和一个边界网格化的 shapefile。 我无法隔离 lapply 生成的列表的“值”列。 目的是使用栅格平均值列作为 shapefile 中的新列 最后一行应该分配 Mean_Raster_Values 列表 object 中的整个第三列，但 [[3]] 只提供第三行如何从 Me ...
我有一个光栅砖（ncell=28536 和 nlayers=181）。 我需要在原始积木上运行数学函数并创建另外两块相同大小的积木。 两个 output 块相互依赖的地方。 inputBrick 有 181 层，每层有 28536 个单元格。 outputBrick1 将通过分析 outputB ...
我正在尝试使用 ggplot 突出显示 DEM 的一部分——本质上，只是一个彩色的、半透明的区域子集。 只要我不向叠加子集添加“alpha”参数，一切看起来都很好——我可以将我的 DEM plot 作为 geom_raster，并在顶部添加我的第二层，作为另一个 geom_raster 或 geo ...
我正在尝试使用 R 中的ranger R实现随机森林 (RF) 回归，但我收到此错误： Error: Missing data in columns: pop （ pop是我的自变量）在运行predict function 时。 作为参考，当使用randomForest package 时，我可以 ...
R 中的库 (rerddap) 中的 griddap function 提供了一个优秀的 function，用于以简单直接的方式从 ERDDAP 服务器访问卫星数据。 下面是一个可重现的示例，它使用ncdcOisst21Agg_LonPM180数据集提取大西洋的最新海面温度 (NOAA OISS ...
在示例栅格r ，如何将像素值84的比例从9%更改为25% ，如您在per_84中看到的那样？ 即我想在以下条件下增加像素值为84的像素数并减少其他值的像素数。 84的比例从9%增加到25%如下：栅格中已经有9% （共84个）。 在剩余的16%中， 6%来自像素值90 ， 5%来自像素值85 ，剩余 ...
我有一个 Sentinel-1 图像和 Sentinel-2 图像，我想让 Sentinel-1 中的像素与 Sentinel-2 重合。 更简单的解决方案是使用S2 <- resample(S2, S1)将两个图像放在相同的几何形状（空间范围和分辨率）中，但我担心光谱信息的完整性。 更好地重 ...
我是 R 中光栅 package 的新手，我试图用 package 光栅打开一个 .nc 文件，但弹出了一些错误。 如果您想尝试，我使用的是来自哥白尼的 2018 年和 2019 年每月海水盐度数据集（网格是魁北克圣劳伦斯河口和加斯佩西海岸）。 以前打开过类似的数据文件，但从来没有出现过这个错误， ...
我每天堆叠 181 个 ndvi 光栅。 我想为每个像素找到峰值，然后将峰值分配给 rest。原始数据曲线： 这就是我要的： 此结果适用于单个像素。 我已经在 excel 中完成了，但为了方便起见，我想在 R 中创建一个代码。 现在每个像素在不同的日期达到峰值（因此堆栈中的不同层）。 这是我试过 ...
我已将我的数据上传到 github： https://github.com/dschmidt42/Benguela_Chla_Danielle_Schmidt 。 设置我的数据： 以下是结果图： 正如您所看到的，尽管它们具有不同的值，但比例都是相同的。 我怎样才能改变它，使一张图表上的 4.0 与 ...
我已将我的数据上传到 github： https://github.com/dschmidt42/Benguela_Chla_Danielle_Schmidt 。 它涉及非洲东海岸的 chl a 数据，时间跨度为 26 年。 我试图将数据分成每个月，然后比较几个月内的差异。 这是我第一次使用 .n ...