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

Added Absolute Curvature colorize

Browse Source
This commit is contained in:
Paolo Cignoni cignoni 2006-01-20 16:25:39 +00:00
parent c26f0d04e9
commit bf56c1f802
3 changed files with 82 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

62
src/meshlabplugins/meshcolorize/color_curvature.h
View File

@ -23,6 +23,9 @@
/****************************************************************************
History
$Log$
Revision 1.3 2006/01/20 16:25:39 vannini
Added Absolute Curvature colorize
Revision 1.2 2006/01/20 14:46:44 vannini
Code refactoring
Added RMS Curvature colorize
@ -46,6 +49,9 @@ Moved gaussian and mean curvature functions into color_curvature.h
#include "../../meshlab/LogStream.h"
#define DEFAULT_HISTO_FRAC 0.1f
#define DEFAULT_HISTO_RANGE 10000
class histoMinMaxQ
{
public:
@ -98,7 +104,7 @@ namespace vcg
}
}
template<class MESH_TYPE> histoMinMaxQ HK(MESH_TYPE &m, bool computeH, bool useHisto, float histo_frac=0.1f, int histo_range=10000)
template<class MESH_TYPE> histoMinMaxQ HK(MESH_TYPE &m, bool computeH, bool useHisto, float histo_frac=DEFAULT_HISTO_FRAC, int histo_range=DEFAULT_HISTO_RANGE)
{
// Calcola la curvatura gaussiana (K) oppure la media (H) in base a computeH
// e salva il risultato in Q
@ -135,7 +141,6 @@ namespace vcg
{
for(fi=m.face.begin();fi!=m.face.end();++fi) if(!(*fi).IsD())
{
angle0 = math::Abs(Angle( (*fi).P(1)-(*fi).P(0),(*fi).P(2)-(*fi).P(0) ));
angle1 = math::Abs(Angle( (*fi).P(0)-(*fi).P(1),(*fi).P(2)-(*fi).P(1) ));
angle2 = M_PI-(angle0+angle1);
@ -235,8 +240,8 @@ namespace vcg
typename MESH_TYPE::VertexIterator vi;
vcg::Histogram<float> histo;
float histo_frac=0.1f;
int histo_range=10000;
float histo_frac=DEFAULT_HISTO_FRAC;
int histo_range=DEFAULT_HISTO_RANGE;
float minQ = std::numeric_limits<float>::max();
float maxQ = -std::numeric_limits<float>::max();
@ -278,6 +283,55 @@ namespace vcg
}
template<class MESH_TYPE> void ColorAbsolute(MESH_TYPE &m, GLLogStream *log)
{
float *H = new float[m.vn];
int i;
typename MESH_TYPE::VertexIterator vi;
vcg::Histogram<float> histo;
float histo_frac=DEFAULT_HISTO_FRAC;
int histo_range=DEFAULT_HISTO_RANGE;
float minQ = std::numeric_limits<float>::max();
float maxQ = -std::numeric_limits<float>::max();
//Compute H
HK(m,true,false);
i=0;
for(vi=m.vert.begin(); vi!=m.vert.end(); ++vi,++i) if(!(*vi).IsD())
H[i] = (*vi).Q();
//Compute K
HK(m,false,false);
//Compute abs(H+sqrt(H*H-K)) + abs(H-sqrt(H*H-K))
i=0;
for(vi=m.vert.begin(); vi!=m.vert.end(); ++vi,++i) if(!(*vi).IsD())
{
(*vi).Q() = math::Abs(H[i] + math::Sqrt(powf(H[i], 2.0f) - (*vi).Q())) + math::Abs(H[i] - math::Sqrt(powf(H[i], 2.0f) - (*vi).Q()));
if ((*vi).Q() < minQ) minQ = (*vi).Q();
if ((*vi).Q() > maxQ) maxQ = (*vi).Q();
}
histo.SetRange(minQ, maxQ, histo_range);
for(vi=m.vert.begin(); vi!=m.vert.end(); ++vi) if(!(*vi).IsD())
histo.Add((*vi).Q());
minQ = histo.Percentile(histo_frac);
maxQ = histo.Percentile(1.0f - histo_frac);
for(vi=m.vert.begin(); vi!=m.vert.end(); ++vi) if(!(*vi).IsD())
(*vi).Q() = math::Clamp((*vi).Q(), minQ, maxQ);
if (log)
log->Log(GLLogStream::Info, "Absolute Curvature: minQ=%f maxQ=%f", minQ, maxQ);
delete[] H;
}
}
#endif

21
src/meshlabplugins/meshcolorize/meshcolorize.cpp
View File

@ -23,6 +23,9 @@
/****************************************************************************
History
$Log$
Revision 1.18 2006/01/20 16:25:39 vannini
Added Absolute Curvature colorize
Revision 1.17 2006/01/20 14:46:44 vannini
Code refactoring
Added RMS Curvature colorize
@ -87,6 +90,7 @@ ExtraMeshColorizePlugin::ExtraMeshColorizePlugin() {
actionList << new QAction(ST(CP_GAUSSIAN), this);
actionList << new QAction(ST(CP_MEAN), this);
actionList << new QAction(ST(CP_RMS), this);
actionList << new QAction(ST(CP_ABSOLUTE), this);
actionList << new QAction(ST(CP_SELFINTERSECT), this);
actionList << new QAction(ST(CP_BORDER), this);
actionList << new QAction(ST(CP_COLORNM), this);
@ -100,6 +104,8 @@ const QString ExtraMeshColorizePlugin::ST(ColorizeType c) {
return QString("Mean Curvature");
case CP_RMS:
return QString("Root mean square Curvature");
case CP_ABSOLUTE:
return QString("Absolute Curvature");
case CP_SELFINTERSECT:
return QString("Self Intersections");
case CP_BORDER:
@ -127,7 +133,12 @@ const ActionInfo &ExtraMeshColorizePlugin::Info(QAction *action)
if( action->text() == ST(CP_RMS) )
{
ai.Help = tr("Colorize vertex and faces depending on root mean square curvature.");
ai.ShortHelp = tr("Colorize by mean curvature");
ai.ShortHelp = tr("Colorize by root mean square curvature");
}
if( action->text() == ST(CP_ABSOLUTE) )
{
ai.Help = tr("Colorize vertex and faces depending on absolute curvature.");
ai.ShortHelp = tr("Colorize by absolute curvature");
}
if( action->text() == ST(CP_SELFINTERSECT) )
{
@ -186,6 +197,14 @@ void ExtraMeshColorizePlugin::Compute(QAction * mode, MeshModel &m, RenderMode &
return;
}
if(mode->text() == ST(CP_ABSOLUTE))
{
ColorAbsolute<CMeshO>(m.cm, log);
vcg::tri::UpdateColor<CMeshO>::VertexQuality(m.cm);
rm.colorMode = GLW::CMPerVert;
return;
}
if(mode->text() == ST(CP_SELFINTERSECT))
{
vector<CFaceO *> IntersFace;

5
src/meshlabplugins/meshcolorize/meshcolorize.h
View File

@ -23,6 +23,9 @@
/****************************************************************************
History
$Log$
Revision 1.14 2006/01/20 16:25:39 vannini
Added Absolute Curvature colorize
Revision 1.13 2006/01/20 14:46:44 vannini
Code refactoring
Added RMS Curvature colorize
@ -57,7 +60,7 @@ class ExtraMeshColorizePlugin : public QObject, public MeshColorizeInterface
public:
enum ColorizeType {CP_GAUSSIAN,CP_MEAN,CP_RMS,CP_SELFINTERSECT,CP_BORDER,CP_COLORNM};
enum ColorizeType {CP_GAUSSIAN,CP_MEAN,CP_RMS,CP_ABSOLUTE,CP_SELFINTERSECT,CP_BORDER,CP_COLORNM};
const QString ST(ColorizeType c);
ExtraMeshColorizePlugin();