first release version
This commit is contained in:
parent
a25d0197a9
commit
fcf300015c
|
|
@ -0,0 +1,73 @@
|
|||
/****************************************************************************
|
||||
* VCGLib o o *
|
||||
* Visual and Computer Graphics Library o o *
|
||||
* _ O _ *
|
||||
* Copyright(C) 2004 \/)\/ *
|
||||
* Visual Computing Lab /\/| *
|
||||
* ISTI - Italian National Research Council | *
|
||||
* \ *
|
||||
* All rights reserved. *
|
||||
* *
|
||||
* This program is free software; you can redistribute it and/or modify *
|
||||
* it under the terms of the GNU General Public License as published by *
|
||||
* the Free Software Foundation; either version 2 of the License, or *
|
||||
* (at your option) any later version. *
|
||||
* *
|
||||
* This program is distributed in the hope that it will be useful, *
|
||||
* but WITHOUT ANY WARRANTY; without even the implied warranty of *
|
||||
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
|
||||
* GNU General Public License (http://www.gnu.org/licenses/gpl.txt) *
|
||||
* for more details. *
|
||||
* *
|
||||
****************************************************************************/
|
||||
#ifndef __VCGLIB_EXPORTERFIELD
|
||||
#define __VCGLIB_EXPORTERFIELD
|
||||
|
||||
namespace vcg {
|
||||
namespace tri {
|
||||
namespace io {
|
||||
|
||||
/**
|
||||
This class encapsulate a filter for saving field formats
|
||||
*/
|
||||
template <class MeshType>
|
||||
class ExporterFIELD
|
||||
{
|
||||
typedef typename MeshType::ScalarType ScalarType;
|
||||
typedef typename MeshType::FaceType FaceType;
|
||||
typedef typename MeshType::VertexType VertexType;
|
||||
typedef typename MeshType::CoordType CoordType;
|
||||
|
||||
public:
|
||||
|
||||
///Save a 4 rosy format file as used by
|
||||
///Interactive Visualization of Rotational Symmetry Fields on Surfaces
|
||||
///Jonathan Palacios and Eugene Zhang
|
||||
static void Save4ROSY(MeshType &mesh,
|
||||
const char *path)
|
||||
{
|
||||
FILE *f = fopen(path,"wt");
|
||||
fprintf(f,"%d\n",mesh.vn);
|
||||
fprintf(f,"4\n");
|
||||
for (unsigned int i=0;i<mesh.vert.size();i++)
|
||||
{
|
||||
float dirX=(float)mesh.vert[i].PD1().X();
|
||||
float dirY=(float)mesh.vert[i].PD1().Y();
|
||||
float dirZ=(float)mesh.vert[i].PD1().Z();
|
||||
fprintf(f,"%f %f %f \n",dirX,dirY,dirZ);
|
||||
|
||||
}
|
||||
fclose(f);
|
||||
}
|
||||
|
||||
|
||||
}; // end class
|
||||
|
||||
|
||||
|
||||
} // end namespace tri
|
||||
} // end namespace io
|
||||
} // end namespace vcg
|
||||
|
||||
#endif
|
||||
|
||||
|
|
@ -0,0 +1,275 @@
|
|||
/****************************************************************************
|
||||
* VCGLib o o *
|
||||
* Visual and Computer Graphics Library o o *
|
||||
* _ O _ *
|
||||
* Copyright(C) 2004 \/)\/ *
|
||||
* Visual Computing Lab /\/| *
|
||||
* ISTI - Italian National Research Council | *
|
||||
* \ *
|
||||
* All rights reserved. *
|
||||
* *
|
||||
* This program is free software; you can redistribute it and/or modify *
|
||||
* it under the terms of the GNU General Public License as published by *
|
||||
* the Free Software Foundation; either version 2 of the License, or *
|
||||
* (at your option) any later version. *
|
||||
* *
|
||||
* This program is distributed in the hope that it will be useful, *
|
||||
* but WITHOUT ANY WARRANTY; without even the implied warranty of *
|
||||
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
|
||||
* GNU General Public License (http://www.gnu.org/licenses/gpl.txt) *
|
||||
* for more details. *
|
||||
* *
|
||||
****************************************************************************/
|
||||
#ifndef __VCGLIB_IMPORTERFIELD
|
||||
#define __VCGLIB_IMPORTERFIELD
|
||||
#include <vcg/complex/algorithms/parametrization/tangent_field_operators.h>
|
||||
|
||||
namespace vcg {
|
||||
namespace tri {
|
||||
namespace io {
|
||||
|
||||
/**
|
||||
This class encapsulate a filter for opening field formats
|
||||
*/
|
||||
template <class MeshType>
|
||||
class ImporterFIELD
|
||||
{
|
||||
typedef typename MeshType::ScalarType ScalarType;
|
||||
typedef typename MeshType::FaceType FaceType;
|
||||
typedef typename MeshType::VertexType VertexType;
|
||||
typedef typename MeshType::CoordType CoordType;
|
||||
|
||||
public:
|
||||
|
||||
static bool LoadGrad(MeshType *mesh,
|
||||
const char *path)
|
||||
{
|
||||
FILE *f = fopen(path,"rt");
|
||||
if (!f)
|
||||
{
|
||||
return false;
|
||||
}
|
||||
int numF;
|
||||
fscanf(f,"%d\n",&numF);
|
||||
assert(numF==mesh->fn);
|
||||
char skipstr[200];
|
||||
//int readed0;
|
||||
for (int i=0;i<9;i++)
|
||||
fscanf(f,"%s",&skipstr[0]);
|
||||
|
||||
for (int i=0;i<mesh->fn;i++)
|
||||
{
|
||||
int i0=-1;
|
||||
int i1=-1;
|
||||
int i2=-1;
|
||||
double u0,v0,u1,v1,u2,v2;
|
||||
int readed1=fscanf(f,"%d %d %d %lf %lf %lf %lf %lf %lf",&i0,&i1,&i2,&u0,&v0,&u1,&v1,&u2,&v2);
|
||||
assert(readed1==9);
|
||||
vcg::Point2<ScalarType> UV[3];
|
||||
UV[0]= vcg::Point2<ScalarType>(u0,v0);
|
||||
UV[1]= vcg::Point2<ScalarType>(u1,v1);
|
||||
UV[2]= vcg::Point2<ScalarType>(u2,v2);
|
||||
CoordType dir1;
|
||||
CoordType dir2;
|
||||
vcg::tri::CrossField<MeshType>::GradientToCross(mesh->face[i],UV[0],UV[1],UV[2],dir1,dir2);
|
||||
dir1.Normalize();
|
||||
dir2.Normalize();
|
||||
mesh->face[i].PD1()=dir1;
|
||||
mesh->face[i].PD2()=dir2;
|
||||
}
|
||||
fclose(f);
|
||||
return true;
|
||||
}
|
||||
|
||||
///load a field on the mesh, it could be a vfield file (per vertex)
|
||||
///or an ffield file (per face)
|
||||
static bool LoadFFIELD(MeshType *mesh,
|
||||
const char *path,
|
||||
bool per_vertex=false)
|
||||
{
|
||||
|
||||
FILE *f = fopen(path,"rt");
|
||||
if (!f)
|
||||
{
|
||||
return false;
|
||||
}
|
||||
{
|
||||
char word[512]; word[0]=0;
|
||||
fscanf(f,"%s",word);
|
||||
char c=0;
|
||||
if (word[0]=='#') {
|
||||
// skip comment line
|
||||
while (fscanf(f,"%c",&c)!=EOF) if (c=='\n') break;
|
||||
}
|
||||
else
|
||||
{
|
||||
return false;
|
||||
}
|
||||
int nnv = -1;
|
||||
if (fscanf(f,"%d",&nnv)!=1)
|
||||
{
|
||||
while (fscanf(f,"%c",&c)!=EOF) if (c=='\n') break; // skip
|
||||
fscanf(f,"%d",&nnv);
|
||||
}
|
||||
int targetnum=mesh->fn;
|
||||
if (per_vertex)
|
||||
targetnum=mesh->vn;
|
||||
if (nnv != (int)targetnum)
|
||||
{
|
||||
//if (errorMsg) sprintf(errorMsg,"Wrong element number. Found: %d. Expected: %d.",nnv,mesh->vn);
|
||||
return false;
|
||||
}
|
||||
while (fscanf(f,"%c",&c)!=EOF) if (c=='\n') break; // skip
|
||||
// skip strange string line
|
||||
while (fscanf(f,"%c",&c)!=EOF) if (c=='\n') break;
|
||||
for (int i=0; i<nnv; i++){
|
||||
vcg::Point3<float> u,v;
|
||||
float a,b;
|
||||
if (fscanf(f,
|
||||
"%f %f %f %f %f %f %f %f",
|
||||
&a,&b,
|
||||
&(v.X()),&(v.Y()),&(v.Z()),
|
||||
&(u.X()),&(u.Y()),&(u.Z())
|
||||
)!=8) {
|
||||
//if (errorMsg) sprintf(errorMsg,"Format error reading vertex n. %d",i);
|
||||
return false;
|
||||
}
|
||||
|
||||
u.Normalize();
|
||||
v.Normalize();
|
||||
|
||||
if (per_vertex)
|
||||
{
|
||||
mesh->vert[i].PD1().X()=(ScalarType) u.X();
|
||||
mesh->vert[i].PD1().Y()=(ScalarType) u.Y();
|
||||
mesh->vert[i].PD1().Z()=(ScalarType) u.Z();
|
||||
mesh->vert[i].PD2().X()=(ScalarType) v.X();
|
||||
mesh->vert[i].PD2().Y()=(ScalarType) v.Y();
|
||||
mesh->vert[i].PD2().Z()=(ScalarType) v.Z();
|
||||
}
|
||||
else
|
||||
{
|
||||
mesh->face[i].PD1().X()=(ScalarType) u.X();
|
||||
mesh->face[i].PD1().Y()=(ScalarType) u.Y();
|
||||
mesh->face[i].PD1().Z()=(ScalarType) u.Z();
|
||||
mesh->face[i].PD2().X()=(ScalarType) v.X();
|
||||
mesh->face[i].PD2().Y()=(ScalarType) v.Y();
|
||||
mesh->face[i].PD2().Z()=(ScalarType) v.Z();
|
||||
mesh->face[i].PD1().Normalize();
|
||||
mesh->face[i].PD2().Normalize();
|
||||
}
|
||||
}
|
||||
}
|
||||
fclose(f);
|
||||
return true;
|
||||
}
|
||||
|
||||
///Load a 4 rosy format file as used by
|
||||
///Interactive Visualization of Rotational Symmetry Fields on Surfaces
|
||||
///Jonathan Palacios and Eugene Zhang
|
||||
static void Load4ROSY(MeshType &mesh,
|
||||
const char *path)
|
||||
{
|
||||
FILE *f = fopen(path,"rt");
|
||||
int num,symm;
|
||||
fscanf(f,"%d",&num);
|
||||
assert(num==mesh.vn);
|
||||
fscanf(f,"%d\n",&symm);
|
||||
assert(symm==4);
|
||||
for (unsigned int i=0;i<num;i++)
|
||||
{
|
||||
float dirX,dirY,dirZ;
|
||||
fscanf(f,"%f %f %f \n",&dirX,&dirY,&dirZ);
|
||||
mesh.vert[i].PD1()=CoordType(dirX,dirY,dirZ);
|
||||
mesh.vert[i].PD2()=mesh.vert[i].PD1()^mesh.vert[i].N();
|
||||
mesh.vert[i].PD1().Normalize();
|
||||
mesh.vert[i].PD2().Normalize();
|
||||
}
|
||||
fclose(f);
|
||||
}
|
||||
|
||||
|
||||
static bool LoadSeamsMMFromOBJ(MeshType &mesh,std::string PathOBJ)
|
||||
{
|
||||
///per face per edge of mmatch in the solver
|
||||
typename MeshType::template PerFaceAttributeHandle<vcg::Point3i> Handle_MMatch;
|
||||
///seam per face
|
||||
typename MeshType::template PerFaceAttributeHandle<vcg::Point3<bool> > Handle_Seams;
|
||||
|
||||
bool HasHandleMMatch=vcg::tri::HasPerFaceAttribute(mesh,std::string("MissMatch"));
|
||||
if (!HasHandleMMatch)
|
||||
Handle_MMatch = vcg::tri::Allocator<MeshType>::template AddPerFaceAttribute<vcg::Point3i>(mesh,std::string("MissMatch"));
|
||||
else
|
||||
Handle_MMatch = vcg::tri::Allocator<MeshType>::template GetPerFaceAttribute<vcg::Point3i>(mesh,std::string("MissMatch"));
|
||||
|
||||
bool HasHandleSeams=vcg::tri::HasPerFaceAttribute(mesh,std::string("Seams"));
|
||||
if (!HasHandleSeams)
|
||||
Handle_Seams=vcg::tri::Allocator<MeshType>::template AddPerFaceAttribute<vcg::Point3<bool> >(mesh,std::string("Seams"));
|
||||
else
|
||||
Handle_Seams=vcg::tri::Allocator<MeshType>::template GetPerFaceAttribute<vcg::Point3<bool> >(mesh,std::string("Seams"));
|
||||
|
||||
FILE *f = fopen(PathOBJ.c_str(),"rt");
|
||||
if (!f)
|
||||
return false;
|
||||
|
||||
for (unsigned int i=0;i<mesh.face.size();i++)
|
||||
{
|
||||
for (int j=0;j<3;j++)
|
||||
{
|
||||
Handle_Seams[i][j]=false;
|
||||
Handle_MMatch[i][j]=0;
|
||||
}
|
||||
}
|
||||
|
||||
while (!feof(f))
|
||||
{
|
||||
|
||||
int f_int,v_int,rot;
|
||||
int readed=fscanf(f,"sm %d %d %d\n",&f_int,&v_int,&rot);
|
||||
///skip lines
|
||||
if (readed==0)
|
||||
{
|
||||
char buff[200];
|
||||
fscanf(f,"%s\n",&buff[0]);
|
||||
}
|
||||
else ///add the actual seams
|
||||
{
|
||||
VertexType *v=&mesh.vert[v_int-1];
|
||||
FaceType *f0=&mesh.face[f_int-1];
|
||||
int e0=-1;
|
||||
if (f0->V(0)==v)e0=0;
|
||||
if (f0->V(1)==v)e0=1;
|
||||
if (f0->V(2)==v)e0=2;
|
||||
e0=(e0+2)%3;
|
||||
assert(e0!=-1);
|
||||
FaceType *f1;
|
||||
int e1;
|
||||
f1=f0->FFp(e0);
|
||||
e1=f0->FFi(e0);
|
||||
Handle_Seams[f0][e0]=true;
|
||||
Handle_Seams[f1][e1]=true;
|
||||
|
||||
Handle_MMatch[f0][e0]=rot;
|
||||
int rot1;
|
||||
if (rot==0)rot1=0;
|
||||
if (rot==1)rot1=3;
|
||||
if (rot==2)rot1=2;
|
||||
if (rot==3)rot1=1;
|
||||
Handle_MMatch[f1][e1]=rot1;
|
||||
}
|
||||
}
|
||||
//printf("NEED %d LINES\n",i);
|
||||
return true;
|
||||
}
|
||||
|
||||
|
||||
}; // end class
|
||||
|
||||
|
||||
|
||||
} // end namespace tri
|
||||
} // end namespace io
|
||||
} // end namespace vcg
|
||||
|
||||
#endif
|
||||
|
||||
Loading…
Reference in New Issue