3

我有几个 netcdf 文件。每个 nc 文件都有几个变量。我只对两个变量“Soil_Moisture”和“Soil_Moisture_Dqx”感兴趣。

我想根据“Soil_Moisture_Dqx”过滤“Soil_Moisture”。每当相应的“Soil_Moisture_Dqx”像素的值大于 0.04 时，我想用 NA 替换“Soil_Moisture”中的值。

:这里是要下载的文件：

1- 我尝试了这个循环，但是当我输入 f[1] 或 f[2] 时，我得到了一些奇怪的东西，这意味着我的循环不正确。我很感激任何帮助我纠正我的循环。

 a<-list.files("C:\\3 nc files", "*.DBL", full.names = TRUE)

for(i in 1:length(a)){
f=open.ncdf(a[i])
A1 = get.var.ncdf(nc=f,varid="Soil_Moisture",verbose=TRUE)
A1* -0.000030518509475997 ## scale factor
 A2 = get.var.ncdf(nc=f,varid="Soil_Moisture_Dqx",verbose=TRUE)
A2*-0.0000152592547379985## scale factor
A1[A2>0.04]=NA ## here is main calculation I need
 }

2- 谁能告诉我再写一遍？

4

2 回答 2

3

缺失值是 netCDF 文件中的特殊值，其值被视为表示数​​据“缺失”。所以你需要使用set.missval.ncdf来设置这个值。

a<-list.files("C:\\3 nc files", "*.DBL", full.names = TRUE)

SM_NAME <- "Soil_Moisture"
SM_SDX_NAME <- "Soil_Moisture_Dqx"
library(ncdf)
lapply(a, function(filename){
  nc <- open.ncdf( filename,write=TRUE )
  SM <- get.var.ncdf(nc=nc,varid=SM_NAME)
  SM_dqx <- get.var.ncdf(nc=nc,varid=SM_SDX_NAME)
  SM[SM_dqx > 0.4] <- NA
  newMissVal <- 999.9
  set.missval.ncdf( nc, SM_NAME, newMissVal )
  put.var.ncdf( nc, SM_NAME, SM )
  close.ncdf(nc)
 })

编辑添加一些检查

计算有多少点将被标记为错过是很有趣的。

在不应用奇数比例因子的情况下，我们有：

lapply(a, function(filename){
  nc <- open.ncdf( filename,write=TRUE )
  SM_dqx <- get.var.ncdf(nc=nc,varid=SM_SDX_NAME)
   table(SM_dqx > 0.4)
  })

[[1]]
[1] 810347     91

[[2]]
[1] 810286    152

[[3]]
[1] 810287    151

[[4]]
[1] 810355     83
于 2013-03-05T15:16:49.860 回答
0

这也可以使用 CDO 从命令行完成。

据我了解，这两个变量都包含在您的输入文件中（我将其称为“datafile.nc”，您可能希望在文件列表的循环中执行以下操作），所以首先我们将提取这两个变量分成两个单独的文件：

cdo selvar,Soil_Moisture     datafile.nc soil_moisture.nc
cdo selvar,Soil_Moisture_Dqx datafile.nc dqx.nc

现在我们将定义一个掩码文件，当 dqx<0.04 时包含 1，但当 dqx>=0.04 时包含 NAN

cdo setctomiss,0 -ltc,0.04 dqx.nc mask.nc

ltc 是“比常数”（您可能希望 lec 代替 <= ）， setctomiss 用 NAN 替换所有零。

现在我们将它们与 CDO 相乘 - NAN*C=NAN 和 1*C=C，因此这会为您提供一个带有所需字段的 netcdf：

cdo mul mask.c soil_moisture.nc masked_soil_moisture.nc

如果您愿意，您实际上可以将最后两行组合在一起，并避免编写掩码文件的 I/O：

cdo mul -setctomiss,0 -ltc,0.04 dqx.nc soil_moisture.nc masked_soil_moisture.nc

但是单独解释这些步骤更容易:-)

您可以在 bash 中轻松地将整个内容放在文件的循环中。

于 2017-09-20T12:09:11.080 回答