+/*# ---------------------------------------------------------------------
+#
+# 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.
+# ------------------------------------------------------------------------ */
+
#include "vtkClipping3DDataViewer.h"
#include "vtkDataSetReader.h"
vtkClipping3DDataViewer::vtkClipping3DDataViewer()
{
int i;
- for (i=0; i<VTKMPRDATA_MAXTISSUE; i++){
+ for (i=0; i<VTKMPRDATA_MAXTISSUE; i++)
+ {
_visibleTissue[i] = false;
_representationType[i] = true;
}
- _visibleVolume=false;
+ _visibleVolume = false;
+ _isRayCasting = true;
+ _isMIP = false;
+ _interpolation = false;
+ _shade = false;
+
+ _volumerendererdata = VolumeRendererData::New();
+ _volumeMapper = vtkVolumeRayCastMapper::New();
+ _newvol = vtkVolume::New();
+ _volumeProperty = vtkVolumeProperty::New();
+ _compositeFunction = vtkVolumeRayCastCompositeFunction::New();
+ _compositeFunctionMIP = vtkVolumeRayCastMIPFunction::New();
+
}
//-------------------------------------------------------------------
vtkClipping3DDataViewer::~vtkClipping3DDataViewer()
{
int i;
for (i=0; i< VTKMPRDATA_MAXTISSUE ; i++){
- _mCubes[i] ->Delete();
- _tissueStripper[i] ->Delete();
- _tissueMapper[i] ->Delete();
- _tissuePlanes[i] ->Delete();
- _tissueClipper[i] ->Delete();
+ _mCubes[i]->Delete();
+ _tissueStripper[i]->Delete();
+ _tissueMapper[i]->Delete();
+ _tissuePlanes[i]->Delete();
+ _tissueClipper[i]->Delete();
}
- _outlineData ->Delete();
- _mapOutline ->Delete();
- _outline ->Delete();
+ _outlineData->Delete();
+ _mapOutline->Delete();
+ _outline->Delete();
// Volume
- _tfun ->Delete();
- _ctfun ->Delete();
- _compositeFunction ->Delete();
- _volumeMapper ->Delete();
- _volumeProperty ->Delete();
- _newvol ->Delete();
- _volumePlanes ->Delete();
+ _tfun->Delete();
+ _ctfun->Delete();
+ _compositeFunction->Delete();
+ _volumeMapper->Delete();
+ _volumeProperty->Delete();
+ _newvol->Delete();
+ _volumePlanes->Delete();
for (i=0; i< VTKMPRDATA_MAXTISSUE ; i++){
- _tissue[i] -> Delete();
- _observerS[i] -> Delete();
+ _tissue[i]-> Delete();
+ _observerS[i]-> Delete();
}
- _observerV -> Delete();
+ _observerV-> Delete();
}
//-------------------------------------------------------------------
vtkActor* vtkClipping3DDataViewer::GetOutlineActor()
// _newvol->VisibilityOn();
-
+
// _volumeMapper = vtkVolumeRayCastMapper::New();
// _volumeMapper->SetInput( this->GetVtkMPRBaseData()->GetImageData() );
// _volumeMapper->SetVolumeRayCastFunction(_compositeFunction);
_visibleVolume = visiblevolume;
}
+bool vtkClipping3DDataViewer::GetVisibleVolumeBox()
+{
+ return _visibleVolumeBox;
+}
+//-------------------------------------------------------------------
+void vtkClipping3DDataViewer::SetVisibleVolumeBox(bool visibleBox)
+{
+ _visibleVolumeBox = visibleBox;
+}
+
//-------------------------------------------------------------------
-void vtkClipping3DDataViewer::Configure_Tissue()
+void vtkClipping3DDataViewer::SetImage()
{
+ int i;
+ vtkImageData *imagedata=_vtkmprbasedata->GetImageData();
+ for (i=0; i< VTKMPRDATA_MAXTISSUE ; i++)
+ {
+ _mCubes[i]->SetInput( imagedata );
+ } // for
+ _volumeMapper->SetInput( imagedata );
+ _volumerendererdata->SetImageData( imagedata );
+}
- double range[2];
+//-------------------------------------------------------------------
+void vtkClipping3DDataViewer::Configure_Tissue()
+{
int i;
- for (i=0; i< VTKMPRDATA_MAXTISSUE ; i++){
+ double range[2];
+ vtkImageData *imagedata=_vtkmprbasedata->GetImageData();
+ imagedata->GetScalarRange( range );
+ for (i=0; i< VTKMPRDATA_MAXTISSUE ; i++)
+ {
// Visualisation - original volume
- vtkImageData *imagedata=_vtkmprbasedata->GetImageData();
- imagedata->GetScalarRange( range );
_mCubes[i] = vtkMarchingCubes::New( );
_mCubes[i]->SetInput( imagedata );
-
_mCubes[i]->SetValue( 0, range[1]*(4+i) / 8 );
// _mCubes[i]->SetValue( 0, 1500 );
-
_tissueStripper[i] = vtkStripper::New();
_tissueStripper[i]->SetInput( _mCubes[i]->GetOutput( ) );
-
-
-
_tissuePlanes[i] = vtkPlanes::New();
-
int x1,x2,y1,y2,z1,z2;
imagedata->GetExtent(x1,x2,y1,y2,z1,z2);
_tissuePlanes[i]->SetBounds (x1,x2,y1,y2,z1,z2);
-
_tissueClipper[i] = vtkClipPolyData::New();
_tissueClipper[i]->SetInput( _tissueStripper[i]->GetOutput() );
_tissueClipper[i]->SetClipFunction( _tissuePlanes[i] );
_tissueClipper[i]->InsideOutOn( );
-
_tissueMapper[i] = vtkPolyDataMapper::New( );
_tissueMapper[i]->SetInput( _tissueClipper[i]->GetOutput() );
_tissueMapper[i]->ScalarVisibilityOff( );
// _tissueMapper[i]->Update();
-
- }
-
-
+ } // for i
// vtkActor tissue
for (i=0; i< VTKMPRDATA_MAXTISSUE ; i++){
{
// Volume
+ /** if(_newvol != NULL)
+ {
+ _newvol->Delete();
+ _volumeMapper->Delete();
+ _volumeProperty->Delete();
+ }*/
+
_tfun = vtkPiecewiseFunction::New();
_ctfun = vtkColorTransferFunction::New();
double range[2];
this->_vtkmprbasedata->GetImageData()->GetScalarRange(range);
double max = range[1];
-
+
/*
adding the poinst of the transference function
*/
//X
- greyValuesTransferenceFVector.push_back(max * 0/2);
- greyValuesTransferenceFVector.push_back(max * 1/2);
+ greyValuesTransferenceFVector.push_back(max * 0/2); // JPR ?? 0/2
+ greyValuesTransferenceFVector.push_back(max * 1/2); // JPR ?? division entiere?
greyValuesTransferenceFVector.push_back(max * 2/2);
//Y
intensityValuesTransferenceFVector.push_back(0.0);
_tfun->AddPoint(max * 2/2 , 1.0);
/*
- Adding the colors to the vectors
+ Adding the colors to the vectors
*/
//RED
redColorsOfColorTransferenceFVector.push_back(0.0);
greenColorsOfColorTransferenceFVector.push_back(0.0);
greenColorsOfColorTransferenceFVector.push_back(1.0);
greenColorsOfColorTransferenceFVector.push_back(0.2);
-
+
//BLUE
blueColorsOfColorTransferenceFVector.push_back(0.0);
blueColorsOfColorTransferenceFVector.push_back(0.0);
blueColorsOfColorTransferenceFVector.push_back(1.0);
blueColorsOfColorTransferenceFVector.push_back(0.0);
blueColorsOfColorTransferenceFVector.push_back(0.0);
- //GREY VALUE
+ //GREY VALUE
greyValueColorsOfColorTransferenceFVector.push_back(max*0/4);
greyValueColorsOfColorTransferenceFVector.push_back(max*1/4);
greyValueColorsOfColorTransferenceFVector.push_back(max*2/4);
// imagedata->GetExtent(x1,x2,y1,y2,z1,z2);
// _volumePlanes->SetBounds(x1,x2,y1,y2,z1,z2);
- _compositeFunction = vtkVolumeRayCastCompositeFunction::New();
+ // EED 13/03/2011
- _volumeMapper = vtkVolumeRayCastMapper::New();
+ // _volumeMapper = vtkVolumeRayCastMapper::New();
_volumeMapper->SetInput( this->GetVtkMPRBaseData()->GetImageData() );
- _volumeMapper->SetVolumeRayCastFunction(_compositeFunction);
+
+ _volumeMapper->SetVolumeRayCastFunction(_compositeFunction);
+
+ // _volumeMapper->SetVolumeRayCastFunction(_compositeFunctionMIP);
+
+//
_volumeMapper->SetClippingPlanes( _volumePlanes );
_volumeMapper->AutoAdjustSampleDistancesOn();
- _volumeProperty = vtkVolumeProperty::New();
+ _volumeProperty->SetInterpolationTypeToNearest();
+ //_volumeProperty = vtkVolumeProperty::New();
_volumeProperty->SetColor(_ctfun);
_volumeProperty->SetScalarOpacity( _tfun );
+
// EED 31/03/2008
- _volumeProperty->SetInterpolationTypeToLinear();
- _volumeProperty->ShadeOn();
_volumeProperty->DisableGradientOpacityOn();
+
// _volumeProperty->SetInterpolationTypeToNearest();
// _volumeProperty->ShadeOff();
// _volumeProperty->SetAmbient(0.3);
// _volumeProperty->SetSpecular(0.8);
// _volumeProperty->DisableGradientOpacityOn();
- _newvol = vtkVolume::New();
+ // _newvol = vtkVolume::New();
_newvol->SetMapper(_volumeMapper );
_newvol->SetProperty(_volumeProperty );
- _observerV = boxSurfaceObserver::New();
+ _observerV = boxSurfaceObserver::New();
_observerV->SetPlanes( _volumePlanes );
_observerV->SetActor( _newvol );
_observerV->SetvtkVolumeRayCastMapper( _volumeMapper );
}
+//-----------------------------------------------------------------------------
+
+void vtkClipping3DDataViewer::Configure_VolumeBox()
+{
+ _volumerendererdata->SetImageData(_vtkmprbasedata->GetImageData());
+
+
+ vector<double> *gv = this->GetGreyValuesTransferenceFVector();
+ vector<double> *op = this->GetIntensityValuesTransferenceFVector();
+
+ _volumerendererdata->SetLookUpTable((vtkLookupTable*)(_ctfun));
+ _volumerendererdata->OpacityChanged(*gv, *op);
+
+ //_wxvtkclipping3Dview->VisibleVolumeActor( false );
+ //_volumerendererdata->ShowVolume(event.IsChecked() );
+ //_volumerendererdata->BoundingBoxChanged(event.IsChecked() );
+}
+
+//-----------------------------------------------------------------------------
+
+void vtkClipping3DDataViewer::SetRenderer(vtkRenderer* renderer)
+{
+ _volumerendererdata->SetRenderer(renderer);
+
+}
+
+//-----------------------------------------------------------------------------
+void vtkClipping3DDataViewer::SetInteractor(vtkRenderWindowInteractor* interactor)
+{
+
+ _volumerendererdata->SetInteractor(interactor);
+}
+
+
+//-----------------------------------------------------------------------------
+void vtkClipping3DDataViewer::BoxActorChanged(bool changed)
+{
+ vector<double> *gv = this->GetGreyValuesTransferenceFVector();
+ vector<double> *op = this->GetIntensityValuesTransferenceFVector();
+
+ _volumerendererdata->ShowVolume(changed);
+ _volumerendererdata->OpacityChanged(*gv, *op);
+ _volumerendererdata->BoundingBoxChanged(changed);
+
+}
+
+//-----------------------------------------------------------------------------
+void vtkClipping3DDataViewer::UpdateVolumeBox(vector<double> gf, vector<double> vf, vtkColorTransferFunction* ctfun)
+{
+ _volumerendererdata->SetLookUpTable((vtkLookupTable*)(ctfun));
+ _volumerendererdata->OpacityChanged(gf, vf);
+}
+
+
+//-----------------------------------------------------------------------------
+void vtkClipping3DDataViewer::updateVolume()
+{
+
+ if(_isRayCasting)
+ {
+ _volumeMapper->SetVolumeRayCastFunction(_compositeFunction);
+ _volumerendererdata->changeCompositeMIPFunction(0);
+ }
+ else if(_isMIP)
+ {
+ _volumeMapper->SetVolumeRayCastFunction(_compositeFunctionMIP);
+ _volumerendererdata->changeCompositeMIPFunction(1);
+ }
+
+ if(_interpolation)
+ {
+ _volumeProperty->SetInterpolationTypeToLinear();
+ }
+ else
+ {
+ _volumeProperty->SetInterpolationTypeToNearest();
+ }
+
+ if(_shade)
+ {
+ _volumeProperty->ShadeOn();
+ }
+ else
+ {
+ _volumeProperty->ShadeOff();
+ }
+
+}
//-------------------------------------------------------------------
void vtkClipping3DDataViewer::Configure()
{
Configure_Tissue();
Configure_Volume();
-
+ Configure_VolumeBox();
// An outline provides context around the data.
//
- _outlineData = vtkOutlineFilter::New();
+ _outlineData = vtkOutlineFilter::New();
_outlineData->SetInput((vtkDataSet *) _vtkmprbasedata->GetImageData() );
- _mapOutline = vtkPolyDataMapper::New();
- _mapOutline->SetInput(_outlineData->GetOutput());
- _outline = vtkActor::New();
+ _mapOutline = vtkPolyDataMapper::New();
+
+//EED 4 nov 2015 Estelle
+// _mapOutline->SetInput(_outlineData->GetOutput());
+
+ _outline = vtkActor::New();
_outline->SetMapper(_mapOutline);
_outline->GetProperty()->SetColor(0,0,0);
}
}
//--------------------------------------------------------------------
std::vector<double>* vtkClipping3DDataViewer::GetBlueColorsOfColorTransferenceFVector()
-{
+{
return &blueColorsOfColorTransferenceFVector;
}
//--------------------------------------------------------------------
/*
int i=0,xi,yi,r,g,b,gValue;
vtkImageData *imagedata = this->_vtkmprbasedata->GetImageData();
-
+
HistogramDialog* hDlg=new HistogramDialog(NULL,_T("Histogram Dialog"),imagedata);
- //
+ //
// put in a method
//
int tfSize=this->greyValuesTransferenceFVector.size();
hDlg->addPointToTransferenceFunction(g,in*100);
i++;
}
-
+
}
int ctfSize=this->redColorsOfColorTransferenceFVector.size();
//setting variables if the user wants to do refresh
hDlg->setCTF(_ctfun);
hDlg->setTF(_tfun);
-
+
//
// when the user had changed the transference Function
//
if(hDlg->ShowModal()== wxID_OK )
- {
+ {
// -- vtkPiecewiseFunction --
this->_tfun->RemoveAllPoints();
greyValuesTransferenceFVector.clear();
intensityValuesTransferenceFVector.clear();
-
+
int nTFPoints=hDlg->getSizeTransferenceFunction();
i=0;
while(i<nTFPoints)
greyValuesTransferenceFVector.push_back(xi);
intensityValuesTransferenceFVector.push_back(yi/100.0);
i++;
- }
+ }
// -- vtkColorTransferFunction --
this->_ctfun->RemoveAllPoints ();
//clean colors
greyValueColorsOfColorTransferenceFVector.clear();
int nCTFpoints=hDlg->getSizeBarColor();
- i=0;
+ i=0;
while(i<nCTFpoints)
{
hDlg->getDataBarColorPoint(i,xi,r,g,b);
greyValueColorsOfColorTransferenceFVector.push_back(xi);
i++;
}
-
+
this->_volumeMapper->Update();
this->_newvol->Update();
}
-
+
else
{
-
+
if(hDlg->getRefreshed())
{
int i=0,size;
double in2=intensityValuesTransferenceFVector[i];
this->_tfun->AddPoint( grey1 , in2 );
}
-
+
// -- vtkColorTransferFunction --
_ctfun->RemoveAllPoints ();
-
+
i=0;
- size=greyValueColorsOfColorTransferenceFVector.size();
+ size=greyValueColorsOfColorTransferenceFVector.size();
for(i=0;i<size;i++)
{
double grey2=(greyValueColorsOfColorTransferenceFVector)[i];
this->_volumeMapper->Update();
this->_newvol->Update();
}
-
+
}
hDlg->Destroy();
*/
// -- vtkPiecewiseFunction --
this->_tfun->RemoveAllPoints();
- fscanf(ff,"%s",tmp);
+ fscanf(ff,"%s",tmp);
while ( strcmp(tmp,"--")!=0 )
{
x=atof(tmp);
- fscanf(ff,"%s",tmp);
+ fscanf(ff,"%s",tmp);
val=atof(tmp);
this->_tfun->AddPoint( x*max , val );
- fscanf(ff,"%s",tmp);
+ fscanf(ff,"%s",tmp);
}
this->_ctfun->RemoveAllPoints ();
// -- vtkColorTransferFunction --
while ( !feof(ff))
{
- fscanf(ff,"%s",tmp);
+ fscanf(ff,"%s",tmp);
x=atof(tmp);
- fscanf(ff,"%s",tmp);
+ fscanf(ff,"%s",tmp);
r=atof(tmp);
- fscanf(ff,"%s",tmp);
+ fscanf(ff,"%s",tmp);
g=atof(tmp);
- fscanf(ff,"%s",tmp);
+ fscanf(ff,"%s",tmp);
b=atof(tmp);
- this->_ctfun->AddRGBPoint( x*max , r,g,b );
+ this->_ctfun->AddRGBPoint( x*max , r,g,b );
}
reader->SetFileName(namefile);
reader->Update();
_tissueStripper[3]->SetInput( reader->GetPolyDataOutput() );
-}
\ No newline at end of file
+}
+
+void vtkClipping3DDataViewer::setColorTransferFunction(vtkColorTransferFunction* colortable){
+
+ _volumeProperty->SetColor(colortable);
+
+}
+
+
+void vtkClipping3DDataViewer::SetRayCasting(bool active)
+{
+ _isRayCasting = active;
+
+}
+
+void vtkClipping3DDataViewer::SetMIPActive(bool active)
+{
+ _isMIP = active;
+}
+
+void vtkClipping3DDataViewer::SetInterpolation(bool active)
+{
+ _interpolation = active;
+}
+
+void vtkClipping3DDataViewer::SetShade(bool active)
+{
+ _shade = active;
+}
+
git://git.creatis.insa-lyon.fr / creaMaracasVisu.git / blobdiff