【已解决】求助-ncl傅里叶分析

是冉冉升起的冉

ncl官网中傅里叶分析用的ezfftf和ezfftb这两个函数，这个例子中分解1波和2波时：
https://www.ncl.ucar.edu/Applications/Scripts/fanal_3.ncl


CF = ezfftf (z)                             ; ezfftf works on right dim  
printVarSummary(CF)                         ; [2] x [12] x [121] x [120]  

cf = CF  
cf(:,:,:,1:)  = 0.0                         ; set waves 2 and higher to 0.0


z_wave1 = ezfftb (cf, 0.0)                  ; [12] x [121] x [240]  
copy_VarMeta(z, z_wave1)                    ; [time | 12] x [lat | 121] x [lon | 240]

这一行cf(:,:,:,1:)  = 0.0                         ; set waves 2 and higher to 0.0
求1波时设置2波以及更高的波数为0，那么为什么是设置最后一维1:为0，不应该是第一维，如果是最后一维那不就是经度么？

是冉冉升起的冉

顶顶我自己，两种设置我都试过了，值也print出来看过，如果是例子中就是经度，第一维就是波，两种画出来的图差异还挺大的。

穿着人字拖

因为ezfftb这个函数经度那一维展开成了2项，展开后的最后一维并不是经度，你可以看下这个函数的说明

是冉冉升起的冉

穿着人字拖 发表于 2022-12-7 15:39
因为ezfftb这个函数经度那一维展开成了2项，展开后的最后一维并不是经度，你可以看下这个函数的说明

感谢！！！

是冉冉升起的冉

ezfftf是对最右边的维度进行转换，如果需要对指定的维度转换，可以使用ezfftf_n。

气一小象

感谢！！

kaiofchina

{:5_275:}