%读取影像头文件和影像文件
% function data=File_read()
%%%%%----------------------------------------------------------------------------------------
%读取ENVI处理的影像头文件
[File_Name, File_Path] = uigetfile('.hdr', 'Select image headerfile to open');
if (File_Name == 0)
File_Name = [];
File_Path = [];
return
else
Fid_header = fopen([File_Path File_Name],'r');
end
elements={'samples ' 'lines ' 'data type ' 'byte order ' 'header offset ' 'bands '};
%对头文件进行分割
%line 影像的行数
%col 影像的列数
%data_type 影像的数据类型
%byte_order 影像的字节顺序，0 为leee-le,1 为leee-be
%datatype_size 每个象元的字节数
%headeroffset 影像的头文件信息字节数
while 1
fline=fgetl(Fid_header);
if ~ischar(fline)
break;
end
[first,second]=strtok(fline,'=');
switch first
case elements(1)
[f,num]=strtok(second);
col=str2num(num);
case elements(2)
[f,num]=strtok(second);
line=str2num(num);
case elements(3)
[f,num]=strtok(second);
data_type=str2num(num);
switch data_type
case 1
datatype_format='uchar';
datatype_size=1;
case 2
datatype_format='int16';
datatype_size=2;
case 3
datatype_format='int32';
datatype_size=2;
case 4
datatype_format='float';
datatype_size=4;
case 5
datatype_format='double';
datatype_size=8;
case 12
datatype_format='uint16';
datatype_size=2;
case 13
datatype_format='uint32';
datatype_size=4;
case 6
datatype_format='float'; %浮点复数数据
datatype_size=4;
end
case elements(4)
[f,num]=strtok(second);
byte_order=str2num(num);
if byte_order==0
byteorder='ieee-le';
else
byteorder='ieee-be';
end
case elements(5)
[f,num]=strtok(second);
headeroffset=str2num(num);
case elements(6)
[f,num]=strtok(second);
bands=str2num(num);
end
end
fclose(Fid_header);
%%%%%-----------------------------------------------------------------------------------
%读取影像文件
[File_Name, File_Path] = uigetfile('.\*.*', 'Select image file to open');
if (File_Name == 0)
File_Name = [];
File_Path = [];
return
end
headeroffset=0;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% if data_type~=6
%读图像
data=multibandread([File_Path File_Name],[line col bands],datatype_format, 0, 'bsq', 'ieee-le');
% for b=1:bands
% %每一波段读数起始位置
% Fid_imag = fopen([File_Path File_Name],'r');
% b_header=(b-1)*col*datatype_size;
% for i=1:line
% pointer=b_header+headeroffset+(i-1)*col*datatype_size; %读取每一行时的读数起始位置
% fseek(Fid_imag,pointer,'bof');
% data(i,:,b)=fread(Fid_imag,col,datatype_format,'ieee-le');%byteorder);
% end
% fclose(Fid_imag);
% end
% else
% %读复图像
% data=zeros(line,col);
% for i=1:line
% pointer=headeroffset+(i-1)*col*2*datatype_size;
% fseek(Fid_imag,pointer,'bof');
% data_temp=fread(Fid_imag,col*2,datatype_format,byteorder);
% slc_real=data_temp(1:2:(col*2));
% slc_imag=data_temp(2:2:(col*2));
% data(i,:)=complex(slc_real,slc_imag);
% end
% clear data_temp
% clear slc_real slc_imag
% fclose(Fid_imag);
% end
评论0