IDL编程学习之一维中值滤波and中文图例

15195775117

数据见附件，只有一列PM10浓度数据，从上往下按时间排列读取数据：

  a=read_ascii('C:\Users\Administrator\Desktop\data.txt',data_start=0)
  pm=a.field1
8位中值滤波，位数太小，突兀的值去不掉，位数太大，滤波后的曲线与原始数据贴合度差
  re=median(pm,8)
画原始数据的曲线，name为图例名，是后面legend函数的目标：
  fig1=plot(pm,color='k',linestyle=6,sym_filled=1,$
    sym=2,/overplot,name='原始数据')
画滤波后数据的曲线，xShowText=0表示x轴数字和标题都不显示：
  fig2=plot(re,color='r',/overplot,thick=5,$
    ytitle='颗粒物浓度',xtitle='时间序列',name='中值滤波',xSHOWTEXT=0)
加载图例，字体选楷体（kaiti）可以显示中文，不会是乱码
      leg=legend(target=[fig1,fig2],font_name='KaiTi',position=[0.3,0.8])