X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;ds=inline;f=packages%2Fvtk%2Fsrc%2FbbvtkSegmentationConnectivity.cxx;h=6ec480d58c641086a2ad5a9f1c0ecfe2cc9f8c8a;hb=2b1ebe5f1c885888b35d860a8c083e5abdc53678;hp=5c33b54c367cfa1aaee0174ac790dd2bc3835348;hpb=17ac9460f5ea5600d1065fbe3ee4ad826d7bb4fc;p=bbtk.git diff --git a/packages/vtk/src/bbvtkSegmentationConnectivity.cxx b/packages/vtk/src/bbvtkSegmentationConnectivity.cxx index 5c33b54..6ec480d 100644 --- a/packages/vtk/src/bbvtkSegmentationConnectivity.cxx +++ b/packages/vtk/src/bbvtkSegmentationConnectivity.cxx @@ -1,54 +1,56 @@ -/*========================================================================= +/* + # --------------------------------------------------------------------- + # + # 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: bbtk Module: $RCSfile: bbvtkSegmentationConnectivity.cxx,v $ Language: C++ - Date: $Date: 2009/05/15 14:58:01 $ - Version: $Revision: 1.5 $ + Date: $Date: 2012/11/16 08:51:58 $ + Version: $Revision: 1.8 $ =========================================================================*/ -/* --------------------------------------------------------------------- - -* Copyright (c) CREATIS-LRMN (Centre de Recherche en Imagerie Medicale) -* Authors : Eduardo Davila, Laurent Guigues, Jean-Pierre Roux -* -* 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. -* ------------------------------------------------------------------------ */ - /** * \file * \brief */ - #ifdef _USE_VTK_ - #include "bbvtkSegmentationConnectivity.h" #include "bbvtkPackage.h" namespace bbvtk { - BBTK_ADD_BLACK_BOX_TO_PACKAGE(vtk,SegmentationConnectivity); BBTK_BLACK_BOX_IMPLEMENTATION(SegmentationConnectivity, bbtk::AtomicBlackBox); - - - // -------------------------------------------------------------- + +// -------------------------------------------------------------- + void SegmentationConnectivity::bbUserSetDefaultValues() { thresh2 = NULL; @@ -69,7 +71,8 @@ namespace bbvtk bbSetOutputOut(NULL); } - // -------------------------------------------------------------- +// -------------------------------------------------------------- + void SegmentationConnectivity::bbUserInitializeProcessing() { thresh2 = vtkImageThreshold::New(); @@ -88,8 +91,9 @@ namespace bbvtk cast4->SetInput(connect2->GetOutput()); cast4->SetOutputScalarTypeToUnsignedShort(); } - - // -------------------------------------------------------------- + +// -------------------------------------------------------------- + void SegmentationConnectivity::bbUserFinalizeProcessing() { if (thresh2!=NULL) @@ -111,11 +115,11 @@ namespace bbvtk { cast4->Delete(); cast4=NULL; - } - + } } - // -------------------------------------------------------------- +// -------------------------------------------------------------- + void SegmentationConnectivity::DoProcess() { vtkImageData *imagedata = bbGetInputIn(); @@ -128,17 +132,32 @@ namespace bbvtk thresh2->Update(); cast2->Update(); connect2->RemoveAllSeeds (); - connect2->AddSeed( bbGetInputPositionXYZ()[0] , bbGetInputPositionXYZ()[1] , bbGetInputPositionXYZ()[2] ); - connect2->Update(); - cast4->Update(); - - bbSetOutputOut(cast4->GetOutput() ); + + int x = bbGetInputPositionXYZ()[0]; + int y = bbGetInputPositionXYZ()[1]; + int z = bbGetInputPositionXYZ()[2]; + + int ext[6]; + imagedata->GetWholeExtent(ext); + int maxx = ext[1]-ext[0]+1; + int maxy = ext[3]-ext[2]+1; + int maxz = ext[5]-ext[4]+1; + if ( x<0 ) { printf("Warnnig bbtk::vtk::SegmentationConnectivity::DoProcess limites of the image out of range x: %d -> 0 \n",x ); x=0; } + if ( y<0 ) { printf("Warnnig bbtk::vtk::SegmentationConnectivity::DoProcess limites of the image out of range y: %d -> 0 \n",y ); y=0;} + if ( z<0 ) { printf("Warnnig bbtk::vtk::SegmentationConnectivity::DoProcess limites of the image out of range z: %d -> 0 \n",z ); z=0;} + + if ( x>=maxx ) { printf("Warnnig bbtk::vtk::SegmentationConnectivity::DoProcess limites of the image out of range x: %d -> %d \n",x, maxx-1 ); x=maxx-1; } + if ( y>=maxy ) { printf("Warnnig bbtk::vtk::SegmentationConnectivity::DoProcess limites of the image out of range y: %d -> %d \n",y, maxy-1 ); y=maxy-1; } + if ( z>=maxz ) { printf("Warnnig bbtk::vtk::SegmentationConnectivity::DoProcess limites of the image out of range z: %d -> %d \n",y, maxz-1 ); z=maxz-1; } + + connect2->AddSeed( x,y,z ); + connect2->Update(); + cast4->Update(); + bbSetOutputOut(cast4->GetOutput() ); } - }// EO namespace bbvtk - #endif //_USE_VTK_