【发布时间】:2020-11-21 02:53:26
【问题描述】:
我有大约 40 个栅格 (tiff)。我平均了栅格,这很简单。但是,作为第二个输出,我想确定用于计算每个像素均值的每个像素的非 NaN 值的 n 个数,并将其显示为栅格。基本上,我需要知道每个像素有多少值用于计算平均值。这可能吗?我正在寻找使用 GIS(ArcGIS 或 QGIS 都可以)、Python 或 R. Cheers 的解决方案。
【问题讨论】:
标签: python r raster arcgis qgis
我有大约 40 个栅格 (tiff)。我平均了栅格,这很简单。但是,作为第二个输出,我想确定用于计算每个像素均值的每个像素的非 NaN 值的 n 个数,并将其显示为栅格。基本上,我需要知道每个像素有多少值用于计算平均值。这可能吗?我正在寻找使用 GIS(ArcGIS 或 QGIS 都可以)、Python 或 R. Cheers 的解决方案。
【问题讨论】:
标签: python r raster arcgis qgis
以下是使用 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))
【讨论】: