ubuntu - 有没有办法计算和保存一个新变量，它是多个 netCDF 或 tif 文件的函数？

Question

我有多个 .nc也在 .tif 中我想用来根据统计模型计算概率的文件。我有 5 个文件(V1.nc、V2.nc、V3.nc、V4.nc、V5.nc)，统计模型为:y = 1 / (1 + exp(-(Intercept + V1 * 1.87 + V2 * 1.53 + V3 * (-4.18) + V4 * 3.58 + V5 * 9.13))) .
我通常使用 Qgis , RStudio和 Ubuntu ，而且我不在乎使用这些程序中的哪一个。
有没有办法让我使用这些文件来计算整体 y ?

Answer 1

您可以尝试将文件分类为一个使用

cdo cat V?.nc combined.nc 

然后使用表达式命令创建一个新变量

intercept=2.3123 # or whatever the value is 
cdo expr,'y=1 / (1 + exp(-(${intercept} + V1 * 1.87 + V2 * 1.53 - 4.18*V3 + V4 * 3.58 + V5 * 9.13)))' combined.nc newvar.nc

我假设 V1.nc 中的变量称为“V1”等，如果没有相应地替换名称。
编辑 :根据注释，每个文件都使用相同的变量名。因此，要么在上述之前使用 NCO 在每个文件中重命名变量，要么可以使用单独的 cdo 函数组合原始文件，如下所示:

cdo -pow,-1 -addc,1 -exp -mulc,-1 -addc,$intercept -add -mulc,1.87 V1.nc -add -mulc,1.53 V2.nc -add -mulc,-4.18 V4.nc -add -mulc,3.58 V4.nc -mulc,9.13 V5.nc mynewvariable.nc 

我已经使用管道将它放在一条线上并避免中间文件，仔细检查我的方程式，但我希望它不包含错误。使用管道有助于从右到左阅读。
文件 mynewvariable.nc 中的变量也将称为 Band1；您可能希望使用 nco 适本地更改元数据，即使文件仅供您自己使用，也始终是一个好习惯......