matlab读取高光谱影像
1.函数multibandread读取读取多波段二进制影像文件(ENVI主菜单file—save file as—envi standard得到的就是二进制影像文件，有时甚至会看到后缀名为bsq、bil、bip等影像)。
im_hyper = multibandread(filename, size, precision, offset, interleave, byteorder)
后附参考代码有实例
precision=\’uint8=>uint8\’;%头文件中datatype=1对应ENVI中数据类型为Byte，对应MATLAB中数据类型为uint8
precision=\’int16=>int16\’;%头文件中datatype=2对应ENVI中数据类型为Integer，对应MATLAB中数据类型为int16
precision=\’uint16=>uint16\’;%头文件中datatype=12对应ENVI中数据类型为Unsighed Int，对应MATLAB中数据类型为uint16
precision=\’int32=>int32\’;%头文件中datatype=3对应ENVI中数据类型为Long Integer，对应MATLAB中数据类型为int32
precision=\’uint32=>uint32\’;%头文件中datatype=13对应ENVI中数据类型为Unsighed Long，对应MATLAB中数据类型为uint32
precision=\’float32=>float32\’;%头文件中datatype=4对应ENVI中数据类型为Floating Point，对应MATLAB中数据类型为float32
precision=\’double=>double\’;%头文件中datatype=5对应ENVI中数据类型为Double Precision，对应MATLAB中数据类型为double
）
function data=read_ENVIimagefile(imgfilename) %本函数读取img格式，前提是img图像显式带有\'.img\'后缀名。 if length(imgfilename)>=4 switch strcmp(imgfilename(length(imgfilename)-3:end), \'.img\') case 0 hdrfilename=strcat(imgfilename, \'.hdr\'); case 1 hdrfilename=strcat(imgfilename(1: (length(imgfilename)-4)), \'.hdr\'); end else hdrfilename=strcat(imgfilename, \'.hdr\'); end %读取ENVI标准格式图像文件 %读取图像头文件 fid = fopen(hdrfilename, \'r\'); info = fread(fid,\'char=>char\'); info=info\';%默认读入列向量，须要转置为行向量才适于显示 fclose(fid); %查找列数 a=strfind(info,\'samples = \'); b=length(\'samples = \'); c=strfind(info,\'lines\'); samples=[]; for i=a+b:c-1 samples=[samples,info(i)]; end samples=str2num(samples); %查找行数 a=strfind(info,\'lines = \'); b=length(\'lines = \'); c=strfind(info,\'bands\'); lines=[]; for i=a+b:c-1 lines=[lines,info(i)]; end lines=str2num(lines); %查找波段数 a=strfind(info,\'bands = \'); b=length(\'bands = \'); c=strfind(info,\'header offset\'); bands=[]; for i=a+b:c-1 bands=[bands,info(i)]; end bands=str2num(bands); %查找数据类型 a=strfind(info,\'data type = \'); b=length(\'data type = \'); c=strfind(info,\'interleave\'); datatype=[]; for i=a+b:c-1 datatype=[datatype,info(i)]; end datatype=str2num(datatype); precision=[]; switch datatype case 1 precision=\'uint8=>uint8\';%头文件中datatype=1对应ENVI中数据类型为Byte，对应MATLAB中数据类型为uint8 case 2 precision=\'int16=>int16\';%头文件中datatype=2对应ENVI中数据类型为Integer，对应MATLAB中数据类型为int16 case 12 precision=\'uint16=>uint16\';%头文件中datatype=12对应ENVI中数据类型为Unsighed Int，对应MATLAB中数据类型为uint16 case 3 precision=\'int32=>int32\';%头文件中datatype=3对应ENVI中数据类型为Long Integer，对应MATLAB中数据类型为int32 case 13 precision=\'uint32=>uint32\';%头文件中datatype=13对应ENVI中数据类型为Unsighed Long，对应MATLAB中数据类型为uint32 case 4 precision=\'float32=>float32\';%头文件中datatype=4对应ENVI中数据类型为Floating Point，对应MATLAB中数据类型为float32 case 5 precision=\'double=>double\';%头文件中datatype=5对应ENVI中数据类型为Double Precision，对应MATLAB中数据类型为double otherwise error(\'invalid datatype\');%除以上几种常见数据类型之外的数据类型视为无效的数据类型 end %查找数据格式 a=strfind(info,\'interleave = \'); b=length(\'interleave = \'); c=strfind(info,\'sensor type\'); interleave=[]; for i=a+b:c-1 interleave=[interleave,info(i)]; end interleave=strtrim(interleave);%删除字符串中的空格 %读取图像文件 fid = fopen(imgfilename, \'r\'); data = multibandread(imgfilename ,[lines, samples, bands],precision,0,interleave,\'ieee-le\'); data= double(data); end