1 /*# ---------------------------------------------------------------------
3 # Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image
5 # Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
6 # Previous Authors : Laurent Guigues, Jean-Pierre Roux
7 # CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil
9 # This software is governed by the CeCILL-B license under French law and
10 # abiding by the rules of distribution of free software. You can use,
11 # modify and/ or redistribute the software under the terms of the CeCILL-B
12 # license as circulated by CEA, CNRS and INRIA at the following URL
13 # http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
14 # or in the file LICENSE.txt.
16 # As a counterpart to the access to the source code and rights to copy,
17 # modify and redistribute granted by the license, users are provided only
18 # with a limited warranty and the software's author, the holder of the
19 # economic rights, and the successive licensors have only limited
22 # The fact that you are presently reading this means that you have had
23 # knowledge of the CeCILL-B license and that you accept its terms.
24 # ------------------------------------------------------------------------ */
26 #include "CutModelFigureSphere.h"
28 //----------------------------------------------------------------------------
29 //----------------------------------------------------------------------------
30 //----------------------------------------------------------------------------
31 CutModelFigureSphere::CutModelFigureSphere()
33 _vtksphere = vtkSphereSource::New();
34 _vtksphere->SetThetaResolution (20);
35 _vtksphere->SetPhiResolution (20);
37 //----------------------------------------------------------------------------
38 CutModelFigureSphere::~CutModelFigureSphere() // virtual
42 //----------------------------------------------------------------------------
43 bool CutModelFigureSphere::IfPointInside(double x, double y, double z) // virtual
50 _inversModel->MultiplyPoint (in, out);
53 if (sqrt( out[0]*out[0] + out[1]*out[1] + out[2]*out[2] )<0.5 )
59 //----------------------------------------------------------------------------
60 double CutModelFigureSphere::GetTheoricVolume() // virtual
62 double piTMP=3.14159265;
63 return (4.0/3.0) * piTMP * (_sx/2)*(_sy/2)*(_sz/2);
65 //----------------------------------------------------------------------------
66 const char *CutModelFigureSphere::GetName() // virtual
71 vtkPolyData* CutModelFigureSphere::getPolyData(){
72 return _vtksphere->GetOutput();