X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;ds=sidebyside;f=packages%2Fvtk%2Fsrc%2FbbvtkSegmentationConnectivity.cxx;h=2deda14b1e2856a14a06bb4298fcd581e09d02b5;hb=4e201b91961b385dfe6fa4625a37f43c37b7d5d5;hp=e7172541089a7452181d7417db5f936cecddda02;hpb=41aba70a65336c651cae6e5211c069d8f33efb51;p=bbtk.git diff --git a/packages/vtk/src/bbvtkSegmentationConnectivity.cxx b/packages/vtk/src/bbvtkSegmentationConnectivity.cxx index e717254..2deda14 100644 --- a/packages/vtk/src/bbvtkSegmentationConnectivity.cxx +++ b/packages/vtk/src/bbvtkSegmentationConnectivity.cxx @@ -1,42 +1,59 @@ -/*========================================================================= - +/*========================================================================= Program: bbtk Module: $RCSfile: bbvtkSegmentationConnectivity.cxx,v $ Language: C++ - Date: $Date: 2008/04/24 14:11:31 $ - Version: $Revision: 1.2 $ - - Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de - l'Image). All rights reserved. See Doc/License.txt or - http://www.creatis.insa-lyon.fr/Public/bbtk/License.html for details. - - This software is distributed WITHOUT ANY WARRANTY; without even - the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR - PURPOSE. See the above copyright notices for more information. - + Date: $Date: 2011/06/16 12:28:43 $ + Version: $Revision: 1.7 $ =========================================================================*/ + +/* --------------------------------------------------------------------- + +* 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::bbUserConstructor() - { - bbSetInputIn(NULL); + BBTK_ADD_BLACK_BOX_TO_PACKAGE(vtk,SegmentationConnectivity); + BBTK_BLACK_BOX_IMPLEMENTATION(SegmentationConnectivity, + bbtk::AtomicBlackBox); + +// -------------------------------------------------------------- + + void SegmentationConnectivity::bbUserSetDefaultValues() + { + thresh2 = NULL; + cast2 = NULL; + connect2 = NULL; + cast4 = NULL; + + bbSetInputIn(NULL); std::vector position; position.push_back(0); position.push_back(0); @@ -47,47 +64,95 @@ namespace bbvtk threshold.push_back(0); bbSetInputThresholdMinMax(threshold); bbSetOutputOut(NULL); - thresh2 = vtkImageThreshold::New(); - thresh2->SetInValue(255); - thresh2->SetOutputScalarTypeToUnsignedShort(); - thresh2->SetOutValue(0); - cast2 = vtkImageCast::New(); - cast2->SetInput(thresh2->GetOutput()); - cast2->SetOutputScalarTypeToUnsignedChar(); - connect2 = vtkImageSeedConnectivity::New(); - connect2->SetInput(cast2->GetOutput()); - connect2->SetInputConnectValue(255); - connect2->SetOutputConnectedValue(255); - connect2->SetOutputUnconnectedValue(0); - cast4 = vtkImageCast::New(); - cast4->SetInput(connect2->GetOutput()); - cast4->SetOutputScalarTypeToUnsignedShort(); - } + } +// -------------------------------------------------------------- + + void SegmentationConnectivity::bbUserInitializeProcessing() + { + thresh2 = vtkImageThreshold::New(); + thresh2->SetInValue(255); + thresh2->SetOutputScalarTypeToUnsignedShort(); + thresh2->SetOutValue(0); + cast2 = vtkImageCast::New(); + cast2->SetInput(thresh2->GetOutput()); + cast2->SetOutputScalarTypeToUnsignedChar(); + connect2 = vtkImageSeedConnectivity::New(); + connect2->SetInput(cast2->GetOutput()); + connect2->SetInputConnectValue(255); + connect2->SetOutputConnectedValue(255); + connect2->SetOutputUnconnectedValue(0); + cast4 = vtkImageCast::New(); + cast4->SetInput(connect2->GetOutput()); + cast4->SetOutputScalarTypeToUnsignedShort(); + } + +// -------------------------------------------------------------- + + void SegmentationConnectivity::bbUserFinalizeProcessing() + { + if (thresh2!=NULL) + { + thresh2->Delete(); + thresh2=NULL; + } + if (cast2!=NULL) + { + cast2->Delete(); + cast2=NULL; + } + if (connect2!=NULL) + { + connect2->Delete(); + connect2=NULL; + } + if (cast4!=NULL) + { + cast4->Delete(); + cast4=NULL; + } + } +// -------------------------------------------------------------- + void SegmentationConnectivity::DoProcess() { - vtkImageData *imagedata = bbGetInputIn(); - imagedata->UpdateInformation(); - imagedata->SetUpdateExtent(imagedata->GetWholeExtent()); - imagedata->Update(); - thresh2->ThresholdBetween(3000, 3001); - thresh2->SetInput(imagedata); - thresh2->ThresholdBetween(bbGetInputThresholdMinMax()[0], bbGetInputThresholdMinMax()[1]); - thresh2->Update(); - cast2->Update(); - connect2->RemoveAllSeeds (); - connect2->AddSeed( bbGetInputPositionXYZ()[0] , bbGetInputPositionXYZ()[1] , bbGetInputPositionXYZ()[2] ); - connect2->Update(); - cast4->Update(); - - bbSetOutputOut(cast4->GetOutput() ); + vtkImageData *imagedata = bbGetInputIn(); + imagedata->UpdateInformation(); + imagedata->SetUpdateExtent(imagedata->GetWholeExtent()); + imagedata->Update(); + thresh2->ThresholdBetween(3000, 3001); + thresh2->SetInput(imagedata); + thresh2->ThresholdBetween(bbGetInputThresholdMinMax()[0], bbGetInputThresholdMinMax()[1]); + thresh2->Update(); + cast2->Update(); + connect2->RemoveAllSeeds (); + + 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_