ModQuadricT.cc
Go to the documentation of this file.00001 /*===========================================================================*\ 00002 * * 00003 * OpenMesh * 00004 * Copyright (C) 2001-2011 by Computer Graphics Group, RWTH Aachen * 00005 * www.openmesh.org * 00006 * * 00007 *---------------------------------------------------------------------------* 00008 * This file is part of OpenMesh. * 00009 * * 00010 * OpenMesh is free software: you can redistribute it and/or modify * 00011 * it under the terms of the GNU Lesser General Public License as * 00012 * published by the Free Software Foundation, either version 3 of * 00013 * the License, or (at your option) any later version with the * 00014 * following exceptions: * 00015 * * 00016 * If other files instantiate templates or use macros * 00017 * or inline functions from this file, or you compile this file and * 00018 * link it with other files to produce an executable, this file does * 00019 * not by itself cause the resulting executable to be covered by the * 00020 * GNU Lesser General Public License. This exception does not however * 00021 * invalidate any other reasons why the executable file might be * 00022 * covered by the GNU Lesser General Public License. * 00023 * * 00024 * OpenMesh is distributed in the hope that it will be useful, * 00025 * but WITHOUT ANY WARRANTY; without even the implied warranty of * 00026 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * 00027 * GNU Lesser General Public License for more details. * 00028 * * 00029 * You should have received a copy of the GNU LesserGeneral Public * 00030 * License along with OpenMesh. If not, * 00031 * see <http://www.gnu.org/licenses/>. * 00032 * * 00033 \*===========================================================================*/ 00034 00035 /*===========================================================================*\ 00036 * * 00037 * $Revision: 362 $ * 00038 * $Date: 2011-01-26 10:21:12 +0100 (Mi, 26 Jan 2011) $ * 00039 * * 00040 \*===========================================================================*/ 00041 00046 //============================================================================= 00047 // 00048 // CLASS ModQuadric - IMPLEMENTATION 00049 // 00050 //============================================================================= 00051 00052 #define OPENMESH_DECIMATER_MODQUADRIC_CC 00053 00054 //== INCLUDES ================================================================= 00055 00056 #include <OpenMesh/Tools/Decimater/ModQuadricT.hh> 00057 00058 00059 //== NAMESPACE =============================================================== 00060 00061 namespace OpenMesh { // BEGIN_NS_OPENMESH 00062 namespace Decimater { // BEGIN_NS_DECIMATER 00063 00064 00065 //== IMPLEMENTATION ========================================================== 00066 00067 00068 template<class DecimaterType> 00069 void 00070 ModQuadricT<DecimaterType>:: 00071 initialize() 00072 { 00073 using Geometry::Quadricd; 00074 // alloc quadrics 00075 if (!quadrics_.is_valid()) 00076 Base::mesh().add_property( quadrics_ ); 00077 00078 // clear quadrics 00079 typename Mesh::VertexIter v_it = Base::mesh().vertices_begin(), 00080 v_end = Base::mesh().vertices_end(); 00081 00082 for (; v_it != v_end; ++v_it) 00083 Base::mesh().property(quadrics_, v_it).clear(); 00084 00085 // calc (normal weighted) quadric 00086 typename Mesh::FaceIter f_it = Base::mesh().faces_begin(), 00087 f_end = Base::mesh().faces_end(); 00088 00089 typename Mesh::FaceVertexIter fv_it; 00090 typename Mesh::VertexHandle vh0, vh1, vh2; 00091 typedef Vec3d Vec3; 00092 double a,b,c,d, area; 00093 00094 for (; f_it != f_end; ++f_it) 00095 { 00096 fv_it = Base::mesh().fv_iter(f_it.handle()); 00097 vh0 = fv_it.handle(); ++fv_it; 00098 vh1 = fv_it.handle(); ++fv_it; 00099 vh2 = fv_it.handle(); 00100 00101 Vec3 v0, v1, v2; 00102 { 00103 using namespace OpenMesh; 00104 00105 v0 = vector_cast<Vec3>(Base::mesh().point(vh0)); 00106 v1 = vector_cast<Vec3>(Base::mesh().point(vh1)); 00107 v2 = vector_cast<Vec3>(Base::mesh().point(vh2)); 00108 } 00109 00110 Vec3 n = (v1-v0) % (v2-v0); 00111 area = n.norm(); 00112 if (area > FLT_MIN) 00113 { 00114 n /= area; 00115 area *= 0.5; 00116 } 00117 00118 a = n[0]; 00119 b = n[1]; 00120 c = n[2]; 00121 d = -(vector_cast<Vec3>(Base::mesh().point(vh0))|n); 00122 00123 Quadricd q(a, b, c, d); 00124 q *= area; 00125 00126 Base::mesh().property(quadrics_, vh0) += q; 00127 Base::mesh().property(quadrics_, vh1) += q; 00128 Base::mesh().property(quadrics_, vh2) += q; 00129 } 00130 } 00131 00132 00133 //============================================================================= 00134 } // END_NS_DECIMATER 00135 } // END_NS_OPENMESH 00136 //=============================================================================