Changeset 10274

Timestamp:

01/03/14 01:17:02 (10 years ago)

Author:

ascheibe

Message:

#1886 added calculation of the centroid of a convex hull and measures for motion characteristics

Location:

branches/HeuristicLab.Analysis.AlgorithmBehavior/HeuristicLab.Analysis.AlgorithmBehavior.Analyzers/3.3

Files:

• ConvexHullMeasures.cs (modified) (1 diff)
• DoubleArrayExtensions.cs (modified) (1 diff)

branches/HeuristicLab.Analysis.AlgorithmBehavior/HeuristicLab.Analysis.AlgorithmBehavior.Analyzers/3.3/ConvexHullMeasures.cs

@@ -89 +89 @@
       return result;
     }
+
+    public static double[] CalculateCentroid(List<double[]> convexHull) {
+      int n = convexHull.Count;
+      int dim = convexHull.First().Length;
+      double[] centroid = new double[dim];
+
+      for (int i = 0; i < dim; i++) {
+        for (int j = 0; j < n; j++) {
+          centroid[i] += convexHull[j][i];
+        }
+        centroid[i] /= n;
+      }
+      return centroid;
+    }
+
+    //measures to calculate: AVG, SUM, STD.DEV
+    public static double[] CalculateMovementDistances(List<double[]> centroids) {
+      double[] distances = new double[centroids.Count - 1];
+
+      // it is assumed that the centroids are sorted chronological
+      for (int i = 0; i < centroids.Count - 1; i++) {
+        distances[i] = centroids[i].EuclideanDistance(centroids[i + 1]);
+      }
+
+      return distances;
+    }
+
+    public static double CalculateOverallMovementDistances(List<double[]> centroids) {
+      // it is assumed that the centroids are sorted chronological
+      return centroids[0].EuclideanDistance(centroids[centroids.Count - 1]);
+    }
+
+    //calculate theta from two vectors (dot product)
+    public static double CalculateTheta(double[] a, double[] b) {
+      return Math.Acos(a.DotProduct(b) / (a.EuclideanNorm() * b.EuclideanNorm()));
+    }
+
+    //calculate theta using law of cosines
+    public static double CalculateTheta(double[] a, double[] b, double[] c) {
+      double ab = a.EuclideanDistance(b);
+      double bc = b.EuclideanDistance(c);
+      double ca = c.EuclideanDistance(a);
+
+      return Math.Acos((ab * ab + bc * bc - ca * ca) / (2 * ab * bc));
+    }
+
+    //measures to calculate: AVG, SUM, STD.DEV
+    public static double[] CalculateCentroidsMotion(List<double[]> centroids) {
+      double[] thetas = new double[centroids.Count - 2];
+
+      for (int i = 1; i < centroids.Count - 1; i++) {
+        thetas[i - 1] = CalculateTheta(centroids[i - 1], centroids[i], centroids[i + 1]);
+      }
+
+      return thetas;
+    }
   }
 }

branches/HeuristicLab.Analysis.AlgorithmBehavior/HeuristicLab.Analysis.AlgorithmBehavior.Analyzers/3.3/DoubleArrayExtensions.cs

@@ -31 +31 @@
       return Math.Sqrt(result);
     }
+
+    public static double EuclideanNorm(this double[] p) {
+      double norm = 0.0;
+      for (int i = 0; i < p.Length; i++) {
+        norm += p[i] * p[i];
+      }
+      return Math.Sqrt(norm);
+    }
+
+    public static double DotProduct(this double[] p1, double[] p2) {
+      double sum = 0.0;
+      for (int i = 0; i < p1.Length; i++) {
+        sum += p1[i] * p2[i];
+      }
+      return sum;
+    }
   }
 }