IDL编程学习之统计实例

15195775117

附件gaga.rar中包括了5月12-16日5个文件夹，每个文件夹中有SO2和NO2在不同时刻的txt文件，每个文件分为3列，第一列是方位角，第二列是半径（我也不知道是什么半径），第三列是浓度。
现在想要得到SO2（和NO2）在各个方位的平均分布。
观察数据，发现方位角都是30,60,90......330,360，浓度值也有9999这样的异常值，所以在统计时需要剔除。
代码：
book=[-99,-99,-99];用来存放合并结果
for date=2,6 do begin
file=findfile('C:\Users\YGDY\Desktop\gaga\051'+strtrim(string(date),2)+'\x_*SO2*.txt');先计算SO2的，然后把S改成N即可
n_file=n_elements(file)
for i=0,n_file-1 do begin
  a=read_ascii(file[i],data_start=1);第一行不要
x=a.field1
book=[[book],[x]];合并
endfor
endfor
help,book
for i=30,360,30 do begin
  p=where(book[0,*] eq i,count);筛选出每个方位角的数值
  if(count gt 0)then begin
    b=x[1,p]
    p2=where(b gt 0,count)
    if(count gt 0)then begin
      r=mean(b[p2])
      endif
    c=x[2,p]
    p3=where( c gt 0 and c ne 9999,count);大于0也就舍去了第一行的-99,9999异常值也舍去
    if(count gt 0)then begin
      s=mean(c[p3])
    endif
    print,i,r,s;输出某方位角下的平均半径和浓度
    endif
endfor
SO2和NO2的结果见附件