卫星图像读取

matlab版本卫星遥感图像数据的读取 与打开

function [ colslinesbandsdata_typeoffsetimaging_way ] = myNewMultiBandRead（ hdrFileNamehdrFilePath ）
%UNtitleD2 此处显示有关此函数的摘要
%   此处显示详细说明
elements={‘samples ‘ ‘lines ‘ ‘bands ‘ ‘header offset ‘ ‘data type ‘ ‘interleave ‘};
way={‘ bsq‘ ‘ bil‘ ‘ bip‘};
%检查相应的头文件是否能够打开
rfid = fopen（fullfile（hdrFilePathhdrFileName）‘r‘）;
if （rfid == -1）
     error（‘Input header files does not exist‘）;
end
cols=2048;lines=2100;bands=150;data_type=‘uint16‘;offset=0;imaging_way=‘bil‘;
%如果能够正确打开
 while 1
    tline = fgetl（rfid）;
    if ~ischar（tline） break end
    [firstsecond]=strtok（tline‘=‘）;
    switch first
        case elements（1）
            [~s]=strtok（second）;
            cols=str2double（s）;
        case elements（2）
            [~s]=strtok（second）;
            lines=str2double（s）;
        case elements（3）
            [~s]=strtok（second）;
            bands=str2double（s）;
        case elements（4）
            [~s]=strtok（second）;
            offset=str2double（s）;
        case elements（5）
            [~s]=strtok（second）;
            t=str2double（s）;
            switch t