X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=packages%2Fvtk%2Fsrc%2FbbvtkFlip.cxx;h=a252888c5a04d9bffdaddef6a76f19e3dab1f7b8;hb=29a64ec37c5dc46cc05277a4f7358a77026bfbda;hp=4f9244fe7a101ad0d78cd51f1f4bdc3e6c7db3f9;hpb=e4f0bedceaed7f7041b846375bf2bcafad4e47c0;p=bbtk.git diff --git a/packages/vtk/src/bbvtkFlip.cxx b/packages/vtk/src/bbvtkFlip.cxx index 4f9244f..a252888 100644 --- a/packages/vtk/src/bbvtkFlip.cxx +++ b/packages/vtk/src/bbvtkFlip.cxx @@ -1,33 +1,38 @@ -/*========================================================================= +/* + # --------------------------------------------------------------------- + # + # 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: bbvtkFlip.cxx,v $ Language: C++ - Date: $Date: 2012/04/11 16:27:00 $ - Version: $Revision: 1.2 $ + Date: $Date: 2012/11/16 08:51:58 $ + Version: $Revision: 1.5 $ =========================================================================*/ -/* --------------------------------------------------------------------- - -* 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 @@ -49,6 +54,7 @@ namespace bbvtk //std::cout << "-------- entree ds Flip::bbUserSetDefaultValues()\n" << std::endl; bbSetInputIn(NULL); + bbSetInputAxis("Y"); mImageOut = NULL; bbSetOutputOut(NULL); } @@ -58,8 +64,6 @@ namespace bbvtk void Flip::bbUserInitializeProcessing() { //std::cout << "-------- entree ds Flip::bbUserInitalizeProcessing()\n" << std::endl; - //bbUserFinalizeProcessing(); - mImageOut = vtkImageData::New(); // Alloc depends on bbGetInputIn().size() } //--------------------------------------------------------------------- @@ -70,8 +74,8 @@ namespace bbvtk // WTF? we never enter here // JPR bbUserFinalizeProcessing() JPR if (mImageOut!=NULL) { - // mImageOut->Delete(); - // mImageOut=NULL; + mImageOut->Delete(); + mImageOut=NULL; } bbSetOutputOut(mImageOut); } @@ -85,83 +89,99 @@ namespace bbvtk /// void Flip::Process() { - vtkImageData* imageIn = bbGetInputIn(); - vtkImageData * mImageOut = flip(imageIn); - - // Devrait etre dans bbUserFinalizeProcessing() ? // JPR - bbSetOutputOut(mImageOut); -} + + if (bbGetInputIn()!=NULL) + { + if (mImageOut!=NULL) + { + mImageOut->Delete(); + mImageOut=NULL; + } // if mImageOut + + mImageOut = vtkImageData::New(); // Alloc depends on bbGetInputIn().size() + mImageOut->Initialize(); + mImageOut->SetScalarType( bbGetInputIn()->GetScalarType() ); + mImageOut->SetSpacing( bbGetInputIn()->GetSpacing() ); + mImageOut->SetDimensions( bbGetInputIn()->GetDimensions() ); + mImageOut->SetNumberOfScalarComponents( bbGetInputIn()->GetNumberOfScalarComponents() ); + mImageOut->AllocateScalars(); + + int inputdims[3]; + // int outputdims[3]; + bbGetInputIn()->GetDimensions (inputdims); + bbGetInputIn()->Update(); + int nbScalComp = bbGetInputIn()->GetNumberOfScalarComponents(); + int scalarSize = bbGetInputIn()->GetScalarSize(); + int lineSize = inputdims[0]*scalarSize*nbScalComp; + int planeSize = inputdims[1]*lineSize; + int volumeSize = inputdims[2]*planeSize; + char *pixelsIn = (char *)bbGetInputIn()->GetScalarPointer(); + char *pixelsOut = (char *)mImageOut->GetScalarPointer(); + + char *lineIn; + char *lineOut; + char *debPlanIn; + char *debPlanOut; + int i,j,k; + + if ( (bbGetInputAxis()=="O") || (bbGetInputAxis()=="o") ) + { + pixelsIn = (char *)bbGetInputIn()->GetScalarPointer(0,0,0); + pixelsOut = (char *)mImageOut->GetScalarPointer(0,0,0); + memcpy(pixelsOut, pixelsIn, volumeSize ); + } else if ( (bbGetInputAxis()=="X") || (bbGetInputAxis()=="x") ) + { + + for(k=0; kGetScalarPointer(i,j,k); + pixelsOut = (char *)mImageOut->GetScalarPointer(inputdims[0]-1-i,j,k); + memcpy(pixelsOut, pixelsIn, scalarSize ); + } // for i + } // for j + } // for k + + } else if ( (bbGetInputAxis()=="Z") || (bbGetInputAxis()=="z") ) { + + for(k=0; kGetScalarPointer(i,j,k); + pixelsOut = (char *)mImageOut->GetScalarPointer(i,j,inputdims[2]-1-k); + memcpy(pixelsOut, pixelsIn, scalarSize ); + } // for i + } // for j + } // for k + + } else { + + for(k=0; kModified(); + bbSetOutputOut( mImageOut ); + } // bbGetInputIn - -vtkImageData * Flip::flip(vtkImageData *imageIn) -{ - int inputdims[3]; -// int outputdims[3]; - imageIn->GetDimensions (inputdims); - imageIn->Update(); - int nbScalComp = imageIn->GetNumberOfScalarComponents(); - int scalarSize = imageIn->GetScalarSize(); - int lineSize = inputdims[0]*scalarSize*nbScalComp; - int planeSize = inputdims[1]*lineSize; - char *pixelsIn = (char *)imageIn->GetScalarPointer(); - -/* - outputdims[0] = inputdims[1]; - outputdims[1] = inputdims[0]; - outputdims[2] = inputdims[2]; - std::cout << " inputdims[0] " << inputdims[0] <<" inputdims[1] " <SetDimensions( outputdims ); - vtkImageOut->SetExtent(0,outputdims[0]-1,0,outputdims[1]-1,0,outputdims[2]-1); - vtkImageOut->SetWholeExtent(0,outputdims[0]-1,0,outputdims[1]-1,0,outputdims[2]-1); - int nbScalComp = imageIn->GetNumberOfScalarComponents(); - vtkImageOut->SetNumberOfScalarComponents(nbScalComp); - vtkImageOut->SetSpacing( imageIn->GetSpacing() ); - vtkImageOut->SetScalarType(imageIn->GetScalarType() ); - vtkImageOut->AllocateScalars(); - vtkImageOut->Update(); - char *pixelsOut = (char *)vtkImageOut->GetScalarPointer(); - unsigned short *output =(unsigned short *)(vtkImageOut->GetScalarPointer()); - */ - - char *temp = (char *)malloc(lineSize); - char *line1; - char *line2; - char *debPlan; - for(int k=0; kSetScalarComponentFromDouble(j,i,k,l,imageIn->GetScalarComponentAsDouble(i,j,k,l)); - } - } - - } - } - return vtkImageOut; -*/ - free (temp); - return imageIn; } + }//namespace bbtk #endif // _USE_VTK_