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之前使用基本语句读取到上述四个值，得到一个三位矩阵（表示高光谱数据尺寸）和一个数据类型。
clcclearfname='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，=后为secondd={'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);endendt='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里显示和修改。
