IDL编程小技巧

1. IDL+ENVI

  compile_opt IDL2
  envi, /restore_base_save_files 
  ENVI_BATCH_INIT

 envi_open_data_file, file, r_fid=fid  
  envi_file_query, fid, ns=ns,nb=nb, nl=nl,dims=dims    
  d= ENVI_GET_DATA(fid=fid, dims=dims, pos=0)/AU
  free_lun,fid

2. 二进制文件的读取

openr,lun,datapath,/get_lun
    dc=fltarr(128,128)
    readu,lun,dc
free_lun,lun

3. ＡＳＣＩＩ文件的读取

  openr,lun,ffpath,/get_lun
    readf,lun,mff    
  free_lun,lun

4. ENVI文件的头文件的添加

 ENVI_SETUP_HEAD, fname=outname, $
        ns=128, nl=nl, nb=32, $
        interleave=1, data_type=4,wl=wl, $
        offset=0, $
        /write, /open

5. 图片生成与输出

    window,0, XSIZE=800, YSIZE=600
    !p.BACKGROUND=!p.color
    a=[0]
    !p.font=1
    device,set_font='微软雅黑*32'
    
    plot,a,color='000000'x,xtitle='Points',ytitle='Radiance', $
        xrange=[0,12],$
        yrange=[0,max([max(m3[i,1:np-1]),max(iim[i,1:np-1])])],$
        title='The Randiance of M3(Red) with IIM(Blue) at'+d2s(iim[i,0]),$
        ySTYLE=0,xthick=2,ythick=2,xstyle=1,$
        ;YTICKs=5,yTICKINTERVAL=5,$
        xTICKINTERVAL=5,xticks=4
    oplot,m3[i,1:np-1],color='0000ff'x,thick=3
    oplot,m3[i,1:np-1],color='0000ff'x,psym=4, SYMSIZE=0.7,thick=6
    oplot,iim[i,1:np-1],color='ff0000'x,thick=3
    oplot,iim[i,1:np-1],color='ff0000'x,psym=7, SYMSIZE=0.7,thick=6
    WRITE_png, path+d2s(iim[i,0])+'.png', TVRD(/TRUE,/order),/order
    wdelete,0

6. 二进制文件的输出

  openw,lun,output,/get_lun
    writeu,lun,o
  free_lun,lun

7. 文件搜索

files=FILE_SEARCH(m3path+'*-m3*.txt',count=nfile)