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

draft version. For the triangle is not tehe smallest enclosing sphere and for the set of spheres works only for two spheres

This commit is contained in:
ganovelli 2006-07-06 12:37:18 +00:00
parent 9374d88dbe
commit 69ddd70cf4
1 changed files with 100 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

100
vcg/space/smallest_enclosing.h Normal file
View File

@ -0,0 +1,100 @@
#pragma once
/****************************************************************************
* 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. *
* *
****************************************************************************/
/****************************************************************************
History
$Log: not supported by cvs2svn $
****************************************************************************/
#include <vcg/space/triangle3.h>
#include <vcg/space/tetra3.h>
#include <assert.h>
namespace vcg{
/** \addtogroup space */
/*@{*/
/**
Class for function computing the smallest enclosing bounding volume
*/
struct SmallestEnclosing {
/// computes the smallest enclosing sphere of a triangle
template <class TriangleType>
static Sphere3<typename TriangleType::ScalarType> SphereOfTriangle(const TriangleType & t);
/// computes the smallest enclosing sphere of a tetrahedron
template <class TetraType>
static Sphere3<typename TetraType::ScalarType> SphereOfTetra(const TetraType & t);
/// computes the smallest enclosing sphere of a container of spheres
template <class SphereContType>
static typename SphereContType::value_type SphereOfSpheres( const SphereContType & t);
};
/*@}*/
template <class TriangleType>
Sphere3<typename TriangleType::ScalarType>
static SmallestEnclosing::SphereOfTriangle(const TriangleType & t){
return Sphere3<typename TriangleType::ScalarType>(t.Barycenter(),(t.Barycenter()-t.cP(0)).Norm() );
}
template <class TetraType>
Sphere3<typename TetraType::ScalarType>
static SmallestEnclosing::SphereOfTetra(const TetraType & t){
return Sphere3<typename TetraType::ScalarType>( Barycenter(t),( Barycenter(t)-t.cP(0)).Norm() );
}
template <class SphereContType>
static typename SphereContType::value_type
SmallestEnclosing::
SphereOfSpheres( const SphereContType & spheres)
{
typename SphereContType::value_type::ScalarType radius;
vcg::Point3f center;
if(spheres.size()==2){
const typename SphereContType::value_type & s0 = spheres[0];
const typename SphereContType::value_type & s1 = spheres[1];
float dst = (s1.Center()-s0.Center()).Norm() ;
radius = (dst+s1.Radius()+s0.Radius())/2;
Point3f a=s0.Center();
Point3f b=s1.Center();
Point3f dir = (b-a).Normalize();
a = a - dir*s0.Radius();
b = b + dir*s1.Radius();
center = (a+b)/2.0;
}
else{
assert(0);
}
return typename SphereContType::value_type(center,radius);
}
}