X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=packages%2Fvtk%2Fsrc%2FbbvtkMIPCreator.cxx;h=c822385d7eaba8a256b5c4148a0cadf43148d8c6;hb=d19c5512748aa8e2dd5c3d01d9ec680ded1f17ce;hp=58ee8ae42ec02c5da30f1e1239a353239079f01f;hpb=a35823bf58c3560e9fe218ecd0f1fab421b6ebed;p=bbtk.git diff --git a/packages/vtk/src/bbvtkMIPCreator.cxx b/packages/vtk/src/bbvtkMIPCreator.cxx index 58ee8ae..c822385 100644 --- a/packages/vtk/src/bbvtkMIPCreator.cxx +++ b/packages/vtk/src/bbvtkMIPCreator.cxx @@ -1,40 +1,48 @@ -/*========================================================================= +/* + # --------------------------------------------------------------------- + # + # 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: bbvtkMIPCreator.cxx,v $ Language: C++ - Date: $Date: 2010/04/08 14:37:59 $ - Version: $Revision: 1.7 $ + Date: $Date: 2012/11/16 08:51:58 $ + Version: $Revision: 1.12 $ =========================================================================*/ -/* --------------------------------------------------------------------- - -* 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 + * \file + * \brief */ #ifdef _USE_VTK_ +#include +#include + #include "bbvtkMIPCreator.h" #include "bbvtkPackage.h" @@ -42,44 +50,48 @@ namespace bbvtk { BBTK_ADD_BLACK_BOX_TO_PACKAGE(vtk,MIPCreator) BBTK_BLACK_BOX_IMPLEMENTATION(MIPCreator,bbtk::AtomicBlackBox); - + //--------------------------------------------------------------------- - void MIPCreator::bbUserSetDefaultValues() - { - bbSetInputShift(0); - bbSetInputScale(1.0); - mCast = NULL; - mMIP = NULL; - mMapper = NULL; - mVolume = NULL; + void MIPCreator::bbUserSetDefaultValues() + { + + + + std::vector opaValues; + std::vector denPoints; + opaValues.push_back(0.0); + opaValues.push_back(1.0); + bbSetInputDensityPoints(denPoints); + bbSetInputOpacityValues(opaValues); + bbSetInputMinColorWindow(0.0); + bbSetInputMaxColorWindow(255.0); + + opacityTransferFunction = vtkPiecewiseFunction::New(); + grayTransferFunction = vtkPiecewiseFunction::New(); + + mCast = NULL; + mMIP = NULL; + mMapper = NULL; + mVolume = NULL; + + } //--------------------------------------------------------------------- - void MIPCreator::bbUserInitializeProcessing() - { + void MIPCreator::bbUserInitializeProcessing() + { + + // Create the pipeline mCast = vtkImageShiftScale::New(); mCast->SetOutputScalarTypeToUnsignedChar(); mCast->ClampOverflowOn(); - - mMIP = vtkVolumeRayCastMIPFunction::New(); - mMIP->SetMaximizeMethodToScalarValue(); - - mMapper = vtkVolumeRayCastMapper::New(); - mMapper->SetVolumeRayCastFunction(mMIP); - mMapper->SetInput(mCast->GetOutput()); // (smoother.GetOutput()) - - mVolume = vtkVolume::New(); - mVolume->SetMapper(mMapper); - - // mMapper->ScalarVisibilityOff(); - // mMapper->ImmediateModeRenderingOn(); - - bbSetOutputOut(mVolume); + + } //--------------------------------------------------------------------- - void MIPCreator::bbUserFinalizeProcessing() - { + void MIPCreator::bbUserFinalizeProcessing() + { if (mCast!=NULL) { mCast->Delete(); @@ -108,15 +120,131 @@ namespace bbvtk void MIPCreator::Process() { - mCast->SetInput( bbGetInputIn() ); - mCast->SetScale( bbGetInputScale() / 100.0 ); - mCast->SetShift( -bbGetInputShift() ); - mMapper->Update(); - + + + double ranges[2]; + vtkImageData *image = bbGetInputIn(); +//EED 2017-01-01 Migration VTK7 +#if (VTK_MAJOR_VERSION <= 5) + image->Update(); +#endif +#if (VTK_MAJOR_VERSION >= 6) + // .. +#endif + image->GetScalarRange(ranges); + + scale = 255/(ranges[1]/100); + +//EED 2017-01-01 Migration VTK7 +#if (VTK_MAJOR_VERSION <= 5) + mCast->SetInput( bbGetInputIn() ); +#endif +#if (VTK_MAJOR_VERSION >= 6) + mCast->SetInputData( bbGetInputIn() ); +#endif + mCast->SetScale( scale / 100.0 ); + //mCast->SetShift( -bbGetInputShift() ); + + + + + + + + + + DensityPoints = bbGetInputDensityPoints(); + + //We assign by default the values of density by regarding at range from the image + if(DensityPoints.size()==0) + { + DensityPoints.push_back(ranges[0]); + DensityPoints.push_back(ranges[1]); + } + + + OpacityValues = bbGetInputOpacityValues(); + + while(DensityPoints.size()>OpacityValues.size()) + { + OpacityValues.push_back(1.0); + } + + + + double factor = scale/ 100.0 ; + + double minCF = bbGetInputMinColorWindow(); + double maxCF = bbGetInputMaxColorWindow(); + + grayTransferFunction->AddSegment( minCF , 0.0 , maxCF , 1.0 ); + + + + + for(int i =0; i< DensityPoints.size(); i++) + { + cout<<"MIPCreator::Process::dPoint"<AddPoint(DensityPoints.at(i)*factor, OpacityValues.at(i)); + } + //opacityTransferFunction = vtkPiecewiseFunction::New(); +// opacityTransferFunction->AddPoint( 0 , 0.0/100.0 ); +// opacityTransferFunction->AddPoint( 85 , 2.0/100.0 ); +// opacityTransferFunction->AddPoint( 160 , 25.0/100.0 ); +// opacityTransferFunction->AddPoint( 210 , 50.0/100.0 ); +// opacityTransferFunction->AddPoint( 255 , 90.0/100.0 ); + + + + // Create a set of properties for mip + vtkVolumeProperty *mipProperty; + mipProperty = vtkVolumeProperty::New(); + mipProperty->SetScalarOpacity( opacityTransferFunction ); + mipProperty->SetColor( grayTransferFunction ); + +// mipProperty->SetInterpolationTypeToLinear(); + mipProperty->SetInterpolationTypeToNearest(); + +// mipProperty->ShadeOff(); + mipProperty->ShadeOn(); + +// mipProperty->SetAmbient(0.9); +// mipProperty->SetDiffuse(0.9); +// mipProperty->SetSpecular(0.9); + + + + + + mMIP = vtkVolumeRayCastMIPFunction::New(); + mMIP->SetMaximizeMethodToScalarValue(); + + mMapper = vtkVolumeRayCastMapper::New(); + mMapper->SetVolumeRayCastFunction(mMIP); + +//EED 2017-01-01 Migration VTK7 +#if (VTK_MAJOR_VERSION <= 5) + mMapper->SetInput(mCast->GetOutput()); // (smoother.GetOutput()) +#endif +#if (VTK_MAJOR_VERSION >= 6) + mMapper->SetInputData(mCast->GetOutput()); // (smoother.GetOutput()) +#endif + + mVolume = vtkVolume::New(); + mVolume->SetMapper(mMapper); + mVolume->SetProperty( mipProperty ); + + + + mMapper->Update(); + + bbSetOutputOut(mVolume); // mVolume->GetProperty()->SetColor( bbGetInputColour()[0], bbGetInputColour()[1], bbGetInputColour()[2] ); // mVolume->GetProperty()->SetOpacity( bbGetInputOpacity() ); - + } }//namespace bbtk