public static void marchingCubesDemo()
throws MeshNotOrientedException, DanglingTriangleException, IOException {
PointCloud pc = ObjReader.readAsPointCloud("objs/dragon_withNormals.obj", true);
// PointCloud pc = PlyReader.readPointCloud("objs/angel_points.ply", true);
pc.normalizeNormals();
HashOctree tree = new HashOctree(pc, 9, 1, 1.3f);
tree.refineTree(2);
LinearSystem system = SSDMatrices.ssdSystem(tree, pc, 1, 0.001f, 10);
// Test Data: create an octree
ArrayList<Float> functionByVertex = new ArrayList<Float>();
SCIPY.solve(system, "whatev", functionByVertex);
// System.out.println(functionByVertex);
MarchingCubes mc = new MarchingCubes(tree);
mc.dualMC(functionByVertex);
WireframeMesh mesh = mc.result;
GLWireframeMesh glMesh = new GLWireframeMesh(mesh);
// And show off...
// visualization of the per vertex values (blue = negative,
// red = positive, green = 0);
MyDisplay d = new MyDisplay();
glMesh.configurePreferredShader(
"shaders/trimesh_flat.vert", "shaders/trimesh_flat.frag", "shaders/trimesh_flat.geom");
d.addToDisplay(glMesh);
// discrete approximation of the zero level set: render all
// tree cubes that have negative values.
GLHashtree gltree = new GLHashtree(tree);
gltree.addFunctionValues(functionByVertex);
gltree.configurePreferredShader(
"shaders/octree_zro.vert", "shaders/octree_zro.frag", "shaders/octree_zro.geom");
d.addToDisplay(gltree);
}
