File_read.rar_MATLAB envi_matlab读取envi_读取envi_读取envi格式

%读取影像头文件和影像文件 % 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