lucaspalomodevelop/meshlab
1
0
Fork 0
mirror of https://github.com/lucaspalomodevelop/meshlab.git synced 2026-03-20 03:16:10 +00:00
Code Issues Packages Projects Releases Wiki Activity

Added geometric measures for quad edges (right angle discrepancy and edge min/max ratio

Browse Source
This commit is contained in:
Paolo Cignoni cignoni 2011-09-08 22:31:37 +00:00
parent 57bcedf29c
commit 4fc5cbcb48
1 changed files with 46 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

48
src/meshlabplugins/filter_measure/filter_measure.cpp
View File

@ -82,7 +82,7 @@ FilterMeasurePlugin::FilterMeasurePlugin()
{
switch(filterId) {
case FP_MEASURE_TOPO : return QString("Compute Topological Measures");
case FP_MEASURE_TOPO_QUAD : return QString("Compute Topological Measures for Quad Meshes");
case FP_MEASURE_TOPO_QUAD : return QString("Compute Measures for Quad Meshes");
case FP_MEASURE_GEOM : return QString("Compute Geometric Measures");
case FP_MEASURE_GAUSSCURV : return QString("Compute Integral of Gaussian Curvature");
case FP_MEASURE_VERTEX_QUALITY_DISTRIBUTION : return QString("Per Vertex Quality Stat");
@ -240,7 +240,51 @@ bool FilterMeasurePlugin::applyFilter(QAction *filter, MeshDocument &md, RichPar
int nTris = tri::Clean<CMeshO>::CountBitTris(m);
int nPolys = tri::Clean<CMeshO>::CountBitPolygons(m);
Log("Mesh has %i tri %i quad and %i polig",nTris,nQuads,nPolys);
Log("Mesh has %i tri %i quad and %i polig",nTris,nQuads,nPolys);
//
// i
//
//
// i+1 i+2
tri::UpdateFlags<CMeshO>::FaceClearV(m);
Distribution<float> ad; // angle distributio
Distribution<float> rd; // ratio distribution
for(CMeshO::FaceIterator fi=m.face.begin();fi!=m.face.end();++fi)
{
// Collect the vertices
Point3f qv[4];
for(int i=0;i<3;++i)
{
if(!(*fi).IsF(i) && !(*fi).IsF((i+1)%3))
{
qv[0]= (*fi).cP0(i);
qv[1]= (*fi).cP1(i);
qv[2]= (*fi).cP2(i);
face::Pos<CFaceO> pp(&*fi,(i+2)%3,fi->V(i));
assert(pp.F()->IsF(pp.E()));
pp.FlipF();pp.FlipE();pp.FlipV();
qv[3]= pp.V()->cP();
break;
}
for(int i=0;i<4;++i)
ad.Add(fabs(90-math::ToDeg(Angle(qv[(i+0)%4] - qv[(i+1)%4], qv[(i+2)%4] - qv[(i+1)%4]))));
float edgeLen[4];
for(int i=0;i<4;++i)
edgeLen[i]=Distance(qv[(i+0)%4],qv[(i+1)%4]);
std::sort(edgeLen,edgeLen+4);
rd.Add(edgeLen[0]/edgeLen[3]);
}
}
Log("Right Angle Discrepancy Avg %4.3f Min %4.3f Max %4.3f StdDev %4.3f Percentile 0.05 %4.3f percentile 95 %4.3f",
ad.Avg(), ad.Min(), ad.Max(),ad.StandardDeviation(),ad.Percentile(0.05),ad.Percentile(0.95));
Log("Quad Ratio Avg %4.3f Min %4.3f Max %4.3f", rd.Avg(), rd.Min(), rd.Max());
}
break;