Method compute_proxy_mesh

code/method/hypothesis_generator.cpp:335–480

Map* HypothesisGenerator::compute_proxy_mesh(Map* bbox_mesh) {
	MapFacetAttribute<Plane3d*>					bbox_mesh_face_supporting_plane(bbox_mesh, "FacetSupportingPlane");
	MapHalfedgeAttribute< std::set<Plane3d*> >	bbox_mesh_edge_source_planes(bbox_mesh, "EdgeSourcePlanes");
	MapVertexAttribute< std::set<Plane3d*> >	bbox_mesh_vertex_source_planes(bbox_mesh, "VertexSourcePlanes");

	Map* mesh = new Map;
	MapBuilder builder(mesh);

	MapFacetAttribute<Color> color(mesh, "color");
	MapFacetAttribute<VertexGroup*> facet_supporting_vertex_group(mesh, Method::facet_attrib_supporting_vertex_group);
	MapFacetAttribute<Plane3d*>		face_supporting_plane(mesh, "FacetSupportingPlane");
	MapHalfedgeAttribute< std::set<Plane3d*> >	edge_source_planes(mesh, "EdgeSourcePlanes");
	MapVertexAttribute< std::set<Plane3d*> >	vertex_source_planes(mesh, "VertexSourcePlanes");

	builder.begin_surface();
	int idx = 0;
	for (std::size_t i = 0; i < plane_segments_.size(); ++i) {
		VertexGroup* g = plane_segments_[i];
		Plane3d* plane = vertex_group_plane_[g];

		std::vector<vec3> points;
		std::vector< std::set<Plane3d*> > point_source_planes;
		FOR_EACH_EDGE_CONST(Map, bbox_mesh, it) {
			const vec3& s = it->prev()->vertex()->point();
			const vec3& t = it->vertex()->point();
			Sign ss = plane->orient(s);
			Sign st = plane->orient(t);
			if ((ss == POSITIVE && st == NEGATIVE) || (ss == NEGATIVE && st == POSITIVE)) {
				vec3 p;
				if (plane->intersection(Line3d::from_two_points(s, t), p)) {
					points.push_back(p);

					std::set<Plane3d*> planes = bbox_mesh_edge_source_planes[it];
					planes.insert(plane);
					point_source_planes.push_back(planes);
				}
				else
					Logger::err("-") << "fatal error. Should have intersection" << std::endl;
			}
			else {
				if (ss == ZERO) {
					points.push_back(s);

					std::set<Plane3d*> planes = bbox_mesh_vertex_source_planes[it->prev()->vertex()];
					point_source_planes.push_back(planes);
				}
				else if (st == ZERO) {
					points.push_back(t);

					std::set<Plane3d*> planes = bbox_mesh_vertex_source_planes[it->vertex()];
					point_source_planes.push_back(planes);
				}
				else {
					// no intersection with the plane
				}
			}
		}