matlab讀取高光譜影像

 

1.函數multibandread讀取讀取多波段二進制影像文件(ENVI主菜單file—save file as—envi standard得到的就是二進制影像文件，有時甚至會看到后綴名為bsq、bil、bip等影像)。

 im_hyper = multibandread(filename, size, precision, offset, interleave, byteorder)

                                                后附參考代碼有實例

       filename： 文件名

       size：圖像尺寸和波段數，size = [ 行數 列數 波段數 ]

       precision：讀取的圖像的數據格式，例如'uint8','uint16','double'等

       offset：偏移（一般為0）

       interleave：存儲的圖像的數據格式，有 bsq,bil,bip三種格式

       byteorder : 數據存儲的字節排列方式，有'ieee-le'(小端)，'ieee-be'(大端)

 

注：precision參數與matlab數據類型相應的關系（FORM  https://malagis.com/use-matlab-read-envi-image.html，辛苦博主整理。）

    （

　　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

）

 

2. 如果是.tif格式的圖像，可先使用 envi將其轉化為envi standard 格式，然后使用enviread讀取

ENVI軟件標准格式的高光譜遙感圖像含有兩部分：

　　一個是高光譜圖像 '*.img '，

　　另一個是遙感圖像頭文件 '*.hdr'，該文件記錄了遙感圖像的信息，如圖像尺寸、波段數、數據類型和大小端等。如果缺少頭文件，將無法對遙感圖像進行讀取

讀取代碼：來源： https://malagis.com/use-matlab-read-envi-image.html， 亦可參考：http://blog.sina.com.cn/s/blog_16ecb33400102ybql.html

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