/// <summary>
/// Computes the center moments for the specified image.
/// </summary>
///
/// <param name="image">The image whose moments should be computed.</param>
/// <param name="area">The region of interest in the image to compute moments for.</param>
///
public override void Compute(float[,] image, Rectangle area)
{
RawMoments raw = new RawMoments(image, area, Order);
CentralMoments center = new CentralMoments(raw);
this.Compute(center);
}
public void ComputeTest()
{
Bitmap image = Resources.hu;
CentralMoments target = new CentralMoments(image, order: 3);
Assert.AreEqual(86424.0 / target.Mu00, 1);
Assert.AreEqual(0, target.Mu01);
Assert.AreEqual(0, target.Mu10);
Assert.AreEqual(5.868206472635379E8 / target.Mu02, 1, 1e-2);
Assert.AreEqual(6348920.945848465 / target.Mu11, 1, 1e-2);
Assert.AreEqual(9.084235762166061E8 / target.Mu20, 1, 1e-3);
Assert.AreEqual(-2.155191E9 / target.Mu12, 1, 1e-5);
Assert.AreEqual(7.125893E8 / target.Mu21, 1, 1e-3);
Assert.AreEqual(-1.26244547E10 / target.Mu30, 1, 1e-8);
Assert.AreEqual(1.71818829E9 / target.Mu03, 1, 1e-8);
SizeF size = target.GetSize();
float angle = target.GetOrientation();
Assert.AreEqual(410.207916f, size.Height);
Assert.AreEqual(329.534637f, size.Width);
Assert.AreEqual(0.0196384024f, angle);
}
/// <summary>
/// Computes the Hu moments from the specified central moments.
/// </summary>
///
/// <param name="moments">The central moments to use as base of calculations.</param>
///
public void Compute(CentralMoments moments)
{
double inv = 1.0 / moments.Mu00;
double inv2 = 1.0 / (moments.Mu00 * moments.Mu00);
double inv5d2 = Math.Sqrt(inv2 * inv2 * inv);
float n20 = (float)(moments.Mu20 * inv2);
float n02 = (float)(moments.Mu02 * inv2);
float n11 = (float)(moments.Mu11 * inv2);
float n21 = (float)(moments.Mu21 * inv5d2);
float n12 = (float)(moments.Mu12 * inv5d2);
float n30 = (float)(moments.Mu30 * inv5d2);
float n03 = (float)(moments.Mu03 * inv5d2);
// (η20 + η02)
I1 = (n20 + n02);
// (η20 − η02)² + 4 η11²
I2 = (n20 - n02) * (n20 - n02) + 4 * (n11 * n11);
// (η30 − 3 η12)²
I3 = (n30 - 3 * n12) * (n30 - 3 * n12)
// + (3 η21 − η03)²
+ (3 * n21 - n03) * (3 * n21 - n03);
// (η30 + η12)² + (η21 + η03)²
I4 = (n30 + n12) * (n30 + n12) + (n21 + n03) * (n21 + n03);
// (η30 − 3 η12) (η30 + η12) [(η30 + η12)² −3 (η21 + η03)² ]
I5 = (n30 - 3 * n12) * (n30 + n12) * ((n30 + n12) * (n30 + n12) - 3 * (n21 + n03) * (n21 + n03))
// (3 η21 − η03) (η21 + η03) [3 (η30 + η12)² − (η21 + η03)² ]
+ (3 * n21 - n03) * (n21 + n03) * (3 * (n30 + n12) * (n30 + n12) - (n21 + n03) * (n21 + n03));
// (η20 − η02) [(η30 + η12)² − (η21 + η03)² ]
I6 = (n20 - n02) * ((n30 + n12) * (n30 + n12) - (n21 + n03) * (n21 + n03))
// + 4 η11 (η30 + η12) (η21 + η03)
+ 4 * n11 * (n30 + n12) * (n21 + n03);
// (3 η21 − η03) (η30 + η12) [(η30 + η12)² − 3 (η21 + η03)² ]
I7 = (3 * n21 - n03) * (n30 + n12) * (n30 + n12) * ((n30 + n12) * (n30 + n12) - 3 * (n21 + n03) * (n21 + n03))
// - (η30 − 3 η12) (η21 + η03) [3 (η30 + η12)² − (η21 + η03)² ]
- (n30 - 3 * n12) * (n21 + n03) * (3 * (n30 + n12) * (n30 + n12) - (n21 + n03) * (n21 + n03));
}
public void ComputeTest2()
{
// 0 and 1 are only translated
var cm0 = new CentralMoments(Resources.hu0, 3);
var cm1 = new CentralMoments(Resources.hu1, 3);
Assert.AreEqual(cm0.Mu00, cm1.Mu00);
Assert.AreEqual(cm0.Mu01, cm1.Mu01);
Assert.AreEqual(cm0.Mu10, cm1.Mu10);
Assert.AreNotEqual(cm0.Mu11, cm1.Mu11);
}
public void ComputeTest2()
{
// 0 and 1 are only translated
var cm0 = new CentralMoments(Accord.Imaging.Image.Clone(Resources.hu0), 3);
var cm1 = new CentralMoments(Accord.Imaging.Image.Clone(Resources.hu1), 3);
Assert.AreEqual(cm0.Mu00, cm1.Mu00);
Assert.AreEqual(cm0.Mu01, cm1.Mu01);
Assert.AreEqual(cm0.Mu10, cm1.Mu10);
Assert.AreNotEqual(cm0.Mu11, cm1.Mu11);
}
/// <summary>
/// Process a new video frame.
/// </summary>
public void ProcessFrame(UnmanagedImage frame)
{
filterImage = filter.Apply(frame);
Blob blob = extractBlob();
if (blob == null)
{
trackingObject.Reset();
return;
}
trackingObject.Rectangle = blob.Rectangle;
trackingObject.Center = (IntPoint)blob.CenterOfGravity;
if (rotation)
{
// Locate moments
CentralMoments moments = new CentralMoments();
moments.Compute(filterImage, blob.Rectangle);
trackingObject.Angle = moments.GetOrientation();
}
if (extract)
{
blobCounter.ExtractBlobsImage(filterImage, blob, false);
trackingObject.Image = blob.Image;
}
}
/// <summary>
/// Computes the center moments for the specified image.
/// </summary>
///
/// <param name="image">The image whose moments should be computed.</param>
/// <param name="area">The region of interest in the image to compute moments for.</param>
///
public override void Compute(float[,] image, Rectangle area)
{
RawMoments raw = new RawMoments(image, area, Order);
CentralMoments center = new CentralMoments(raw);
this.Compute(center);
}
/// <summary>
/// Computes the Hu moments from the specified central moments.
/// </summary>
///
/// <param name="moments">The central moments to use as base of calculations.</param>
///
public void Compute(CentralMoments moments)
{
double inv = 1.0 / moments.Mu00;
double inv2 = 1.0 / (moments.Mu00 * moments.Mu00);
double inv5d2 = Math.Sqrt(inv2 * inv2 * inv);
float n20 = (float)(moments.Mu20 * inv2);
float n02 = (float)(moments.Mu02 * inv2);
float n11 = (float)(moments.Mu11 * inv2);
float n21 = (float)(moments.Mu21 * inv5d2);
float n12 = (float)(moments.Mu12 * inv5d2);
float n30 = (float)(moments.Mu30 * inv5d2);
float n03 = (float)(moments.Mu03 * inv5d2);
// (η20 + η02)
I1 = (n20 + n02);
// (η20 − η02)² + 4 η11²
I2 = (n20 - n02) * (n20 - n02) + 4 * (n11 * n11);
// (η30 − 3 η12)²
I3 = (n30 - 3 * n12) * (n30 - 3 * n12)
// + (3 η21 − η03)²
+ (3 * n21 - n03) * (3 * n21 - n03);
// (η30 + η12)² + (η21 + η03)²
I4 = (n30 + n12) * (n30 + n12) + (n21 + n03) * (n21 + n03);
// (η30 − 3 η12) (η30 + η12) [(η30 + η12)² −3 (η21 + η03)² ]
I5 = (n30 - 3 * n12) * (n30 + n12) * ((n30 + n12) * (n30 + n12) - 3 * (n21 + n03) * (n21 + n03))
// (3 η21 − η03) (η21 + η03) [3 (η30 + η12)² − (η21 + η03)² ]
+ (3 * n21 - n03) * (n21 + n03) * (3 * (n30 + n12) * (n30 + n12) - (n21 + n03) * (n21 + n03));
// (η20 − η02) [(η30 + η12)² − (η21 + η03)² ]
I6 = (n20 - n02) * ((n30 + n12) * (n30 + n12) - (n21 + n03) * (n21 + n03))
// + 4 η11 (η30 + η12) (η21 + η03)
+ 4 * n11 * (n30 + n12) * (n21 + n03);
// (3 η21 − η03) (η30 + η12) [(η30 + η12)² − 3 (η21 + η03)² ]
I7 = (3 * n21 - n03) * (n30 + n12) * (n30 + n12) * ((n30 + n12) * (n30 + n12) - 3 * (n21 + n03) * (n21 + n03))
// - (η30 − 3 η12) (η21 + η03) [3 (η30 + η12)² − (η21 + η03)² ]
- (n30 - 3 * n12) * (n21 + n03) * (3 * (n30 + n12) * (n30 + n12) - (n21 + n03) * (n21 + n03));
}
