ENVI: Layer Stacking

最新推荐文章于 2023-11-03 17:13:44 发布

l_iu_xing

最新推荐文章于 2023-11-03 17:13:44 发布

阅读量5.7k

点赞数 3

分类专栏：遥感图像

本文链接：https://blog.csdn.net/l_iu_xing/article/details/84990939

版权

本文介绍如何使用 ENVI 的 envi_layer_stacking_doit 函数进行波段融合操作，探讨了采用面向对象编程重构代码的可能性。

摘要由CSDN通过智能技术生成

利用 envi_layer_stacking_doit 实现波段融合。可考虑用面向对象的编程方式去重新编码（暂时复制以前的代码）。

compile_opt idl2
  
  envi, /restore_base_save_files 
  
  start_path = 'F:\Temporary\data\LC8 Results\ENVI_clip\temp\ROI\'
  cd,start_path
  out_path = 'F:\Temporary\data\LC8 Results\ENVI_clip\temp\ROI\'
  
  Files = FILE_Search(start_path + 'L*roi.tif', count = Filecount)
  print, files
  If FileCount Eq 0 Then Begin
    print, 'File not exist'
    Return
  Endif
  
  test_file = FILE_Search(out_path + '*Layer_stacking*', count = test_count)
  if test_count ne 0 then begin
    file_delete, out_path +"Layer_stacking"
    file_delete, out_path +"Layer_stacking.hdr"
  endif
  
  nb = Filecount
  fid = lonarr(nb)
  pos = lonarr(nb)
  dims = lonarr(5,nb)
  out_bname = strarr(nb) ;为输出的每个波段命名
  
  For i = 0, Filecount - 1 do begin

    A_Name = Files[i]
    envi_open_file