基于matlab读取envi格式高光谱图像数据

1关于光谱：

电磁波按波长大小分为射线，可见光，近红外（780nm-2526nm），无线电等。

波长：一个振动周期内传播的距离。

什么是光谱：复色光色散后按波长大小排列的图案。

什么是高光谱：1光谱信息丰富的光谱信息2光谱分辨率足够高。

什么叫光谱信息丰富：区分光波段范围足够细。

光谱分辨率足够高：1在一固定波长内波段越多的光谱，分辨率越高。2波段越小，分辨率越高。

什么是半高全宽：看分辨率有多高，波段有多宽。

高光谱遥感数据光谱特征的三层体系

高光谱遥感影像是一个空间-光谱数据立方体。

每一层为一个波段，每一个像元各波段的属性值都构成一个光谱向量

通过对光谱向量的分析，可以提取许多对分类、信息提取都用的信息

高光谱数据格式之一：标准envi格式是两个文件，一个HDR一个DAT。

2基于matlab简单读取函数

function [image]=freadhyperspectraldata(fname)fid=fopen(fname);%打开冰箱
image=fread(fid);%放进大象
fclose(fid);     %关上冰箱门
%%简单的读取函数，得到一个一维结果。1*161200000.

将高光谱文件每一个数值进行读取，输出一个一维数组。

3按头文件，将上述数组reshape

我们无法处理上述一维数据，所以需要把高光谱数据变成三维矩阵，以便矩阵运算。

我们需要处理矩阵行数，列数，通道数（波段数），数据类型，这四个值在头文件有标明。

所以reshape之前使用基本语句读取到上述四个值，得到一个三位矩阵（表示高光谱数据尺寸）和一个数据类型。

clc
clear
fname='G:\轻工大学\1高光谱\高光谱数据\T1\1\capture\1-1'
rfid = fopen(strcat(fname,'.hdr'),'r'); %打开头文件%strcat是连接字符串，rfid是指定头文件
while 1tline = fgetl(rfid);%fget是读取每一行并去掉换行符if ~ischar(tline), break, end[first,second]=strtok(tline,'=');%将=作为分割节点，%=前为first，=后为second
d={'bit8' 'int16' 'int32' 'float32' 'float64' 'uint16' 'uint32' 'int64' 'uint64'};switch first%获取头文件中samples，lines，bands的值，即行数、列数、波段数信息case 'samples '[f,s]=strtok(second);p(1)=str2num(s);first% a=1case 'lines '[f,s]=strtok(second);p(2)=str2num(s);case 'bands '[f,s]=strtok(second);p(3)=str2num(s);endend
t='uint16'%本次实验为了加快进度未设置如何搜索数据类型，直接赋值
fname='G:\轻工大学\1高光谱\高光谱数据\T1\1\capture\1-1.dat'
fid=fopen(fname);
image=fread(fid,t);
image=reshape(image,[p(1),p(2),p(3)]);
fclose(fid);

运行结果：image：1024*948*224

p指光谱数据矩阵尺寸，image指reshape后的光谱数据已经被读取并转换为三维矩阵。

ps：envi标准高光谱数据格式分为两个文件：1头文件 hdr文件2数据文件dat（有一点不懂为什么记事本打开dat文件全是乱码而不是数字矩阵，按头文件说的应该uint16。如果不是数字那么matlab里的image三维矩阵又是什么呢。）

今天学习了元数据的读取和数据整理，下一次我们可以看看这些数据如何在matlab里显示和修改。