function data=readimg(imgfilename)

%本函数读取img格式，前提是img图像显式带有'.img'后缀名。

if length(imgfilename)>=4

    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

%读取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

        error('invalid datatype');%除以上几种常见数据类型之外的数据类型视为无效的数据类型

end

%查找数据格式

a=strfind(info,'interleave = ');

b=length('interleave = ');

c=strfind(info,'byte order');

interleave=[];

for i=a+b:c-1

    interleave=[interleave,info(i)];

end

    interleave=strtrim(interleave);%删除字符串中的空格

%读取图像文件

fid = fopen(imgfilename, 'r');
%data=fread(fid,[lines, samples],precision,'ieee-le');
data = multibandread(imgfilename ,[lines, samples, bands],precision,0,interleave,'ieee-be');
data=single(data);
fclose(fid);
end

