IDL读取ENVI标准文件源码

; PURPOSE:
; IDL program, which reads standard ENVI image files (*.img).
;
; CALLING SEQUENCE:
; read_envi_file, infile, img, xs, ys, type,offset
;
; INPUTS:
; infile - 传入数据文件名
;
; OPTIONAL INPUTS:
; None
;
; KEYWORD PARAMETERS:
; None
;
; OUTPUTS:--这些均是返回值，注意是位置参数
; img - ENVI的图像文件；
; xs - 列号；
; ys - 行号；
; type - 数据类型代码
; offset - 头文件偏移量；
; mapinfo - map及地理坐标信息；
;
; EXAMPLE:read_envi_image, 'C:Program FilesITTIDL708productsenvi46dataddd', data;;;data即读取的数据值
;
; MODIFICATION HISTORY:
; Written by: Carsten Pathe, cp@ipf.tuwien.ac.at
; Date: 25.08.2003
;
; Modified By DYQ
; 2009-8-18 修正：数据无扩展名会读取错误
; 简化数组创建代码
; 2009-12-22 修正：添加BIP、BIL多波段数据支持
; 
; SUXIN,2020.6.11
; 修改头文件读取代码
; 头文件是否存在判断
; Mapinfo支持两个方向的分辨率
;-
PRO Read_envi_image, infile, img, bandsnum,xs, ys, type, offset, mapinfo

  COMPILE_OPT IDL2

  image = infile
  
  header = file_dirname(image)+PATH_SEP()+FILE_BASENAME(image,'.dat')+'.hdr'
  
  if ~file_test(header) then begin
    d = DIALOG_MESSAGE('头文件不存在:'+header,/ERROR)
    return 
  endif

  ;打开解析头文件
  OPENR, unit, header, /get_lun

  header_line = ''
  ;文件未读取完之前一直循环
  WHILE NOT EOF(unit) DO BEGIN

    READF, unit, header_line
    tmp = STRSPLIT(header_line[0], '=', /extract)
    header_keyword = STRSPLIT(tmp[0], ' ', /extract)
    ;解析头文件中的信息
    IF header_keyword[0] EQ 'samples' THEN xs = LONG(tmp[1])
    IF header_keyword[0] EQ 'lines' THEN ys = LONG(tmp[1])
    IF header_keyword[0] EQ 'header' THEN offset = LONG(tmp[1])
    IF header_keyword[0] EQ 'bands' THEN bandsnum = LONG(tmp[1])
    IF header_keyword[0] EQ 'interleave' THEN dataInter = STRUPCASE(STRTRIM(tmp[1],2))
    IF header_keyword[0] EQ 'data' THEN type = LONG(tmp[1])
    ;如读取map信息则解析map信息
    IF header_keyword[0] EQ 'map' THEN BEGIN

      mapinfo_tmp=STRSPLIT(tmp[1],'{',/extract)
      mapinfo_tmp=STRSPLIT(mapinfo_tmp[1],',',/extract)

      mapinfo={ULX:0.,ULY:0.,XResolution:0.,YResolution:0.}
      mapinfo.ULX=mapinfo_tmp[3]
      mapinfo.ULY=mapinfo_tmp[4]
      mapinfo.XResolution=mapinfo_tmp[5]
      mapinfo.YResolution=mapinfo_tmp[6]
    ENDIF
  ENDWHILE

  ;关闭头文件
  CLOSE,unit & FREE_LUN, unit

  ;modified by dyq

  ;打开数据文件定位到数据位置
  OPENR, unit,IMAGE, /get_lun
  POINT_LUN, unit, offset

  ;判断数据存储类型
  ;完全按照ENVI的数据存储格式定义
  CASE dataInter OF
    'BSQ': BEGIN

      img = MAKE_ARRAY(xs,ys,bandsnum,type = type)
      ;读取数据文件中的数据
      READU,unit,img
    END
    ;按行保存
    'BIL': BEGIN

      img = MAKE_ARRAY(xs,bandsnum, ys,type = type)
      ;读取数据文件中的数据
      READU,unit,img
      ;转换为BSQ
      img= TRANSPOSE(img,[0,2,1])
    END
    'BIP': BEGIN
      ; IF bandsnum GT 1 THENimg = MAKE_ARRAY(bandsnum,xs,ys,type = type) $
      ; ELSE img = MAKE_ARRAY(xs,ys,type = type)
      img = MAKE_ARRAY(bandsnum,xs,ys,type = type)
      ;读取数据文件中的数据
      READU,unit,img
      ;转换为BSQ
      img= TRANSPOSE(img,[1,2,0])
    END
    ELSE:
  ENDCASE
  FREE_LUN, unit
  img = REFORM(img)
end

喜欢 (1)赏

支付宝[853060844@qq.com]

分享 (0)