Fork me on GitHub

大话IDL编程之函数功能调用(envi_doit、ENVIRaster、ENVITask)

  2020年2月1日。好长时间没更新博客,还真有点不习惯。受新型冠性病毒的影响,平时街上熙熙攘攘的人流了无踪影,2020的春节竟然来的如此冷清。为响应“呆在家里就是做贡献的号召”,在家一宅就是十多天。闲来无事,就总结一下ENVI-IDL二次开发的学习经验吧。

  ENVI二次开发主要有两种编程方式:IDL二次开发、ENVI-IDL交互式开发。

  其中,IDL二次开发方式较为繁琐,但是也的确算是基本功,需要熟练掌握。而ENVI-IDL交互开发的方式简洁高效,是现在的主流开发方式。

  其中,对于某些功能函数,ENVI已经给出了新的调用方法。即envi_doit、ENVIRaster虚拟栅格、ENVITask这三种方法。其中虚拟栅格的方法将数据在内存中运行,运行完

之后释放。而ENVITASK也同样简洁高效、易懂。

  •   首先讲ENVITask的实现方式(这里以ENVIRadiometricCalibration为例)
    e = ENVI()
    file = 'a.dat'
    raster = e.OpenRaster(file)
    Task = ENVITask('RadiometricCalibration')
    Task.INPUT_RASTER = raster
    Task.OUTPUT_TYPE = 'Double'
    Task.OUTPUT_RASTER_UI = 'C:/a.dat'    ;第一种保存方法,设置输出路径,保存数据
    Task.Execute

    ;第二种ENVITask的文件保存方式
    outputRaster = TASK.OUTPUT_RASTER
    outputRaster.export, 'C:/a.dat', 'envi'

    一般来说,ENVITask方法只要设置好该方法的输入输出,就可以得到结果,较为方便。此外,对于自定义task, 在有.pro文件后,可以编写.task文件用来调用。

  • ui = e.UI   r = ui.SelectTaskParameters,添加该代码显示参数窗口。

  •  
  •      其次讲ENVIRaster的实现方式,即虚拟栅格(ENVIRadiometricCalibrationTask)
    e = ENVI()
    file = 'LC80410302013213LGN00_MTL.txt'
    raster = e.OpenRaster(file)
    OLIBands = raster[0]    ;Landsat8数据被存储在5-element数组中,第一个数组阵列是波段数据
    refRaster = ENVICalibrateRaster(OLIBands,CALIBRATION = 'Top-of-Atmophere Reflectance')
    refRaster.export,'c://new.dat' ;data.export输出文件并保存 e.DATA.add, refRaster ;数据添加到交互界面 view
    = e.getview() ;创建视图 layer = view.createlayer(refRaster) ;添加图层
  •     最后是 envi_doit方法/IDL方法
    fn = 'a.dat'
    envi_open_file, fn, r_fid=fid
    envi_file_query, fid, ns=ns, nl=nl, nb=nb, dims=dims,pos=pos
    map_info = envi_get_map_info(fid=fid)
    
    ;接下来有两种方法, envi_doit和普通的IDL方法
    pos=[0,1,2]
    envi_doit,'envi_avhrr_calibrate_doit',dims=dims,fid=fid,pos=pos,out_name='b.dat'
    
    ;普通的IDL方法
    fn_calib = dialog_pickfile(filter='*.txt', title='calibrate')
    openr,lun,fn,fn_calib,/get_lun   ;openr 打开txt文件
    data=fltarr(2,6)
    readf,lun,data  ;readf 读取到data
    free_lun, lun
    gain = data[0,*]
    bais = data[1,*]
    
    L = fltarr(ns,nl,nb)
    for i=0,nb-1 do begin
        data_band=envi_get_data(fid=fid,dims=dims,pos=i)   
        ;envi_get_data 读取数据
        L[*,*,i]=gain[i]*data_band+bias[i]
    endfor
    
    outfile='new.dat'
    envi_write_envi_file,L,out_name=outfile,ns=ns,nl=nl,nb=nb,$ data_type=4,interleave=interleave,offset=offset,bnames=bnames,map_info=map_info
    #保存文件 envi_write_envi_file

     如有错误,还请指摘。

posted @ 2020-02-01 18:41  Rser_ljw  阅读(4526)  评论(0编辑  收藏  举报