NCL小波分析的一些问题，求教求讨论哇~~~

Sunshining

借鉴的是官网上http://www.ncl.ucar.edu/Applications/wavelet.shtml的第三个例子的脚本，我的时间序列是1948-2015年春季的海温资料，使用减去11年滑动平均的方法来滤去年代际变化，有用的序列长度变成1953-2010，以这个序列输入，用官网给出的脚本画出了下面的图，我有个问题，我看了下脚本，显著性检验的部分是这么写的：
; compute significance ( >= 1 is significant)                                       
  SIG              = power                            ; transfer meta data         
  SIG              = power/conform (power,w@signif,0)                              
  SIG@long_name    = "Significance"                                                
  SIG@units        = " "         
官网上给的两个包括显著性检验的部分都是这么写的，我想知道这个进行小波分析，是所有的都是SIG>=1就是显著的吗？ 画图画SIG部分的代码是这样的 ：
opt                      = True                                 
opt@gsnShadeHigh         = 6           ; cross hatch pattern     
opt@gsnShadeFillType     = "pattern"                             
iplot_hatch_def = gsn_contour_shade(iplot_hatch_def,0, 0.8, opt)
这是画线的部分不是小于0和大于0.8的吗？ 为什么跟之前的>=1是通过显著性检验不一样呢？                                                

Sunshining

{:eb302:}{:eb302:}求教求教哇~~

Sunshining

wavelet函数中这些参数的设置，是根据什么来的呢？mother  = 0                                                                           
     param   = 6.0                                                                          
     dt      = 0.25    ;timesteps in units of years                                         
     s0      = dt                                                                           
     dj      = 0.25                                                                        
     jtot    = 1+floattointeger(((log10(N*dt/s0))/dj)/log10(2.))                           
     npad    = N                                                                           
     nadof   = 0                                                                           
     noise   = 1  
     isigtest= 0                                                                          
     siglvl  = .05                                                                          
     nadof   = 0                                          
比如我的58年的数据，这些参数设置要改动吗？

liyf

Sunshining 发表于 2017-4-30 13:18
wavelet函数中这些参数的设置，是根据什么来的呢？mother  = 0                                           ...

ncl官网有函数wavelet参数的解释，另外ncl这个函数和matlab里面的wavelet函数是一样的，看懂一个就行

yeah...

liyf 发表于 2017-5-2 11:47
ncl官网有函数wavelet参数的解释，另外ncl这个函数和matlab里面的wavelet函数是一样的，看懂一个就行

你好 我想请教下 我在看官网的例子的时候发现
ninoseas = asciiread("sst_nino3.dat",-1,"float")
  ninoseas!0 = "time"
  ntime = dimsizes(ninoseas)
  timeo = fspan(1871.25,1996.,ntime)
  ninoseas&time = timeo

  ninomam = dim_avg_Wrap(ninoseas(0::4))
  ninojja = dim_avg_Wrap(ninoseas(1::4))
  ninoson = dim_avg_Wrap(ninoseas(2::4))
  ninodjf = dim_avg_Wrap(ninoseas(3::4))

  ninoseas(0::4) = ninoseas(0::4) - ninomam
  ninoseas(1::4) = ninoseas(1::4) - ninojja
  ninoseas(2::4) = ninoseas(2::4) - ninoson
  ninoseas(3::4) = ninoseas(3::4) - ninodjf
这个是先对数据进行了一个处理？请问一下在做小波分析前需要对数据进行什么处理呢？我是想通过小波分析看数据存在多少年的周期性
谢谢谢谢

liyf

yeah... 发表于 2017-5-4 10:56
你好 我想请教下 我在看官网的例子的时候发现
ninoseas = asciiread("sst_nino3.dat",-1,"float")
  n ...

你可以将数据直接放进程序里面进行小波分析，之前不用进行处理

流浪的心

{:eb348:}

林泓锦

{:lxm_23:}