]> Creatis software - creaContours.git/blobdiff - lib/Interface_ManagerContour_NDimensions/wxContourMainFrame.cxx
#3327 creaContours Feature New Normal - Jump Slice
[creaContours.git] / lib / Interface_ManagerContour_NDimensions / wxContourMainFrame.cxx
index 760f5431b8e2618120e0b19db779df4059a8f100..33db75f8652b799e2e366ebc8186d9058709bb5e 100644 (file)
@@ -94,6 +94,8 @@ char wxContourMainFrame::COPY = 'C';
                _viewThresholdImage                     = NULL;
                _viewThresholdImagePanel        = NULL;
                _viewColorLayerImagePanel       = NULL;
+
+               _frameShowResultImages          = NULL;
        }
 
        wxContourMainFrame::wxContourMainFrame(wxWindow* parent, wxWindowID id,const wxString& title,const wxPoint& pos,const wxSize& size,std::vector<vtkImageData*> images, long style,std::string datadir)
@@ -103,7 +105,6 @@ char wxContourMainFrame::COPY = 'C';
                _contourGroup                           = 1;
                _creatingContoursActive         = false;
                _theViewPanel                           = NULL;
-
                _instantPanel                           = NULL;
                _refLineControl                         = NULL;
                _refLineModel                           = NULL;
@@ -123,7 +124,7 @@ char wxContourMainFrame::COPY = 'C';
                inredo                                          = 0;
                inundo                                          = 0;
                _pannew = interfMainPanel::getInstance(this,datadir+"/data/Icons");//, eventHandler);
-               _performingOperation            = new PerformingOperation();
+               _performingOperation        = new PerformingOperation();
                if(images.size() > 0)
                {
                        this->setVectImages(images);
@@ -134,15 +135,15 @@ char wxContourMainFrame::COPY = 'C';
 void wxContourMainFrame::setVectImages(std::vector<vtkImageData*> imgs)
 {
        _images = imgs;
-       #if defined(__GNUC__)
-               std::string str_home(getenv("HOME"));
-       #elif defined(_WIN32)
+       #if defined(_WIN32)
                std::string str_home(getenv("USERPROFILE"));
+       #elif defined(__GNUC__)
+               std::string str_home(getenv("HOME"));
        #endif
-       std::string strCreaContourDataTmp = str_home + "/.creaContourDataTemp/";
-       std::vector<std::string> conceptNameVect;
-       std::vector<int> conceptSizeVect;
-       notebook = this->createNotebook();
+       std::string                             strCreaContourDataTmp = str_home + "/.creaContourDataTemp/";
+       std::vector<std::string>        conceptNameVect;
+       std::vector<int>                        conceptSizeVect;
+       notebook                = this->createNotebook();
        kernelManager   = new KernelManagerContour( imgs , _datadir+"/data/" , strCreaContourDataTmp );
        _instantPanel   = new wxInstantChooserPanel( notebook, "Instant Chooser", true, false, "c" );
        _theViewPanel   = new wxContourViewPanel( kernelManager->getSourceImage(), notebook );
@@ -320,20 +321,25 @@ void wxContourMainFrame::onCreateContourSpline( ){
        createContour( 1 );
 
 }
-void wxContourMainFrame::onCreateContourRectangle( ){
+
+void wxContourMainFrame::onCreateContourRectangle( )
+{
     //JCP 20-10-08 Undo redo implementation
        saveState();
        //JCP 20-10-08 Undo redo implementation
        createContour( 2 );
 }
-void wxContourMainFrame::onCreateContourCircle( ){
+
+void wxContourMainFrame::onCreateContourCircle( )
+{
     //JCP 20-10-08 Undo redo implementation
        saveState();
        //JCP 20-10-08 Undo redo implementation
        createContour( 3 );
 }
 
-void wxContourMainFrame::onCreateContourLine( ){
+void wxContourMainFrame::onCreateContourLine( )
+{
     //JCP 20-10-08 Undo redo implementation
        saveState();
        //JCP 20-10-08 Undo redo implementation
@@ -341,7 +347,8 @@ void wxContourMainFrame::onCreateContourLine( ){
 }
 
 // RaC 09-09 --------------------------------------
-void wxContourMainFrame::onCreateContourPoints( ){
+void wxContourMainFrame::onCreateContourPoints( )
+{
     //JCP 20-10-08 Undo redo implementation
        saveState();
        //JCP 20-10-08 Undo redo implementation
@@ -350,7 +357,8 @@ void wxContourMainFrame::onCreateContourPoints( ){
 // RaC 09-09 --------------------------------------
 
 // RaC 10-09 --------------------------------------
-void wxContourMainFrame::onCreateContourPolygon( ){
+void wxContourMainFrame::onCreateContourPolygon( )
+{
     //JCP 20-10-08 Undo redo implementation
        saveState();
        //JCP 20-10-08 Undo redo implementation
@@ -385,15 +393,19 @@ void wxContourMainFrame :: onCreateContourBullEye(wxPanel* panel )
                _theViewPanel->getSceneManager()->configureViewControlTo( theName, manModelContour,spc , 4) ;
        }
 }
-void wxContourMainFrame :: createContour( int typeContour )
+
+void wxContourMainFrame::createContour( int typeContour )
 {
        //Creating the manualContourModel and including in the model
        manualBaseModel * manModelContour = kernelManager->factoryManualContourModel(typeContour);
-
        std::vector<int> instantVector;
        _instantPanel->getInstant( instantVector );
-       std::string theName = kernelManager->createOutline(manModelContour, instantVector);
 
+       std::vector<int> tempVector;
+       getInstantVector( tempVector );
+       manModelContour->SetLabel2( std::to_string(tempVector[1]) );
+
+       std::string theName = kernelManager->createOutline(manModelContour, instantVector);
        /*std::vector<int> instantVector;
        _instantPanel->getInstant( instantVector );
        std::string theName;
@@ -407,9 +419,9 @@ void wxContourMainFrame :: createContour( int typeContour )
                //Adding the manualViewContour to interface objects structure
                //_sceneManager->setControlActiveStateOfALL( false );//This call is being done here because if the ROI is created underneath the previously created ROIS will still be active.
                _theViewPanel->getSceneManager()->configureViewControlTo( theName, manModelContour,spc , typeContour) ;
-       }
-
+       } // if addedModel
 }
+
 manualBaseModel * wxContourMainFrame::factoryManualContourModel(wxPanel* panel){
 
        manualContourModelBullEye *manModelContourBullEye = new manualContourModelBullEye();
@@ -541,11 +553,17 @@ void wxContourMainFrame :: deleteContour( std::string theKeyName )
 //JCP 21 - 11 - 08
 }
 
-void wxContourMainFrame::onDeleteContoursActSlice(){
+void wxContourMainFrame::onDeleteContoursActSlice()
+{
        //JCP 20-10-08 Undo redo implementation
        saveState();
        //JCP 20-10-08 Undo redo implementation
-       deleteContours( (std::vector<std::string>)_theViewPanel->getSceneManager()->GetlstContoursNameActualSlice() );
+
+       std::vector<int> tempVector;
+       getInstantVector( tempVector );
+       deleteContours( getOutlinesName( tempVector[1] ) );
+//     deleteContours( (std::vector<std::string>)_theViewPanel->getSceneManager()->GetlstContoursNameActualSlice() );
+
 }
 void wxContourMainFrame::onDeleteAllContours(){
        //JCP 20-10-08 Undo redo implementation
@@ -632,21 +650,27 @@ void wxContourMainFrame::ShowResultImages(int typeContourGroup, int selection, i
        nTypeView.push_back(2);
        nTypeView.push_back(0);
 
-       wxFrame *frame = new wxFrame(this, -1,_T("Result image extracted.."),wxDefaultPosition,wxSize(600,600), wxDEFAULT_FRAME_STYLE|wxSTAY_ON_TOP);
-       wxAuiManager *m_mgr_noteBookViewerResult = new wxAuiManager();
-       m_mgr_noteBookViewerResult->SetManagedWindow(frame);
+       if (_frameShowResultImages==NULL)
+       {
+//             wxFrame *_frameShowResultImages = new wxFrame(this, -1,_T("Result image extracted.."),wxDefaultPosition,wxSize(600,600), wxDEFAULT_FRAME_STYLE|wxSTAY_ON_TOP);
+               _frameShowResultImages = new wxFrame(this, -1,_T("Result image extracted.."),wxDefaultPosition,wxSize(600,600), 
+                    wxFRAME_TOOL_WINDOW|wxSYSTEM_MENU | wxCAPTION |  wxCLIP_CHILDREN |wxFRAME_FLOAT_ON_PARENT | wxRESIZE_BORDER  );
+
+               wxAuiManager *m_mgr_noteBookViewerResult = new wxAuiManager();
+               m_mgr_noteBookViewerResult->SetManagedWindow(_frameShowResultImages);
 
-       wxAuiNotebook* noteBookViewerResult = new wxAuiNotebook(frame, -1, wxPoint(600, 600), wxSize(430,200), m_notebook_style);
-       m_mgr_noteBookViewerResult->AddPane(noteBookViewerResult, wxAuiPaneInfo().Name(wxT("notebook_content")).CenterPane().PaneBorder(false));
+               wxAuiNotebook* noteBookViewerResult = new wxAuiNotebook(_frameShowResultImages, -1, wxPoint(600, 600), wxSize(430,200), m_notebook_style);
+               m_mgr_noteBookViewerResult->AddPane(noteBookViewerResult, wxAuiPaneInfo().Name(wxT("notebook_content")).CenterPane().PaneBorder(false));
 
-       wxPanel *panelViewer1 = new wxMaracas_N_ViewersWidget(noteBookViewerResult, value,&nTypeView);
-       wxPanel *panelViewer2 = new wxMaracas_N_ViewersWidget(noteBookViewerResult, mask,&nTypeView);
+               wxPanel *panelViewer1 = new wxMaracas_N_ViewersWidget(noteBookViewerResult, value,&nTypeView);
+               wxPanel *panelViewer2 = new wxMaracas_N_ViewersWidget(noteBookViewerResult, mask,&nTypeView);
 
 
-       noteBookViewerResult->AddPage( panelViewer1, wxT("Image result") );
-       noteBookViewerResult->AddPage( panelViewer2, wxT("Mask result") );
-       m_mgr_noteBookViewerResult->Update();
-       frame->Show();
+               noteBookViewerResult->AddPage( panelViewer1, wxT("Image result") );
+               noteBookViewerResult->AddPage( panelViewer2, wxT("Mask result") );
+               m_mgr_noteBookViewerResult->Update();
+               _frameShowResultImages->Show();
+       } // if _frameShowResultImages
 }
 
 
@@ -712,13 +736,14 @@ void wxContourMainFrame::updateInstantImageData()
        _instantPanel->getInstant(inst);
 
        vtkImageData* img = kernelManager->getImageAtInstant(inst);
-       if(img!=NULL){
+       if(img!=NULL)
+       {
                _theViewPanel->changeImage(img);
                showAxis(false);
        }else{
                int z = _instantPanel->getConceptDataOf("Axe Depth")->getActualValue();
                _theViewPanel->setImageSlice(z);
-       }
+       } // if img
 }
 
 void wxContourMainFrame::updateInstantAxes()
@@ -1703,19 +1728,18 @@ printf("EED wxContourMainFrame::SegmentationOneSliceITK  Skypping this code.  No
                {
                        vecXo.push_back(p[0]);
                        vecYo.push_back(p[1]);
-               }
-               else
-               {
-               vecX.push_back( p[0] );
-               myfile <<p[0]<<","<<p[1]<<"\n";
-               //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;
-               vecY.push_back( p[1] );
-               vecZ.push_back( 900 );
-               xAct=x;
-               yAct=y;
-               }
+               } else {
+                       vecX.push_back( p[0] );
+                       myfile <<p[0]<<","<<p[1]<<"\n";
+                       //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;
+                       vecY.push_back( p[1] );
+//EED 2017-12-16
+                       vecZ.push_back( -900 );
+                       xAct = x;
+                       yAct = y;
+               } // if xAct
 
 
        }
@@ -1725,7 +1749,8 @@ printf("EED wxContourMainFrame::SegmentationOneSliceITK  Skypping this code.  No
                vecX.push_back(vecXo.back());
                //std::cout<<" x Siguiente "<<vecXo.back();
                vecXo.pop_back();
-               vecZ.push_back( 900 );
+//EED 2017-12-16
+               vecZ.push_back( -900 );
        }
        while(!vecYo.empty())
        {
@@ -1876,18 +1901,15 @@ void wxContourMainFrame::SegmentationOneSlice( int x, int y, int z, int isovalue
 {
        int typeofcontour = 1;
        //--Extracting Contour
-       //vtkImageData  *imagedata      = _theViewPanel->getSceneManager()->GetImageData();
        vtkImageData    *imagedata      = getImageData();
-//             double                  *range          = imagedata->GetScalarRange();
-//             double                  thr                     = 1;
 
        vtkImageReslice *imageReslice = vtkImageReslice::New();
 //EED
-//             double spc[3];
-//             imagedata->GetSpacing(spc);
-//             x = x*spc[0];
-//             y = y*spc[1];
-//             z = z*spc[3];
+               double spc[3];
+               imagedata->GetSpacing(spc);
+               x = x*spc[0];
+               y = y*spc[1];
+               z = z;
 
 //EED 2017-01-01 Migration VTK7
 #if VTK_MAJOR_VERSION <= 5
@@ -1897,7 +1919,7 @@ void wxContourMainFrame::SegmentationOneSlice( int x, int y, int z, int isovalue
 #endif
        imageReslice->SetInformationInput(imagedata);
        imageReslice->SetResliceAxesDirectionCosines(1,0,0, 0,1,0 ,0,0,1);
-       imageReslice->SetResliceAxesOrigin(0,0,z);
+       imageReslice->SetResliceAxesOrigin(0,0,z*spc[2]);
        imageReslice->SetOutputDimensionality(2);
        imageReslice->SetInterpolationModeToLinear();
 
@@ -1945,7 +1967,15 @@ void wxContourMainFrame::SegmentationOneSlice( int x, int y, int z, int isovalue
        cpd->UpdateInformation();
 
        vtkPolyDataConnectivityFilter* conn = vtkPolyDataConnectivityFilter::New( );
-       conn->SetExtractionModeToClosestPointRegion( );
+
+//     conn->SetExtractionModeToPointSeededRegions();
+//     conn->SetExtractionModeToCellSeededRegions();
+//     conn->SetExtractionModeToSpecifiedRegions();
+//     conn->SetExtractionModeToLargestRegion();
+//     conn->SetExtractionModeToAllRegions();
+
+       conn->SetExtractionModeToClosestPointRegion();
+
        //conn->SetMaxRecursionDepth( 3000 );
 
 //EED 2017-01-01 Migration VTK7
@@ -1970,11 +2000,15 @@ void wxContourMainFrame::SegmentationOneSlice( int x, int y, int z, int isovalue
        cpd2->UpdateInformation();
 
        vtkStripper* vtkstripper = vtkStripper::New( );
+vtkstripper->SetJoinContiguousSegments(true);
+       vtkstripper->SetMaximumLength(20000);
 //EED 2017-01-01 Migration VTK7
 #if VTK_MAJOR_VERSION <= 5
        vtkstripper->SetInput( cpd2->GetOutput() );
 #else
        vtkstripper->SetInputData( cpd2->GetOutput() );
+//     vtkstripper->SetInputData( cntVTK->GetOutput() );
+
 #endif
        vtkstripper->Update();
        vtkstripper->UpdateInformation();
@@ -1999,116 +2033,144 @@ myfile << "-------------------------------------\n";
 polyDataResult->GetLines()->Print(myfile);
 myfile.close();
 */
-
-       cntVTK          -> Delete();
-       cpd2            -> Delete();
-       cpd                     -> Delete();
-       conn            -> Delete();
-
+       cntVTK  -> Delete();
+       cpd2    -> Delete();
+       cpd             -> Delete();
+       conn    -> Delete();
 
 //--Calculating control points
-
        std::vector<double> vecX;
        std::vector<double> vecY;
        std::vector<double> vecZ;
-
        std::vector<double> vecCtrlPointX;
        std::vector<double> vecCtrlPointY;
        std::vector<double> vecCtrlPointZ;
 
+printf("EED wxContourMainFrame::SegmentationOneSlice GetNumberOfLines=%d \n",polyDataResult->GetNumberOfLines() );
 
-       double *p;
-       int ii,size=polyDataResult->GetNumberOfPoints();
-       int id;
-       for (ii=1;ii<=size;ii++)
-       {
-               id      = polyDataResult->GetLines()->GetData()->GetValue(ii);
-               p       = polyDataResult->GetPoint(id);
-//             double x=p[0];
-//             double y=p[1];
-               vecX.push_back( p[0] );
-               vecY.push_back( p[1] );
-               vecZ.push_back( 900 );
-//                     vecZ.push_back( p[2] );
-       }
+int numberOfPointsInCell;
+int numberOfcells=polyDataResult->GetNumberOfCells();
+for (int jj=0; jj<numberOfcells ; jj++)
+{
+       vtkCell* vtkcell = polyDataResult->GetCell( jj );
+       numberOfPointsInCell=vtkcell->GetNumberOfPoints();
+       printf("EED wxContourMainFrame::SegmentationOneSlice cell=%d  points=%d \n",jj,numberOfPointsInCell );
+}// for jj
 
+printf("EED wxContourMainFrame::SegmentationOneSlice PrintSelf \n");
 
-       ExtractControlPoints2D *extractcontrolpoints2d = new ExtractControlPoints2D();
-       extractcontrolpoints2d->SetContour( &vecX , &vecY , &vecZ );
+polyDataResult->PrintSelf( std::cout, vtkIndent(0) );
 
-//PROOFS
-       /*if (methodRadiobox->GetSelection()==0){
-               extractcontrolpoints2d->GetInitialControlPoints( &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
-       }
-       if (methodRadiobox->GetSelection()==1){
-               extractcontrolpoints2d->GetControlPoints(  &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
-       }
-       if (methodRadiobox->GetSelection()==2){
-               extractcontrolpoints2d->SetSamplingControlPoints( sampling );
-               extractcontrolpoints2d->GetSamplingControlPoints(  &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
-       }*/
 
-       if (method==0){
-               extractcontrolpoints2d->GetInitialControlPoints( &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
-       }
-       else if (method==1){
-               extractcontrolpoints2d->GetControlPoints(  &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
-       }
-       else if (method==2){
-               extractcontrolpoints2d->SetSamplingControlPoints( sampling );
-               extractcontrolpoints2d->GetSamplingControlPoints(  &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
-       }
+       if (polyDataResult->GetNumberOfCells()>0)
+       {
+               double *p;
+       printf("EED wxContourMainFrame::SegmentationOneSlice A \n");
+               vtkCell* vtkcell = polyDataResult->GetCell( 0 );
+       printf("EED wxContourMainFrame::SegmentationOneSlice B \n");
+       //      int ii,size=polyDataResult->GetNumberOfPoints();
+               int ii,size=vtkcell->GetNumberOfPoints();
+       printf("EED wxContourMainFrame::SegmentationOneSlice C \n");
+               int id;
+               for (ii=1;ii<=size;ii++)
+               {
+       //              id      = polyDataResult->GetLines()->GetData()->GetValue(ii);
+       //              p       = polyDataResult->GetPoint(id);
+                       p       = vtkcell->GetPoints()->GetPoint(ii);
 
-       //--Adding contour to the system
+       printf("EED wxContourMainFrame::SegmentationOneSlice (%d) id=%d  px=%f py=%f\n",ii,id, p[0],p[1]);
+       //              double x=p[0];
+       //              double y=p[1];
+                       vecX.push_back( p[0] );
+                       vecY.push_back( p[1] );
+       //EED 2017-12-16
+                       vecZ.push_back( -900 );
+       //                      vecZ.push_back( p[2] );
+               }
 
-       std::vector<int> actualInstantVector;
-       _instantPanel->getInstant( actualInstantVector );
-       actualInstantVector[1]=z;
+               ExtractControlPoints2D *extractcontrolpoints2d = new ExtractControlPoints2D();
+               extractcontrolpoints2d->SetContour( &vecX , &vecY , &vecZ );
 
-       int j,sizeCtrPt = vecCtrlPointX.size();
-       manualBaseModel *manModelContour =  kernelManager->factoryManualContourModel( typeofcontour );
-       manModelContour->SetNumberOfPointsSpline( ((sizeCtrPt/100)+1)*100 );
-       if (sizeCtrPt>=3){
-               for (j=0 ; j<sizeCtrPt ; j++)
+       //PROOFS
+               /*if (methodRadiobox->GetSelection()==0){
+                       extractcontrolpoints2d->GetInitialControlPoints( &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
+               }
+               if (methodRadiobox->GetSelection()==1){
+                       extractcontrolpoints2d->GetControlPoints(  &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
+               }
+               if (methodRadiobox->GetSelection()==2){
+                       extractcontrolpoints2d->SetSamplingControlPoints( sampling );
+                       extractcontrolpoints2d->GetSamplingControlPoints(  &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
+               }*/
+
+               if (method==0)
                {
-                       manModelContour->AddPoint( vecCtrlPointX[j] , vecCtrlPointY[j] , vecCtrlPointZ[j]  );
-               } // for
-               std::string theName;
-               //theName = _modelManager->createOutline( manModelContour, actualInstantVector );
-               theName = kernelManager->createOutline( manModelContour, actualInstantVector );
-               bool addedModel = theName.compare("") != 0;
-               if( addedModel )
+                       extractcontrolpoints2d->GetInitialControlPoints( &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
+               } else if (method==1){
+                       extractcontrolpoints2d->GetControlPoints(  &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
+               } else if (method==2){
+                       extractcontrolpoints2d->SetSamplingControlPoints( sampling );
+                       extractcontrolpoints2d->GetSamplingControlPoints(  &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
+               }
+
+               //--Adding contour to the system
+               std::vector<int> actualInstantVector;
+               _instantPanel->getInstant( actualInstantVector );
+               actualInstantVector[1]=z;
+
+               int j,sizeCtrPt = vecCtrlPointX.size();
+               manualBaseModel *manModelContour =  kernelManager->factoryManualContourModel( typeofcontour );
+               manModelContour->SetNumberOfPointsSpline( ((sizeCtrPt/100)+1)*100 *2 );
+               manModelContour->SetLabel2( std::to_string(z) );
+               if (sizeCtrPt>=3)
                {
-                       double spc[3];//Si no hay imagen pero hay contornos que spacing se pone por default
-                       _theViewPanel->getSpacing(spc);
-                       //Adding the manualContourControler to interface objects structure
-                       //Adding the manualViewContour to interface objects structure
-                       //_theViewPanel->getSceneManager()->setControlActiveStateOfALL( false );//This call is being done here because if the ROI is created underneath the previously created ROIS will still be active.
-                       _theViewPanel->configureViewControlTo(theName, manModelContour, spc, typeofcontour);
-                       //_theViewPanel->getSceneManager()->configureViewControlTo( theName, manModelContour,spc, typeofcontour ) ;
-               }       // if addedModel
-       } // if sizeCtrPt
+                       for (j=0 ; j<sizeCtrPt ; j++)
+                       {
+                               manModelContour->AddPoint( vecCtrlPointX[j] / spc[0], vecCtrlPointY[j] / spc[1], vecCtrlPointZ[j] );
+                       } // for
+                       std::string theName;
+                       //theName = _modelManager->createOutline( manModelContour, actualInstantVector );
+                       theName = kernelManager->createOutline( manModelContour, actualInstantVector );
+                       bool addedModel = theName.compare("") != 0;
+                       if( addedModel )
+                       {
+                               double spc[3];//Si no hay imagen pero hay contornos que spacing se pone por default
+                               _theViewPanel->getSpacing(spc);
+                               //Adding the manualContourControler to interface objects structure
+                               //Adding the manualViewContour to interface objects structure
+                               //_theViewPanel->getSceneManager()->setControlActiveStateOfALL( false );//This call is being done here because if the ROI is created underneath the previously created ROIS will still be active.
+                               _theViewPanel->configureViewControlTo(theName, manModelContour, spc, typeofcontour);
+                               //_theViewPanel->getSceneManager()->configureViewControlTo( theName, manModelContour,spc, typeofcontour ) ;
+                       }       // if addedModel
+               } // if sizeCtrPt
+       } // if NumberOfCells > 0
+
+printf("EED wxContourMainFrame::SegmentationOneSlice  End\n");
 }
-int wxContourMainFrame::GetImageDataSizeZ(){
+
+int wxContourMainFrame::GetImageDataSizeZ()
+{
        return _theViewPanel->GetImageDataSizeZ();
 }
 
-void wxContourMainFrame::GetImageDataRange(double *range){
+void wxContourMainFrame::GetImageDataRange(double *range)
+{
        _theViewPanel->GetImageDataRange(range);
 }
 
-void wxContourMainFrame::onSegmentationAllSlice(int minZ,int maxZ,int isovalue,int sampling,int method){
+void wxContourMainFrame::onSegmentationAllSlice(int minZ,int maxZ,int isovalue,int sampling,int method)
+{
        //JCP 20-10-08 Undo redo implementation
        saveState();
        //JCP 20-10-08 Undo redo implementation
 
        wxBusyCursor wait;
-       int                                     x                                       = _theViewPanel->GetX();
-       int                                     y                                       = _theViewPanel->GetY();
-       int z;
-       double porcent;
-       wxString tmpString;
-       double totalZ = maxZ-minZ+1;
+       int                                     x = _theViewPanel->GetX();
+       int                                     y = _theViewPanel->GetY();
+       int                             z;
+       double                          porcent;
+       wxString                        tmpString;
+       double                          totalZ = maxZ-minZ+1;
        for( z=minZ ; z<=maxZ ; z++ )
        {
                porcent = 100.0* (z-minZ)/totalZ;
@@ -2347,18 +2409,44 @@ void wxContourMainFrame::onInterpolation(bool interpolate)
        RefreshInterface();
 }
 
+void wxContourMainFrame::onBeforeAfterContour(bool before,bool after)
+{
+       _theViewPanel->onBeforeAfterContour(before,after);
+       RefreshInterface();
+}
+
+void wxContourMainFrame::onJumpSlice(int step, int shift)
+{
+       std::vector<int> tempVector;
+       _instantPanel->getInstant( tempVector );
+       int ss  = abs(step);
+       int pos         = tempVector[1];
+       int maxZ        = interfMainPanel::getInstance()->GetImageDataSizeZ();
+       pos=(pos/ss)*ss+step+shift;
+       if ((pos>=0) && (pos<maxZ)) 
+       {
+               tempVector[1]=pos;
+               _instantPanel->setInstant( tempVector );
+               RefreshInterface();
+       } // if pos     
+}
+
+void wxContourMainFrame::onShowTextContour(bool show)
+{
+       _theViewPanel->onShowTextContour(show);
+       RefreshInterface();
+}
+
 void wxContourMainFrame::onChangeInstant(std::string name,int actual)
 {
        _instantPanel->setConceptValue( name, actual );
 }
 
-
 void wxContourMainFrame::resetAppend()
 {
        kernelManager->resetAppend();
 }
 
-
 void wxContourMainFrame::onSpreadAdd()
 {
        std::vector<double> vecX;
@@ -2469,6 +2557,7 @@ void wxContourMainFrame::onSpreadGo(int type)
 
 void wxContourMainFrame::getInstantVector(std::vector<int>& tempVector)
 {
+       _instantPanel->getInstant( tempVector );
 }
 
 std::vector<std::string> wxContourMainFrame::getOutlinesName(int slide)
@@ -2945,7 +3034,6 @@ void wxContourMainFrame::getMaskValue(vtkImageData ** mask,vtkImageData ** value
        int z;
        int ii,sizeLstContourThings;
 
-
        double totalZ;
        double porcent;
 
@@ -3026,13 +3114,41 @@ void wxContourMainFrame::getMaskValue(vtkImageData ** mask,vtkImageData ** value
 
 }
 
-void wxContourMainFrame::SaveImageResult(std::string directory,std::string namefile, vtkImageData* mask,vtkImageData* value)
+std::string wxContourMainFrame::GetStringTypeImage( vtkImageData* image)
 {
+       std::string result="";
+       if ( image->GetScalarType()==VTK_CHAR )
+       {
+               result="CHAR";
+       } else if ( image->GetScalarType()==VTK_UNSIGNED_CHAR ){
+               result="UCHAR";
+       } else if ( image->GetScalarType()==VTK_FLOAT ){
+               result="FLOAT";
+       } else if ( image->GetScalarType()==VTK_DOUBLE ){
+               result="DOUBLE";
+       } else if ( image->GetScalarType()==VTK_INT ){
+               result="INT";
+       } else if ( image->GetScalarType()==VTK_UNSIGNED_INT ){
+               result="UINT";
+       } else if ( image->GetScalarType()==VTK_SHORT ){
+               result="SHORT";
+       } else if ( image->GetScalarType()==VTK_UNSIGNED_SHORT ){
+               result="USHORT";
+       } else if ( image->GetScalarType()==VTK_LONG ){
+               result="LONG";
+       } else if ( image->GetScalarType()==VTK_UNSIGNED_LONG ){
+               result="ULONG";
+       }
+       return result;
+}
 
+void wxContourMainFrame::SaveImageResult(std::string directory,std::string namefile, vtkImageData* mask,vtkImageData* value)
+{
+       std::string stringType;
        double spc[3];
        kernelManager->GetSpacing(spc,0);   // returns the spacing of the first image
-       
-       std::string temp = directory + "/" + namefile + "-image-Value.mhd";
+       stringType              = GetStringTypeImage(value);
+       std::string temp = directory + "/" + namefile + "-image-Value-"+stringType+".mhd";
 
        wxString filename = crea::std2wx(temp);
 
@@ -3046,6 +3162,7 @@ void wxContourMainFrame::SaveImageResult(std::string directory,std::string namef
        valueChInfo->SetInputData(value);
 #endif
        valueChInfo->SetSpacingScale(spc);      
+       valueChInfo->Update();  
        vtkMetaImageWriter *writerValueImage = vtkMetaImageWriter::New( );
 
 //EED 2017-01-01 Migration VTK7
@@ -3060,7 +3177,8 @@ void wxContourMainFrame::SaveImageResult(std::string directory,std::string namef
        writerValueImage->Write( );
 
 // Image Mask
-       temp = directory + "/" + namefile + "-image-Mask.mhd";
+       stringType              = GetStringTypeImage(mask);
+       temp = directory + "/" + namefile + "-image-Mask-"+stringType+".mhd";
 
        // LG 14/01/09 : using crea
        filename = crea::std2wx(temp);
@@ -3075,6 +3193,7 @@ void wxContourMainFrame::SaveImageResult(std::string directory,std::string namef
        maskChInfo->SetInputData(mask);
 #endif
        maskChInfo->SetSpacingScale(spc);       
+       maskChInfo->Update();   
        vtkMetaImageWriter *writerMaskImage = vtkMetaImageWriter::New( );
 
 //EED 2017-01-01 Migration VTK7
@@ -3153,7 +3272,7 @@ void wxContourMainFrame::SetContourGroup(int contourGroup)
 
 wxPanel *wxContourMainFrame::getMaskImageViewPanel(wxWindow *parent)
 {
-
        if      (_viewMaskImagePanel==NULL)
        {
                _contourextractdata = new ContourExtractData(true);
@@ -3162,15 +3281,29 @@ wxPanel *wxContourMainFrame::getMaskImageViewPanel(wxWindow *parent)
 //             _theViewPanel->getImageData()->GetScalarRange(range);
                //int min = (int)floor (range[0]);
                //int max = (int)ceil (range[1]);
+
+//EED 2017-12-17
                int min = 254;
                int max = 256;
-               _viewMaskImagePanel = new ThresholdImageViewPanel( parent, -1, -1 , 0);
+
+               _viewMaskImagePanel = new ThresholdImageViewPanel( parent, -1, -1 , -1,-1,  0);
                _viewMaskImage = new ThresholdImageView();
-               _viewMaskImage->SetImage(_contourextractdata->GetVtkImageMaskResult() );
+               _viewMaskImage->SetImage( _contourextractdata->GetVtkImageMaskResult() );
                _viewMaskImage->SetminMaxValue( min, max);
                _viewMaskImage->SetBaseColor( 0,0,1 );
+
+
+               std::vector<int> colorbarposition;
+               colorbarposition.push_back( -999 );
+               colorbarposition.push_back( -999 );
+               _viewMaskImage->SetColorBarPosition( colorbarposition );
+               double spc[3];
+               _viewMaskImage->SetSpcBase( _theViewPanel->getImageData()->GetSpacing() );
+               _viewMaskImage->SetNewSpacingLayer( _theViewPanel->getImageData()->GetSpacing() );
+
                _viewMaskImagePanel->SetThresholdImageView( _viewMaskImage );
                _viewMaskImage->SetwxVtkBaseView( _theViewPanel->getWxVtkBaseView() );
+               _viewMaskImage->onThresholdChangeOpacity( 25 );
        }
        return _viewMaskImagePanel;
 }
@@ -3186,6 +3319,21 @@ wxPanel *wxContourMainFrame::getColorLayerImageViewPanel(wxWindow *parent)
                int max = (int)ceil (range[1]);
                _viewColorLayerImagePanel = new ColorLayerImageViewPanel( parent, min, max, 100, 1);
                _viewColorLayerImagePanel->GetColorLayerImageViewManager()->SetwxVtkBaseView( 0 , _theViewPanel->getWxVtkBaseView() );
+       //    std::vector<double> range = bbGetInputlstRangeForColorBar();
+       //    clivp->GetColorLayerImageViewManager()->SetRangeForColorBar( range );
+               std::vector<int> colorbarposition;
+               colorbarposition.push_back( 50 );
+               colorbarposition.push_back( 105 );
+               _viewColorLayerImagePanel->GetColorLayerImageViewManager()->SetColorBarPosition( colorbarposition );
+       //    std::vector<double> base_color;
+       //    _viewColorLayerImagePanel->GetColorLayerImageViewManager()->SetBaseColors( base_color );
+       //    std::vector<double> grey_level_boundaries;
+       //    clivp->GetColorLayerImageViewManager()->SetGreyLevelBoundaries( grey_level_boundaries );
+               std::vector<double> transparence_level_boundaries;
+               transparence_level_boundaries.push_back(0);
+               _viewColorLayerImagePanel->GetColorLayerImageViewManager()->SetBaseTransparence( transparence_level_boundaries );
+               _viewColorLayerImagePanel->GetColorLayerImageViewManager()->SetPlainOrGradientColor( false );
+//             _viewColorLayerImagePanel->SetFittingMode( 3 );
        }
        return _viewColorLayerImagePanel;
 }
@@ -3199,13 +3347,25 @@ wxPanel *wxContourMainFrame::getThresholdImageViewPanel(wxWindow *parent)
                _theViewPanel->getImageData()->GetScalarRange(range);
                int min = (int)floor (range[0]);
                int max = (int)ceil (range[1]);
-               _viewThresholdImagePanel = new ThresholdImageViewPanel( parent, min, max , 1);
+               int start= (max+min)/2 - (max-min)*0.25; 
+               int end= (max+min)/2 + (max-min)*0.25; 
+               _viewThresholdImagePanel = new ThresholdImageViewPanel( parent, min, max, start, end , 1);
                _viewThresholdImage = new ThresholdImageView();
+               _viewThresholdImage->SetminMaxValue(start,end);
                _viewThresholdImage->SetImage( _theViewPanel->getImageData() );
-               _viewThresholdImage->SetminMaxValue( min, max);
+               std::vector<int> colorbarposition;
+               colorbarposition.push_back( 10 );
+               colorbarposition.push_back( 105 );
+               _viewThresholdImage->SetColorBarPosition( colorbarposition );
+               double spc[3];
+               _viewThresholdImage->SetSpcBase( _theViewPanel->getImageData()->GetSpacing() );
+               _viewThresholdImage->SetNewSpacingLayer( _theViewPanel->getImageData()->GetSpacing() );
+
                _viewThresholdImage->SetBaseColor( 1,0,0 );
                _viewThresholdImagePanel->SetThresholdImageView( _viewThresholdImage );
                _viewThresholdImage->SetwxVtkBaseView( _theViewPanel->getWxVtkBaseView() );
+               _viewThresholdImage->onThresholdChangeOpacity( 25 );
+               RefreshInterface();
        }
        return _viewThresholdImagePanel;
 }
@@ -3315,7 +3475,6 @@ void wxContourMainFrame::SetFileLocation(std::string newFileLocation)
 
 void wxContourMainFrame::FlipSelectedContours(int flipDirection)
 {
-       printf("EED wxContourMainFrame::FlipSelectedContours  flipDirection=%d  \n",  flipDirection ); 
        FlipLstOfContours( 
                            (std::vector<std::string>)_theViewPanel->getSceneManager()->getSelectedObjects() , 
                            flipDirection
@@ -3325,7 +3484,6 @@ void wxContourMainFrame::FlipSelectedContours(int flipDirection)
 
 void wxContourMainFrame::FlipActualSliceContours(int flipDirection)
 {
-       printf("EED wxContourMainFrame::FlipActualSliceContours  flipDirection=%d  \n",  flipDirection );       
        FlipLstOfContours( 
                            (std::vector<std::string>)_theViewPanel->getSceneManager()->GetlstContoursNameActualSlice() , 
                            flipDirection
@@ -3334,7 +3492,6 @@ void wxContourMainFrame::FlipActualSliceContours(int flipDirection)
 
 void wxContourMainFrame::FlipAllContours(int flipDirection)
 {
-       printf("EED wxContourMainFrame::FlipAllContours  flipDirection=%d  \n",  flipDirection );
        FlipLstOfContours( 
                            kernelManager->GetLstNameThings() , 
                            flipDirection
@@ -3405,6 +3562,11 @@ void wxContourMainFrame::FlipContours(int typeContourFlip, int flipDirection)
 }
 
 
+void wxContourMainFrame::SetXY(int x, int y)
+{
+       _theViewPanel->SetXY(x,y);
+}
+