EOF分解平均值和距平值的结果一样

刘勇的下铺

这两天做EOF分解发现一个问题。先用原始数据做了一次EOF分解，再用dim_rmvmean_n_Wrap这个函数求数据的距平值再做一次EOF分解，两次的结果一样，有点不理解怎么会出现这样的结果。想问下有人遇到过这样的情况吗？
下面是处理数据和EOF的NCL脚本
f0 = addfile("precip.mon.mean.2.5x2.5.nc","r")
    vname = getfilevarnames(f0)
  print(vname)
  ;exit

  yrStrt = 1958
  yrLast = 2013
  TIME   = f0->time
  YYYY   = cd_calendar(TIME,-1)/100                 ; entire file
  iYYYY  = ind(YYYY.ge.yrStrt .and. YYYY.le.yrLast)
  prep    = f0->precip(iYYYY,:,:)
  printVarSummary(prep)
   
;prep = dim_rmvmean_n_Wrap(prep,0)
Prep = month_to_season(prep(:,{0:60},{100:160}),"JJA")
nyrs   = dimsizes(Prep&time)

  w = sqrt(cos(0.01745329*Prep&lat))
  wp = Prep*conform(Prep, w, 1)
  copy_VarCoords(Prep, wp)
  x = wp(lat|:,lon|:,time|:)

  optEOF = True      
  optEOF@jopt = 0  
  eof1 = eofunc_Wrap(x, 2, optEOF)
  eof = eof1*10
  copy_VarCoords(eof1, eof)
  ;optEOF=True
  ;optEOF@jopt = 1
  sig_pcv = eofunc_north(eof@pcvar, nyrs, False)
  eof_ts = eofunc_ts_Wrap(x, eof, False)
  eof_ts = dim_standardize_n(eof_ts, 1, 1)

刘勇的下铺

言深深 发表于 2018-11-14 22:45
你怎么做eof的？

eofunc_Wrap用ncl的这个函数。昨天做过的同学告诉我ncl函数会自动把元数据计算距平再做eof分解，如果元数据是距平就直接做分解。

刘勇的下铺

有人遇到过这个问题吗

刘勇的下铺

@Aires @言深深 麻烦两位大大帮忙看看，谢谢啦

言深深

你怎么做eof的？

言深深

刘勇的下铺 发表于 2018-11-15 10:59
eofunc_Wrap用ncl的这个函数。昨天做过的同学告诉我ncl函数会自动把元数据计算距平再做eof分解，如果元数 ...

是的，据我所知也是，ncl的函数默认是先距平，再分解的

言深深

刘勇的下铺 发表于 2018-11-15 10:59
eofunc_Wrap用ncl的这个函数。昨天做过的同学告诉我ncl函数会自动把元数据计算距平再做eof分解，如果元数 ...

是的，据我所知也是，ncl的函数默认是先距平，再分解的

笨拙的土豆

我用站点数据做eof，距平和原数据的结果是一样的

刘勇的下铺

言深深 发表于 2018-11-16 13:16
是的，据我所知也是，ncl的函数默认是先距平，再分解的

那应该是这样了。还想问下版主有没有做过矢量的EOF分解，我现在分别对U,V做EOF分解，然后把分解的结果用画流线的方法进行合成，不知道对不对。

刘勇的下铺

笨拙的土豆 发表于 2018-11-17 15:12
我用站点数据做eof，距平和原数据的结果是一样的

用NCL的EOF函数，得到的结果是一样的。如果是距平数据函数会直接进行分解，是原始数据函数会先计算距平在进行分解。