关于滤波

易小凯

subtropical 发表于 2016-3-9 18:20
非常感谢您的回复，我其实一直都有一个疑问就是怎么把原来改了顺序的变量参数再变回来，求大神指教。谢谢 ...

额(⊙o⊙)…
s3=wgt_runave(s2(lat|:,lon|:,time|:)
s4=s3(time|:,lat|:,lon|:)就Ok了啊

subtropical

易小凯 发表于 2016-3-10 11:03
额(⊙o⊙)…
s3=wgt_runave(s2(lat|:,lon|:,time|:)
s4=s3(time|:,lat|:,lon|:)就Ok了啊

嗯我后来也这样写了，但是出不来图，麻烦您帮忙看一下



只出来一个图，我的脚本哪里有问题吗谢谢您

subtropical

易小凯 发表于 2016-3-10 11:03
额(⊙o⊙)…
s3=wgt_runave(s2(lat|:,lon|:,time|:)
s4=s3(time|:,lat|:,lon|:)就Ok了啊

是不是与设置缺省值有关哪？我看官网上

但是我写了
s3=wgt_runave_Wrap(s2(lat|:,lon|:,time|:),wgts,0)
  ;s3 =wgt_runave_leftdim(s2,wgts,0)   ; 10 year
; copy_VarMeta(s2,s3)
s3@_FillValue  =  -999.0
还是出现这些 请大神帮忙看一下，非常感谢

易小凯

subtropical 发表于 2016-3-10 16:03
嗯我后来也这样写了，但是出不来图，麻烦您帮忙看一下

看了下你的代码，不知道你s2的时间长度有多少？
因为你用来120个月来做低通滤波，这说明在滤波完后s3的前面60个月和最后60个月都是缺省
你可以检查下，是不是这个问题。

subtropical

易小凯 发表于 2016-3-11 07:37
看了下你的代码，不知道你s2的时间长度有多少？
因为你用来120个月来做低通滤波，这说明在滤波完后s3的 ...

哦哦，我的s2是130年的海温去过线性趋势后的值，我想研究年际变化，所以我用原来的场s2把低频波s3去掉,之后再把我需要的年份画出来的,那为什么低通滤波之后s3的前面60个月和最后60个月都是缺省呢？非常感谢您的指导，我想应该就是这个原因，麻烦您您再教一下该怎么设置才能正确出图？非常非常感谢您的回复{:eb511:}

易小凯

subtropical 发表于 2016-3-11 14:25
哦哦，我的s2是130年的海温去过线性趋势后的值，我想研究年际变化，所以我用原来的场s2把低频波s3去掉,之 ...

滤波的话就会造成信息缺失的问题，这是无法避免的，否则你只能通过增加数据量来达到。
不知道你需要滤除多少的波，一般31点的滤波就差不多了。
换言之，如果你需要1980-2000之间的高频信息，那么原始数据建议向前向后多15年，即1965-2015的数据，滤波得到1980-2000的高频信息。

subtropical

易小凯 发表于 2016-3-12 07:54
滤波的话就会造成信息缺失的问题，这是无法避免的，否则你只能通过增加数据量来达到。
不知道你需要滤除 ...

哦哦好的我懂您的意思了，非常感谢

zzli

thanks~

qq469015280

易小凯 发表于 2012-7-29 17:02
用低通滤波器提取出年代际变化，然后用原始值减去年代际变化
11a滑动平均是低通滤波器的一种，可以试着选 ...

请教大神，我用ncl的傅里叶滤波对5年(1825天)滤出8-20天的波
是应该对365天（每一年分开）滤吃8-20天波后再合起来，
还是可以直接对1825天（5年）这整时间序列滤？？？？？？？

易小凯

qq469015280 发表于 2016-9-1 14:52
请教大神，我用ncl的傅里叶滤波对5年(1825天)滤出8-20天的波
是应该对365天（每一年分开）滤吃8-20天波 ...

第二种方法