example01.m

来自「用matlab实现图像融合」· M 代码 · 共 29 行

% Small technical example.
% Performs a check on the computation of certain central moments.
%
disp('Small technical example.');
disp('Performs a check on the computation of certain central moments.');
disp('FOR MORE INFORMATION:  help mupq');
disp(' ');
disp('See also the report http://ftp.cwi.nl/CWIreports/PNA/PNA-R0224.pdf');
disp('Dr. Paul M. de Zeeuw <[e-mail]>');
disp('Copyright (C) 1998-2006 Stichting CWI, Amsterdam, The Netherlands');
disp(' ');
%---INSERT YOUR IMAGE HERE-----------------------------------------------------
if exist('imread','file') == 2
  Orig = double(imread('zenithgray.TIF','tiff'));
else
  load zenithgray; Orig = zenithgray; clear zenithgray;
end
%
disp([' Dimensions ' int2str( size(Orig) )]);
%------------------------------------------------------------------------------
disp([' Mass  ' num2str( mupq(Orig, 0, 0), '%+12.4e')]);
%
disp(' NICE CHECK: First order central moment (x-dir) should vanish');
disp([' reads ' num2str( mupq(Orig, 1, 0), '%+12.4e')]);
%
disp(' NICE CHECK: First order central moment (y-dir) should vanish');
disp([' reads ' num2str( mupq(Orig, 0, 1), '%+12.4e')]);
%