function?[CV1CV2rhov1v2s11s22s12]?=?cancorr（fname1fname2?...
????nrowsncolsnvar1nvar2varargin）

%?[CV1CV2rhov1v2s11s22s12]?=?cancorr（fname1fname2?...
%???nrowsncolsnvar1nvar2varargin）
%
%?CANCORR?-?canonical?correlation?analysis
%
%?Input
%?fname1??-?file?name?of?multivariate?band?sequential?byte?float32?or?int16?input?image
%???????????number?one
%?fname2??-?file?name?of?multivariate?band?sequential?byte?float32?or?int16?input?image
%???????????number?two
%?nrows???-?number?of?rows?in?input?image?number?one
%?ncols???-?number?of?columns?in?input?image?number?one
%?nvars1??-?number?of?variables?or?bands?or?channels?in?input?image?number?one
%?nvars2??-?number?of?variables?or?bands?or?channels?in?input?image?number?two
%?fnameo1?-?outout?file?name?for?CV1?BSQ?float32?image??（optio-
%?fnameo2?-?outout?file?name?for?CV2?BSQ?float32?image???nal）
%
%?If?fnameo1?occurs?so?must?fnameo2.
%?
%?Output
%?the?canonical?variates?（CV1?and?CV2）
%?the?canonical?correlations?（rho）
%?the?eigenvectors?（v1?and?v2）?normed?to?give?CVs?unit?variance?and
%?the?relevant?（variance-）covariance?matrices.
%
%?Output?arrays?CV1?and?CV2?consist?of?transposed?images?must?be?viewed?with?e.g.
%???????imshow（reshape（CV1（::1）ncolsnrows）‘[-3?3]）
%
%?If?disk?output?is?requested?a?primitive?.hdr?file?for?the?output?file?is?written;
%?（a??full?ENVI?or?another?header?file?must?be?constructed?manually）.

%?（c）?Copyright?2005
%?Allan?Aasbjerg?Nielsen
%?aa@imm.dtu.dk?www.imm.dtu.dk/~aa
%?15?Feb?2005

if?nargin<6?error（‘Not?enough?input?arguments.‘）;?end
if?nargin>8?error（‘Too?many?input?arguments.‘）;?end
if?~ischar（fname1）?error（‘fname1?should?be?a?char?string‘）;?end
if?~ischar（fname2）?error（‘fname2?should?be?a?char?string‘）;?end
if?nargin==8
????fnameo1?=?varargin{1};
????if?~ischar（fnameo1）?error（‘fnameo1?should?be?a?char?string‘）;?end
????fnameo2?=?varargin{2};
????if?~ischar（fnameo2）?error（‘fnameo2?should?be?a?char?string‘）;?end
end
if?nvar2>nvar1
????error（‘input?with?highest?number?of?variables?must?be?first?set‘）;
end

%?open?as?byte?（uint8）?image?if?unsuccesful?open?as?float32?or?int16
fid1?=?fopen（fname1‘r‘）;
if?fid1==-1?error（strcat（fname1‘?not?found‘））;?end
[xcount1]?=?fread（fid1‘uint8‘）;
fclose（fid1）;
if?count1~=（nrows*ncols*nvar1）
????warning（‘data?in?fname1?do?not?match?nrows?ncols?nvars?for?uint8?try?float32‘）;
????fid1?=?fopen（fname1‘r‘）;
????if?fid1==-1?error（strcat（fname1‘?not?found‘））;?end
????[xcount1]?=?fread（fid1‘float32‘）;
????fclose（fid1）;
????if?count1~=（nrows*ncols*nvar1）
????????warning（‘data?in?fname1?do?not?match?nrows?ncols?nvars?for?float32?try?int16‘）;
????????fid1?=?fopen（fname1‘r‘）;
????????if?fid1==-1?error（strcat（fname1‘?not?found‘））;?end
????????[xcount1]?=?fread（fid1‘int16‘）;
????????fclose（fid1）;
????????if?count1~=（nrows*ncols*nvar1）
????????????error（‘data?in?fname1?do?not?match?nrows?ncols?nvars?for?int16?either‘）;
????????end
????end
end

fid2?=?fopen（fname2‘r‘）;
if?fid2==-1?er