00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057
00058
00059
00060 #ifndef __VCGLIB_TRISUBSET
00061 #define __VCGLIB_TRISUBSET
00062
00063 #include <vcg/complex/trimesh/update/flag.h>
00064
00065 namespace vcg {
00066 namespace tri {
00067
00068
00069 template <class I_MESH_TYPE>
00070 struct InsertedV
00071 {
00072 typedef I_MESH_TYPE IMeshType;
00073 typedef typename IMeshType::VertexPointer VertexPointer;
00074 typedef typename IMeshType::FacePointer FacePointer;
00075
00076 InsertedV(VertexPointer _v, FacePointer _f, int _z)
00077 : v(_v), f(_f), z(_z)
00078 {}
00079
00080 VertexPointer v;
00081 FacePointer f;
00082 int z;
00083
00084 bool operator < (const InsertedV & o) const
00085 {
00086 return (v<o.v);
00087 }
00088
00089 bool operator ==(const InsertedV & o)
00090 {
00091 return (v==o.v);
00092 }
00093
00094 bool operator !=(const InsertedV & o)
00095 {
00096 return (v!=o.v);
00097 }
00098 };
00099
00100
00101
00102
00103
00104
00105
00106 template <class S_MESH_TYPE, class STL_CONT>
00107 void SubSet(S_MESH_TYPE & m, STL_CONT & subSet)
00108 {
00109 std::vector< InsertedV<S_MESH_TYPE> > newVertices;
00110 typename STL_CONT::const_iterator pfi;
00111 typename S_MESH_TYPE::VertexIterator vi;
00112 typename S_MESH_TYPE::FaceIterator fi;
00113 typedef typename S_MESH_TYPE::VertexType S_VertexType;
00114 std::vector<typename S_MESH_TYPE::VertexPointer> redirect;
00115
00116
00117 fi = vcg::tri::Allocator<S_MESH_TYPE>::AddFaces(m,subSet.size());
00118 for(pfi=subSet.begin(); pfi!=subSet.end(); ++pfi)
00119 {
00120 assert(!(*pfi)->IsD());
00121 (*fi).ImportData(**pfi);
00122 for(int ii = 0 ; ii < (*fi).VN(); ++ii)
00123 (*fi).V(ii) = (S_VertexType*)(void*)(*pfi)->V(ii);
00124 ++fi;
00125 }
00126
00127
00128 for(fi=m.face.begin(); fi!=m.face.end(); ++fi)
00129 for(int ii = 0 ; ii < (*fi).VN(); ++ii)
00130 newVertices.push_back(InsertedV<S_MESH_TYPE>((*fi).V(ii), &(*fi),ii));
00131
00132 sort(newVertices.begin(), newVertices.end());
00133
00134 typename std::vector< InsertedV<S_MESH_TYPE> >::iterator curr, next;
00135 int pos=0;
00136 curr=next=newVertices.begin();
00137 while(next!=newVertices.end())
00138 {
00139 if((*curr)!=(*next))
00140 pos++;
00141 (*next).f->V((*next).z)=(typename S_MESH_TYPE::VertexPointer)pos;
00142 curr=next;
00143 next++;
00144 }
00145
00146 typename std::vector< InsertedV<S_MESH_TYPE> >::iterator newE=unique(newVertices.begin(), newVertices.end());
00147
00148 vi = vcg::tri::Allocator<S_MESH_TYPE>::AddVertices(m,newE-newVertices.begin());
00149 for(curr=newVertices.begin(); curr!=newE; ++curr,++vi)
00150 (*vi).ImportData(*((*curr).v));
00151
00152 for(vi=m.vert.begin(); vi!=m.vert.end(); ++vi)
00153 redirect.push_back(&(*vi));
00154
00155 for(fi=m.face.begin(); fi!=m.face.end(); ++fi)
00156 {
00157 for(int ii = 0 ; ii < (*fi).VN(); ++ii)
00158 (*fi).V(ii)=redirect[(size_t)(*fi).V(ii)];
00159 }
00160 m.vn=(int)m.vert.size();
00161 m.fn=(int)m.face.size();
00162 }
00163
00164
00165 }
00166 }
00167
00168
00169 #endif
00170
00171