【问题标题】:When averaging rasters, is there a way to calculate the n number of non-NaN values per pixel used to calculate the pixel means?在平均栅格时,有没有办法计算用于计算像素均值的每个像素的 n 个非 NaN 值?
【发布时间】:2020-11-21 02:53:26
【问题描述】:

我有大约 40 个栅格 (tiff)。我平均了栅格,这很简单。但是,作为第二个输出,我想确定用于计算每个像素均值的每个像素的非 NaN 值的 n 个数,并将其显示为栅格。基本上,我需要知道每个像素有多少值用于计算平均值。这可能吗?我正在寻找使用 GIS(ArcGIS 或 QGIS 都可以)、Python 或 R. Cheers 的解决方案。

【问题讨论】:

    标签: python r raster arcgis qgis


    【解决方案1】:

    以下是使用 R 使用 terra 包的方法

    library(terra)
    s <- rast(system.file("ex/logo.tif", package="terra"))   
    
    m <- mean(s)
    n <- sum(!is.na(s))
    

    raster

    library(raster)
    b <- brick(system.file("external/rlogo.grd", package="raster"))
    mb <- mean(b)
    nb <- sum(!is.na(b))
    

    【讨论】:

    • 抱歉,我只是修改了我的评论,说我使用了光栅包说明并且它有效。
    猜你喜欢
    • 2019-05-03
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    相关资源
    最近更新 更多