matlab读取并显示遥感img图像，包含单波段显示与多波段显示

%freadenvi('06newmul01') function data = freadenvi(imgfilename) %freadenvi(imgfilename) %本函数读取img格式，前提是img图像显式带有'.img'后缀名 % if length(imgfilename)>=4 %freadenvi('06newmul01') % 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 hdrfilename = strcat(imgfilename, '.hdr'); imgfilename = strcat(imgfilename,'.img');%'.img' dat %读取ENVI标准格式图像文件 %读取图像头文件 fid = fopen(hdrfilename, 'r'); info = fread(fid,'char=>char');%从二进制文件中读取数据 fread(fid,1,'char=>char') info = deblank(info); info=info';%默认读入列向量，须要转置为行向量才适于显示 if fid == -1 error('Input header file does not exist'); end %r to display fclose(fid); %查找列数 a=strfind(info,'samples = '); b=length('samples = '); c=strfind(info,'lines'); samples=[]; for i=a+b:c-1 %读取hdr文件中samples的数值，除去换行和空行，所以c-3，或者c-1 samples=[samples,info(i)]; end %samples2 = strrep(samples, ' ', '') %删除字符串中间的空格 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');%对于多波段 %c=strfind(info,'data type');%对于单多波段; 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');%pp='invalid datatype';%除以上几种常见数据类型之外的数据类型视为无效的数据类型 end %查找数据格式 a=strfind(info,'interleave = '); b=length('interleave = '); %c=strfind(info,'byte order');%对于单波段 file type c=strfind(info,'sensor type');%对于多波段 interleave=[]; for i=a+b:c-1 interleave=[interleave,info(i)]; end interleave=strtrim(interleave);%删除字符串前面的空格 %读取图像文件 data = multibandread(imgfilename ,[lines, samples, bands],precision,0,interleave,'ieee-le'); data = im2double(data);%data=double(data); end % if(bands==1) % imagesc(data),colormap(cool) % colorbar,title('show by imagesc in figure1--band1_cool');axis off; % end % if(bands==6) % datar=data(:,:,1);%取第一页矩阵 % datag=data(:,:,2);%取第二页矩阵 % datab=data(:,:,3);%取第三页矩阵 % data4=data(:,:,4);%取第四页矩阵 % data5=data(:,:,5);%取第五页矩阵 % data6=data(:,:,6);%取第六页矩阵 % datamul=cat(3,datar,datag,datab);%合成3维矩阵 % datamul_2=cat(3,data4,data5,data6);%合成3维矩阵 % datamul_3=cat(3,datar,zeros(size(datar)),datab);%合成3维矩阵 % % subplot(2,2,1),imagesc(datar),colormap(cool)%整个figure中有2个图是排成一行的，一共2行，在第一块显示datar % colorbar,title('show by imagesc in figure1--band1_cool');axis off; % subplot(2,2,2),imagesc(datag),colormap(colorcube) % colorbar,title('show by imagesc in figure2--band2_colorcube');axis off; % subplot(2,2,3),imagesc(datab),colormap(colorcube) % colorbar,title('show by imagesc in figure3--band3_colorcube');axis off; % subplot(2,2,4),imagesc(data4),colormap(gray)%colormap(colorcube) % colorbar,title('show by imagesc in figure4--band4_gray');axis off; % figure(5),imshow(datamul),colormap(colorcube)%figure(4),imagesc(datar);colormap(colorcube) % colorbar,title('show by imagesc in figure5--band_1+2+3合成图_colorcube');axis off; % % figure(6),image(datamul_2); colormap(colorcube) % colorbar,title('show by image in figure6--band_4+5+6合成图_colorcube');axis off; % % figure(7),imshow(datamul_3),colormap(prism) % colorbar,title('show by imshow(I) in figure7--band_1+3合成图_prism');axis off; % % figure(8),imagesc(data5);colormap(copper) % colorbar,title('show by imagesc in figure8--band5_copper');axis off; % % figure(9),imagesc(data6);colormap(prism); % colorbar,title('show by imagesc in figure9--band6_prism');axis off; % end % end