raster是一个用于地理数据分析和建模的R包。
我需要从光栅堆栈数据中替换或插入每个单元格向量中的 NA。我有两个函数 fun_sub 用于替换 NA 和 fun_interp 用于插值 NA。 我发现 fun_sub wo...
我正在尝试平均我的降水数据,它是一个光栅砖形式的对象(该对象称为“Prec”)。 “Prec”有 95 层。然而,这个想法是
我正在使用 raster 和 terra 包中的 rasterize 函数将一些区域 shapefile 数据转换为栅格。 然而,我发现所有栅格最终都具有与
我对光栅包中的“stackApply”函数有问题。首先,我想堆叠三个栅格层(每个层有一个波段) - 这是可行的。然后我想创建一个光栅-
我下载了世界气候数据并将其更改为栅格数据。 直到 19 为止都有像 wc2.1_5m_bio_1 这样的名称,我想使用 substr func 将这些变量重命名为 bio_1 (start = 10, stop = 16)...
#####首先将单个文件夹中的所有文件作为列表导入##### 列表 <- list.files(path = ".", pattern='.tif$', all.files=TRUE, full.names=FALSE) ######load them in a stack ######
我的项目工作流程的一部分需要从栅格创建轮廓。我能够使用 raster::rasterToContour 生成轮廓,并可以使用 R 中的 sf 包转换为多边形...
我正在尝试使用 SLGA 包来查找澳大利亚土壤的沙子含量。 我在工作时受到防火墙保护,因此之前无法从 github 安装该软件包,但现在可以下载...
我生成了一个堆栈“class.stack”,其中包含 175 个栅格。每个栅格中的每个像元可以是四个类别之一(用整数 0-3 表示)或 NA。 NA 应该出现在完全相同的位置
我正在使用随机森林生成物种分布模型。这些模型试图根据各种环境属性来预测某个物种发生的概率。对于我...
我有一组栅格,我想使用简单的栅格数学从另一组栅格中减去它们,但在将正确的栅格与其伙伴匹配时遇到困难。 图书馆(土地) 斯...
我正在尝试使用光栅包中的堆栈功能将多个光栅文件(.tif)放入一个堆栈中。当我这样做时,我遇到了两个错误,即: > covStack <- stack(paste0(getwd(),&q...
如何使用来自多个 csv 的传单在 R 中创建热图?合并/重新采样问题
我一直以此为灵感。我有 c.400 个 GPS 文件,其中的 GPS 位置位于 csv 中,每个 csv 都有一个文件属性,我想保留该文件属性以根据该属性创建热图。 https://gis.
terra::clamp() 比 ´<- NA´ on a SpatRaster?
我有一个 **SpatRaster **并且想要确保 NA 标志(在我的例子中为 -9999)正确设置为 NA。 我应该使用 r1 <- terra::rast() r1[r1 < 0 ] <- NA or using r1 <- terra::rast() t...
基本上,我以 ASCII 的形式计算了一个全局分布概率模型,比如: 国内生产总值。 gdpm的值都在0到1之间。 然后我从形状文件导入了本地地图: 形状 <- file.
大型栅格的并行处理(拆分为多个部分,运行函数,然后重新组合)
总体目标是在堆叠栅格上运行相关分析,其中我使用移动窗口并将焦点像素堆栈与移动窗口中每个其他像素的堆栈相关联。最后
我想使用 bbox 或已知范围(即行和列 10 个像素)裁剪栅格。 下面您可以看到一个可重现的示例: 图书馆(土地) <- terra::set.ext(rast(volcano), terra::ext(0, ...
我导入了一个包含不同森林补丁(称为“补丁”)的多多边形 shp 文件: 修补 然后,我在名为“SenBel”的栅格中为同一区域砖砌了不同的 Sentinel-2 波段 d...
terra 包中的 rast 函数在加载一个 tiff 或列表时将图层名称更改为 bigfile[,varname]
所以这是我的第一个问题。令人兴奋!并对我可能犯的任何错误提前表示歉意。我当然可以进行编辑以帮助解决我的问题。 我正在处理来自 ClimateNA 的 tif 文件栅格...
我已经开始使用 elsa 包(https://cran.r-project.org/web/packages/elsa/vignettes/elsa.html)来测量分类数据的基于熵的空间关联,但它旋转很多