IDL 打开多波段图像

转自：http://blog.csdn.net/rrrrssss00/article/details/4843841

方法1：

pro testEnviImg,path

;  compile_opt idl2

;  envi,/restore_base_save_files

;  envi_batch_init

  path = ENVI_PICKFILE(title='pick file')

  envi_open_file, path, r_fid=fid

 

  if (fid eq -1) then return

  ;ENVI_SELECT, fid=fid ,pos=pos

  envi_file_query, fid, dims=dims, nb=nb

  pos = lindgen(nb)

  num_cols = dims[2]-dims[1]+1

  num_rows = dims[4]-dims[3]+1

  image = fltarr(nb,num_cols,num_rows)

  for i=0,nb-1 do image[i,*,*]=$  

  envi_get_data(fid=fid,dims=dims,pos=pos[i])

  tv,image,/true

end

方法2：

pro testenviselect

  envi_select,fid=fid, dims=dims,pos=pos

 

  num_cols = dims[2]-dims[1]+1

  num_rows = dims[4]-dims[3]+1

  num_bands = n_elements(pos)

 

  ; BIP array

  image = fltarr(num_bands,num_cols,num_rows)

 

  for i=0,num_bands-1 do image[i,*,*] = envi_get_data(fid=fid,dims=dims,pos=pos[i])

  tv,image,/true

end