IDL时间序列图如何平滑?

飞天木脑壳

T T求大神相助..
  用IDL画出来的线画图是这样的

以为用smooth平滑以后是像参考文献里这样华丽的..

哪知道....T T..丑哭了...就像用左脚画出来的一样....求大神指导~~改怎么平滑或者拟合成漂亮的曲线呢?

步绝尘_南山居士

@爱听雨的青蛙来看看什么情况

飞天木脑壳

大神们能不能支援一下阿~~T T好苦阿~~

mali.

这个是要插值 不是平滑。
建议你用loop_interpol.pro没有的话的google一下
假如x=[0,1,2,3,4,5,6,7,8,9]
      y=[11,33,45,22,55,22,44,66,99,9]
      plot,x,y后效果很差曲曲折折的
插值后弄就好了
  loop_interpol,1,x,y,x1,y1（x1,y1 是插值后的数据）
plots,x1,y1
看下效果把

飞天木脑壳

mali. 发表于 2014-3-14 16:42

                               
登录/注册后可看大图

这个是要插值 不是平滑。
建议你用loop_interpol.pro没有的话的google一下
假如x=[0,1,2,3,4,5,6,7,8,9]
...

阿克塞pet.jpgo..(≧v≦)o~~谢谢呐~~真的平滑好咯..

飞天木脑壳

mali. 发表于 2014-3-14 16:42
这个是要插值 不是平滑。
建议你用loop_interpol.pro没有的话的google一下
假如x=[0,1,2,3,4,5,6,7,8,9]
...

你好,我想再麻烦问问为什么我用这个函数插值会出现这样的情况呐?
360截图20140416214856825. ...

mali.

飞天木脑壳 发表于 2014-4-16 21:51
你好,我想再麻烦问问为什么我用这个函数插值会出现这样的情况呐?
360截图20140416214856825. ...

你这是什么数据额  数据本身没问题么 不插值之前是什么效果  光贴图。。。。

飞天木脑壳

mali. 发表于 2014-4-17 20:35
你这是什么数据额  数据本身没问题么 不插值之前是什么效果  光贴图。。。。

恩恩..数据本身没问题的..插值之前是这样哒~

                               
登录/注册后可看大图

下面是我的程序..

pro petdrawing
openr,lun,'f:\data\After\choosen\cru.pet.txt',/get_lun
;pet0为读取原始月潜在蒸发量,pet为年潜在蒸发量,pet1为平滑后的潜在蒸发量
pet0=make_array(1,239,/float)
readf,lun,pet0
pet=fltarr(20)
x=indgen(20)+1974
for j=0,10 do begin
  pet[0]=pet[0]+pet0[j]
endfor
for i=1,19 do begin
  for j=0,11 do begin
   pet=pet0[(i*11+j)]+pet
  endfor  
endfor
free_lun,lun
loop_interpol,0.02,x,pet,x1,pet1
plot,x1,pet

end

mali.

飞天木脑壳 发表于 2014-4-19 11:34
恩恩..数据本身没问题的..插值之前是这样哒~下面是我的程序..

pro petdrawing

loop_interpol,0.02,x,pet,x1,pet1 loop_interpol源码中 这个dpix使用来计算插值的线上的像素点

飞天木脑壳

mali. 发表于 2014-4-19 15:14
loop_interpol,0.02,x,pet,x1,pet1 loop_interpol源码中 这个dpix使用来计算插值的线上的像素点个数ns的 ...

谢谢了,我原先太粗心只是大致看了看公式,看到貌似用的是样条插值法所以把这个值当成了步长了,而且我的X和Y都是浮点型数据,X内插了对应的Y值也内插了小数~但是即使我把dpix改成了大于1的值,就算到了10也是一样的绕起来的图呢~