使用R省略netcdf文件中的某些数据点

时间:2015-02-11 04:50:47

标签: r

我一直试图用6个变量操作netcdf文件(让他们将它们命名为1,2,3,4,5,6),每个变量都有尺寸(x,y,z,t)。我想省略其中一个变量(比如变量4)大于某个值(比如10)的所有数据行,然后使用所有变量执行一些算术函数。

有关最佳方法的任何输入吗?

我已经能够创建一个新变量,其中省略了变量4的不需要的值,但是这样做会丢失网格信息,并且我无法使用变量4对此文件执行进一步的算术运算。是否存在一种方法来重建索引"这个新的变量到它的旧网格位置?

我应该替换变量4的所有值吗? 10用NA或其他东西,然后操纵数据?

任何建议/想法都会非常有用。

我不知道如何发布此数据的可重现示例,但这里是open.ncdf(filename)的输出:

[1] "file clouds.ncf has 6 dimensions:"
[1] "TSTEP   Size: 240"
[1] "DATE-TIME   Size: 2"
[1] "LAY   Size: 1"
[1] "VAR   Size: 6"
[1] "ROW   Size: 299"
[1] "COL   Size: 459"
[1] "------------------------"
[1] "file clouds.ncf has 7 variables:"
[1] "int TFLAG[DATE-TIME,VAR,TSTEP]  Longname:TFLAG            Missval:NA"
[1] "float Augm3[COL,ROW,LAY,TSTEP]  Longname:Augm3          Missval:1e+30"
[1] "float Atry [COL,ROW,LAY,TSTEP]  Longname:Atry     Missval:1e+30"
[1] "float H2OT[COL,ROW,LAY,TSTEP]  Longname:H2OT         Missval:1e+30"
[1] "float AS[COL,ROW,LAY,TSTEP]  Longname:AS         Missval:1e+30"
[1] "float salt[COL,ROW,LAY,TSTEP]  Longname:salt            Missval:1e+30"
[1] "float rain[COL,ROW,LAY,TSTEP]  Longname:rain   Missval:1e+30"

谢谢!

0 个答案:

没有答案