Velocity Moments Explained

In the field of computer vision, velocity moments are weighted averages of the intensities of pixels in a sequence of images, similar to image moments but in addition to describing an object's shape also describe its motion through the sequence of images. Velocity moments can be used to aid automated identification of a shape in an image when information about the motion is significant in its description. There are currently two established versions of velocity moments: Cartesian^[1] and Zernike.^[2]

Cartesian velocity moments

Cartesian moments for single images

A Cartesian moment of a single image is calculated by

m _{pq} = \sum ^{M} _{x =1} \sum ^{N} _{y =1} x ^{p} y ^{q} P _{xy}

where

M and N are the dimensions of the image, P _{xy} is the intensity of the pixel at the point (x,y) in the image, and x ^{p} y ^{q} is the basis function.

Cartesian velocity moments for sequences of images

Cartesian velocity moments are based on these Cartesian moments. A Cartesian velocity moment

vm _{pq\mu\gamma} is defined byvm _{pq\mu\gamma} = \sum ^{images} _{i =2} \sum ^{M} _{x =1} \sum ^{N} _{y =1} U(i,\mu,\gamma) C(i,p,g) P _{i _{xy}}

where

M and N are again the dimensions of the image, images is the number of images in the sequence, and P _{i _{xy}} is the intensity of the pixel at the point (x,y) in image i .C(i,p,q) is taken from Central moments, added so the equation is translation invariant, defined asC(i,p,q) = (x - \overline{x ^{p} _{i})} (y - \overline{y ^{q} _{i})}

where

\overline{x _{i}} is the x coordinate of the centre of mass for image i , and similarly for y .U(i,\mu,\gamma) introduces velocity into the equation asU(i,\mu,\gamma) = (\overline{x _{i}-\overline{x} _{i-1} })^\mu (\overline-\overline)^\gamma

where

\overline{x _{i-1} } is the x coordinate of the centre of mass for the previous image, i -1 , and again similarly for y .

After the Cartesian velocity moment is calculated, it can be normalised by

\overline{vm _{pq\mu\gamma} } = \frac

where

A is the average area of the object, in pixels, and I is the number of images. Now the value is not affected by the number of images in the sequence or the size of the object.

As Cartesian moments are non-orthogonal, so are Cartesian velocity moments, so different moments can be closely correlated. These velocity moments do however provide translation and scale invariance (unless the scale changes within the sequence of images).

Zernike velocity moments

Zernike moments for single images

A Zernike moment of a single image is calculated by

A _{mn} = (m +1)/(\pi) \sum _{x} \sum _{y} [V _{mn} (r,\theta)] ^{*} P _{xy}

where

^{*} denotes the complex conjugate, m is an integer between 0 and inf ty , and n is an integer such that m - |n| is even and |n| < m . For calculating Zernike moments, the image, or section of the image which is of interest is mapped to the unit disc, then P _{xy} is the intensity of the pixel at the point (x,y) on the disc and x ^{2} + y ^{2} \le 1 is a restriction on values of x and y . The coordinates are then mapped to polar coordinates, and r and \theta are the polar coordinates of the point (x,y) on the unit disc map.V _{mn} (r,\theta) is derived from Zernike polynomials and is defined byV _{mn} (r,\theta) = R _{mn} (r)e ^{jn\theta} R _{mn} (r) = \sum ^{(m - |n|)/(2)} _{s =0} ( -1 ) ^{s} F(m,n,s,r) F(m,n,s,r) = ((m - s)!)/(s! ( (m + |n|)/(2- s)!) ( (m - |n|)/(2- s)!)) r ^{m-2s}

Zernike velocity moments for sequences of images

Zernike velocity moments are based on these Zernike moments. A Zernike velocity moment

A _{mn\mu\gamma} is defined byA _{mn\mu\gamma} = (m +1)/(\pi) \sum ^{images} _{i =2} \sum _{x=1} \sum _{y=1} U(i,\mu,\gamma) [V _{mn} (r,\theta)] ^{*} P _{i _{xy}}

where

images is again the number of images in the sequence, and P _{i _{xy}} is the intensity of the pixel at the point (x,y) on the unit disc mapped from image i .U(i,\mu,\gamma) introduces velocity into the equation in the same way as in the Cartesian velocity moments and [V _{mn} (r,\theta)] ^{*} is from the Zernike moments equation above.

Like the Cartesian velocity moments, Zernike velocity moments can be normalised by

\overline{A _{mn\mu\gamma} } = \frac

where

A is the average area of the object, in pixels, and I is the number of images.

As Zernike velocity moments are based on the orthogonal Zernike moments, they produce less correlated and more compact descriptions than Cartesian velocity moments. Zernike velocity moments also provide translation and scale invariance (even when the scale changes within the sequence).

Comparison of methods

Velocity moment type	Translation invariance	Scale invariance	Orthogonal
Cartesian	Yes	Yes (if the object does not change scale within the sequence)	No
Zernike	Yes	Yes	Yes

References

1. J. D. Shutler, M. S. Nixon, C. J. Harris, "Statistical Gait Description via Temporal Moments", Proc. SSIAI 2000 - Austin, Texas, :pp. 291-295, 2000
2. J. D. Shutler and M. S. Nixon, "Zernike Velocity Moments for Description and Recognition of Moving Shapes", Proc. BMVC 2001, Manchester, UK, 2:pp. 705-714, 2001

External links

• CVonline Velocity Moments page