Added tutorial 701

2016-06-21 11:54:53 +02:00
parent c0e04685e2
commit 54e2fe0935
7 changed files with 121 additions and 0 deletions
@@ -0,0 +1,59 @@
+import sys, os
+
+# Add the igl library to the modules search path
+import math
+
+sys.path.insert(0, os.getcwd() + "/../")
+import pyigl as igl
+
+from shared import TUTORIAL_SHARED_PATH, check_dependencies
+
+dependencies = []
+check_dependencies(dependencies)
+
+if __name__ == "__main__":
+    V = igl.eigen.MatrixXd()
+    F = igl.eigen.MatrixXi()
+
+    # Load meshes in OFF format
+    igl.readOBJ(TUTORIAL_SHARED_PATH + "horse_quad.obj", V, F)
+
+    # Count the number of irregular vertices, the border is ignored
+    irregular = igl.is_irregular_vertex(V, F)
+    vertex_count = V.rows()
+    irregular_vertex_count = sum(irregular)
+    irregular_ratio = irregular_vertex_count / vertex_count
+
+    print("Irregular vertices: \n%d/%d (%.2f%%)\n" % (
+        irregular_vertex_count, vertex_count, irregular_ratio * 100))
+
+    # Compute areas, min, max and standard deviation
+    area = igl.eigen.MatrixXd()
+    igl.doublearea(V, F, area)
+    area /= 2.0
+
+    area_avg = area.mean()
+    area_min = area.minCoeff() / area_avg
+    area_max = area.maxCoeff() / area_avg
+    area_ns = (area - area_avg) / area_avg
+    print(area_ns[0])
+    # area_ns *= area_ns
+    print(area_ns[0])
+    area_sigma = math.sqrt(area_ns.mean())
+
+    print("Areas (Min/Max)/Avg_Area Sigma: \n%.2f/%.2f (%.2f)\n" % (
+        area_min, area_max, area_sigma))
+
+    # Compute per face angles, min, max and standard deviation
+    angles = igl.eigen.MatrixXd()
+    igl.internal_angles(V, F, angles)
+    angles = 360.0 * (angles / (2 * math.pi))
+
+    angle_avg = angles.mean()
+    angle_min = angles.minCoeff()
+    angle_max = angles.maxCoeff()
+    angle_ns = (angles - angle_avg) * (angles - angle_avg)
+    angle_sigma = math.sqrt(angle_ns.mean())
+
+    print("Angles in degrees (Min/Max) Sigma: \n%.2f/%.2f (%.2f)\n" % (
+        angle_min, angle_max, angle_sigma))