manolas
/
vcglib
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

aggiunto vettore delle face sel, per poi ripulirle....

This commit is contained in:
Paolo Cignoni 2004-11-19 16:12:57 +00:00
parent 75c9f2e194
commit 191168cc47
1 changed files with 10 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
vcg/complex/intersection.h
View File

@ -36,10 +36,10 @@ bool Intersect( GridType & grid,Plane3<ScalarType> plane, vector<typename Grid
for(i = 0 ; i <= grid.siz[axis0]; ++i){ for(i = 0 ; i <= grid.siz[axis0]; ++i){
for(j = 0 ; j <= grid.siz[axis1]; ++j) for(j = 0 ; j <= grid.siz[axis1]; ++j)
{ {
seg.P0()[axis0] = grid.bbox.min[axis0]+ (i+0.1) * grid.voxel[axis0] ; seg.P0()[axis0] = grid.bbox.min[axis0]+ (i+0.01) * grid.voxel[axis0] ;
seg.P1()[axis0] = grid.bbox.min[axis0]+ (i+0.1) * grid.voxel[axis0]; seg.P1()[axis0] = grid.bbox.min[axis0]+ (i+0.01) * grid.voxel[axis0];
seg.P0()[axis1] = grid.bbox.min[axis1]+ (j+0.1) * grid.voxel[axis1]; seg.P0()[axis1] = grid.bbox.min[axis1]+ (j+0.01) * grid.voxel[axis1];
seg.P1()[axis1] = grid.bbox.min[axis1]+ (j+0.1) * grid.voxel[axis1]; seg.P1()[axis1] = grid.bbox.min[axis1]+ (j+0.01) * grid.voxel[axis1];
if ( Intersection(pl,seg,p)) if ( Intersection(pl,seg,p))
{ {
pi[axis] = min(max(0,floor((p[axis ]-grid.bbox.min[axis])/grid.voxel[axis])),grid.siz[axis]); pi[axis] = min(max(0,floor((p[axis ]-grid.bbox.min[axis])/grid.voxel[axis])),grid.siz[axis]);
@ -76,7 +76,8 @@ bool Intersection( TriMeshType & m, Plane3<ScalarType> pl,EdgeMeshType & em,dou
typedef GridStaticPtr<FaceContainer> GridType; typedef GridStaticPtr<FaceContainer> GridType;
EdgeMeshType::VertexIterator vi; EdgeMeshType::VertexIterator vi;
TriMeshType::FaceIterator fi; TriMeshType::FaceIterator fi;
vector<TriMeshType::FaceType*> v;
v.clear();
Intersect(*grid,pl,cells); Intersect(*grid,pl,cells);
Segment3<ScalarType> seg; Segment3<ScalarType> seg;
ave_length = 0.0; ave_length = 0.0;
@ -91,8 +92,10 @@ bool Intersection( TriMeshType & m, Plane3<ScalarType> pl,EdgeMeshType & em,dou
if(!face.IsS()) if(!face.IsS())
{ {
face.SetS(); face.SetS();
v.push_back(&face);
if(vcg::Intersection(pl,face,seg))// intersezione piano triangolo if(vcg::Intersection(pl,face,seg))// intersezione piano triangolo
{ {
face.SetS();
// add to em // add to em
ave_length+=seg.Length(); ave_length+=seg.Length();
vcg::edge::Allocator<EdgeMeshType>::AddEdges(em,1); vcg::edge::Allocator<EdgeMeshType>::AddEdges(em,1);
@ -108,6 +111,8 @@ bool Intersection( TriMeshType & m, Plane3<ScalarType> pl,EdgeMeshType & em,dou
}//end while }//end while
} }
ave_length/=em.en; ave_length/=em.en;
vector<TriMeshType::FaceType*>::iterator v_i;
for(v_i=v.begin(); v_i!=v.end(); ++v_i) (*v_i)->ClearS();
return true; return true;
} }