function?[varargout]?=?enviread（varargin）
%?ENVIREAD?Easily?import?ENVI?raster?files?（BSQBILBIP）?from?header?info.
%???Z?=?ENVIREAD（FILENAME）;?Reads?an?ENVI?binary?file?（BSQBILBIP）?into
%???????an?array?using?the?information?from?the?corresponding?header
%???????file?FILENAME.hdr.?The?output?array?will?be?of?dimensions
%???????（mnb）?where?b?is?the?number?of?bands.
%???Z?=?ENVIREAD（FILENAMEHEADERFILE）;?Uses?the?header?information?in
%???????headerfile.
%???[ZXY]?=?ENVIREAD（....）;?Returns?the?map?coordinate?vectors?for?geo-
%???????registered?data.
%???[ZXYinfo]?=?ENVIREAD（....）;?Returns?the?header?information?as?a
%???????structure.
%
%???NOTES:??-Requires?READ_ENVIHDR?to?read?header?data.
%??????????-Geo-registration?does?not?currently?support?rotated?images.
%
%?Ian?M.?Howat?Applied?Physics?Lab?University?of?Washington
%?ihowat@apl.washington.edu
%?Version?1:?11-Jul-2007?15:11:13

%%?check?for?header?reader
if?exist（‘read_envihdr.m‘‘file‘）?==?0
????error（‘This?function?requires?READ_ENVIHDR.m‘）
end
%%?READ?HEADER?INFO
%?Read?Filename
file?=?varargin{1};
hdrfile?=?[file‘.hdr‘];
if?nargin?==?2
????hdrfile?=?varargin{2};
end
%?Get?image?size?and?map?data?from?header
info?=?read_envihdr（hdrfile）;
%%?Make?geo-location?vectors
if?isfield（info.map_info‘mapx‘）?&&?isfield（info.map_info‘mapy‘）
????xi?=?info.map_info.image_coords（1）;
????yi?=?info.map_info.image_coords（2）;
????xm?=?info.map_info.mapx;
????ym?=?info.map_info.mapy;
????%adjust?points?to?corner?（1.51.5）
????if?yi?>?1.5?
???????ym?=??ym?+?（（yi*info.map_info.dy）-info.map_info.dy）;
????end
????if?xi?>?1.5?
????????xm?=?xm?-?（（xi*info.map_info.dy）-info.map_info.dx）
????end
????????
????varargout{2}?=?xm?+?（（0:info.samples-1）.*info.map_info.dx）;
????varargout{3}?=?fliplr（ym?-?（（0:info.lines-1）.*info.map_info.dy））;
end
%%?Set?binary?format?parameters
switch?info.byte_order
????case?{0}
????????machine?=?‘ieee-le‘;
????case?{1}
????????machine?=?‘ieee-be‘;
????otherwise
????????machine?=?‘n‘;
end
switch?info.data_type
????case?{1}
????????format?=?‘int8‘;
????case?{2}
????????format=?‘int16‘;
????case{3}
????????format=?‘int32‘;
????case?{4}
????????format=?‘float‘;
????case?{5}
????????format=?‘double‘;
????case?{6}
????????disp（‘>>?Sorry?Complex?（2x32?bits）data?currently?not?supported‘）;
????????disp（‘>>?Importing?as?double-precision?instead‘）;
????????format=?‘double‘;
case?{9}
????????error（‘Sorry?double-precision?complex?（2x64?bits）?data?currently?not?supported‘）;
case?{12}
????????format=?‘uint16‘;
case?{13}
????????format=?‘uint32‘;
case?{14}
?????????format=?‘int64‘;
case?{15}
????????format=?‘uint64‘;
????otherwise
????????error（[‘File?type?number:?‘num2str（dtype）‘?not?supported‘]）;
end
%%?Read?File
Z?=?fread（fopen（file）info.samples*info.lines*info.bandsformat0machine）;?fclose?all;
switch?lower（info.interleave）
????case?{‘bsq‘}
????????Z?=?reshape（Z[info.samplesinfo.lines

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3556??2007-08-19?10:53??myenviread.m

-----------?---------??----------?-----??----

?????????????????3556????????????????????1