gmpTetGenBuildGroups.h

/************************************************************************
**
** Copyright (C) 2014 by Carlos Augusto Teixera Mendes
** All rights reserved.
**
** This file is part of the "GeMA" software. It's use should respect
** the terms in the license agreement that can be found together
** with this source code.
** It is provided AS IS, with NO WARRANTY OF ANY KIND,
** INCLUDING THE WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR
** A PARTICULAR PURPOSE.
**
************************************************************************/
 
#ifndef _GEMA_PLUGIN_TETGEN_BUILDGROUPS_H_
#define _GEMA_PLUGIN_TETGEN_BUILDGROUPS_H_
 
#include <vector>
class GmCellMesh;
class GmValueAccessor;
class GmCellAccessor;
 
/*
*   Group cells into volumes based on the surfaces received.
*   It uses a flood fill approach, so surfaces must be 'closed enough' to form
*   volumes in the given grid discretization and bounding box.
*   Also tries to group similar volumes together
*   based on their intersection histogram if groupThreshold is set to a value > 0 
* 
        Given a mesh:
          - mpoints: list of mesh points (x1,y1,z1,x2,y2,z2,...)
          - mcells: list of mesh cell connectivity (n1_1, n1_2, n1_3, n1_4, n2_1, n2_2, n2_3, n2_4,...)
          - mcellSize: list of cell sizes (number of nodes per cell, e.g. 4 for tetrahedra)
        Set of triangle surfaces forming volumes:
          - surfPoints: list of surface points (S1_points, S2_points, ...), points = (x1,y1,z1,x2,y2,z2,...)
          - surfTri: list of surface triangle connectivity (S1_triangles, S2_triangles, ...),  triangles = (n1_1, n1_2, n1_3, n2_1, n2_2, n2_3,...)
          - surfPCount: list of surface points count (number of points in each surfPoints entry)
          - surfTCount: list of surface triangle count (number of triangles in each surfTri entry)
        Bounding box:
          - bmin: minimum coordinates of the bounding box (x_min, y_min, z_min)
          - bmax: maximum coordinates of the bounding box (x_max, y_max, z_max)
        Segmentation/Grouping parameters:
          - groupThreshold: minimum percentage of grid intersections in a volume to be not ignored (remove outliers)
          - gridRes: resolution of the grid used for flood fill (e.g. number of cells along the longest dimension of the bounding box)
 
         Returns a tuple with:
          - Element classes (segmentation) 
          - Grouping; for each class 'i', group[i] is what other class it groups to
*/
std::tuple < std::vector<int>, std::vector<int>> 
        buildGroups(const std::vector<double>& mpoints, const std::vector<int>& mcells, const std::vector<int>& mcellSize,
        const std::vector<std::vector<double>>& surfPoints, const std::vector<std::vector<int>>& surfTri,
        double bmin[3], double bmax[3], double groupThreshold, int gridRes[3]);
 
/*
        Build groups and set directly as a mesh property
        - Mesh: mesh to build groups
        - propAc: cell property accessor to set the group property
*/
std::tuple < std::vector<int>, std::vector<int>> 
        buildGroups(GmCellMesh* mesh, GmCellAccessor* propAc, GmValueAccessor* coordAc,
        const std::vector<std::vector<double>>& surfPoints, const std::vector<std::vector<int>>& surfTri,
        double bmin[3], double bmax[3], double groupThreshold, int gridRes[3]);
 
#endif