c1驾照年审需要多少钱:IDL直接读取ENVI的img数据格式源码(终结版) - ENVI/IDL空间 - ALL...

来源：百度文库 编辑：偶看新闻 时间：2024/07/05 14:42:08

IDL直接读取ENVI的img数据格式源码(终结版)

上一篇 / 下一篇  2009-12-22 18:35:03 / 个人分类：IDL技术相关

查看( 55 ) / 评论( 4 ) [i=s] 本帖最后由 lianyi_dyq 于 2009-12-22 21:05 编辑

IDL直接解析ENVI的img格式源码终结版，在原来代码的基础上加入了多波段、多文件类型（BIP、BIL）的支持。
   原帖地址

CODE:

; 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 Files\ITT\IDL708\products\envi46\data\ddd', 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多波段数据支持
;-
PRO Read_envi_image, infile, img, bandsnum,xs, ys, type, offset, mapinfo
  COMPILE_OPT idl2
  
  image = infile
  ;Add by DYQ-for数据文件无扩展名的bug
  pointExist = STRPOS(infile,'.')
  IF pointExist[0] NE -1 THEN BEGIN
    header = Strsplit(infile,'.',/extract)
    header = header[N_ELEMENTS(header)-2]+'.hdr'
  END ELSE header = infile + '.hdr'
  
  ;打开解析头文件
  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.,spacing:0.}
      mapinfo.ulx=mapinfo_tmp[3]
      mapinfo.uly=mapinfo_tmp[4]
      mapinfo.spacing=mapinfo_tmp[5]
    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 THEN  img = 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)
  
img文件时干什么的，在dos下怎样读取？ 数据格式的问题！ 怎么用gis的软件ENVI？谁会？ ASP中数据格式的问题 我的delphi7中设计IDL时没有corba object 有谁有envi 或pci 的实用手册（电子版） envi软件 如何用excel进行数据格式的查询 请问SP2数据格式是什么？如何表示的？ 这是什么布线软件的数据格式？ 有没有可以直接把ISO文件直接转换成IMG文件的工具？ 如何直接读取QQ的聊天记录文件MsgEx.db？？ IMG是什么文件的 img格式的特点 浮点数的数据格式是：（ ），阶码，（ ），尾数． 磁带上的歌曲怎么录制成数据格式？ IDL接口和CORBA对象系统是什么来的????请详细点.谢谢 遥感图像处理软件ENVI的教程和资料谁有？谁会？ 我的光驱可以直接播发音乐 但是就是不能读取任何光盘 总是提示磁盘损坏切无法读取 IDL是什么东西？ img是什么意思啊~~论坛的IMG开关指什么啊~? .img文件的打开方式? 有关 img 代码的问题!! img文件的打开方法