Egalware/EgtGeomKernel
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

EgtGeomKernel 2.7f2 :

Browse Source 
- Aggiunte funzioni per calcolo di Offset per superfici chiuse TriMesh
- Piccola miglioria alla triangolazione (con SaraP)
- Migliorie per rimozioni TJunction, calcolo delle normali dei triangoli e creazione di una TriMesh a partire da uno ZMap (con SaraP)
- Aggiunte funzioni di SubtractMap e piccole modifiche per estensione dei Box di creazione per gli Zmap.
This commit is contained in:
Riccardo Elitropi
2025-06-16 11:34:23 +02:00
parent 0bffa0039c
commit 626d5b0e51
14 changed files with 811 additions and 81 deletions
Download Patch File Download Diff File Expand all files Collapse all files
SurfTriMeshUtilities.cpp
+37 -31
View File
@@ -138,7 +138,7 @@ SurfTriMesh::FlipTriangles( int nTA, int nTB)
//----------------------------------------------------------------------------
bool
SurfTriMesh::RemoveTJunctions( bool& bModified)
SurfTriMesh::RemoveTJunctions( bool& bModified, double dMinSqDist)
{
bModified = false ;
@@ -147,6 +147,11 @@ SurfTriMesh::RemoveTJunctions( bool& bModified)
// Ciclo sui triangoli della superficie per determinare gli altri vertici sul loro perimetro
for ( int nT = 0 ; nT < int( m_vTria.size()) ; ++ nT) {
// se adiacenze tutte valide, passo al successivo
if ( m_vTria[nT].nIdAdjac[0] != SVT_DEL && m_vTria[nT].nIdAdjac[0] != SVT_NULL &&
m_vTria[nT].nIdAdjac[1] != SVT_DEL && m_vTria[nT].nIdAdjac[1] != SVT_NULL &&
m_vTria[nT].nIdAdjac[2] != SVT_DEL && m_vTria[nT].nIdAdjac[2] != SVT_NULL)
continue ;
// Se il triangolo non è valido, passo al successivo
Triangle3d trTria ;
if ( ! GetTriangle( nT, trTria) || ! trTria.Validate( true))
@@ -190,9 +195,9 @@ SurfTriMesh::RemoveTJunctions( bool& bModified)
if ( ! GetVertex( m_vTria[vNearTria[nI]].nIdVert[nVert], ptVert))
continue ;
double dProj = ( ptVert - ptSegSt) * vtSeg ;
double dOrt = ( ( ptVert - ptSegSt) - dProj * vtSeg).SqLen() ;
if ( dProj > EPS_SMALL && dProj < dSegLen - EPS_SMALL && dOrt < SQ_EPS_TRIA_H)
vVertOtl.emplace_back( m_vTria[vNearTria[nI]].nIdVert[nVert]) ;
double dOrt = ( ( ptVert - ptSegSt) - dProj * vtSeg).SqLen() ;
if ( dProj > EPS_SMALL && dProj < dSegLen - EPS_SMALL && dOrt < dMinSqDist)
vVertOtl.emplace_back( m_vTria[vNearTria[nI]].nIdVert[nVert]) ;
}
}
// Riordino i vertici sul segmento
@@ -338,8 +343,10 @@ ChooseGoodStartPoint( PNTULIST& PointList)
// -------------------------------------------------------------
bool
SurfTriMesh::AdjustLoop( PNTULIST& PointList, double dMaxEdgeLen, int nF, bool& bModif) const
SurfTriMesh::AdjustLoop( PNTULIST& PointList, double dMaxEdgeLen, double dTolAlign, bool& bModif) const
{
// vettore dei loop della faccia adiacente
POLYLINEVECTOR LoopVec ;
// Ciclo sui punti del loop
auto itLast = PointList.begin() ;
for ( auto it = next( itLast) ; it != PointList.end() ; ++ it) {
@@ -347,30 +354,26 @@ SurfTriMesh::AdjustLoop( PNTULIST& PointList, double dMaxEdgeLen, int nF, bool&
// bisogna fermarsi per analizzare il tratto corrente alla ricerca di punti allineati se dal punto corrente
// inizia un tratto adiacente ad un'altra faccia oppure se il punto corrente non verrà eliminato dal loop
// della faccia adiacente
bool bAnalyze = ( itLast->second != it->second) ;
if ( bAnalyze)
LoopVec.clear() ;
if ( ! bAnalyze && itLast->second != - 1) {
// recupero i loop della faccia adiacente
POLYLINEVECTOR LoopVec ;
GetFacetLoops( int( itLast->second), LoopVec) ;
bool bFound = false ;
for ( int i = 0 ; i < int( LoopVec.size()) && ! bFound ; i ++) {
if ( LoopVec.empty())
GetFacetLoops( int( itLast->second), LoopVec) ;
for ( int i = 0 ; i < int( LoopVec.size()) && ! bAnalyze ; i ++) {
const PNTULIST& PointListAdj = LoopVec[i].GetUPointList() ;
for ( auto itAdj = PointListAdj.begin() ; itAdj != prev( PointListAdj.end()) && ! bFound ; ++ itAdj) {
int nSamePoints = 0 ;
for ( auto itAdj = PointListAdj.begin() ; itAdj != prev( PointListAdj.end()) ; ++ itAdj) {
// cerco il punto corrente sul loop della faccia adiacente
if ( itAdj->second == nF && AreSamePointApprox( it->first, itAdj->first)) {
bFound = true ;
auto itPrev = ( itAdj == PointListAdj.begin() ? prev( prev( PointListAdj.end())) : prev( itAdj)) ;
auto itNext = next( itAdj) ;
DistPointLine PointLineDistCalc( itAdj->first, itPrev->first, itNext->first) ;
double dDist ;
if ( PointLineDistCalc.GetDist( dDist) && dDist > 50 * EPS_SMALL)
// se il punto non è allineato al precedente e al successivo non verrà rimosso dal loop della faccia adiacente,
// quindi non deve essere rimosso dal loop corrente altrimenti si genererebbe una T-Junction. Quindi posso
// interrompere e analizzare il tratto trovato fino ad ora
bAnalyze = true ;
if ( AreSamePointApprox( it->first, itAdj->first))
++ nSamePoints ;
if ( ( nSamePoints == 1 && int( PointListAdj.size()) <= 4) || nSamePoints > 1) {
bAnalyze = true ;
break ;
}
}
}
}
}
}
if ( bAnalyze) {
@@ -383,7 +386,7 @@ SurfTriMesh::AdjustLoop( PNTULIST& PointList, double dMaxEdgeLen, int nF, bool&
}
PL.AddUPoint( ++nPar, it->first) ;
// Provo ad eliminare i punti allineati
PL.RemoveAlignedPoints( 50 * EPS_SMALL) ;
PL.RemoveAlignedPoints( dTolAlign) ;
if ( PL.GetPointNbr() < nPar + 1) {
// rimuovo dalla lista dei punti gli eliminati (salto gli estremi)
int nUCurr = 1 ;
@@ -475,7 +478,7 @@ SurfTriMesh::AdjustLoop( PNTULIST& PointList, double dMaxEdgeLen, int nF, bool&
//----------------------------------------------------------------------------
bool
SurfTriMesh::SimplifyFacets( double dMaxEdgeLen, bool bForced)
SurfTriMesh::SimplifyFacets( double dMaxEdgeLen, bool bForced, double dTolAlign)
{
// La trimesh deve essere valida
if ( ! IsValid())
@@ -502,13 +505,17 @@ SurfTriMesh::SimplifyFacets( double dMaxEdgeLen, bool bForced)
// Lista dei punti del loop
PNTULIST& PointList = LoopVec[nL].GetUPointList() ;
// Se il loop è un triangolo, non va modificato
if ( int( PointList.size()) <= 4)
continue ;
// Mi assicuro che il punto iniziale/finale non sia all'interno di un possibile segmento
if ( ! ChooseGoodStartPoint( PointList))
continue ;
// Sistemo il loop
bool bModif = false ;
if ( ! AdjustLoop( PointList, dMaxEdgeLen, nF, bModif))
if ( ! AdjustLoop( PointList, dMaxEdgeLen, dTolAlign, bModif))
return false ;
if ( bModif)
bToRetriangulate = true ;
@@ -536,12 +543,11 @@ SurfTriMesh::SimplifyFacets( double dMaxEdgeLen, bool bForced)
// Eseguo la ritriangolazione della faccia
PNTVECTOR vPt ;
INTVECTOR vTr ;
if ( Triangulate().Make( LoopVec, vPt, vTr)) {
if ( Triangulate().Make( LoopVec, vPt, vTr) && ! vTr.empty())
FacetMap.emplace( nF, make_pair( vPt, vTr)) ;
}
// Se non riesco a triangolare anche solo questa faccia, interrompo tutto
else
return false ;
else
return false ;
}
}