IDL数据平滑的两种方法及其案例

15195775117

附件HUMID_LWC_R1.txt文件给出了两列数据，
第一列：比湿（g/kg）；第二列：含水量（mg/m3）
作图要求：画含水量与比湿的关系。
检查数据：比湿从0-0.0135866，都是小数，可转换以kg/kg为单位，即都乘1000，则比湿范围0-13.5866
以每0.1kg/kg比湿为一段统计含水量的均值：
a=read_ascii('C:\Users\Administrator\Desktop\HUMID_LWC_R1.txt',data_start=0)
x=a.field1
help,x
x[0,*]=x[0,*]*1000.0
b=[-99,-99]
for i=0.0,20.0,0.1 do begin
  p=where(x[0,*] gt i and x[0,*] lt i+0.1,count)
  if(count gt 10)then begin
    b=[,[i,mean(x[1,p])]]
    endif
endfor
help,b
b=b[*,1:-1]
help,b
fig=plot(b[1,*],b[0,*],color='r',sym=24,yrange=[0,14],yminor=0,xminor=1)
得：

该图波动较大，考虑使用滑动平均：
5步滑动平均：

10步滑动平均：

20步滑动平均：

滑动平均的效果是有效的，想要平滑还可以扩大比湿间隔：
以0.5kg/kg为间隔统计含水量均值：

以1kg/kg为间隔统计含水量均值：

总结：
平滑方法：1，用滑动平均的smooth函数；2，扩大统计间隔

meichh

太感谢了，正需要它

15195775117

meichh 发表于 2017-12-11 18:25
太感谢了，正需要它

荣幸！荣幸！荣幸！荣幸！