主页 / coding / 问题 / 79528685
Accepted
Breeze
Asked: 2025-03-23 16:55:14 +0800 CST

对于相同的经度或纬度栅格数据使用基于 terra 的函数

  • 772

假设我需要计算同一经度的栅格数据的方差，使用该sapp函数是正确的方法吗？但我失败了。

我现在可以使用更复杂的方法来转换为 data.frame 并分组来计算方差：

library(terra)
library(tidyverse)

r <- rast(nrows = 5, ncols = 10, vals = c(1:50))
plot(r)


r |> 
  terra::as.data.frame(xy= TRUE) |> 
  group_by(x) |> 
  summarise(lonsd = sd(lyr.1))

有没有更直接的函数terra？如果你知道的话请告诉我。

还有更复杂的场景，需要同时对多层栅格数据进行此操作，该怎么做？

1 个回答

  1. Best Answer

    您可以使用聚合来计算栅格中每列（或每行）的内容

    示例数据

    library(terra)
set.seed(1)
r <- rast(nrows = 5, ncols = 10, vals = sample(100), nlyr=2)

    解决方案

    a <- aggregate(r, c(nrow(r), 1), var, na.rm=TRUE)
data.frame(lon=xFromCol(r), values(a))
#    lon  lyr.1  lyr.2
#1  -162  468.2 1364.3
#2  -126 1007.3  493.0
#3   -90  557.5  822.2
#4   -54  358.3 1089.7
#5   -18 1224.3  665.7
#6    18 1125.3  260.3
#7    54 1180.3  207.2
#8    90  631.3  433.0
#9   126 1564.7 1669.7
#10  162  881.7  752.3

    或者

    a <- as.data.frame(r, xy=TRUE)
aggregate(a[,3:4], a["x"], var, na.rm=TRUE)
    • 1

相关问题