+/*# ---------------------------------------------------------------------
+#
+# Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image
+# pour la Sant�)
+# Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
+# Previous Authors : Laurent Guigues, Jean-Pierre Roux
+# CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil
+#
+# This software is governed by the CeCILL-B license under French law and
+# abiding by the rules of distribution of free software. You can use,
+# modify and/ or redistribute the software under the terms of the CeCILL-B
+# license as circulated by CEA, CNRS and INRIA at the following URL
+# http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
+# or in the file LICENSE.txt.
+#
+# As a counterpart to the access to the source code and rights to copy,
+# modify and redistribute granted by the license, users are provided only
+# with a limited warranty and the software's author, the holder of the
+# economic rights, and the successive licensors have only limited
+# liability.
+#
+# The fact that you are presently reading this means that you have had
+# knowledge of the CeCILL-B license and that you accept its terms.
+# ------------------------------------------------------------------------ */
+
/*=========================================================================
Program: wxMaracas
Module: $RCSfile: CutModelData.cxx,v $
Language: C++
- Date: $Date: 2010/03/31 12:41:58 $
- Version: $Revision: 1.6 $
+ Date: $Date: 2012/11/15 14:15:48 $
+ Version: $Revision: 1.10 $
Copyright: (c) 2002, 2003
License:
=========================================================================*/
#include "CutModelData.h"
-
+#include "vtkRenderWindow.h"
/**
** Start of the manager class
**/
createActor();
createShapes();
ChangeShape(0);
+ interactor->GetRenderWindow ()->Render();
checkInvariant();
}
inversModel = vtkTransform::New();
vtkMatrix4x4* matrix = vtkMatrix4x4::New();
- matrix->Identity();
+ matrix->Identity();
+ double *ori = img->GetOrigin();
int *ext = img->GetExtent();
double *spc = img->GetSpacing();
matrix->SetElement(1,1,(ext[3]-ext[2])/4*spc[1]);
matrix->SetElement(2,2,(ext[5]-ext[4])/4*spc[2]);
- double orgx = (ext[1]-ext[0])/2*spc[0];
- double orgy = (ext[3]-ext[2])/2*spc[1];
- double orgz = (ext[5]-ext[4])/2*spc[2];
+ double orgx = (ext[1]+ext[0])/2.*spc[0];
+ double orgy = (ext[3]+ext[2])/2.*spc[1];
+ double orgz = (ext[5]+ext[4])/2.*spc[2];
matrix->SetElement(0,3,orgx);
matrix->SetElement(1,3,orgy);
copyimage->GetExtent(ext);
- for (xx=ext[0];xx<ext[1]; xx++)
+ for (xx=ext[0];xx<=ext[1]; xx++)
{
- for (yy=ext[2];yy<ext[3]; yy++)
+ for (yy=ext[2];yy<=ext[3]; yy++)
{
- for (zz=ext[4];zz<ext[5];zz++)
+ for (zz=ext[4];zz<=ext[5];zz++)
{
inside=actualCuttingModel->IfPointInside(xx,yy,zz);
if ( ((inside==true)&&(isinside==true)) || ((!inside==true)&&(!isinside)) )
copyimage->GetExtent(ext);
- for (xx=ext[0];xx<ext[1]; xx++)
+ for (xx=ext[0];xx<=ext[1]; xx++)
{
- for (yy=ext[2];yy<ext[3]; yy++)
+ for (yy=ext[2];yy<=ext[3]; yy++)
{
- for (zz=ext[4];zz<ext[5];zz++)
+ for (zz=ext[4];zz<=ext[5];zz++)
{
inside=actualCuttingModel->IfPointInside(xx,yy,zz);
if ( ((inside==true)&&(isinside==true)) || ((!inside==true)&&(!isinside)) )
void CutModelData::setTransform(vtkTransform* transform, vtkImageData* img)throw( CutModelException){
- double *spc = img->GetSpacing(); // spc unused // JPR
+// double *spc = img->GetSpacing(); // spc unused // JPR
currentmatrix->Identity();