快捷導(dǎo)航

C/C++下讀取ENVI柵格文件格式的示例代碼

更新時(shí)間：2024年10月22日 11:37:52 作者：ENVI-IDL技術(shù)殿堂

ENVI使用的是通用柵格數(shù)據(jù)格式,包含一個(gè)簡(jiǎn)單的二進(jìn)制文件（ a simple flat binary ）和一個(gè)相關(guān)的ASCII（文本）的頭文件,下面我們就來(lái)看看如何使用C++讀取ENVI柵格文件格式吧

ENVI使用的是通用柵格數(shù)據(jù)格式，包含一個(gè)簡(jiǎn)單的二進(jìn)制文件（ a simple flat binary ）和一個(gè)相關(guān)的ASCII（文本）的頭文件。

利用其他語(yǔ)言如C/C++等直接讀取ENVI的數(shù)據(jù)，則可以先對(duì)hdr文件進(jìn)行解析，獲取數(shù)據(jù)類型。

hdr的文件結(jié)構(gòu)如下

ENVI
description = {
  Canon City, Colorado, Landsat TM, Calibrated to Reflectance }
samples = 640
lines   = 400
bands   = 6
header offset = 0
file type = ENVI Standard
data type = 1
interleave = bsq
sensor type = Landsat TM
wavelength units = Micrometers
z plot range = {0.00, 100.00}
z plot titles = {Wavelength, Reflectance}
band names = {
 TM Band 1, TM Band 2, TM Band 3, TM Band 4, TM Band 5, TM Band 7}
wavelength = {
 0.48500, 0.56000, 0.66000, 0.83000, 1.65000, 2.21500}

解析的關(guān)鍵信息有samples：640（列），lines：400（行），header offset：0（頭信息偏移量-單位為字節(jié)），data type=1（數(shù)據(jù)類型代碼，見下表）。

數(shù)據(jù)類型	代碼
字節(jié)型	1
16位有符號(hào)整型	2
32位有符號(hào)長(zhǎng)整型	3
32位無(wú)符號(hào)長(zhǎng)整型	13
浮點(diǎn)型	4
雙精度浮點(diǎn)型	5

對(duì)常用數(shù)據(jù)類型文件進(jìn)行了讀寫的測(cè)試，值完全一致。

利用IDL進(jìn)行文件寫出：

/*
 C++讀取ENVI格式技術(shù)測(cè)試代碼
 輸出不同數(shù)據(jù)類型的二進(jìn)制文件
 
 Author: DYQ 2011年6月2日
 
  BBS： http://bbs.esrichina-bj.cn/ESRI/forum-28-1.html
  E-Mail: dongyq@esrichina-bj.cn
  Blog: http://hi.baidu.com/dyqwrp
 */
PRO test_out_bin
  outdir = 'c:\temp\'
  if file_test(outdir,/directory) ne 1 then file_mkdir,outdir
  //字節(jié)byte
 OPENW,lun,outdir+'a.dat',/get_lun
 WRITEU,lun,BINDGEN(10)
  FREE_LUN,lun
  //整型int
 OPENW,lun,outdir+'b.dat',/get_lun
 WRITEU,lun,INDGEN(10)
  FREE_LUN,lun
  //浮點(diǎn)float
 OPENW,lun,outdir+'c.dat',/get_lun
 WRITEU,lun,FINDGEN(10)
  FREE_LUN,lun
  //長(zhǎng)整型long
 OPENW,lun,outdir+'d.dat',/get_lun
 WRITEU,lun,LINDGEN(10)
  FREE_LUN,lun
  //雙精度double
 OPENW,lun,outdir+'e.dat',/get_lun
 WRITEU,lun,DINDGEN(10)
  FREE_LUN,lun
END

C++下讀取文件：

//C++讀取ENVI格式技術(shù)測(cè)試代碼
#include "stdafx.h"
#include "iostream.h"
int main(int argc, char* argv[])
{
      printf("Hello ! Successful Using C++! ^_^ \n");
      int i,n;
      FILE*fp;
      //二進(jìn)制字節(jié)型
      char *bdata=new char[10];
      fp=fopen("c:\\temp\\a.dat","rb");
      n=fread(bdata,1,10,fp);
      fclose(fp);
      for(i=0;i<10;i++)
      {
             cout<<"二進(jìn)制";
             cout<<i<<":"<<short(bdata[i])<<endl;
      }
   //二進(jìn)制整型文件
      short int *idata=new short int[10];
      fp=fopen("c:\\temp\\b.dat","rb");
      n=fread(idata,2,10,fp);
      fclose(fp);
      for(i=0;i<10;i++)
      {
             cout<<"整型";
             cout<<i<<":"<<idata[i]<<endl;
      }
      //二進(jìn)制浮點(diǎn)文件
      float *fdata=new float[10];
      fp=fopen("c:\\temp\\c.dat","rb");
      n=fread(fdata,4,10,fp);
      fclose(fp);
      for(i=0;i<10;i++)
      {
             cout<<"浮點(diǎn)";
             cout<<i<<":"<<fdata[i]<<endl;
      }
      //二進(jìn)制長(zhǎng)整型文件
      long *ldata=new long[10];
      fp=fopen("c:\\temp\\d.dat","rb");
      n=fread(ldata,4,10,fp);
      fclose(fp);
      for(i=0;i<10;i++)
      {
             cout<<"長(zhǎng)整型";
             cout<<i<<":"<<ldata[i]<<endl;
      }
   //雙精度double
      double *ddata=new double[10];
      fp=fopen("c:\\temp\\e.dat","rb");
      n=fread(ddata,8,10,fp);
      fclose(fp);
      for(i=0;i<10;i++)
      {
             cout<<"雙精度型";
             cout<<i<<":"<<ddata[i]<<endl;
      }
      return 0;
}

最后輸出：