X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=lib%2FmaracasVisuLib%2Fsrc%2Finterface%2FwxWindows%2Fwidgets%2FvtkClipping3DDataViewer.cxx;h=0e301e2d79a1cbfa8590a6bd525682c5cb480a42;hb=3d1c0689f4e1844c8f170d9c20f8ad38b34d98b6;hp=d1104dff66740b5c9d6c5a473427dff93b4d1e90;hpb=cea19159becbb13fffac6d87a914480c385f8026;p=creaMaracasVisu.git diff --git a/lib/maracasVisuLib/src/interface/wxWindows/widgets/vtkClipping3DDataViewer.cxx b/lib/maracasVisuLib/src/interface/wxWindows/widgets/vtkClipping3DDataViewer.cxx index d1104df..0e301e2 100644 --- a/lib/maracasVisuLib/src/interface/wxWindows/widgets/vtkClipping3DDataViewer.cxx +++ b/lib/maracasVisuLib/src/interface/wxWindows/widgets/vtkClipping3DDataViewer.cxx @@ -1,3 +1,28 @@ +/*# --------------------------------------------------------------------- +# +# 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" @@ -11,35 +36,47 @@ vtkClipping3DDataViewer::vtkClipping3DDataViewer() _representationType[i] = true; } _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() @@ -76,7 +113,7 @@ void vtkClipping3DDataViewer::Refresh() // _newvol->VisibilityOn(); - + // _volumeMapper = vtkVolumeRayCastMapper::New(); // _volumeMapper->SetInput( this->GetVtkMPRBaseData()->GetImageData() ); // _volumeMapper->SetVolumeRayCastFunction(_compositeFunction); @@ -154,6 +191,15 @@ void vtkClipping3DDataViewer::SetVisibleVolume(bool visiblevolume) _visibleVolume = visiblevolume; } +bool vtkClipping3DDataViewer::GetVisibleVolumeBox() +{ + return _visibleVolumeBox; +} +//------------------------------------------------------------------- +void vtkClipping3DDataViewer::SetVisibleVolumeBox(bool visibleBox) +{ + _visibleVolumeBox = visibleBox; +} //------------------------------------------------------------------- void vtkClipping3DDataViewer::Configure_Tissue() @@ -177,11 +223,11 @@ void vtkClipping3DDataViewer::Configure_Tissue() _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] ); @@ -193,7 +239,7 @@ void vtkClipping3DDataViewer::Configure_Tissue() // _tissueMapper[i]->Update(); } - + // vtkActor tissue @@ -227,19 +273,26 @@ void vtkClipping3DDataViewer::Configure_Volume() { // 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); @@ -251,7 +304,7 @@ void vtkClipping3DDataViewer::Configure_Volume() _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); @@ -266,14 +319,14 @@ void vtkClipping3DDataViewer::Configure_Volume() 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); @@ -292,22 +345,28 @@ void vtkClipping3DDataViewer::Configure_Volume() // 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); @@ -315,23 +374,104 @@ void vtkClipping3DDataViewer::Configure_Volume() // _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 *gv = this->GetGreyValuesTransferenceFVector(); + vector *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 *gv = this->GetGreyValuesTransferenceFVector(); + vector *op = this->GetIntensityValuesTransferenceFVector(); + + _volumerendererdata->ShowVolume(changed); + _volumerendererdata->OpacityChanged(*gv, *op); + _volumerendererdata->BoundingBoxChanged(changed); + +} + +void vtkClipping3DDataViewer::UpdateVolumeBox(vector gf, vector 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(); @@ -400,7 +540,7 @@ std::vector* vtkClipping3DDataViewer::GetGreenColorsOfColorTransferenceF } //-------------------------------------------------------------------- std::vector* vtkClipping3DDataViewer::GetBlueColorsOfColorTransferenceFVector() -{ +{ return &blueColorsOfColorTransferenceFVector; } //-------------------------------------------------------------------- @@ -430,9 +570,9 @@ void vtkClipping3DDataViewer::ReadVolumeFunctions() /* 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(); @@ -447,7 +587,7 @@ void vtkClipping3DDataViewer::ReadVolumeFunctions() hDlg->addPointToTransferenceFunction(g,in*100); i++; } - + } int ctfSize=this->redColorsOfColorTransferenceFVector.size(); @@ -469,17 +609,17 @@ void vtkClipping3DDataViewer::ReadVolumeFunctions() //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_ctfun->RemoveAllPoints (); //clean colors @@ -499,7 +639,7 @@ void vtkClipping3DDataViewer::ReadVolumeFunctions() greyValueColorsOfColorTransferenceFVector.clear(); int nCTFpoints=hDlg->getSizeBarColor(); - i=0; + i=0; while(igetDataBarColorPoint(i,xi,r,g,b); @@ -510,14 +650,14 @@ void vtkClipping3DDataViewer::ReadVolumeFunctions() greyValueColorsOfColorTransferenceFVector.push_back(xi); i++; } - + this->_volumeMapper->Update(); this->_newvol->Update(); } - + else { - + if(hDlg->getRefreshed()) { int i=0,size; @@ -531,12 +671,12 @@ void vtkClipping3DDataViewer::ReadVolumeFunctions() double in2=intensityValuesTransferenceFVector[i]; this->_tfun->AddPoint( grey1 , in2 ); } - + // -- vtkColorTransferFunction -- _ctfun->RemoveAllPoints (); - + i=0; - size=greyValueColorsOfColorTransferenceFVector.size(); + size=greyValueColorsOfColorTransferenceFVector.size(); for(i=0;i_volumeMapper->Update(); this->_newvol->Update(); } - + } hDlg->Destroy(); */ @@ -578,29 +718,29 @@ void vtkClipping3DDataViewer::ReadVolumeFunctions(char *namefile) // -- 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 ); } @@ -620,4 +760,33 @@ void vtkClipping3DDataViewer::ReadMeshVTK(char *namefile) 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; +} +