1. 程式人生 > >matlab讀取img格式檔案

matlab讀取img格式檔案



hdrfilename='a.hdr';
imgfilename='a.img';
%讀取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
            pp='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);


  
   datar=data(:,:,1);
   datag=data(:,:,2);
   datab=data(:,:,3);
  % subplot(2,2,1),imagesc(datar),colormap(gray);
   %subplot(2,2,2),imagesc(datag),colormap(gray);
   %subplot(2,2,3),imagesc(datab),colormap(gray);


   

剛開始用這個程式可能會遇到一些問題,比如提示什麼str2num的問題,最好開啟hdr檔案看一下是,有可能多一個空格或者少一個空格就會引起報錯。