Update tangent_field_operators.h
added functions InitBorderField , SmoothIterative and PropagateFromSelF
This commit is contained in:
nico
parent
ce75b4e68f
commit
ed63e3f7ee
|
|
@ -153,6 +153,8 @@ class CrossField
|
||||||
|
|
||||||
private:
|
private:
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
static void SubDivideDir(const CoordType &Edge0,
|
static void SubDivideDir(const CoordType &Edge0,
|
||||||
const CoordType &Edge1,
|
const CoordType &Edge1,
|
||||||
std::vector<CoordType> &SubDEdges,
|
std::vector<CoordType> &SubDEdges,
|
||||||
|
|
@ -219,9 +221,9 @@ private:
|
||||||
}
|
}
|
||||||
|
|
||||||
static void FindSubDir(vcg::Triangle3<ScalarType> T3,
|
static void FindSubDir(vcg::Triangle3<ScalarType> T3,
|
||||||
size_t Nvert,
|
size_t Nvert,
|
||||||
std::vector<CoordType> &SubDEdges,
|
std::vector<CoordType> &SubDEdges,
|
||||||
int Nsub)
|
int Nsub)
|
||||||
{
|
{
|
||||||
CoordType P0=T3.P0(Nvert);
|
CoordType P0=T3.P0(Nvert);
|
||||||
CoordType P1=T3.P1(Nvert);
|
CoordType P1=T3.P1(Nvert);
|
||||||
|
|
@ -234,9 +236,9 @@ private:
|
||||||
}
|
}
|
||||||
|
|
||||||
static void SubdivideTris(vcg::Triangle3<ScalarType> T3,
|
static void SubdivideTris(vcg::Triangle3<ScalarType> T3,
|
||||||
size_t Nvert,
|
size_t Nvert,
|
||||||
std::vector<vcg::Triangle3<ScalarType> > &SubTris,
|
std::vector<vcg::Triangle3<ScalarType> > &SubTris,
|
||||||
int Nsub)
|
int Nsub)
|
||||||
{
|
{
|
||||||
std::vector<CoordType> SplittedPos;
|
std::vector<CoordType> SplittedPos;
|
||||||
FindSubDir(T3,Nvert,SplittedPos,Nsub);
|
FindSubDir(T3,Nvert,SplittedPos,Nsub);
|
||||||
|
|
@ -244,8 +246,8 @@ private:
|
||||||
//then create the triangles
|
//then create the triangles
|
||||||
for (size_t j=0;j<SplittedPos.size()-1;j++)
|
for (size_t j=0;j<SplittedPos.size()-1;j++)
|
||||||
{
|
{
|
||||||
TriangleType T(P0,SplittedPos[j+1],SplittedPos[j]);
|
TriangleType T(P0,SplittedPos[j+1],SplittedPos[j]);
|
||||||
SubTris.push_back(T);
|
SubTris.push_back(T);
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
@ -328,18 +330,18 @@ private:
|
||||||
alpha=1-(ScalarType)(sub_int+MiddlePos)/(ScalarType)SubDivisionSize;
|
alpha=1-(ScalarType)(sub_int+MiddlePos)/(ScalarType)SubDivisionSize;
|
||||||
}
|
}
|
||||||
else
|
else
|
||||||
if (sub_int>MiddlePos)
|
if (sub_int>MiddlePos)
|
||||||
{
|
{
|
||||||
IndexF0=interval;
|
IndexF0=interval;
|
||||||
IndexF1=(interval+1) % OriginalFace.size();
|
IndexF1=(interval+1) % OriginalFace.size();
|
||||||
alpha=1-(sub_int-MiddlePos)/(ScalarType)SubDivisionSize;
|
alpha=1-(sub_int-MiddlePos)/(ScalarType)SubDivisionSize;
|
||||||
}
|
}
|
||||||
else //sub_int==MiddlePos
|
else //sub_int==MiddlePos
|
||||||
{
|
{
|
||||||
IndexF0=interval;
|
IndexF0=interval;
|
||||||
IndexF1=(interval+1) % OriginalFace.size();
|
IndexF1=(interval+1) % OriginalFace.size();
|
||||||
alpha=1;
|
alpha=1;
|
||||||
}
|
}
|
||||||
assert((IndexF0>=0)&&(IndexF0<(int)OriginalFace.size()));
|
assert((IndexF0>=0)&&(IndexF0<(int)OriginalFace.size()));
|
||||||
assert((IndexF1>=0)&&(IndexF1<(int)OriginalFace.size()));
|
assert((IndexF1>=0)&&(IndexF1<(int)OriginalFace.size()));
|
||||||
|
|
||||||
|
|
@ -396,80 +398,80 @@ private:
|
||||||
CoordType &Interpolated,
|
CoordType &Interpolated,
|
||||||
int &Face)
|
int &Face)
|
||||||
{
|
{
|
||||||
//find smallest edge
|
//find smallest edge
|
||||||
ScalarType smallestE=std::numeric_limits<ScalarType>::max();
|
ScalarType smallestE=std::numeric_limits<ScalarType>::max();
|
||||||
for (int j=0;j<3;j++)
|
for (int j=0;j<3;j++)
|
||||||
{
|
{
|
||||||
ScalarType L0=(t0.P0(j)-t0.P1(j)).Norm();
|
ScalarType L0=(t0.P0(j)-t0.P1(j)).Norm();
|
||||||
ScalarType L1=(t1.P0(j)-t1.P1(j)).Norm();
|
ScalarType L1=(t1.P0(j)-t1.P1(j)).Norm();
|
||||||
if (L0<smallestE) smallestE=L0;
|
if (L0<smallestE) smallestE=L0;
|
||||||
if (L1<smallestE) smallestE=L1;
|
if (L1<smallestE) smallestE=L1;
|
||||||
}
|
}
|
||||||
|
|
||||||
//safety check
|
//safety check
|
||||||
assert(t0.P(0)==t1.P(0));
|
assert(t0.P(0)==t1.P(0));
|
||||||
|
|
||||||
CoordType Origin=t0.P(0);
|
CoordType Origin=t0.P(0);
|
||||||
TriangleType T0Rot(CoordType(0,0,0),t0.P(1)-Origin,t0.P(2)-Origin);
|
TriangleType T0Rot(CoordType(0,0,0),t0.P(1)-Origin,t0.P(2)-Origin);
|
||||||
TriangleType T1Rot(CoordType(0,0,0),t1.P(1)-Origin,t1.P(2)-Origin);
|
TriangleType T1Rot(CoordType(0,0,0),t1.P(1)-Origin,t1.P(2)-Origin);
|
||||||
|
|
||||||
//then rotate normal of T0 to match with normal of T1
|
//then rotate normal of T0 to match with normal of T1
|
||||||
CoordType N0=vcg::Normal(T0Rot.cP(0),T0Rot.cP(1),T0Rot.P(2));
|
CoordType N0=vcg::Normal(T0Rot.cP(0),T0Rot.cP(1),T0Rot.P(2));
|
||||||
CoordType N1=vcg::Normal(T1Rot.cP(0),T1Rot.cP(1),T1Rot.cP(2));
|
CoordType N1=vcg::Normal(T1Rot.cP(0),T1Rot.cP(1),T1Rot.cP(2));
|
||||||
N0.Normalize();
|
N0.Normalize();
|
||||||
N1.Normalize();
|
N1.Normalize();
|
||||||
vcg::Matrix33<ScalarType> rotation=vcg::RotationMatrix(N0,N1);
|
vcg::Matrix33<ScalarType> rotation=vcg::RotationMatrix(N0,N1);
|
||||||
|
|
||||||
//transform the first triangle
|
//transform the first triangle
|
||||||
T0Rot.P(1)=rotation*T0Rot.P(1);
|
T0Rot.P(1)=rotation*T0Rot.P(1);
|
||||||
T0Rot.P(2)=rotation*T0Rot.P(2);
|
T0Rot.P(2)=rotation*T0Rot.P(2);
|
||||||
|
|
||||||
//also rotate the directions
|
//also rotate the directions
|
||||||
CoordType Dir0Rot=rotation*Dir0;
|
CoordType Dir0Rot=rotation*Dir0;
|
||||||
CoordType Dir1Rot=Dir1;
|
CoordType Dir1Rot=Dir1;
|
||||||
CoordType Sep0Rot=rotation*Sep0;
|
CoordType Sep0Rot=rotation*Sep0;
|
||||||
CoordType Sep1Rot=Sep1;
|
CoordType Sep1Rot=Sep1;
|
||||||
|
|
||||||
//find the interpolation angles
|
//find the interpolation angles
|
||||||
ScalarType Angle0=vcg::Angle(Dir0Rot,Sep0Rot);
|
ScalarType Angle0=vcg::Angle(Dir0Rot,Sep0Rot);
|
||||||
ScalarType Angle1=vcg::Angle(Dir1Rot,Sep1Rot);
|
ScalarType Angle1=vcg::Angle(Dir1Rot,Sep1Rot);
|
||||||
assert(Angle0>=0);
|
assert(Angle0>=0);
|
||||||
assert(Angle1>=0);
|
assert(Angle1>=0);
|
||||||
assert(Angle0<=(M_PI/2));
|
assert(Angle0<=(M_PI/2));
|
||||||
assert(Angle1<=(M_PI/2));
|
assert(Angle1<=(M_PI/2));
|
||||||
ScalarType alpha=0.5;//Angle0/(Angle0+Angle1);
|
ScalarType alpha=0.5;//Angle0/(Angle0+Angle1);
|
||||||
|
|
||||||
//then interpolate the direction
|
//then interpolate the direction
|
||||||
//CoordType Interp=Dir0Rot*(1-alpha)+Dir1Rot*alpha;
|
//CoordType Interp=Dir0Rot*(1-alpha)+Dir1Rot*alpha;
|
||||||
Interpolated=Sep0Rot*(1-alpha)+Sep1Rot*alpha;
|
Interpolated=Sep0Rot*(1-alpha)+Sep1Rot*alpha;
|
||||||
Interpolated.Normalize();
|
Interpolated.Normalize();
|
||||||
Interpolated*=smallestE;
|
Interpolated*=smallestE;
|
||||||
|
|
||||||
//then find the triangle which falls into
|
//then find the triangle which falls into
|
||||||
CoordType bary0,bary1;
|
CoordType bary0,bary1;
|
||||||
bool Inside0=vcg::InterpolationParameters(T0Rot,Interpolated,bary0);
|
bool Inside0=vcg::InterpolationParameters(T0Rot,Interpolated,bary0);
|
||||||
bool Inside1=vcg::InterpolationParameters(T1Rot,Interpolated,bary1);
|
bool Inside1=vcg::InterpolationParameters(T1Rot,Interpolated,bary1);
|
||||||
assert(Inside0 || Inside1);
|
assert(Inside0 || Inside1);
|
||||||
// if (!(Inside0 || Inside1))
|
// if (!(Inside0 || Inside1))
|
||||||
// {
|
// {
|
||||||
// std::cout << "Not Inside " << Interpolated.X() << "," << Interpolated.Y() << "," << Interpolated.Z() << std::endl;
|
// std::cout << "Not Inside " << Interpolated.X() << "," << Interpolated.Y() << "," << Interpolated.Z() << std::endl;
|
||||||
// std::cout << "bary0 " << bary0.X() << "," << bary0.Y() << "," << bary0.Z() << std::endl;
|
// std::cout << "bary0 " << bary0.X() << "," << bary0.Y() << "," << bary0.Z() << std::endl;
|
||||||
// std::cout << "bary1 " << bary1.X() << "," << bary1.Y() << "," << bary1.Z() << std::endl;
|
// std::cout << "bary1 " << bary1.X() << "," << bary1.Y() << "," << bary1.Z() << std::endl;
|
||||||
// std::cout << "Diff0 " << fabs(bary0.Norm() - 1) << std::endl;
|
// std::cout << "Diff0 " << fabs(bary0.Norm() - 1) << std::endl;
|
||||||
// std::cout << "Diff1 " << fabs(bary1.Norm() - 1) << std::endl;
|
// std::cout << "Diff1 " << fabs(bary1.Norm() - 1) << std::endl;
|
||||||
// }
|
// }
|
||||||
|
|
||||||
if (Inside0)
|
if (Inside0)
|
||||||
{
|
{
|
||||||
Interpolated=t0.P(0)*bary0.X()+t0.P(1)*bary0.Y()+t0.P(2)*bary0.Z();
|
Interpolated=t0.P(0)*bary0.X()+t0.P(1)*bary0.Y()+t0.P(2)*bary0.Z();
|
||||||
Interpolated-=Origin;
|
Interpolated-=Origin;
|
||||||
Face=0;
|
Face=0;
|
||||||
}
|
}
|
||||||
else
|
else
|
||||||
Face=1;
|
Face=1;
|
||||||
|
|
||||||
//otherwise is already in the tangent space of t0
|
//otherwise is already in the tangent space of t0
|
||||||
Interpolated.Normalize();
|
Interpolated.Normalize();
|
||||||
}
|
}
|
||||||
|
|
||||||
static void ReduceOneDirectionField(std::vector<CoordType> &directions,
|
static void ReduceOneDirectionField(std::vector<CoordType> &directions,
|
||||||
|
|
@ -521,12 +523,122 @@ private:
|
||||||
|
|
||||||
public:
|
public:
|
||||||
|
|
||||||
|
static void InitBorderField(MeshType & mesh)
|
||||||
|
{
|
||||||
|
typedef typename MeshType::FaceType FaceType;
|
||||||
|
typedef typename MeshType::VertexType VertexType;
|
||||||
|
typedef typename MeshType::CoordType CoordType;
|
||||||
|
typedef typename MeshType::ScalarType ScalarType;
|
||||||
|
|
||||||
|
vcg::tri::UpdateTopology<MeshType>::FaceFace(mesh);
|
||||||
|
for (size_t i=0;i<mesh.face.size();i++)
|
||||||
|
for (int j=0;j<mesh.face[i].VN();j++)
|
||||||
|
{
|
||||||
|
FaceType *f0=&mesh.face[i];
|
||||||
|
FaceType *f1=f0->FFp(j);
|
||||||
|
assert(f1!=NULL);
|
||||||
|
if (f0!=f1)continue;
|
||||||
|
|
||||||
|
CoordType dir=f0->P0(j)-f0->P1(j);
|
||||||
|
dir.Normalize();
|
||||||
|
f0->PD1()=dir;
|
||||||
|
f0->PD2()=f0->N()^dir;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
static void SmoothIterative(MeshType &mesh,int NDir=4,
|
||||||
|
int NSteps=3,bool FixSelected=false)
|
||||||
|
{
|
||||||
|
|
||||||
|
typedef typename MeshType::FaceType FaceType;
|
||||||
|
typedef typename MeshType::VertexType VertexType;
|
||||||
|
typedef typename MeshType::CoordType CoordType;
|
||||||
|
typedef typename MeshType::ScalarType ScalarType;
|
||||||
|
|
||||||
|
vcg::tri::UpdateTopology<MeshType>::FaceFace(mesh);
|
||||||
|
|
||||||
|
for (size_t s=0;s<NSteps;s++)
|
||||||
|
{
|
||||||
|
std::vector<CoordType> NewPD1(mesh.face.size(),CoordType(0,0,0));
|
||||||
|
for (size_t i=0;i<mesh.face.size();i++)
|
||||||
|
{
|
||||||
|
if ((FixSelected)&&(mesh.face[i].IsS()))continue;
|
||||||
|
std::vector<CoordType> TangVect;
|
||||||
|
std::vector<CoordType> Norms;
|
||||||
|
FaceType *f0=&mesh.face[i];
|
||||||
|
for (int j=0;j<f0->VN();j++)
|
||||||
|
{
|
||||||
|
FaceType *f1=f0->FFp(j);
|
||||||
|
assert(f1!=NULL);
|
||||||
|
if (f0==f1)continue;
|
||||||
|
TangVect.push_back(f1->PD1());
|
||||||
|
Norms.push_back(f1->N());
|
||||||
|
}
|
||||||
|
assert(Norms.size()>0);
|
||||||
|
std::vector<ScalarType> Weights;
|
||||||
|
Weights.resize(Norms.size(),1/(ScalarType)Norms.size());
|
||||||
|
NewPD1[i]=InterpolateCrossField(TangVect,Weights,Norms,f0->N(),NDir);
|
||||||
|
}
|
||||||
|
for (size_t i=0;i<mesh.face.size();i++)
|
||||||
|
{
|
||||||
|
if ((FixSelected)&&(mesh.face[i].IsS()))continue;
|
||||||
|
assert(NewPD1[i]!=CoordType(0,0,0));
|
||||||
|
mesh.face[i].PD1()=NewPD1[i];
|
||||||
|
mesh.face[i].PD2()=mesh.face[i].N()^mesh.face[i].PD1();
|
||||||
|
}
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
static void PropagateFromSelF(MeshType &mesh)
|
||||||
|
{
|
||||||
|
vcg::tri::UpdateTopology<MeshType>::FaceFace(mesh);
|
||||||
|
|
||||||
|
//typedef typename std::pair<FaceType*,FaceType*> FacePair;
|
||||||
|
std::vector<int> queue;
|
||||||
|
std::vector<int> Sel0;
|
||||||
|
//initialize the queue
|
||||||
|
for (int i=0; i<(int)mesh.face.size(); i++)
|
||||||
|
{
|
||||||
|
FaceType *f=&(mesh.face[i]);
|
||||||
|
if (f->IsD())continue;
|
||||||
|
if (!f->IsS())continue;
|
||||||
|
queue.push_back(i);
|
||||||
|
}
|
||||||
|
Sel0=queue;
|
||||||
|
assert(queue.size()>0);
|
||||||
|
do
|
||||||
|
{
|
||||||
|
std::vector<int> new_queue;
|
||||||
|
for (int i=0; i<queue.size(); i++)
|
||||||
|
{
|
||||||
|
FaceType *f0=&(mesh.face[queue[i]]);
|
||||||
|
assert(!f0->IsD());
|
||||||
|
for (int i=0;i<f0->VN();i++)
|
||||||
|
{
|
||||||
|
FaceType *f1=f0->FFp(i);
|
||||||
|
if (f1==f0)continue;
|
||||||
|
if (f1->IsS())continue;
|
||||||
|
f1->PD1()=Rotate(*f0,*f1,f0->PD1());
|
||||||
|
f1->PD2()=f1->PD1()^f1->N();
|
||||||
|
f1->SetS();
|
||||||
|
new_queue.push_back(vcg::tri::Index(mesh,f1));
|
||||||
|
}
|
||||||
|
}
|
||||||
|
queue=new_queue;
|
||||||
|
}while (queue.size()>0);
|
||||||
|
|
||||||
|
//restore selected flag
|
||||||
|
vcg::tri::UpdateFlags<MeshType>::FaceClearS(mesh);
|
||||||
|
for (int i=0; i<(int)Sel0.size(); i++)
|
||||||
|
mesh.face[Sel0[i]].SetS();
|
||||||
|
}
|
||||||
|
|
||||||
static size_t FindSeparatrices(const typename vcg::face::Pos<FaceType> &vPos,
|
static size_t FindSeparatrices(const typename vcg::face::Pos<FaceType> &vPos,
|
||||||
std::vector<CoordType> &directions,
|
std::vector<CoordType> &directions,
|
||||||
std::vector<FaceType*> &faces,
|
std::vector<FaceType*> &faces,
|
||||||
std::vector<TriangleType> &WrongTris,
|
std::vector<TriangleType> &WrongTris,
|
||||||
int expVal=-1,
|
int expVal=-1,
|
||||||
int numSub=3)
|
int numSub=3)
|
||||||
{
|
{
|
||||||
|
|
||||||
directions.clear();
|
directions.clear();
|
||||||
|
|
@ -627,7 +739,7 @@ public:
|
||||||
else
|
else
|
||||||
{
|
{
|
||||||
if ((versef0D2 * versef1D2 )< ScalarType(0))
|
if ((versef0D2 * versef1D2 )< ScalarType(0))
|
||||||
InterpolateDir(Dir2F0,Dir2F1,Bary0,Bary1,t0,t1,InterpDir,tri_Index);
|
InterpolateDir(Dir2F0,Dir2F1,Bary0,Bary1,t0,t1,InterpDir,tri_Index);
|
||||||
}
|
}
|
||||||
//no separatrix found continue
|
//no separatrix found continue
|
||||||
if (tri_Index==-1)continue;
|
if (tri_Index==-1)continue;
|
||||||
|
|
@ -842,22 +954,22 @@ public:
|
||||||
const ScalarType &alpha2,
|
const ScalarType &alpha2,
|
||||||
int RefEdge=1)
|
int RefEdge=1)
|
||||||
{
|
{
|
||||||
CoordType axis0=f.cP1(RefEdge)-f.cP0(RefEdge);
|
CoordType axis0=f.cP1(RefEdge)-f.cP0(RefEdge);
|
||||||
axis0.Normalize();
|
axis0.Normalize();
|
||||||
CoordType axis2=f.cN();
|
CoordType axis2=f.cN();
|
||||||
axis2.Normalize();
|
axis2.Normalize();
|
||||||
CoordType axis1=axis2^axis0;
|
CoordType axis1=axis2^axis0;
|
||||||
axis1.Normalize();
|
axis1.Normalize();
|
||||||
|
|
||||||
vcg::Matrix33<ScalarType> Trans=vcg::TransformationMatrix(axis0,axis1,axis2);
|
vcg::Matrix33<ScalarType> Trans=vcg::TransformationMatrix(axis0,axis1,axis2);
|
||||||
vcg::Matrix33<ScalarType> InvTrans=Inverse(Trans);
|
vcg::Matrix33<ScalarType> InvTrans=Inverse(Trans);
|
||||||
|
|
||||||
CoordType PD1=CoordType(cos(alpha1),sin(alpha1),0);
|
CoordType PD1=CoordType(cos(alpha1),sin(alpha1),0);
|
||||||
CoordType PD2=CoordType(cos(alpha2),sin(alpha2),0);
|
CoordType PD2=CoordType(cos(alpha2),sin(alpha2),0);
|
||||||
|
|
||||||
//then transform and store in the face
|
//then transform and store in the face
|
||||||
f.PD1()=(InvTrans*PD1);
|
f.PD1()=(InvTrans*PD1);
|
||||||
f.PD2()=(InvTrans*PD2);
|
f.PD2()=(InvTrans*PD2);
|
||||||
}
|
}
|
||||||
|
|
||||||
///return the 4 directiona of the cross field in 3D
|
///return the 4 directiona of the cross field in 3D
|
||||||
|
|
@ -953,19 +1065,19 @@ public:
|
||||||
|
|
||||||
for (int i=0;i<3;i++)
|
for (int i=0;i<3;i++)
|
||||||
{
|
{
|
||||||
vcg::Matrix33<ScalarType> rotN=vcg::RotationMatrix(f.V(i)->N(),f.N());
|
vcg::Matrix33<ScalarType> rotN=vcg::RotationMatrix(f.V(i)->N(),f.N());
|
||||||
CoordType rotatedDir=rotN*f.V(i)->PD1();
|
CoordType rotatedDir=rotN*f.V(i)->PD1();
|
||||||
|
|
||||||
if (fabs(rotatedDir*tF0)>fabs(rotatedDir*tF1))
|
if (fabs(rotatedDir*tF0)>fabs(rotatedDir*tF1))
|
||||||
{
|
{
|
||||||
mag1+=fabs(f.V(i)->K1());
|
mag1+=fabs(f.V(i)->K1());
|
||||||
mag2+=fabs(f.V(i)->K2());
|
mag2+=fabs(f.V(i)->K2());
|
||||||
}
|
}
|
||||||
else
|
else
|
||||||
{
|
{
|
||||||
mag1+=fabs(f.V(i)->K2());
|
mag1+=fabs(f.V(i)->K2());
|
||||||
mag2+=fabs(f.V(i)->K1());
|
mag2+=fabs(f.V(i)->K1());
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
f.K1()=mag1/(ScalarType)3;
|
f.K1()=mag1/(ScalarType)3;
|
||||||
|
|
@ -1094,10 +1206,11 @@ public:
|
||||||
static CoordType InterpolateCrossField(const std::vector<CoordType> &TangVect,
|
static CoordType InterpolateCrossField(const std::vector<CoordType> &TangVect,
|
||||||
const std::vector<ScalarType> &Weight,
|
const std::vector<ScalarType> &Weight,
|
||||||
const std::vector<CoordType> &Norms,
|
const std::vector<CoordType> &Norms,
|
||||||
const CoordType &BaseNorm)
|
const CoordType &BaseNorm,
|
||||||
|
int NDir=4)
|
||||||
{
|
{
|
||||||
|
|
||||||
CoordType sum=InterpolateNRosy3D(TangVect,Norms,Weight,4,BaseNorm);
|
CoordType sum=InterpolateNRosy3D(TangVect,Norms,Weight,NDir,BaseNorm);
|
||||||
return sum;
|
return sum;
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
@ -1153,23 +1266,23 @@ public:
|
||||||
CoordType t03D=CoordType(t0.X(),t0.Y(),0);
|
CoordType t03D=CoordType(t0.X(),t0.Y(),0);
|
||||||
CoordType t13D=CoordType(t1.X(),t1.Y(),0);
|
CoordType t13D=CoordType(t1.X(),t1.Y(),0);
|
||||||
CoordType Norm=CoordType(0,0,1);
|
CoordType Norm=CoordType(0,0,1);
|
||||||
// CoordType n=CoordType(0,0,1);
|
// CoordType n=CoordType(0,0,1);
|
||||||
// CoordType trans1=K_PI(t13D,t03D,n);
|
// CoordType trans1=K_PI(t13D,t03D,n);
|
||||||
// ScalarType diff=vcg::AngleN(trans0,trans1)/(M_PI_4);
|
// ScalarType diff=vcg::AngleN(trans0,trans1)/(M_PI_4);
|
||||||
//ScalarType diff = 1-fabs(trans0*trans1);
|
//ScalarType diff = 1-fabs(trans0*trans1);
|
||||||
return DifferenceCrossField(t03D,t13D,Norm);
|
return DifferenceCrossField(t03D,t13D,Norm);
|
||||||
}
|
}
|
||||||
|
|
||||||
///return the difference of two cross field, values between [0,1]
|
///return the difference of two cross field, values between [0,1]
|
||||||
static typename FaceType::ScalarType DifferenceLineField(const typename vcg::Point2<ScalarType> &t0,
|
static typename FaceType::ScalarType DifferenceLineField(const typename vcg::Point2<ScalarType> &t0,
|
||||||
const typename vcg::Point2<ScalarType> &t1)
|
const typename vcg::Point2<ScalarType> &t1)
|
||||||
{
|
{
|
||||||
CoordType t03D=CoordType(t0.X(),t0.Y(),0);
|
CoordType t03D=CoordType(t0.X(),t0.Y(),0);
|
||||||
CoordType t13D=CoordType(t1.X(),t1.Y(),0);
|
CoordType t13D=CoordType(t1.X(),t1.Y(),0);
|
||||||
CoordType Norm=CoordType(0,0,1);
|
CoordType Norm=CoordType(0,0,1);
|
||||||
// CoordType n=CoordType(0,0,1);
|
// CoordType n=CoordType(0,0,1);
|
||||||
// CoordType trans1=K_PI(t13D,t03D,n);
|
// CoordType trans1=K_PI(t13D,t03D,n);
|
||||||
// ScalarType diff=vcg::AngleN(trans0,trans1)/(M_PI_4);
|
// ScalarType diff=vcg::AngleN(trans0,trans1)/(M_PI_4);
|
||||||
//ScalarType diff = 1-fabs(trans0*trans1);
|
//ScalarType diff = 1-fabs(trans0*trans1);
|
||||||
return DifferenceLineField(t03D,t13D,Norm);
|
return DifferenceLineField(t03D,t13D,Norm);
|
||||||
}
|
}
|
||||||
|
|
@ -1288,10 +1401,10 @@ public:
|
||||||
|
|
||||||
static bool IsSingular(MeshType &mesh,const VertexType &v)
|
static bool IsSingular(MeshType &mesh,const VertexType &v)
|
||||||
{
|
{
|
||||||
assert(vcg::tri::HasPerVertexAttribute(mesh,std::string("Singular")));
|
assert(vcg::tri::HasPerVertexAttribute(mesh,std::string("Singular")));
|
||||||
typename MeshType::template PerVertexAttributeHandle<bool> Handle_Singular;
|
typename MeshType::template PerVertexAttributeHandle<bool> Handle_Singular;
|
||||||
Handle_Singular = vcg::tri::Allocator<MeshType>::template GetPerVertexAttribute<bool>(mesh,std::string("Singular"));
|
Handle_Singular = vcg::tri::Allocator<MeshType>::template GetPerVertexAttribute<bool>(mesh,std::string("Singular"));
|
||||||
return (Handle_Singular[v]);
|
return (Handle_Singular[v]);
|
||||||
}
|
}
|
||||||
|
|
||||||
static void GradientToCross(const FaceType &f,
|
static void GradientToCross(const FaceType &f,
|
||||||
|
|
@ -1320,8 +1433,8 @@ public:
|
||||||
dirU=PosT0*Bary0.X()+PosT1*Bary0.Y()+PosT2*Bary0.Z();
|
dirU=PosT0*Bary0.X()+PosT1*Bary0.Y()+PosT2*Bary0.Z();
|
||||||
dirV=PosT0*Bary1.X()+PosT1*Bary1.Y()+PosT2*Bary1.Z();
|
dirV=PosT0*Bary1.X()+PosT1*Bary1.Y()+PosT2*Bary1.Z();
|
||||||
|
|
||||||
// dirU-=Origin3D;
|
// dirU-=Origin3D;
|
||||||
// dirV-=Origin3D;
|
// dirV-=Origin3D;
|
||||||
dirU.Normalize();
|
dirU.Normalize();
|
||||||
dirV.Normalize();
|
dirV.Normalize();
|
||||||
//orient coherently
|
//orient coherently
|
||||||
|
|
@ -1329,36 +1442,36 @@ public:
|
||||||
CoordType NTarget=vcg::Normal(f.cP(0),f.cP(1),f.cP(2));
|
CoordType NTarget=vcg::Normal(f.cP(0),f.cP(1),f.cP(2));
|
||||||
if ((Ntest*NTarget)<0)dirV=-dirV;
|
if ((Ntest*NTarget)<0)dirV=-dirV;
|
||||||
|
|
||||||
// //then make them orthogonal
|
// //then make them orthogonal
|
||||||
// CoordType dirAvg=dirU^dirV;
|
// CoordType dirAvg=dirU^dirV;
|
||||||
CoordType dirVTarget=NTarget^dirU;
|
CoordType dirVTarget=NTarget^dirU;
|
||||||
CoordType dirUTarget=NTarget^dirV;
|
CoordType dirUTarget=NTarget^dirV;
|
||||||
|
|
||||||
dirUTarget.Normalize();
|
dirUTarget.Normalize();
|
||||||
dirVTarget.Normalize();
|
dirVTarget.Normalize();
|
||||||
if ((dirUTarget*dirU)<0)dirUTarget=-dirUTarget;
|
if ((dirUTarget*dirU)<0)dirUTarget=-dirUTarget;
|
||||||
if ((dirVTarget*dirV)<0)dirVTarget=-dirVTarget;
|
if ((dirVTarget*dirV)<0)dirVTarget=-dirVTarget;
|
||||||
|
|
||||||
dirU=(dirU+dirUTarget)/2;
|
dirU=(dirU+dirUTarget)/2;
|
||||||
dirV=(dirV+dirVTarget)/2;
|
dirV=(dirV+dirVTarget)/2;
|
||||||
|
|
||||||
dirU.Normalize();
|
dirU.Normalize();
|
||||||
dirV.Normalize();
|
dirV.Normalize();
|
||||||
|
|
||||||
// ///compute non normalized normal
|
// ///compute non normalized normal
|
||||||
// CoordType n = f.cN();
|
// CoordType n = f.cN();
|
||||||
|
|
||||||
// CoordType p0 =f.cP(1) - f.cP(0);
|
// CoordType p0 =f.cP(1) - f.cP(0);
|
||||||
// CoordType p1 =f.cP(2) - f.cP(1);
|
// CoordType p1 =f.cP(2) - f.cP(1);
|
||||||
// CoordType p2 =f.cP(0) - f.cP(2);
|
// CoordType p2 =f.cP(0) - f.cP(2);
|
||||||
|
|
||||||
// CoordType t[3];
|
// CoordType t[3];
|
||||||
// t[0] = -(p0 ^ n);
|
// t[0] = -(p0 ^ n);
|
||||||
// t[1] = -(p1 ^ n);
|
// t[1] = -(p1 ^ n);
|
||||||
// t[2] = -(p2 ^ n);
|
// t[2] = -(p2 ^ n);
|
||||||
|
|
||||||
// dirU = t[1]*UV0.X() + t[2]*UV1.X() + t[0]*UV2.X();
|
// dirU = t[1]*UV0.X() + t[2]*UV1.X() + t[0]*UV2.X();
|
||||||
// dirV = t[1]*UV0.Y() + t[2]*UV1.Y() + t[0]*UV2.Y();
|
// dirV = t[1]*UV0.Y() + t[2]*UV1.Y() + t[0]*UV2.Y();
|
||||||
}
|
}
|
||||||
|
|
||||||
static void MakeDirectionFaceCoherent(FaceType *f0,
|
static void MakeDirectionFaceCoherent(FaceType *f0,
|
||||||
|
|
|
||||||
Loading…
Reference in New Issue