X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=packages%2Fvtk%2Fsrc%2FbbvtkMIPCreator.cxx;h=b65e51d22893d90ba8a461a1a5044b6052be254f;hb=36a131a7d47f451563f40b28047200b038905c92;hp=d45db49bf1ecb040145116bfa8c476cf1c759cdd;hpb=00045104d4fbcd4f71774a9634739bf954259e4e;p=bbtk.git diff --git a/packages/vtk/src/bbvtkMIPCreator.cxx b/packages/vtk/src/bbvtkMIPCreator.cxx index d45db49..b65e51d 100644 --- a/packages/vtk/src/bbvtkMIPCreator.cxx +++ b/packages/vtk/src/bbvtkMIPCreator.cxx @@ -1,36 +1,41 @@ -/*========================================================================= +/* + # --------------------------------------------------------------------- + # + # 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: 2011/03/17 15:49:14 $ - Version: $Revision: 1.8 $ + 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_ @@ -45,79 +50,38 @@ 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(); - - - - - //define range of values - vtkPiecewiseFunction *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 transfer function mapping scalar value to color (grey) - vtkPiecewiseFunction *grayTransferFunction = vtkPiecewiseFunction::New(); - grayTransferFunction->AddSegment( 0 , 0.0 , 255 , 1.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(); -// mMIP->SetMaximizeMethodToOpacity(); - - mMapper = vtkVolumeRayCastMapper::New(); - mMapper->SetVolumeRayCastFunction(mMIP); - mMapper->SetInput(mCast->GetOutput()); // (smoother.GetOutput()) - - mVolume = vtkVolume::New(); - mVolume->SetMapper(mMapper); - mVolume->SetProperty( mipProperty ); - - // mMapper->ScalarVisibilityOff(); - // mMapper->ImmediateModeRenderingOn(); - - bbSetOutputOut(mVolume); } //--------------------------------------------------------------------- - void MIPCreator::bbUserFinalizeProcessing() - { + void MIPCreator::bbUserFinalizeProcessing() + { if (mCast!=NULL) { mCast->Delete(); @@ -146,15 +110,89 @@ namespace bbvtk void MIPCreator::Process() { - mCast->SetInput( bbGetInputIn() ); - mCast->SetScale( bbGetInputScale() / 100.0 ); - mCast->SetShift( -bbGetInputShift() ); + 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 ); - mMapper->Update(); - + // 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); + mCast->Update(); + +//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