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

first release version

Browse Source
This commit is contained in:
Nico Pietroni nicopietroni 2010-05-26 09:03:19 +00:00
parent 22c81d3dda
commit 4a7d71e7e7
1 changed files with 135 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

135
src/meshlabplugins/filter_isoparametrization/iso_transfer.h Normal file
View File

@ -0,0 +1,135 @@
#include <iso_parametrization.h>
#include<vcg/simplex/edge/base.h>
#include<vcg/simplex/vertex/base.h>
#include<vcg/simplex/face/base.h>
#include <vcg/complex/trimesh/base.h>
#include <vcg/complex/trimesh/update/topology.h>
#include <vcg/complex/trimesh/update/edges.h>
#include <vcg/complex/trimesh/update/bounding.h>
#include <vcg/complex/trimesh/update/flag.h>
#include <vcg/space/index/grid_static_ptr.h>
#include <vcg/complex/trimesh/closest.h>
#ifndef _ISO_TRANSFER
#define _ISO_TRANSFER
class IsoTransfer
{
typedef vcg::GridStaticPtr<ParamMesh::FaceType, ParamMesh::ScalarType> TriMeshGrid;
typedef ParamMesh::CoordType CoordType;
typedef ParamMesh::ScalarType ScalarType;
TriMeshGrid TRGrid;
void Clamp(CoordType &bary)
{
/* float eps=0.01;*/
float sum=0;
int bigger=0;
int lower=0;
/* for (int i=0;i<3;i++)
{
if ((bary.V(i)<eps)&&(bary.V(i)>-eps))
bary.V(i)=0;
if ((bary.V(i)<1+eps)&&(bary.V(i)>1-eps))
bary.V(i)=1;
sum+=bary.V(i);
if (bary.V(i)>bary.V(bigger))
bigger=i;
if (bary.V(i)<bary.V(lower))
lower=i;
}
assert(bigger!=lower);
if (sum>(1.0+eps))
{
float diff=sum-1.0;
bary.V(bigger)-=diff;
}
else
if (sum<(1.0-eps))
{
float diff=1.0-sum;
bary.V(lower)+=diff;
}*/
for (int i=0;i<3;i++)
{
if (bary.V(i)<0)
bary.V(i)=0;
if (bary.V(i)>1)
bary.V(i)=1;
sum+=bary.V(i);
if (bary.V(i)>bary.V(bigger))
bigger=i;
if (bary.V(i)<bary.V(lower))
lower=i;
}
//assert(bigger!=lower);
if (sum>(1.0))
{
float diff=sum-1.0;
bary.V(bigger)-=diff;
}
else
if (sum<(1.0))
{
float diff=1.0-sum;
bary.V(lower)+=diff;
}
}
public:
template <class MeshType>
void Transfer(IsoParametrization &IsoParam,
MeshType &to_assing)
{
///put the mesh in the grid
typedef typename MeshType::ScalarType ScalarType;
vcg::tri::UpdateBounding<ParamMesh>::Box(*IsoParam.ParaMesh());
vcg::tri::UpdateNormals<ParamMesh>::PerFaceNormalized(*IsoParam.ParaMesh());
vcg::tri::UpdateNormals<ParamMesh>::PerVertexAngleWeighted(*IsoParam.ParaMesh());
vcg::tri::UpdateNormals<ParamMesh>::NormalizeVertex(*IsoParam.ParaMesh());
vcg::tri::UpdateEdges<ParamMesh>::Set(*IsoParam.ParaMesh());
vcg::tri::UpdateFlags<ParamMesh>::FaceProjection(*IsoParam.ParaMesh());
TRGrid.Set(IsoParam.ParaMesh()->face.begin(),IsoParam.ParaMesh()->face.end());
ScalarType maxDist=IsoParam.ParaMesh()->bbox.Diag();
///then for each vertex find the closest
for (int i=0;i<to_assing.vert.size();i++)
{
MeshType::VertexType *vert=&to_assing.vert[i];
if (!vert->IsD())
{
ScalarType dist;
CoordType closest,norm,bary;
ParamMesh::FaceType * f=NULL;
f=GetClosestFace(*IsoParam.ParaMesh(),TRGrid,vert->P(), maxDist,dist,closest,norm,bary);
assert(f!=NULL);
///then find back the coordinates
if (!((bary.X()>=0)&&(bary.X()<=1)&&
(bary.Y()>=0)&&(bary.Y()<=1)&&
(bary.Z()>=0)&&(bary.Z()<=1)))
{
printf("%d,%3.3f,%3.3f,%3.3f",i,bary.X(),bary.Y(),bary.Z());
system("pause");
}
Clamp(bary);
int I;
vcg::Point2<ScalarType> UV;
IsoParam.Phi(f,bary,I,UV);
///and finally set to the vertex
assert(I>=0);
vert->T().P()=UV;
vert->T().N()=I;
vert->Q()=(typename MeshType::ScalarType)I;
}
}
}
};
#endif