读取栅格列表时出现错误:[rast] 无法打开文件:Apr 2017_masked_idw2.asc 是 terra 包?
我有一个包含 20 年的每月栅格。我只想将 2017 年的栅格作为文件列表读取。但我收到此错误 错误:[rast] 无法打开文件:Apr 2017_masked_idw2.asc 另…
Terra 两个栅格之间的空间相关性
我正在尝试计算两个栅格之间的空间相关性。我有两个具有相同范围、分辨率等的大栅格 class : RasterLayer dimensions : 45598, 53241, 2427683118 (nr…
R terra:: subst 未正确重新分类?
使用数字向量时,terra::subst 函数似乎返回错误地重新分类的栅格值: require(terra) # dummy SpatRast set.seed(234) r1 <- rast(xmin=0, xmax=100,…
如何在 R 中将 SpatVector 列转换为日期时间格式?
我想将 SpatVector 列从字符转换为 R 中的 POSIXct。使用 as.POSIXct() 不会引发错误或警告,但该列后面仍然是字符。 f <- system.file("ex/lux.shp",…
terra::rasterize 包含多个属性的点
我有一个具有 10 个点和 2 个属性 (i_1, i_2) 的 SpatVector。我需要在一个网格单元上总结我的 SpatVector 覆盖这 10 个点,并认为 terra::rasterize …
如何使用类别名称对 terra 中的分类栅格进行重新分类
我希望能够合并分类栅格中的两个类别。到目前为止我找到的唯一解决方案使用级别索引号,而不是类别名称。我如何使用类别名称来做到这一点? library(t…
如何从 terra 包中获取 SpatRasterDataset 的子集?
我正在尝试使用 R 中的 terra 包从 SpatRasterDataset 中选择多个栅格。但是,当我尝试一次选择 x 个栅格时,名称仅显示前 x 个栅格名称,无论我尝试…
使用 terra 函数 rast 创建栅格时的 nrow 限制
我正在尝试使用以下数据框(没有 NA)创建一个栅格,其中包含从 SeaDAS 生成的 HDF4 文件中提取的纬度、经度和值: > dim(df) [1] 10538622 3 > head(…
terra::classify 保留切点并将其用作 terra::as.polygon 中的类名
下面是一些显示问题的代码: 1 -- 重新分类栅格 建立类限制向量 class.limits #[1] 45.16490 50.12986 51.40991 52.68021 54.00000 55.42731 57.71864…