Q1001momentsupto3(F10, F01)

function [mass, mus, sims, orthos, simorthos] = Q1001momentsupto3(F10, F01)
%------------------------------------------------------------------------------
%
% Computation of central moments (upto 3rd order) is called for, followed by 
% call for similitude invariants, orthogonal invariants and similitude with 
% orthogonal invariants combined.
%
% Input:
% F10       = Gridfunction that with F01 constitutes a quincunx gridfunction
% F01       = Gridfunction that with F10 constitutes a quincunx gridfunction
%
% Output:
% mass      = mass of quincunx gridfunction (F10 U F01)
% mus       = vector of length 7, containing the 2nd & 3rd order moments 
%             mu20 mu11 mu02 mu30 mu21 mu12 mu03
% sims      = Central moments made invariant to similitude transforms
% orthos    = Central moments made orthogonally invariant
% simorthos = Invariants w.r.t. both similitude and orthogonal transformations
%
% See pages 180, 181, 185 in:
% Ming-Kuei Hu  Visual Pattern Recognition by Moment Invariants,
% IRE Transactions on Information Theory, pp. 179--187 (1962).
%
% See also: masscenter, mupq, HUinvariants, momentsupto3, Q0011momentsupto3
%
% Design and implementation by:
% Dr. Paul M. de Zeeuw <[e-mail]>  http://homepages.cwi.nl/~pauldz/
% Last Revision: February 5, 2001.
% (c) 1999-2001 Stichting CWI, Amsterdam
%------------------------------------------------------------------------------
% Mass and center
[c, mass] = Q1001masscenter(F10, F01);
%
% Central moments
mu20 = Q1001mupq(F10, F01, 2, 0, c);
mu11 = Q1001mupq(F10, F01, 1, 1, c);
mu02 = Q1001mupq(F10, F01, 0, 2, c);
%
mu30 = Q1001mupq(F10, F01, 3, 0, c);
mu21 = Q1001mupq(F10, F01, 2, 1, c);
mu12 = Q1001mupq(F10, F01, 1, 2, c);
mu03 = Q1001mupq(F10, F01, 0, 3, c);
%
mus = [mu20 mu11 mu02 mu30 mu21 mu12 mu03];
%
% Invariants
[sims, orthos, simorthos] = HUinvariants(mass, mus);
%
%------------------------------------------------------------------------------