wc.m

图像对齐（inage registration）工具箱

function [code,msg ] = wc( filename,matrix )
%WC  write an image in pgm format
% syntax:
%     [ code,msg ] = writecon( filename, matrix  )
%
% writecon write an matrix into file filename.pgm  
%
% code is zero if operation was succesfull, nonzero otherwise.
% msg contains error message.
%

%     Copyright (c) P. Turcaj, R. Turcajova, J. Kautsky, 
%     Flinders University of South Australia.
%     Use of this software without permission of authors is prohibited.
%     Modified by Jan Flusser, UTIA

code=0;
if (nargin<2 | nargin>2 ) 
   code=1;
   msg='writecon: wrong number of parameters';
   return;
end;
[fid, mess] = fopen([filename,'.pgm'],'w');
if fid==-1
   code=2;
   msg=['writecon:can not open file ' filename '.pgm:' mess];
   return;
end;

count = fprintf(fid,'%s\n','P5');
if count~=3
   code=3;
   msg='writecon: sorry, can not write file';
   return;
end;

[n,m]=size(matrix);
matrix=matrix';
col=255;
count = fprintf(fid,'%d ',m);
count = fprintf(fid,'%d\n',n);
count = fprintf(fid,'%d \n',col+1);

prec = ['uint',int2str(8)];
%prec = ['uchar'];

mi=min(min(matrix));
ma=max(max(matrix));
if (mi<0 | ma>255)
 for i=1:m,
  for j=1:n,
   matrix(i,j)=fix((matrix(i,j)-mi)/ma*col);
  end;
 end;
else
matrix=(round(matrix));
end;
count = fwrite(fid,matrix,prec);
if count~=m*n
  code=4;
  msg=['writecon: warning! file ' filename ...
        '.pgm was not written completely!'];
end;

fclose(fid);