]> Creatis software - creaContours.git/blobdiff - lib/Interface_ManagerContour_NDimensions/wxContourMainFrame.cxx
OK
[creaContours.git] / lib / Interface_ManagerContour_NDimensions / wxContourMainFrame.cxx
index 34182bdfeca0d59f525d589d6bfb145e2f30dcd4..365cd88bf99817ae26d9000b2435a21b3e6d6215 100644 (file)
@@ -58,6 +58,8 @@ char wxContourMainFrame::COPY = 'C';
                _refLineControl                         = NULL;
                _refLineModel                           = NULL;
                _refLineView                            = NULL;
+               _thresholdActor                         = NULL;
+               _actorPresent                           = false;
 
                // set up default notebook style
 //             m_notebook_style =wxAUI_NB_TAB_SPLIT | wxAUI_NB_TAB_EXTERNAL_MOVE | wxNO_BORDER;
@@ -77,6 +79,8 @@ char wxContourMainFrame::COPY = 'C';
                _refLineControl                         = NULL;
                _refLineModel                           = NULL;
                _refLineView                            = NULL;
+               _thresholdActor                         = NULL;
+               _actorPresent                           = false;
 
                // set up default notebook style
                m_notebook_style =wxAUI_NB_TAB_SPLIT | wxAUI_NB_TAB_EXTERNAL_MOVE | wxNO_BORDER;
@@ -112,7 +116,7 @@ void wxContourMainFrame::setVectImages( std::vector<vtkImageData*> images ){
        std::vector<std::string> conceptNameVect;
        std::vector<int> conceptSizeVect;
 
-       wxAuiNotebook * notebook = this->createNotebook();              
+       notebook = this->createNotebook();              
 
 
        kernelManager = new KernelManagerContour( images , _datadir+"/data/" , strCreaContourDataTmp );
@@ -133,7 +137,7 @@ void wxContourMainFrame::setVectImages( std::vector<vtkImageData*> images ){
        //*******************changeInstant();
        _theViewPanel->setVerticalConcept( "Axe Depth", data->getMinValue(), data->getMaxValue(), data->getMinShowedValue(),  data->getMaxShowedValue(), data->getActualValue() );
 
-       this->configurePanels( notebook );
+       this->configurePanels( );
 
 }
 
@@ -199,45 +203,44 @@ void wxContourMainFrame::setVectImages( std::vector<vtkImageData*> images ){
        //------------------------------------------------------------------------------------------------------------
        // Creational and initialization methods 
        //------------------------------------------------------------------------------------------------------------
-       bool wxContourMainFrame :: configurePanels(wxAuiNotebook* theNoteBook)
+       bool wxContourMainFrame :: configurePanels( )
        {
                bool configured = _theViewPanel!=NULL;
                                        
                configured &= _theViewPanel!=NULL;
                if( _theViewPanel!=NULL )
                {
-                       theNoteBook->AddPage( _theViewPanel, wxT("       View       ") );               
+                       notebook->AddPage( _theViewPanel, wxT("       View       ") );          
                        m_mgr.Update();
                }
 
 
                if( configured )
                {
-                       theNoteBook->AddPage( _instantPanel, wxT("Instant Page") );
+                       notebook->AddPage( _instantPanel, wxT("Instant Page") );
                        m_mgr.Update();
                }
 
-               m_mgr.AddPane(theNoteBook, wxAuiPaneInfo().Name(wxT("notebook_content")).CenterPane().PaneBorder(false));
+               m_mgr.AddPane(notebook, wxAuiPaneInfo().Name(wxT("notebook_content")).CenterPane().PaneBorder(false));
                m_mgr.Update();
-/*             configured &= _buttonsBar!=NULL;
-               if( _buttonsBar!=NULL )
-               {                       
-                       m_mgr.AddPane(_buttonsBar, wxAuiPaneInfo().
-                                       Name(wxT("TB")).Caption(wxT("Buttons Bar")).
-                                       ToolbarPane().Top().
-                                       LeftDockable(  ).RightDockable( false ).CloseButton(false));
-                       m_mgr.Update();
-               }*/
-               //CreateStatusBar();
-                        
+                
                SetMinSize(wxSize(300,300));
                m_mgr.Update();
                return configured;
        }
 
-       void wxContourMainFrame :: setNotebook( wxAuiNotebook * noteBook )
-       {
 
+       bool wxContourMainFrame::addNewPanel(wxPanel* panel)
+               {
+               bool configured = panel!=NULL;
+               configured &= panel!=NULL;
+               
+               if( configured )
+               {
+                       notebook->AddPage( panel, wxT("       Viewer       "), true );
+                       m_mgr.Update();
+               }
+               return configured;
        }
 
        //------------------------------------------------------------------------------------------------------------
@@ -593,11 +596,12 @@ void wxContourMainFrame::changeInstant()
                updateInstantOutlines();
                updateInstantImageData();
                updateInstantAxes();
-
-               
+               updateActor();
 }
-void wxContourMainFrame :: updateInstantOutlines(){
-       
+
+
+void wxContourMainFrame :: updateInstantOutlines()
+{
        _theViewPanel->removeSceneContours();
        _theViewPanel->addNameWrapperToScene();
 
@@ -614,7 +618,8 @@ std::string wxContourMainFrame::getNameWrapping(int i){
        return kernelManager->getNameWrapping(i);
 }
 
-void wxContourMainFrame::updateInstantImageData(){
+void wxContourMainFrame::updateInstantImageData()
+{
        std::vector<int> inst;
        _instantPanel->getInstant(inst);
 
@@ -628,7 +633,84 @@ void wxContourMainFrame::updateInstantImageData(){
        }
 }
 
-void wxContourMainFrame::updateInstantAxes(){
+void wxContourMainFrame::updateInstantAxes()
+{
+}
+
+void wxContourMainFrame :: updateActor()
+{
+               if (_actorPresent)
+               {
+                       int minMax[2];
+                       _pannew->onThresholdInstantChange(minMax);
+
+                       int     z = _theViewPanel->GetZ();
+
+                       vtkImageData * img = getImageData();
+                       vtkImageReslice *imageReslice = vtkImageReslice::New();
+
+                       imageReslice->SetInput( img );
+                       imageReslice->SetInformationInput(img);
+                       imageReslice->SetResliceAxesDirectionCosines(1,0,0, 0,1,0 ,0,0,1);
+                       imageReslice->SetResliceAxesOrigin(0,0,z);
+                       imageReslice->SetOutputDimensionality(2);
+                       imageReslice->SetInterpolationModeToLinear();
+
+                       img = imageReslice->GetOutput();
+                       img->Update();
+                       img->UpdateInformation();
+
+                       wxVtkBaseView * baseView = _theViewPanel->getWxVtkBaseView();
+                       
+                       double range[2];
+                       img->GetScalarRange(range);
+
+                       int min = floor (range[0]);
+                       int max = ceil (range[1]);
+
+                       //Lookup Table
+                       vtkLookupTable *lookup = vtkLookupTable::New();
+                       lookup->SetNumberOfTableValues(max+1);
+                       lookup->SetTableRange(range); 
+                       lookup->SetAlphaRange(0, 1);
+                       lookup->SetValueRange(0, 1);
+                       lookup->SetSaturationRange(0, 0); 
+                       lookup->SetRampToLinear( );
+
+                       //Assign a fake color for the upper image, and set the white as transparent
+                       int i;
+                       int minVal = minMax[0];
+                       int maxVal = minMax[1];
+
+                       for(i = min; i <= max; i++)
+                       {
+                               if( i >= minVal && i <= maxVal )
+                               {
+                                       lookup->SetTableValue(i, 1.0, 0.0, 0.0, 1);
+                               }
+                               else if( i >= min && i < minVal )
+                               {
+                                       lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                               }
+                               else if( i > maxVal && i < max )
+                               {
+                                       lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                               }
+                               else
+                               {
+                                       lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                               }
+                       }
+                               
+                       lookup->Build( );
+
+                       vtkImageMapToColors *mapperImage = vtkImageMapToColors::New( );
+                       mapperImage->SetLookupTable( lookup );
+                       mapperImage->SetInput( img );
+                       //mapperImage->SetOutputFormatToRGBA( );
+
+                       _thresholdActor->SetInput( mapperImage->GetOutput() );
+               }
 }
 
 void wxContourMainFrame::onChangeDeep(int val){
@@ -942,25 +1024,25 @@ void wxContourMainFrame::SegmentationOneSliceITK(int x, int y, int z, wxString d
   typedef unsigned short OutputPixelType2;
   typedef itk::Image< OutputPixelType2, Dimension > OutputImageType2;
 
-  //Definición del thresholder
+  //DefiniciÛn del thresholder
   typedef itk::BinaryThresholdImageFilter< 
                         InternalImageType, 
                         OutputImageType    >    ThresholdingFilterType;
   
-  //Definición del primer filtro de conversión de pixeles
+  //DefiniciÛn del primer filtro de conversiÛn de pixeles
   typedef itk::CastImageFilter<
                OutputImageType, OutputImageType2 >  CastFilterType;
 
-  //Definición del segundo tipo de conversión de pixeles
+  //DefiniciÛn del segundo tipo de conversiÛn de pixeles
   typedef itk::CastImageFilter<
                OutputImageType2, InternalImageType >  CastFilterType2;
 
-  //Tercer tipo de conversión
+  //Tercer tipo de conversiÛn
   typedef itk::RescaleIntensityImageFilter< 
                                InternalImageType, 
                                OutputImageType >   CastFilterType3;
 
-  //Cuarto tipo de conversión
+  //Cuarto tipo de conversiÛn
   typedef itk::RescaleIntensityImageFilter< 
                                OutputImageType, 
                                OutputImageType >   CastFilterType4;
@@ -973,7 +1055,7 @@ void wxContourMainFrame::SegmentationOneSliceITK(int x, int y, int z, wxString d
   thresholder->SetOutsideValue(  255  );
   thresholder->SetInsideValue(  0 );
 
-  //Definción de conexiónes entre VTK e ITK y el writer
+  //DefinciÛn de conexiÛnes entre VTK e ITK y el writer
   typedef itk::VTKImageToImageFilter<OutputImageType2> ConnectorType;
   typedef itk::ImageToVTKImageFilter<OutputImageType> ConnectorType2;
   typedef  itk::ImageFileWriter<  OutputImageType  > WriterType;
@@ -1026,23 +1108,23 @@ void wxContourMainFrame::SegmentationOneSliceITK(int x, int y, int z, wxString d
                               InternalImageType, 
                               InternalImageType >    ZeroCrossingFilterType;
 ZeroCrossingFilterType::Pointer zeroCrossing =
-                                                                       ZeroCrossingFilterType::New();
+                                               ZeroCrossingFilterType::New();
 
-const double propagationScaling = atof( prop );
+const double propagationScaling = atof( crea::wx2std(prop).c_str() );
 
   geodesicActiveContour->SetPropagationScaling( propagationScaling );
   geodesicActiveContour->SetCurvatureScaling( 1.0 );
   geodesicActiveContour->SetAdvectionScaling( 1.0 );
 
   geodesicActiveContour->SetMaximumRMSError( 0.02 );
-  int it=atoi( iter );
+  int it=atoi(  crea::wx2std(iter).c_str() );
   geodesicActiveContour->SetNumberOfIterations( it );
 
   smoothing->SetInput( filter2->GetOutput() );
   gradientMagnitude->SetInput( smoothing->GetOutput() );
   sigmoid->SetInput( gradientMagnitude->GetOutput() );
   fastMarching->SetInput( sigmoid->GetOutput() );
-  geodesicActiveContour->SetInput(  fastMarching->GetOutput() );
+  geodesicActiveContour->SetInput( fastMarching->GetOutput() );
   geodesicActiveContour->SetFeatureImage( sigmoid->GetOutput() );
   
   zeroCrossing->SetInput( geodesicActiveContour->GetOutput() );
@@ -1056,11 +1138,11 @@ const double propagationScaling = atof( prop );
   smoothing->SetConductanceParameter( 9.0 );
 
 
-  const double sigma = atof( sigm );
+  const double sigma = atof(  crea::wx2std(sigm).c_str() );
   gradientMagnitude->SetSigma(  sigma  );
 
-  const double alpha =  atof( alf );
-  const double beta  =  atof( bet );
+  const double alpha =  atof(  crea::wx2std(alf).c_str() );
+  const double beta  =  atof(  crea::wx2std(bet).c_str() );
 
   sigmoid->SetAlpha( alpha );
   sigmoid->SetBeta(  beta  );
@@ -1074,7 +1156,7 @@ const double propagationScaling = atof( prop );
   seedPosition[0] = x;
   seedPosition[1] = y;
 
-  const double initialDistance = atof( distanc );
+  const double initialDistance = atof( crea::wx2std(distanc).c_str() );
 
   NodeType node;
 
@@ -1086,7 +1168,7 @@ const double propagationScaling = atof( prop );
   seeds->Initialize();
   seeds->InsertElement( 0, node );
 
-  fastMarching->SetTrialPoints(  seeds  );
+  fastMarching->SetTrialPoints( seeds );
 
   fastMarching->SetSpeedConstant( 1.0 );
   
@@ -1189,7 +1271,7 @@ myfile.open ("C:/Creatis/example2.txt");
                        std::cout<<" x Anterior "<<xAct<<" x actual "<<x<<std::endl;
                std::cout<<" y Anterior "<<yAct<<" y actual "<<y<<std::endl;
                std::cout<<" x "<<p[0]<<" y "<<p[1]<<std::endl;
-                       vecZ.push_back( -900 );
+                       vecZ.push_back( 900 );
                        xAct=x;
                        yAct=y;
                        }
@@ -1213,7 +1295,7 @@ myfile.open ("C:/Creatis/example2.txt");
                std::cout<<" y Anterior "<<yAct<<" y actual "<<y<<std::endl;
                std::cout<<" x "<<p[0]<<" y "<<p[1]<<std::endl;
                vecY.push_back( p[1] );
-               vecZ.push_back( -900 );
+               vecZ.push_back( 900 );
                xAct=x;
                yAct=y;
                }
@@ -1226,7 +1308,7 @@ myfile.open ("C:/Creatis/example2.txt");
                vecX.push_back(vecXo.back());
                std::cout<<" x Siguiente "<<vecXo.back();
                vecXo.pop_back();
-               vecZ.push_back( -900 );
+               vecZ.push_back( 900 );
        }
        while(!vecYo.empty())
        {
@@ -1480,7 +1562,7 @@ myfile.close();
 //             double y=p[1];
                vecX.push_back( p[0] );
                vecY.push_back( p[1] );
-               vecZ.push_back( -900 );
+               vecZ.push_back( 900 );
 //                     vecZ.push_back( p[2] );
        }
 
@@ -1577,49 +1659,52 @@ void wxContourMainFrame::onSegmentationAllSlice(int minZ,int maxZ,int isovalue,i
 //Creates a reference line to be used as the axis of the mirroring
 void wxContourMainFrame::referenceLine()
 {
-    wxBusyCursor wait;
-       double spc[3];
+       if (_refLineControl==NULL)
+       {
+               wxBusyCursor wait;
+               double spc[3];
 
-       vtkImageData *vtkimagedata = _theViewPanel->getImageData();
-       vtkimagedata->GetSpacing(spc);
+               vtkImageData *vtkimagedata = _theViewPanel->getImageData();
+               vtkimagedata->GetSpacing(spc);
 
-       wxVtkBaseView * viewer2D = _theViewPanel->getWxVtkBaseView();
+               wxVtkBaseView * viewer2D = _theViewPanel->getWxVtkBaseView();
 
-       _refLineControl = new manualLineControler();
-       _refLineModel   = new manualContourModelLine();
-       _refLineView = new manualViewLine();
-       _refLineView->SetModel( _refLineModel );
-       _refLineView->SetWxVtkBaseView( viewer2D );
-       _refLineView->SetRange( 2 );
-       _refLineView->SetZ( 1200 );
+               _refLineControl = new manualLineControler();
+               _refLineModel   = new manualContourModelLine();
+               _refLineView = new manualViewLine();
+               _refLineView->SetModel( _refLineModel );
+               _refLineView->SetWxVtkBaseView( viewer2D );
+               _refLineView->SetRange( 2 );
+               _refLineView->SetZ( 1200 );
 
-       _refLineView->SetSpacing(spc);
+               _refLineView->SetSpacing(spc);
 
-       _refLineView->SetColorNormalContour(0, 0, 1);
-       _refLineView->SetColorEditContour(0, 0.5, 0.5);
-       _refLineView->SetWidthLine(4);
+               _refLineView->SetColorNormalContour(0, 0, 1);
+               _refLineView->SetColorEditContour(0, 0.5, 0.5);
+               _refLineView->SetWidthLine(4);
 
-       _refLineControl->SetModelView( _refLineModel , _refLineView );
-       viewer2D->GetInteractorStyleBaseView()->AddInteractorStyleMaracas( _refLineControl );
-       _refLineModel->SetCloseContour(false);
-       _refLineControl->CreateNewManualContour();
+               _refLineControl->SetModelView( _refLineModel , _refLineView );
+               viewer2D->GetInteractorStyleBaseView()->AddInteractorStyleMaracas( _refLineControl );
+               _refLineModel->SetCloseContour(false);
+               _refLineControl->CreateNewManualContour();
 
-       double z = _refLineControl->GetZ();
-/*
-       _refLineControl->AddPoint(50,100,z);
-       _refLineControl->AddPoint(80,40,z);
-       _refLineControl->SetCompleteCreation(true);
+               double z = _refLineControl->GetZ();
+       /*
+               _refLineControl->AddPoint(50,100,z);
+               _refLineControl->AddPoint(80,40,z);
+               _refLineControl->SetCompleteCreation(true);
 
-       _refLineModel->AddPoint(50,100,z);
-       _refLineModel->AddPoint(80,40,z);
-       _refLineModel->UpdateSpline();
+               _refLineModel->AddPoint(50,100,z);
+               _refLineModel->AddPoint(80,40,z);
+               _refLineModel->UpdateSpline();
 
-       _refLineView->UpdateViewPoint(0);
-       _refLineView->UpdateViewPoint(1);
-*/
+               _refLineView->UpdateViewPoint(0);
+               _refLineView->UpdateViewPoint(1);
+       */
 
-       _refLineControl->SetActive(true);
-       _refLineView->RefreshContour();
+               _refLineControl->SetActive(true);
+               _refLineView->RefreshContour();
+       }
 }
 
 //Hides the referenceLine
@@ -1627,9 +1712,14 @@ void wxContourMainFrame::refLineHide()
 {
        if (_refLineControl!=NULL)
        {
-               _refLineView->RemoveCompleteContourActor();
-               _refLineControl->SetEditable(false);
-               _refLineControl->SetActive(false);
+               if (_refLineControl->GetActive())
+               {
+                       _refLineView->RemoveCompleteContourActor();
+                       _refLineControl->SetEditable(false);
+                       _refLineControl->SetActive(false);
+
+                       RefreshInterface();
+               }
        }
 }
 
@@ -1638,18 +1728,23 @@ void wxContourMainFrame::refLineShow()
 {
        if (_refLineControl!=NULL)
        {
-               wxVtkBaseView * viewer2D = _theViewPanel->getWxVtkBaseView();
+               if (_refLineControl->GetActive()==false)
+               {
+                       wxVtkBaseView * viewer2D = _theViewPanel->getWxVtkBaseView();
 
-               _refLineView->SetModel( _refLineModel );
-               _refLineView->SetWxVtkBaseView( viewer2D );
-               _refLineView->SetRange( 2 );
-               _refLineView->SetZ( 1200 );
+                       _refLineView->SetModel( _refLineModel );
+                       _refLineView->SetWxVtkBaseView( viewer2D );
+                       _refLineView->SetRange( 2 );
+                       _refLineView->SetZ( 1200 );
 
-               _refLineControl->SetModelView( _refLineModel , _refLineView );
-               viewer2D->GetInteractorStyleBaseView()->AddInteractorStyleMaracas( _refLineControl );
-               _refLineControl->CreateNewManualContour();
-               _refLineControl->SetActive(true);
-               _refLineView->RefreshContour();
+                       _refLineControl->SetModelView( _refLineModel , _refLineView );
+                       viewer2D->GetInteractorStyleBaseView()->AddInteractorStyleMaracas( _refLineControl );
+                       _refLineControl->CreateNewManualContour();
+                       _refLineControl->SetActive(true);
+                       _refLineView->RefreshContour();
+
+                       RefreshInterface();
+               }
        }
 }
 
@@ -1698,57 +1793,208 @@ void wxContourMainFrame::onMirror()
 
 void wxContourMainFrame::onThreshold(int minVal, int maxVal)
 {
-       
-       vtkImageData * img = _theViewPanel->getImageData();
-       
-       double range[2];
-       img->GetScalarRange(range);
+       if (!_actorPresent)
+       {
+               int     z = _theViewPanel->GetZ();
 
-       //vtkImageData * fakeImage = new vtkImageData();
+               vtkImageData * img = getImageData();
+               vtkImageReslice *imageReslice = vtkImageReslice::New();
 
-       //Lookup Table
-       vtkLookupTable *lookup = vtkLookupTable::New();
-       lookup->SetNumberOfTableValues(range[1]);
-       lookup->SetTableRange(range); 
-       lookup->SetAlphaRange(0, 1);
-       lookup->SetValueRange(0, 1);
-       lookup->SetSaturationRange(0, 0); 
-       lookup->SetRampToLinear( );
+               imageReslice->SetInput( img );
+               imageReslice->SetInformationInput(img);
+               imageReslice->SetResliceAxesDirectionCosines(1,0,0, 0,1,0 ,0,0,1);
+               imageReslice->SetResliceAxesOrigin(0,0,z);
+               imageReslice->SetOutputDimensionality(2);
+               imageReslice->SetInterpolationModeToLinear();
 
+               img = imageReslice->GetOutput();
+               img->Update();
+               img->UpdateInformation();
 
-       //Assign a fake color for the upper image, and set the white as transparent
-       for(int i = range[0]; i < range[1]; i++)
-       {
-               if( i >= minVal && i <= maxVal )
+               wxVtkBaseView * baseView = _theViewPanel->getWxVtkBaseView();
+               
+               double range[2];
+               img->GetScalarRange(range);
+
+               int min = floor (range[0]);
+               int max = ceil (range[1]);
+
+               //Lookup Table
+               vtkLookupTable *lookup = vtkLookupTable::New();
+               lookup->SetNumberOfTableValues(max+1);
+               lookup->SetTableRange(range); 
+               lookup->SetAlphaRange(0, 1);
+               lookup->SetValueRange(0, 1);
+               lookup->SetSaturationRange(0, 0); 
+               lookup->SetRampToLinear( );
+
+               //Assign a fake color for the upper image, and set the white as transparent
+               int i;
+               for(i = min; i <= max; i++)
                {
-                       lookup->SetTableValue(i, 1.0, 0.0, 0.0, 1);
+                       if( i >= minVal && i <= maxVal )
+                       {
+                               lookup->SetTableValue(i, 1.0, 0.0, 0.0, 1);
+                       }
+                       else if( i >= min && i < minVal )
+                       {
+                               lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                       }
+                       else if( i > maxVal && i < max )
+                       {
+                               lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                       }
+                       else
+                       {
+                               lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                       }
                }
-               else if( i >= range[0] && i < minVal )
+                       
+               lookup->Build( );
+
+               vtkImageMapToColors *mapperImage = vtkImageMapToColors::New( );
+               mapperImage->SetLookupTable( lookup );
+               mapperImage->SetInput( img );
+               //mapperImage->SetOutputFormatToRGBA( );
+
+               if (_thresholdActor==NULL)
                {
-                       lookup->SetTableValue(i, 0.0, 0.0, 0.0, 0); //transparent
+                       _thresholdActor = vtkImageActor::New( );
+                       _thresholdActor->SetInput( mapperImage->GetOutput() );
+                       _thresholdActor->SetOpacity( 0.6 );
+                       _thresholdActor->InterpolateOn( ); 
                }
-               else if( i > maxVal && i < range[1] )
+               else
                {
-                       lookup->SetTableValue(i, 0.0, 0.0, 0.0, 0); //transparent
+                       _thresholdActor->SetInput( mapperImage->GetOutput() );
                }
+
+               baseView->GetRenderer()->AddActor( _thresholdActor );
+               _actorPresent = true;
        }
+
+       RefreshInterface();
+/*
+       vtkImageViewer2* viewer = vtkImageViewer2::New();        
+       viewer->SetInput( upperImageActor->GetInput() ); 
+       viewer->SetColorLevel((range[1]-range[0])/2);
+       viewer->SetColorWindow(range[1]);
+       viewer->GetRenderer()->AddActor( upperImageActor );
+       viewer->Render();
+*/
+
+}
+
+void wxContourMainFrame::onThresholdChange(int minVal, int maxVal)
+{
+       if (_actorPresent)
+       {
+               int     z = _theViewPanel->GetZ();
+
+               vtkImageData * img = getImageData();
+               vtkImageReslice *imageReslice = vtkImageReslice::New();
+
+               imageReslice->SetInput( img );
+               imageReslice->SetInformationInput(img);
+               imageReslice->SetResliceAxesDirectionCosines(1,0,0, 0,1,0 ,0,0,1);
+               imageReslice->SetResliceAxesOrigin(0,0,z);
+               imageReslice->SetOutputDimensionality(2);
+               imageReslice->SetInterpolationModeToLinear();
+
+               img = imageReslice->GetOutput();
+               img->Update();
+               img->UpdateInformation();
+
+               wxVtkBaseView * baseView = _theViewPanel->getWxVtkBaseView();
                
+               double range[2];
+               img->GetScalarRange(range);
+
+               int min = floor (range[0]);
+               int max = ceil (range[1]);
+
+               //Lookup Table
+               vtkLookupTable *lookup = vtkLookupTable::New();
+               lookup->SetNumberOfTableValues(max+1);
+               lookup->SetTableRange(range); 
+               lookup->SetAlphaRange(0, 1);
+               lookup->SetValueRange(0, 1);
+               lookup->SetSaturationRange(0, 0); 
+               lookup->SetRampToLinear( );
+
+               //Assign a fake color for the upper image, and set the white as transparent
+               int i;
+               for(i = min; i <= max; i++)
+               {
+                       if( i >= minVal && i <= maxVal )
+                       {
+                               lookup->SetTableValue(i, 1.0, 0.0, 0.0, 1);
+                       }
+                       else if( i >= min && i < minVal )
+                       {
+                               lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                       }
+                       else if( i > maxVal && i < max )
+                       {
+                               lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                       }
+                       else
+                       {
+                               lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                       }
+               }
+                       
                lookup->Build( );
 
                vtkImageMapToColors *mapperImage = vtkImageMapToColors::New( );
-               mapperImage->SetInput( img );
                mapperImage->SetLookupTable( lookup );
+               mapperImage->SetInput( img );
+               //mapperImage->SetOutputFormatToRGBA( );
 
-               vtkImageActor * upperImageActor = vtkImageActor::New( );
-               upperImageActor->SetInput( mapperImage->GetOutput() );
-               upperImageActor->SetOpacity( 1 );
+               _thresholdActor->SetInput( mapperImage->GetOutput() );
+       }
 
-               //viewer->SetInput( myImage ); viewer->GetRenderer()->AddActor( myActor );ACTOR!!!
+       RefreshInterface();
+}
 
-               //baseImageActor = vtkImageActor::New( );
-               //baseImageActor->SetInput( mBox->bbGetInputImage1() );
+void wxContourMainFrame::onThresholdInterpolation(bool interpolate)
+{
+       if (_thresholdActor!=NULL)
+       {
+               if (interpolate)
+               {
+                       _thresholdActor->InterpolateOn( );
+               }
+
+               else
+               {
+                       _thresholdActor->InterpolateOff( );
+               }
 
                RefreshInterface();
+       }
+}
+
+void wxContourMainFrame::onThresholdChangeOpacity (int opacity)
+{
+       if (_actorPresent)
+       {
+               _thresholdActor->SetOpacity(opacity*0.1);
+       }
+
+       RefreshInterface();
+}
+
+void wxContourMainFrame::onThresholdRemove()
+{
+       if (_actorPresent)
+       {
+               wxVtkBaseView * baseView = _theViewPanel->getWxVtkBaseView();
+               baseView->GetRenderer()->RemoveActor( _thresholdActor );
+               _actorPresent = false;
+       }
+
+       RefreshInterface();
 }
 
 void wxContourMainFrame::showAxis(bool show)
@@ -1798,7 +2044,6 @@ int wxContourMainFrame::getWindowLevel(){
 void wxContourMainFrame::onBrigthnessColorWindowLevel(int colorwindow,int windowlevel){
        _theViewPanel->onBrigthnessColorWindowLevel(colorwindow, windowlevel);
 
-       
        RefreshInterface();
 }
 
@@ -2093,7 +2338,10 @@ void wxContourMainFrame::onInformationContourLabels(int typeContourGroup,int sel
 
 void wxContourMainFrame::onSaveResults(std::string directory,std::string namefile, std::string filename, int typeContourGroup){
     SaveValuesXYZ( directory , namefile, typeContourGroup );
-       SaveImageResult( directory , namefile, typeContourGroup);
+
+       vtkImageData* mask, *value;
+       this->getMaskValue(mask, value, typeContourGroup);
+       this->SaveImageResult( directory , namefile, mask, value);
 
        
        FILE *pFile=fopen(filename.c_str(),"w+");
@@ -2253,8 +2501,7 @@ void wxContourMainFrame::SaveValuesXYZ(std::string directory,std::string namefil
        delete contourextractdata;
 }
 
-void wxContourMainFrame::SaveImageResult(std::string directory,std::string namefile, int typeContourGroup)
-{
+void wxContourMainFrame::getMaskValue(vtkImageData* mask,vtkImageData* value,int typeContourGroup){
        wxBusyCursor wait;
        ContourExtractData      *contourextractdata = new ContourExtractData( true );
        //int typeContourGroup = this->_radiolstboxContourGroup->GetSelection();
@@ -2333,20 +2580,19 @@ void wxContourMainFrame::SaveImageResult(std::string directory,std::string namef
                        }
 
                        contourextractdata->SetLstManualContourModel( lstManConModTmp );
-
-
-//                             for (ii=0 ; ii<sizeLstContourThings ; ii++)
-//                             {
-//                                     ContourThing **contourthing = lstContourThings[ii];
-//                                     lstManConMod.push_back( (*contourthing)->getModel() );
-//                             }
-//                             contourextractdata->SetLstManualContourModel( lstManConMod );
-
                        contourextractdata->CalculateImageResult(); // with actual Z
 
                } // for  iContourGroup
        } // for z
 
+       value = contourextractdata->GetVtkImageValueResult();
+       mask = contourextractdata->GetVtkImageMaskResult();
+       delete contourextractdata;
+}
+
+void wxContourMainFrame::SaveImageResult(std::string directory,std::string namefile, vtkImageData* mask,vtkImageData* value)
+{
+       
        std::string temp = directory + "/" + namefile + "-Value.mhd";
 
        wxString filename = crea::std2wx(temp);
@@ -2360,7 +2606,7 @@ void wxContourMainFrame::SaveImageResult(std::string directory,std::string namef
 
 // Image Value
        vtkMetaImageWriter *writerValueImage = vtkMetaImageWriter::New( );
-       writerValueImage->SetInput( contourextractdata->GetVtkImageValueResult() );
+       writerValueImage->SetInput( value );
        writerValueImage->SetFileName( (const char *)filename.mb_str() );
        writerValueImage->SetFileDimensionality( 3 );
        writerValueImage->Write( );
@@ -2379,15 +2625,17 @@ void wxContourMainFrame::SaveImageResult(std::string directory,std::string namef
        filename.append(_T("-Mask.mhd"));*/
 
        vtkMetaImageWriter *writerMaskImage = vtkMetaImageWriter::New( );
-       writerMaskImage->SetInput( contourextractdata->GetVtkImageMaskResult() );
+       writerMaskImage->SetInput( mask );
        writerMaskImage->SetFileName( (const char *)filename.mb_str() );
        writerMaskImage->SetFileDimensionality( 3 );
        writerMaskImage->Write( );
 
        interfMainPanel::getInstance()->setStringInfoPanel( _T("") );
-       delete contourextractdata;
+       
 }
 
+
+
 void wxContourMainFrame::onChangeInstantInstantPanel(std::string name,int value, int minshow, int maxshow){
        if( name.compare( _theViewPanel->getVerticalConceptName() ) == 0){
                _theViewPanel->setActualVertical( value);
@@ -2431,4 +2679,5 @@ void wxContourMainFrame::saveFileWithContoursAutomatique(){
        }else{
                onSave();
        }
-}
\ No newline at end of file
+}
+