計算netcdf文件中缺失值的數量-R

Question

有沒有一種快速的方法來知道netcdf文件中有多少個缺失值？ 可能使用R。

目前我必須

hum<-nc_open("rhum.sig995.2008.nc")
rhum<-ncvar_get(hum, "rhum")

然后通過鍵入'hum'手動查找缺失值並將其復制到此操作中

sum(abs(rhum - 9.96920996838687e+36) < -9.96920996838687e+36)

有沒有更直接的方法，特別是如果我必須處理數百個文件時？ 我想避免復制和粘貼缺失的值，而且我不確定應該以哪種精度處理數字。

Answer 1

我的建議是使用出色的柵格數據包：

install.packages(raster)
library(raster)

r <- raster("rhum.sig995.2008.nc", var="rhum")

NAnum <- summary(r)[6]

Answer 2

變量名稱“ var”的丟失數據點總數可以使用以下命令存儲在新的附加變量中：

ncap2 -s "nmiss=number_miss(var)" in.nc out.nc

要么

ncap2 -s "nmiss=var.number_miss()" in.nc out.nc

如果您的數據具有時間維度，並且您想查看空間維度上總計的缺失點總數，則可以使用

cdo info in.nc