123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990 |
- //
- // removeDuplicates.h
- // Preview3D
- //
- // Created by Olga Diamanti on 17/11/11.
- // Copyright (c) 2011 __MyCompanyName__. All rights reserved.
- //
- #ifndef RemoveDuplicates_h
- #define RemoveDuplicates_h
- #include <Eigen/Core>
- namespace igl
- {
- // [ NV, NF ] = removeDuplicates( V,F,epsilon )
- // Merge the duplicate vertices from V, fixing the topology accordingly
- //
- // Input:
- // V,F: mesh description
- // epsilon: minimal distance to consider two vertices identical
- //
- // Output:
- // NV, NF: new mesh without duplicate vertices
-
- template <typename T>
- inline void removeDuplicates(const Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic> &V, const Eigen::MatrixXi &F, Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic> &NV, Eigen::MatrixXi &NF, Eigen::VectorXi &I, const double epsilon = 2.2204e-15);
-
- }
- // Implementation
- template <typename T>
- inline void igl::removeDuplicates(const Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic> &V, const Eigen::MatrixXi &F, Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic> &NV, Eigen::MatrixXi &NF, Eigen::VectorXi &I, const double epsilon)
- {
- assert (F.cols() == 3);
-
- //// build collapse map
- int n = V.rows();
-
- I = Eigen::VectorXi (n);
- I[0] = 0;
-
- bool *VISITED = new bool[n];
- for (int i =0; i <n; ++i)
- VISITED[i] = false;
-
- NV = Eigen::Matrix<T, Eigen::Dynamic, Eigen::Dynamic>(n,V.cols());
- int count = 0;
- Eigen::VectorXd d(n);
- for (int i =0; i <n; ++i)
- {
- if(!VISITED[i])
- {
- NV.row(count) = V.row(i);
- I[i] = count;
- VISITED[i] = true;
- for (int j = i+1; j <n; ++j)
- {
- if((V.row(j) - V.row(i)).norm() < epsilon)
- {
- VISITED[j] = true;
- I[j] = count;
- }
- }
- count ++;
- }
- }
-
- NV.conservativeResize ( count , Eigen::NoChange );
- count = 0;
- NF = Eigen::MatrixXi(F.rows(),3);
- for (int i =0; i <F.rows(); ++i)
- {
- int v0 = I[F(i,0)];
- int v1 = I[F(i,1)];
- int v2 = I[F(i,2)];
- if ( (v0 != v1) && (v2 != v1) && (v0 != v2) )
- NF.row(count++) << v0, v1, v2;
- }
- NF.conservativeResize ( count , Eigen::NoChange );
-
- delete [] VISITED;
- }
- #endif
|