]> Creatis software - creaContours.git/blobdiff - lib/Interface_ManagerContour_NDimensions/wxContourMainFrame.cxx
GetViewer Box
[creaContours.git] / lib / Interface_ManagerContour_NDimensions / wxContourMainFrame.cxx
index 2c37a7f5ec068f3181376cc97f5c5a4da3cbcab7..37cc8d33f88d2aa9c4fc62b5071fb2bbab335de2 100644 (file)
@@ -52,7 +52,7 @@ using namespace std;
        // Constructors & Destructors
        //------------------------------------------------------------------------------------------------------------
 
-wxContourMainFrame* wxContourMainFrame:: instance = NULL;
+wxContourMainFrame* wxContourMainFrame::instance=NULL;
 char wxContourMainFrame::COPY = 'C';
 
        wxContourMainFrame :: wxContourMainFrame(wxWindow* parent, wxWindowID id,const wxString& title,const wxPoint& pos,const wxSize& size,long style)
@@ -60,7 +60,6 @@ char wxContourMainFrame::COPY = 'C';
                //: wxWindow(parent, id, pos, size, style)
                : wxPanel(parent, id, pos, size, style)
        {
-printf("EED wxContourMainFrame :: wxContourMainFrame Start 1\n");
                m_mgr.SetManagedWindow(this);
                _contourGroup                           = 1;
                _creatingContoursActive         = false;
@@ -95,18 +94,16 @@ printf("EED wxContourMainFrame :: wxContourMainFrame Start 1\n");
                _viewThresholdImage                     = NULL;
                _viewThresholdImagePanel        = NULL;
                _viewColorLayerImagePanel       = NULL;
-printf("EED wxContourMainFrame :: wxContourMainFrame End 1\n");
+               _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)
+       wxContourMainFrame::wxContourMainFrame(wxWindow* parent, wxWindowID id,const wxString& title,const wxPoint& pos,const wxSize& size,std::vector<vtkImageData*> images, long style,std::string datadir)
                : wxPanel(parent, id, pos, size, style)
        {
-printf("EED wxContourMainFrame :: wxContourMainFrame Start 2\n");
                m_mgr.SetManagedWindow(this);
                _contourGroup                           = 1;
                _creatingContoursActive         = false;
                _theViewPanel                           = NULL;
-
                _instantPanel                           = NULL;
                _refLineControl                         = NULL;
                _refLineModel                           = NULL;
@@ -117,7 +114,6 @@ printf("EED wxContourMainFrame :: wxContourMainFrame Start 2\n");
                _viewThresholdImage                     = NULL;
                _viewThresholdImagePanel        = NULL;
                _viewColorLayerImagePanel       = NULL;
-printf("EED wxContourMainFrame :: wxContourMainFrame A\n");
                // set up default notebook style
                m_notebook_style                        = wxAUI_NB_TAB_SPLIT | wxAUI_NB_TAB_EXTERNAL_MOVE | wxNO_BORDER;
                m_notebook_theme                        = 0;
@@ -126,47 +122,33 @@ printf("EED wxContourMainFrame :: wxContourMainFrame A\n");
                _datadir                                        = datadir;
                inredo                                          = 0;
                inundo                                          = 0;
-printf("EED wxContourMainFrame :: wxContourMainFrame B\n");
                _pannew = interfMainPanel::getInstance(this,datadir+"/data/Icons");//, eventHandler);
-printf("EED wxContourMainFrame :: wxContourMainFrame C\n");
-               _performingOperation            = new PerformingOperation();
-printf("EED wxContourMainFrame :: wxContourMainFrame C.1\n");
+               _performingOperation        = new PerformingOperation();
                if(images.size() > 0)
                {
-printf("EED wxContourMainFrame :: wxContourMainFrame C.2\n");
                        this->setVectImages(images);
-printf("EED wxContourMainFrame :: wxContourMainFrame C.3\n");
                }
                _contourextractdata                     = NULL;
-printf("EED wxContourMainFrame :: wxContourMainFrame End 2\n");
 }
 
 void wxContourMainFrame::setVectImages(std::vector<vtkImageData*> imgs)
 {
-printf("EED wxContourMainFrame::setVectImages Start\n");
        _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;
-printf("EED wxContourMainFrame::setVectImages Start  01\n");
-       notebook = this->createNotebook();
-printf("EED wxContourMainFrame::setVectImages Start  02\n");
+       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 );
-printf("EED wxContourMainFrame::setVectImages Start  03\n");
        _instantPanel   = new wxInstantChooserPanel( notebook, "Instant Chooser", true, false, "c" );
-printf("EED wxContourMainFrame::setVectImages Start  04\n");
        _theViewPanel   = new wxContourViewPanel( kernelManager->getSourceImage(), notebook );
-printf("EED wxContourMainFrame::setVectImages Start   05\n");
 //-------------------------------------------------------------
        kernelManager->getConceptsInformation(conceptNameVect, conceptSizeVect);
-printf("EED wxContourMainFrame::setVectImages 1\n");
        _instantPanel->addConcepts(conceptNameVect, conceptSizeVect);
-printf("EED wxContourMainFrame::setVectImages 2\n");
        double val = _theViewPanel->getCurrentDeep();
        _instantPanel->setConceptValue( "Axe Depth", (int)val );
        //JCP _instantPanel->setConceptValue("time", num de imagenes);
@@ -174,34 +156,34 @@ printf("EED wxContourMainFrame::setVectImages 2\n");
        //*******************changeInstant();
        _theViewPanel->setVerticalConcept( "Axe Depth", data->getMinValue(), data->getMaxValue(), data->getMinShowedValue(),  data->getMaxShowedValue(), data->getActualValue() );
        this->configurePanels( );
-printf("EED wxContourMainFrame::setVectImages End\n");
 }
 
 
        wxContourMainFrame* wxContourMainFrame :: getInstance(wxWindow* parent, wxWindowID id,const wxString& title,const wxPoint& pos,const wxSize& size,std::vector<vtkImageData*> images, long style,std::string datadir )
        {
-               if(instance == NULL){
-                       instance = new wxContourMainFrame(parent, id, title, pos, size, images, style,datadir);
+               if(instance == NULL)
+               {
+                       instance=new wxContourMainFrame(parent, id, title, pos, size, images, style,datadir);
                }
                return instance;
        }
 
-       std::vector<vtkImageData*> wxContourMainFrame ::getVectImages()
+       std::vector<vtkImageData*> wxContourMainFrame::getVectImages()
        {
                return kernelManager->getVectImages();
        }
 
-       wxContourMainFrame* wxContourMainFrame :: getInstance()
+       wxContourMainFrame* wxContourMainFrame::getInstance()
        {
                return instance;
        }
 
-       void wxContourMainFrame :: resetInstance()
+       void wxContourMainFrame::resetInstance()
        {
                instance->Destroy();
        }
 
-       wxContourMainFrame :: ~wxContourMainFrame()
+       wxContourMainFrame::~wxContourMainFrame()
        {
                delete _contourextractdata;
                deleteAllContours();
@@ -230,7 +212,7 @@ printf("EED wxContourMainFrame::setVectImages End\n");
        // Creational and initialization methods using WxAui
        //------------------------------------------------------------------------------------------------------------
        //
-       wxAuiNotebook * wxContourMainFrame :: createNotebook()
+       wxAuiNotebook * wxContourMainFrame::createNotebook()
        {
                wxSize client_size = GetClientSize();
                wxAuiNotebook* noteBook = new wxAuiNotebook(this, -1, wxPoint(client_size.x, client_size.y), wxSize(430,200), m_notebook_style);
@@ -242,7 +224,7 @@ printf("EED wxContourMainFrame::setVectImages End\n");
        //------------------------------------------------------------------------------------------------------------
        // Creational and initialization methods
        //------------------------------------------------------------------------------------------------------------
-       bool wxContourMainFrame :: configurePanels( )
+       bool wxContourMainFrame::configurePanels( )
        {
                bool configured = _theViewPanel!=NULL;
 
@@ -252,17 +234,13 @@ printf("EED wxContourMainFrame::setVectImages End\n");
                        notebook->AddPage( _theViewPanel, wxT("       View       ") );
                        m_mgr.Update();
                }
-
-
                if( configured )
                {
                        notebook->AddPage( _instantPanel, wxT("Instant Page") );
                        m_mgr.Update();
                }
-
                m_mgr.AddPane(notebook, wxAuiPaneInfo().Name(wxT("notebook_content")).CenterPane().PaneBorder(false));
                m_mgr.Update();
-
                SetMinSize(wxSize(300,300));
                m_mgr.Update();
                return configured;
@@ -343,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
@@ -364,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
@@ -373,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
@@ -408,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;
@@ -430,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();
@@ -465,13 +454,15 @@ manualBaseModel * wxContourMainFrame::factoryManualContourModel(wxPanel* panel){
 
        return manModelContour;
 }
-void wxContourMainFrame :: saveState(){
 
+void wxContourMainFrame :: saveState()
+{
        std::string temp = kernelManager->saveState();
        saveFileWithContours(temp);
 }
 
-void wxContourMainFrame::onDeleteContour(){
+void wxContourMainFrame::onDeleteContour()
+{
        //JCP 20-10-08 Undo redo implementation
                saveState();
                //JCP 20-10-08 Undo redo implementation
@@ -532,7 +523,8 @@ void wxContourMainFrame :: deleteContours( std::vector<std::string>  keyNamesVec
        }
 }
 
-void wxContourMainFrame :: deleteContour( std::string theKeyName ){
+void wxContourMainFrame :: deleteContour( std::string theKeyName )
+{
        /*
        manualContourModel                      * cModel;
        manualViewBaseContour           * cViewer;
@@ -561,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
@@ -621,17 +619,20 @@ void wxContourMainFrame::deleteAllContours(){
 JCP --08-09-2008 */
 }
 
-void wxContourMainFrame::setConceptValue( std::string name, int value ){
+void wxContourMainFrame::setConceptValue( std::string name, int value )
+{
        _instantPanel->setConceptValue(name, value);
 }
 
 
-ConceptDataWrap* wxContourMainFrame::getLastConceptData(){
+ConceptDataWrap* wxContourMainFrame::getLastConceptData()
+{
        return _instantPanel->getLastConceptData();
 }
 
 
-bool wxContourMainFrame::getIfConceptCheckedAt( std::string name, int pos ){
+bool wxContourMainFrame::getIfConceptCheckedAt( std::string name, int pos )
+{
        return _instantPanel->getIfConceptCheckedAt( name, pos );
 }
 
@@ -649,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
 }
 
 
@@ -674,11 +681,10 @@ void wxContourMainFrame::changeInstant()
                std::vector<int> instantVect;
                _instantPanel->getInstant( instantVect );
                int actualSlice = instantVect[1];
-
                kernelManager->setInstant(instantVect);
-
            // Refresh Mask image
-               if(_viewMaskImage!=NULL){
+               if(_viewMaskImage!=NULL)
+               {
                        _viewMaskImage->SetZ(actualSlice);
                        if (_viewMaskImagePanel->IsVisible()==true)
                        {
@@ -686,17 +692,16 @@ void wxContourMainFrame::changeInstant()
                                getMaskValue(&mask,&value, _contourGroup , 0, -1, -1);
                                _viewMaskImage->onThreshold();
                        }
-               }
-
+               } // if _viewMaskImage
        // Refresh Threshold image
-               if(_viewThresholdImage!=NULL){
+               if(_viewThresholdImage!=NULL)
+               {
                        _viewThresholdImage->SetZ(actualSlice);
-                       if (_viewThresholdImagePanel->IsVisible()==true){
+                       if (_viewThresholdImagePanel->IsVisible()==true)
+                       {
                                _viewThresholdImage->onThreshold();
                        }
-               }
-
-
+               } // if _viewThresholdImage
                updateInstantOutlines();
                updateInstantImageData();
                updateInstantAxes();
@@ -704,7 +709,7 @@ void wxContourMainFrame::changeInstant()
 
 
 
-void wxContourMainFrame :: updateInstantOutlines()
+void wxContourMainFrame::updateInstantOutlines()
 {
        _theViewPanel->removeSceneContours();
        _theViewPanel->addNameWrapperToScene();
@@ -715,10 +720,13 @@ void wxContourMainFrame :: updateInstantOutlines()
        //}
 }
 
-int wxContourMainFrame::getNamesWrappingSize(){
+int wxContourMainFrame::getNamesWrappingSize()
+{
        return kernelManager->getNamesWrappingSize();
 }
-std::string wxContourMainFrame::getNameWrapping(int i){
+
+std::string wxContourMainFrame::getNameWrapping(int i)
+{
        return kernelManager->getNameWrapping(i);
 }
 
@@ -728,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()
@@ -750,7 +759,6 @@ void wxContourMainFrame::onChangeDeep(int val){
 
 void wxContourMainFrame::onCopy(){
        std::vector<std::string> currentSelection = _theViewPanel->getSelectedObjects();
-
        std::vector<int> tempVector;
        _instantPanel->getInstant( tempVector );
     _performingOperation->reset();
@@ -759,9 +767,8 @@ void wxContourMainFrame::onCopy(){
        _performingOperation->setKeyNamesOperationElems( currentSelection );
 }
 
-void wxContourMainFrame::onPaste(){
-
-
+void wxContourMainFrame::onPaste()
+{
        char theStartCommand = _performingOperation->getStartCommand();
        if (  theStartCommand == COPY )
        {
@@ -779,7 +786,9 @@ void wxContourMainFrame::onPaste(){
                }
        }
 }
-void wxContourMainFrame::onUndo(){
+
+void wxContourMainFrame::onUndo()
+{
        std::string filename;
        if(kernelManager->onUndoSaveFile(filename)){
                saveFileWithContours(filename);
@@ -843,116 +852,17 @@ void wxContourMainFrame :: createMirrorContourOf ( std::string anExistingKName,
                mp->SetPointX( vecOut[0] );
                mp->SetPointY( vecOut[1] );
                mp->SetPointZ( vecOut[2] );
-       }
+       } // for i
 
        _theViewPanel->getSceneManager()->createCopyContourOf( anExistingKName, cloneName, manualModel , append );
 }
 
-//EED04
-void wxContourMainFrame ::loadState(std::string filename){
-
+void wxContourMainFrame ::loadState(std::string filename)
+{
        deleteAllContours();
        onLoadContours(filename,false);
-
-       /* EED Borrame
-       char tmp[255];
-       FILE *pFile=fopen(filename.c_str(),"r+");
-
-       fscanf(pFile,"%s",tmp); // --CreaContour--
-
-       fscanf(pFile,"%s",tmp); // Version
-       fscanf(pFile,"%s",tmp); // 1.0.3 || 1.0.2 || 1.0.1 || 1.0.0
-
-       deleteAllContours();
-
-       loadContours(pFile, false);
-       loadContours(pFile, true);
-
-       fclose(pFile);
-       */
 }
 
-
-//EED03
-/*EED Borrame
-void wxContourMainFrame::loadContours( FILE *pFile, bool staticContour )
-{
-       char tmp[255];
-
-       if (staticContour==false)
-       {
-               fscanf(pFile,"%s",tmp); // ImageDimensions
-               fscanf(pFile,"%s",tmp); // X
-               fscanf(pFile,"%s",tmp); // Y
-               fscanf(pFile,"%s",tmp); // Z
-
-               fscanf(pFile,"%s",tmp); // ImageSpacing
-               fscanf(pFile,"%s",tmp); // X
-               fscanf(pFile,"%s",tmp); // Y
-               fscanf(pFile,"%s",tmp); // Z
-       }
-
-       fscanf(pFile,"%s",tmp); // NumberOfContours
-       fscanf(pFile,"%s",tmp); // ##
-       int numberOfContours = atoi(tmp);
-
-       std::vector<int> instantVector;
-       int typeContourModel;
-       manualBaseModel *manModelContour;
-       int typeView;
-
-       int i;
-       for (i=0;i<numberOfContours;i++)
-       {
-               instantVector.clear();
-               fscanf(pFile,"%s",tmp); // Instant
-               fscanf(pFile,"%s",tmp); // 1
-               instantVector.push_back( atoi(tmp) );
-               fscanf(pFile,"%s",tmp); // 2
-               instantVector.push_back( atoi(tmp) );
-               fscanf(pFile,"%s",tmp); // 3
-               instantVector.push_back( atoi(tmp) );
-               fscanf(pFile,"%s",tmp); // 4
-               instantVector.push_back( atoi(tmp) );
-               fscanf(pFile,"%s",tmp); // 5
-               instantVector.push_back( atoi(tmp) );
-               fscanf(pFile,"%s",tmp); // 6
-               instantVector.push_back( atoi(tmp) );
-
-               fscanf(pFile,"%s",tmp); // TypeContourModel
-               fscanf(pFile,"%s",tmp); // ##
-               typeContourModel = atoi(tmp);
-
-               manModelContour =  kernelManager->factoryManualContourModel(typeContourModel);
-               manModelContour->Open(pFile);
-
-               fscanf(pFile,"%s",tmp); // TypeView
-               fscanf(pFile,"%s",tmp); // ##
-               typeView = atoi(tmp);
-
-               std::string theName;
-               theName = kernelManager->createOutline( manModelContour, instantVector );
-               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->getSceneManager()->configureViewControlTo( theName, manModelContour,spc , typeView) ;
-               }
-
-               if (staticContour==true)
-               {
-                       Instant instant(&instantVector);
-                       kernelManager->changeContourOfManager( theName , &instant );
-               }
-       }// for  numberOfContours
-}
-*/
-
-
 void wxContourMainFrame::onLoad()
 {
        std::string fileNameContourROI = GetFileLocation();
@@ -969,23 +879,17 @@ void wxContourMainFrame::onLoad()
                if (dialog.ShowModal() == wxID_OK)
                {
                        fileNameContourROI = (const char *)(dialog.GetPath().mb_str());
-               }
-       }
-
+               } // if dialog
+       } // if FileLocation
        onLoadContours(fileNameContourROI,true);
 }
 
-
 void  wxContourMainFrame::SetZForAllContours(int pz)
 {
-       printf("EED wxContourMainFrame::SetZForAllContours\n");
-
        //-- Normal Contours
-
        std::vector< std::string > lstNameThings;       
        int i,sizeLstNameThings;
-       int ii,sizeLstPoints;
-       
+       int ii,sizeLstPoints;   
        lstNameThings           = kernelManager->GetLstNameThings();
        sizeLstNameThings       = lstNameThings.size(); 
        for (i=0 ; i<sizeLstNameThings ; i++)
@@ -1037,17 +941,10 @@ void wxContourMainFrame::onLoadContours(std::string fileNameContourROI, bool int
                        fscanf(pFileData,"%s",tmpD); // NumberOfContours
                        fscanf(pFileData,"%s",tmpD); // #
                }
-
                fscanf(pFile,"%s",tmp); // --CreaContour--
-
                fscanf(pFile,"%s",tmp); // Version
                fscanf(pFile,"%s",tmp); // 1.0.3 || 1.0.2 || 1.0.1 || 1.0.0
                std::string version(tmp);
-
-
-printf("EED wxContourMainFrame::onLoadContours  version%s \n", version.c_str() );
-
-
                //AD:02-06-09
         _tmpReadFileTypeOfTransformation=-1;
                if (version=="1.0.3")
@@ -1056,7 +953,6 @@ printf("EED wxContourMainFrame::onLoadContours  version%s \n", version.c_str() )
                        openContours(pFile,pFileData,false);
                        openContours(pFile,pFileData,true);     //Load StaticContours
                }       
-
                if (version=="1.0.2")
                {
 //EED001
@@ -1065,28 +961,18 @@ printf("EED wxContourMainFrame::onLoadContours  version%s \n", version.c_str() )
                        readDimSpc(pFile,false);                // DimSpc interactive OFF
                        openContours(pFile,pFileData,true);     // Load StaticContours
                        SetZForAllContours(-900);
-               }
-
-               //AD:02-06-09
-               else if (version=="1.0.1")
-               {
+               } else if (version=="1.0.1") { //AD:02-06-09
                        openContours(pFile,pFileData,false); //Skips ImageDimensions and ImageSpacing
                        openContours(pFile,pFileData,true); //Load StaticContours, skips ImageDimensions and ImageSpacing
-               }
-
-               //AD:02-06-09
-               else if (version=="1.0.0")
-               {
+               } else if (version=="1.0.0") { //AD:02-06-09
                        openContours(pFile,pFileData,false); //Skips ImageDimensions and ImageSpacing
                }
-
                fclose(pFile);
-
                if(pFileData != NULL)
                {
                        fclose(pFileData);
                } // if
-       }
+       } // pFile
        _theViewPanel->getSceneManager()->removeSceneContours();
        changeInstant();
 }
@@ -1106,7 +992,8 @@ void wxContourMainFrame::onImport()
                filename = dialog.GetPath().mb_str();
                filenamecontours = kernelManager->parseOsirixFile(filename);
 
-               if(filenamecontours.compare("") != 0){
+               if(filenamecontours.compare("") != 0)
+               {
                        onLoadContours(filenamecontours,true);
                } // if
        } // if
@@ -1121,7 +1008,6 @@ void wxContourMainFrame::onTest()
 #else
        wxFileDialog dialog(NULL, _T("Choose a File"), _T(""), _T(""), _T("*.*"), wxFD_OPEN );
 #endif
-
        if (dialog.ShowModal() == wxID_OK)
        {
                filename = dialog.GetPath().mb_str();
@@ -1129,9 +1015,9 @@ void wxContourMainFrame::onTest()
        } // if
 }
 
-
 //AD: 29-05-09
-void wxContourMainFrame::onSave(){
+void wxContourMainFrame::onSave()
+{
 //CMRU 03-09-09-----------------------------------------------------------------------------------------------
        if(GetFileLocation().empty())
        {
@@ -1148,10 +1034,7 @@ void wxContourMainFrame::onSave(){
                        kernelManager->setCurrentFileName(fileNameContourROI);
                        saveFileWithContours( fileNameContourROI );
                }
-       }
-//CMRU 03-09-09-----------------------------------------------------------------------------------------------
-       else
-       {
+       } else {//CMRU 03-09-09-----------------------------------------------------------------------------------------------
                std::cout<<"Nombre del archivo: "<<GetFileLocation()<<std::endl;
                std::string fileNameContourROI = GetFileLocation();
                std::cout<<"Nombre del archivo: "<<fileNameContourROI<<std::endl;
@@ -1291,12 +1174,13 @@ void wxContourMainFrame::readDimSpc(FILE *pFile, bool interactiveOption)
                {
                        wxDialog* dial = new wxDialog (this,-1,_T("Tools"),wxDefaultPosition, wxSize(460,360));
                        wxSizer* buttonsSizer = dial->CreateSeparatedButtonSizer(wxOK|wxCANCEL);
-                       wxBoxSizer *dialSizer = new wxBoxSizer(wxVERTICAL);
+//                     wxBoxSizer              *dialSizer = new wxBoxSizer(wxVERTICAL);
+                       wxFlexGridSizer *dialSizer = new wxFlexGridSizer(1);
 
                        wxString lstOptOperation[2];
                        lstOptOperation[0]=_T("By Dimension");
                        lstOptOperation[1]=_T("By Spacing");
-                       wxRadioBox * radioOpts = new wxRadioBox(dial, -1, _T("Type of Transformation"), wxDefaultPosition, wxSize(270,45), 2 , lstOptOperation,  2, wxRA_SPECIFY_COLS);
+                       wxRadioBox * radioOpts = new wxRadioBox(dial, -1, _T("Type of Transformation"), wxDefaultPosition, wxSize(270,90), 2 , lstOptOperation,  2, wxRA_SPECIFY_COLS);
 
                        dialSizer->Add( new wxStaticText(dial,-1,_T("The image resolution is not compatible with the "))  , 0, wxGROW );
                        dialSizer->Add( new wxStaticText(dial,-1,_T("selected Roi file. Do you want to tranform it?"))  , 0, wxGROW );
@@ -1482,14 +1366,16 @@ void wxContourMainFrame::openContours( FILE *pFile, FILE *pFileData, bool static
 
 }
 
-void wxContourMainFrame::RefreshInterface(){
+void wxContourMainFrame::RefreshInterface()
+{
        changeInstant();
        _theViewPanel->RefreshInterface();
        //wxVtk2DBaseView *wxvtk2dbaseview = (wxVtk2DBaseView*)this->_theViewPanel->getWxVtkBaseView();
        //wxvtk2dbaseview->Refresh();
 }
 
-vtkImageData* wxContourMainFrame::getImageData(){
+vtkImageData* wxContourMainFrame::getImageData()
+{
        return _theViewPanel->getImageData();
 }
 
@@ -1512,473 +1398,321 @@ void wxContourMainFrame::SegmentationOneSliceITK(int x, int y, int z, wxString d
        int typeofcontour = 1;
        //Image Data
        vtkImageData    *imagedata      = getImageData();
-
        //Tipo de pixeles a utilizar internamente en ITK
        typedef   float  InternalPixelType;
        const     unsigned int    Dimension = 2;
        typedef itk::Image< InternalPixelType, Dimension >  InternalImageType;
-
        //Tipo de pixeles de salida 1
        typedef unsigned char OutputPixelType;
        typedef itk::Image< OutputPixelType, Dimension > OutputImageType;
-
        //Tipo de pixeles de salida 2
        typedef unsigned short OutputPixelType2;
        typedef itk::Image< OutputPixelType2, Dimension > OutputImageType2;
-
        //DefiniciÛn del thresholder
        typedef itk::BinaryThresholdImageFilter<
                                                        InternalImageType,
                                                        OutputImageType    >    ThresholdingFilterType;
-
        //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
        typedef itk::CastImageFilter<
                                OutputImageType2, InternalImageType >  CastFilterType2;
-
        //Tercer tipo de conversiÛn
        typedef itk::RescaleIntensityImageFilter<
                                                                InternalImageType,
                                                                OutputImageType >   CastFilterType3;
-
        //Cuarto tipo de conversiÛn
        typedef itk::RescaleIntensityImageFilter<
                                                                OutputImageType,
                                                                OutputImageType >   CastFilterType4;
-
        ThresholdingFilterType::Pointer thresholder = ThresholdingFilterType::New();
-
        thresholder->SetLowerThreshold( 0.0 );
        thresholder->SetUpperThreshold( 128 );
-
        thresholder->SetOutsideValue(  255  );
        thresholder->SetInsideValue(  0 );
-
        //DefinciÛn de conexiÛnes entre VTK e ITK y el writer
        typedef itk::VTKImageToImageFilter<OutputImageType2> ConnectorType;
        typedef itk::ImageToVTKImageFilter<OutputImageType> ConnectorType2;
 //EED10JUIN2011        typedef  itk::ImageFileWriter<  OutputImageType  > WriterType;
-
-       ConnectorType::Pointer connector= ConnectorType::New();
-       ConnectorType2::Pointer connector2= ConnectorType2::New();
-
-
-       CastFilterType::Pointer filter=CastFilterType::New();
-       CastFilterType2::Pointer filter2=CastFilterType2::New();
-
+       ConnectorType::Pointer connector        = ConnectorType::New();
+       ConnectorType2::Pointer connector2      = ConnectorType2::New();
+       CastFilterType::Pointer filter   = CastFilterType::New();
+       CastFilterType2::Pointer filter2 = CastFilterType2::New();
        connector->SetInput( imagedata );
        filter2->SetInput(connector->GetOutput());
-
        typedef   itk::CurvatureAnisotropicDiffusionImageFilter<
                                                                InternalImageType,
                                                                InternalImageType >  SmoothingFilterType;
-
        SmoothingFilterType::Pointer smoothing = SmoothingFilterType::New();
-
        typedef   itk::GradientMagnitudeRecursiveGaussianImageFilter<
                                                                InternalImageType,
                                                                InternalImageType >  GradientFilterType;
-
        typedef   itk::SigmoidImageFilter<
                                                                InternalImageType,
                                                                InternalImageType >  SigmoidFilterType;
-
        GradientFilterType::Pointer  gradientMagnitude = GradientFilterType::New();
-
        SigmoidFilterType::Pointer sigmoid = SigmoidFilterType::New();
-
        sigmoid->SetOutputMinimum(  0.0  );
        sigmoid->SetOutputMaximum(  255.0  );
-
-
        typedef  itk::FastMarchingImageFilter<
                                                                InternalImageType,
                                                                InternalImageType >    FastMarchingFilterType;
-
-
        FastMarchingFilterType::Pointer  fastMarching = FastMarchingFilterType::New();
-       
-printf("EED wxContourMainFrame::SegmentationOneSliceITK  Skypping this code.  Not compiling   GeodesicActiveContourLevelSetImageFilter     in fedora 17 with ITK 3.20.1 \n ");
-/*EED 17 sept 2010
- > - > - > - > - > - > - > - > - > - > - > - > - > - > - > - > -   Not compiling   GeodesicActiveContourLevelSetImageFilter     in fedora 17 with ITK 3.20.1
        typedef  itk::GeodesicActiveContourLevelSetImageFilter< InternalImageType,
                                        InternalImageType >    GeodesicActiveContourFilterType;
        GeodesicActiveContourFilterType::Pointer geodesicActiveContour =
                                                                                GeodesicActiveContourFilterType::New();
-
        typedef  itk::ZeroCrossingImageFilter<
                                                                InternalImageType,
                                                                InternalImageType >    ZeroCrossingFilterType;
        ZeroCrossingFilterType::Pointer zeroCrossing =
                                                        ZeroCrossingFilterType::New();
-
        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(  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->SetFeatureImage( sigmoid->GetOutput() );
-
        zeroCrossing->SetInput( geodesicActiveContour->GetOutput() );
        //thresholder->SetInput( zeroCrossing->GetOutput() );
        thresholder->SetInput( geodesicActiveContour->GetOutput() );
        connector2->SetInput( thresholder->GetOutput()  );
-< - < - < - < - < - < - < - < - < - < - < - < - < - < - < - < - < - < - < - < - 
-*/
-       
-       
-       
-
        smoothing->SetTimeStep( 0.125 );
        smoothing->SetNumberOfIterations(  5 );
        smoothing->SetConductanceParameter( 9.0 );
-
-
-       const double sigma = atof(  crea::wx2std(sigm).c_str() );
+       const double sigma      = atof(  crea::wx2std(sigm).c_str() );
        gradientMagnitude->SetSigma(  sigma  );
-
-       const double alpha =  atof(  crea::wx2std(alf).c_str() );
-       const double beta  =  atof(  crea::wx2std(bet).c_str() );
-
+       const double alpha      =  atof(  crea::wx2std(alf).c_str() );
+       const double beta       =  atof(  crea::wx2std(bet).c_str() );
        sigmoid->SetAlpha( alpha );
        sigmoid->SetBeta(  beta  );
-
        typedef FastMarchingFilterType::NodeContainer  NodeContainer;
        typedef FastMarchingFilterType::NodeType       NodeType;
-
        NodeContainer::Pointer seeds = NodeContainer::New();
-
        InternalImageType::IndexType  seedPosition;
-       seedPosition[0] = x;
-       seedPosition[1] = y;
-
+       seedPosition[0]         = x;
+       seedPosition[1]         = y;
        const double initialDistance = atof( crea::wx2std(distanc).c_str() );
-
        NodeType node;
-
        const double seedValue = - initialDistance;
-
        node.SetValue( seedValue );
        node.SetIndex( seedPosition );
-
        seeds->Initialize();
        seeds->InsertElement( 0, node );
-
        fastMarching->SetTrialPoints( seeds );
-
        fastMarching->SetSpeedConstant( 1.0 );
-
-       fastMarching->SetOutputSize(
-                       connector->GetOutput()->GetBufferedRegion().GetSize() );
-
+       fastMarching->SetOutputSize( connector->GetOutput()->GetBufferedRegion().GetSize() );
        fastMarching->SetStoppingValue( 800 );
        try
     {
-
-       connector2->Update();
-       vtkImageData *idata = connector2->GetOutput();
-
-       vtkMarchingContourFilter* cntVTK = vtkMarchingContourFilter::New( );
-
-       cntVTK->SetNumberOfContours( 1 );
-       cntVTK->SetValue( 0, 255 );
+//             connector2->Update();
+               vtkImageData *idata = NULL;
+               idata                           = connector2->GetOutput();
+               vtkMarchingContourFilter* cntVTK = vtkMarchingContourFilter::New( );
+               cntVTK->SetNumberOfContours( 1 );
+               cntVTK->SetValue( 0, 255 );
 //EED 2017-01-01 Migration VTK7
 #if VTK_MAJOR_VERSION <= 5
-       cntVTK->SetInput( idata );
-       cntVTK->Update( );
-       cntVTK->UpdateInformation();
+               cntVTK->SetInput( idata );
+               cntVTK->Update( );
+               cntVTK->UpdateInformation();
 #else
-       cntVTK->SetInputData( idata );
+               cntVTK->SetInputData( idata );
 #endif
-
-       vtkCleanPolyData* cpd = vtkCleanPolyData::New( );
-
+               vtkCleanPolyData* cpd = vtkCleanPolyData::New( );
 //EED 2017-01-01 Migration VTK7
 #if VTK_MAJOR_VERSION <= 5
-       cpd->SetInput( cntVTK->GetOutput( ) );
-       cpd->Update( );
-       cpd->UpdateInformation();
+               cpd->SetInput( cntVTK->GetOutput( ) );
+               cpd->Update( );
+               cpd->UpdateInformation();
 #else
-       cpd->SetInputData( cntVTK->GetOutput( ) );
+               cpd->SetInputData( cntVTK->GetOutput( ) );
 #endif
-
-
-       vtkPolyDataConnectivityFilter* conn = vtkPolyDataConnectivityFilter::New( );
-       conn->SetExtractionModeToLargestRegion( );
+               vtkPolyDataConnectivityFilter* conn = vtkPolyDataConnectivityFilter::New( );
+               conn->SetExtractionModeToLargestRegion( );
 //EED 2017-01-01 Migration VTK7
 #if VTK_MAJOR_VERSION <= 5
-       conn->SetInput( cpd->GetOutput( ) );
-       conn->Update( );
-       conn->UpdateInformation();
+               conn->SetInput( cpd->GetOutput( ) );
+               conn->Update( );
+               conn->UpdateInformation();
 #else
-       conn->SetInputData( cpd->GetOutput( ) );
+               conn->SetInputData( cpd->GetOutput( ) );
 #endif
-
-
-       vtkStripper* vtkstripper = vtkStripper::New( );
-
+               vtkStripper* vtkstripper = vtkStripper::New( );
 //EED 2017-01-01 Migration VTK7
 #if VTK_MAJOR_VERSION <= 5
-       vtkstripper->SetInput( conn->GetOutput() );
-       vtkstripper->Update();
-       vtkstripper->UpdateInformation();
+               vtkstripper->SetInput( conn->GetOutput() );
+               vtkstripper->Update();
+               vtkstripper->UpdateInformation();
 #else
-       vtkstripper->SetInputData( conn->GetOutput() );
-       cntVTK->Update();
+               vtkstripper->SetInputData( conn->GetOutput() );
+               cntVTK->Update();
 #endif
-
-       vtkPolyData* polyDataResult =  cntVTK->GetOutput();
-       //std::cout<<"Points "<<polyDataResult->GetNumberOfPoints()<<std::endl;
+               vtkPolyData* polyDataResult =  cntVTK->GetOutput();
+               //std::cout<<"Points "<<polyDataResult->GetNumberOfPoints()<<std::endl;
 
 //EED 2017-01-01 Migration VTK7
 #if VTK_MAJOR_VERSION <= 5
-       polyDataResult->Update( );
-       polyDataResult->UpdateInformation();
+               polyDataResult->Update( );
+               polyDataResult->UpdateInformation();
 #else
        // ..
 #endif
-
-       //EED
-               /*
-       ofstream myfile;
-       myfile.open ("C:/Creatis/example.txt");
-       myfile << "\n";
-       polyDataResult->Print(myfile);
-       myfile << "-------------------------------------\n";
-       polyDataResult->GetLines()->Print(myfile);
-       myfile.close();*/
-
-       cntVTK          -> Delete();
-       cpd                     -> Delete();
-       conn            -> Delete();
-
-
-       //--Calculating control points
-
-       std::vector<double> vecX;
-       std::vector<double> vecY;
-       std::vector<double> vecXo;
-       std::vector<double> vecYo;
-       std::vector<double>::iterator vecXoi;
-       std::vector<double>::iterator vecYoi;
-       std::vector<double> vecZ;
-
-       std::vector<double> vecCtrlPointX;
-       std::vector<double> vecCtrlPointY;
-       std::vector<double> vecCtrlPointZ;
-
-
-       double *p;
-       double xAct=0;
-       double yAct=0;
-       int ii,size=polyDataResult->GetNumberOfPoints();
-       ofstream myfile;
-       myfile.open ("C:/Creatis/example2.txt");
-
-       size=polyDataResult->GetNumberOfPoints();
-       for (ii=0;ii<size;ii++)
-       {
-               if(ii==0)
+               cntVTK          -> Delete();
+               cpd                     -> Delete();
+               conn            -> Delete();
+               //--Calculating control points
+               std::vector<double> vecX;
+               std::vector<double> vecY;
+               std::vector<double> vecXo;
+               std::vector<double> vecYo;
+               std::vector<double>::iterator vecXoi;
+               std::vector<double>::iterator vecYoi;
+               std::vector<double> vecZ;
+               std::vector<double> vecCtrlPointX;
+               std::vector<double> vecCtrlPointY;
+               std::vector<double> vecCtrlPointZ;
+               double *p;
+               double xAct=0;
+               double yAct=0;
+               int ii,size=polyDataResult->GetNumberOfPoints();
+               ofstream myfile;
+               myfile.open ("C:/Creatis/example2.txt");
+               size=polyDataResult->GetNumberOfPoints();
+               for (ii=0;ii<size;ii++)
                {
-                       xAct=x;
-                       yAct=y;
-               }
-               p       = polyDataResult->GetPoint(ii);
-               double x=p[0];
-               double y=p[1];
-               /*if(fabs(yAct-y)>20)
-               {
-                       if((xAct-x)>1 || (xAct-x)<-1)
+                       if(ii==0)
                        {
-                       vecX.push_back( p[0] );
-                       vecY.push_back( p[1] );
-                       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;
-                       vecZ.push_back( 900 );
-                       xAct=x;
-                       yAct=y;
+                               xAct=x;
+                               yAct=y;
                        }
-                       else
+                       p       = polyDataResult->GetPoint(ii);
+                       double x=p[0];
+                       double y=p[1];
+//                     if(fabs(yAct-y)>20)
+//                     {
+//                             if((xAct-x)>1 || (xAct-x)<-1)
+//                             {
+//                                     vecX.push_back( p[0] );
+//                                     vecY.push_back( p[1] );
+//                                     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;
+//                                     vecZ.push_back( 900 );
+//                                     xAct=x;
+//                                     yAct=y;
+//                             } else {
+//                                     vecXo.push_back(p[0]);
+//                                     vecYo.push_back(p[1]);
+//                             }
+//                     } else
+
+                       if(fabs(xAct-x)>11)
                        {
                                vecXo.push_back(p[0]);
                                vecYo.push_back(p[1]);
-                       }
-
-               }
-               else*/ if(fabs(xAct-x)>11)
-               {
-                       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
 
-       }
 
-       while(!vecXo.empty())
-       {
-               vecX.push_back(vecXo.back());
-               //std::cout<<" x Siguiente "<<vecXo.back();
-               vecXo.pop_back();
-               vecZ.push_back( 900 );
-       }
-       while(!vecYo.empty())
-       {
-               vecY.push_back(vecYo.back());
-                       vecYo.pop_back();
-       }
-       myfile.close();
+               } // for
 
-       /*for(int l=0;l<vecX.size();l++)
-       {
-               if(l==0)
+               while(!vecXo.empty())
                {
-            vecXo.push_back(p[0]);
-                       vecYo.push_back(p[1]);
+                       vecX.push_back(vecXo.back());
+                       //std::cout<<" x Siguiente "<<vecXo.back();
+                       vecXo.pop_back();
+       //EED 2017-12-16
+                       vecZ.push_back( -900 );
                }
-               else
+               while(!vecYo.empty())
                {
-                       if(vecXoi[l-1]==)
-                       {
-                       }
+                       vecY.push_back(vecYo.back());
+                               vecYo.pop_back();
                }
+               myfile.close();
+
+//             for(int l=0;l<vecX.size();l++)
+//             {
+//                     if(l==0)
+//                     {
+//               vecXo.push_back(p[0]);
+//                             vecYo.push_back(p[1]);
+//                     } else {
+//                             if(vecXoi[l-1]==)
+//                             {
+//                             }
+//                     } // if
+//             } // for
+
+               ExtractControlPoints2D *extractcontrolpoints2d = new ExtractControlPoints2D();
+               extractcontrolpoints2d->SetContour( &vecX , &vecY , &vecZ );
+
+               int method=2;
+               if (method==0)
+               {
+                       extractcontrolpoints2d->GetInitialControlPoints( &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
+               } else if (method==1){
+                       extractcontrolpoints2d->GetControlPoints(  &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
+               } else if (method==2){
+                       extractcontrolpoints2d->SetSamplingControlPoints( 15 );
+                       extractcontrolpoints2d->GetSamplingControlPoints(  &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
+               }
+               //--Adding contour to the system
 
-       }*/
-
-       ExtractControlPoints2D *extractcontrolpoints2d = new ExtractControlPoints2D();
-
-       extractcontrolpoints2d->SetContour( &vecX , &vecY , &vecZ );
-
-       int method=2;
-       if (method==0){
-               extractcontrolpoints2d->GetInitialControlPoints( &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
-       }
-       else if (method==1){
-               extractcontrolpoints2d->GetControlPoints(  &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
-       }
-       else if (method==2){
-               extractcontrolpoints2d->SetSamplingControlPoints( 15 );
-               extractcontrolpoints2d->GetSamplingControlPoints(  &vecCtrlPointX , &vecCtrlPointY , &vecCtrlPointZ );
-       }
-       //--Adding contour to the system
-
-       std::vector<int> actualInstantVector;
-       _instantPanel->getInstant( actualInstantVector );
-       actualInstantVector[1]=z;
+               std::vector<int> actualInstantVector;
+               _instantPanel->getInstant( actualInstantVector );
+               actualInstantVector[1]=z;
 
-       int j,sizeCtrPt = vecCtrlPointX.size();
+               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++)
-               {
-                       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 )
-               {
-                       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
+               manualBaseModel *manModelContour =  kernelManager->factoryManualContourModel( typeofcontour );
+               manModelContour->SetNumberOfPointsSpline( ((sizeCtrPt/100)+1)*100 );
+               if (sizeCtrPt>=3){
+                       for (j=0 ; j<sizeCtrPt ; j++)
+                       {
+                               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 )
+                       {
+                               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
 
+    }   catch( itk::ExceptionObject & excep ) {
+               std::cerr << "Exception caught !" << std::endl;
+               std::cerr << excep << std::endl;
+    } // try
 
-/*
-//EED10JUIN2011
-                WriterType::Pointer writer = WriterType::New();
-         CastFilterType3::Pointer caster = CastFilterType3::New();
-
-                caster->SetInput( gradientMagnitude->GetOutput() );
-                writer->SetInput( caster->GetOutput() );
-                writer->SetFileName("Gradient Magnitude.png");
-                caster->SetOutputMinimum(   0 );
-                caster->SetOutputMaximum( 255 );
-                writer->Update();
-
-                CastFilterType3::Pointer caster2 = CastFilterType3::New();
-                WriterType::Pointer writer2 = WriterType::New();
-
-                caster2->SetInput( sigmoid->GetOutput() );
-                writer2->SetInput( caster2->GetOutput() );
-                writer2->SetFileName("Sigmoid.png");
-                caster2->SetOutputMinimum(   0 );
-                caster2->SetOutputMaximum( 255 );
-                writer2->Update();
-
-                CastFilterType3::Pointer caster3 = CastFilterType3::New();
-                WriterType::Pointer writer3 = WriterType::New();
-
-                caster3->SetInput( fastMarching->GetOutput() );
-                writer3->SetInput( caster3->GetOutput() );
-                writer3->SetFileName("FastMarching.bmp");
-                caster3->SetOutputMinimum(   0 );
-                caster3->SetOutputMaximum( 255 );
-                writer3->Update();
-
-                CastFilterType3::Pointer caster4 = CastFilterType3::New();
-                WriterType::Pointer writer4 = WriterType::New();
-
-                caster4->SetInput( geodesicActiveContour->GetOutput() );
-                writer4->SetInput( caster4->GetOutput() );
-                writer4->SetFileName("GeodesicActiveContour.png");
-                caster4->SetOutputMinimum(   0 );
-                caster4->SetOutputMaximum( 255 );
-                writer4->Update();
-
-                CastFilterType3::Pointer caster5 = CastFilterType3::New();
-                WriterType::Pointer writer5 = WriterType::New();
-
-                caster5->SetInput( zeroCrossing->GetOutput() );
-                writer5->SetInput( caster5->GetOutput() );
-                writer5->SetFileName("ZeroCrossing.bmp");
-                caster5->SetOutputMinimum(   0 );
-                caster5->SetOutputMaximum( 255 );
-                writer5->Update();
-*/
-    }
-  catch( itk::ExceptionObject & excep )
-    {
-    std::cerr << "Exception caught !" << std::endl;
-    std::cerr << excep << std::endl;
-    }
 }
 
 void wxContourMainFrame::onSegmentationOneSlice(int isovalue,int sampling,int method){
@@ -1998,20 +1732,21 @@ void wxContourMainFrame::onSegmentationOneSlice(int isovalue,int sampling,int me
 
 void wxContourMainFrame::SegmentationOneSlice( int x, int y, int z, int isovalue, int sampling, int method )
 {
+
+printf("EED wxContourMainFrame::SegmentationOneSlice Start=%d  %d\n",sampling,method );
+printf("EED wxContourMainFrame::SegmentationOneSlice %d %d %d  \n",x,y,z );
+
        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
@@ -2021,7 +1756,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();
 
@@ -2069,7 +1804,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
@@ -2094,11 +1837,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();
@@ -2123,116 +1870,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=0;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 );
+
+       //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 );
+               }*/
 
-       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++)
+               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;
@@ -2458,29 +2233,59 @@ int wxContourMainFrame::getWindowLevel()
        return _theViewPanel->getWindowLevel();
 }
 
-void wxContourMainFrame::onBrigthnessColorWindowLevel(int colorwindow,int windowlevel){
+void wxContourMainFrame::onBrigthnessColorWindowLevel(int colorwindow,int windowlevel)
+{
        _theViewPanel->onBrigthnessColorWindowLevel(colorwindow, windowlevel);
 
        RefreshInterface();
 }
 
-void wxContourMainFrame::onInterpolation(bool interpolate){
+void wxContourMainFrame::onInterpolation(bool interpolate)
+{
        _theViewPanel->onInterpolation(interpolate);
        RefreshInterface();
 }
 
-void wxContourMainFrame::onChangeInstant(std::string name,int actual){
-       _instantPanel->setConceptValue( name, actual );
+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::resetAppend(){
-       kernelManager->resetAppend();
+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(){
+void wxContourMainFrame::onSpreadAdd()
+{
        std::vector<double> vecX;
        std::vector<double> vecY;
        std::vector<double> vecZ;
@@ -2492,7 +2297,8 @@ void wxContourMainFrame::onSpreadAdd(){
 }
 
 
-void wxContourMainFrame::onSpreadAddAll(){
+void wxContourMainFrame::onSpreadAddAll()
+{
 //EED02
        std::vector<int> tempVector;
        _instantPanel->getInstant( tempVector );
@@ -2530,7 +2336,8 @@ void wxContourMainFrame::onSpreadAddAll(){
 }
 
 
-void wxContourMainFrame::onSpreadGo(int type){
+void wxContourMainFrame::onSpreadGo(int type)
+{
        //JCP 20-10-08 Undo redo implementation
        saveState();
        //JCP 20-10-08 Undo redo implementation
@@ -2585,11 +2392,13 @@ void wxContourMainFrame::onSpreadGo(int type){
        //RefreshInterface();
 }
 
-void wxContourMainFrame::getInstantVector(std::vector<int>& tempVector){
-
+void wxContourMainFrame::getInstantVector(std::vector<int>& tempVector)
+{
+       _instantPanel->getInstant( tempVector );
 }
 
-std::vector<std::string> wxContourMainFrame::getOutlinesName(int slide){
+std::vector<std::string> wxContourMainFrame::getOutlinesName(int slide)
+{
        std::vector<int> tempVector;
        _instantPanel->getInstant(tempVector);
        //Asignation of slide number should be different ex by name
@@ -3062,7 +2871,6 @@ void wxContourMainFrame::getMaskValue(vtkImageData ** mask,vtkImageData ** value
        int z;
        int ii,sizeLstContourThings;
 
-
        double totalZ;
        double porcent;
 
@@ -3143,13 +2951,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);
 
@@ -3163,6 +2999,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
@@ -3177,7 +3014,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);
@@ -3192,6 +3030,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
@@ -3248,7 +3087,9 @@ void wxContourMainFrame::onSnakePressed(){
 
 }
 
-void wxContourMainFrame::saveFileWithContoursAutomatique(){
+void wxContourMainFrame::saveFileWithContoursAutomatique()
+{
+printf("EED wxContourMainFrame::saveFileWithContoursAutomatique \n");
        std::string filename = kernelManager->getCurrentFileName();
        if(filename.compare("")!=0){
                saveFileWithContours(filename);
@@ -3256,7 +3097,8 @@ void wxContourMainFrame::saveFileWithContoursAutomatique(){
                onSave();
        }
 }
-void wxContourMainFrame::ShowToolsPanel(bool show){
+void wxContourMainFrame::ShowToolsPanel(bool show)
+{
        _pannew->Show(show);
 }
 
@@ -3270,7 +3112,7 @@ void wxContourMainFrame::SetContourGroup(int contourGroup)
 
 wxPanel *wxContourMainFrame::getMaskImageViewPanel(wxWindow *parent)
 {
-
        if      (_viewMaskImagePanel==NULL)
        {
                _contourextractdata = new ContourExtractData(true);
@@ -3279,15 +3121,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;
 }
@@ -3303,10 +3159,37 @@ 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;
 }
 
+void wxContourMainFrame::SetColorLayerImage(vtkImageData* img)
+{
+       if (_viewColorLayerImagePanel!=NULL) 
+       {
+               if (img!=NULL)
+               {
+                       _viewColorLayerImagePanel->SetActive(true);
+                       _viewColorLayerImagePanel->SetImage( img );
+//                     RefreshInterface();
+               } // img
+       } // _viewColorLayerImagePanel
+}
 
 wxPanel *wxContourMainFrame::getThresholdImageViewPanel(wxWindow *parent)
 {
@@ -3314,15 +3197,28 @@ wxPanel *wxContourMainFrame::getThresholdImageViewPanel(wxWindow *parent)
        {
                double range[2];
                _theViewPanel->getImageData()->GetScalarRange(range);
-               int min = (int)floor (range[0]);
-               int max = (int)ceil (range[1]);
-               _viewThresholdImagePanel = new ThresholdImageViewPanel( parent, min, max , 1);
-               _viewThresholdImage = new ThresholdImageView();
+               int min                                         = (int)floor (range[0]);
+               int max                                         = (int)ceil (range[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->SetminMaxValue(range[0],range[1]);
                _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;
 }
@@ -3420,6 +3316,9 @@ void wxContourMainFrame::SetFileLocation(std::string newFileLocation)
 {
        _fileLocation = newFileLocation;
        std::cout<<"GETFILELOCATION:"<<_fileLocation<<std::endl;
+
+//EED 2019-12-17
+       kernelManager->setCurrentFileName(_fileLocation);
 }
 //------------------------------------------------------------------------------------------------------------
 
@@ -3432,7 +3331,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
@@ -3442,7 +3340,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
@@ -3451,7 +3348,6 @@ void wxContourMainFrame::FlipActualSliceContours(int flipDirection)
 
 void wxContourMainFrame::FlipAllContours(int flipDirection)
 {
-       printf("EED wxContourMainFrame::FlipAllContours  flipDirection=%d  \n",  flipDirection );
        FlipLstOfContours( 
                            kernelManager->GetLstNameThings() , 
                            flipDirection
@@ -3472,7 +3368,6 @@ void wxContourMainFrame::FlipLstOfContours(std::vector<std::string> lstNameThing
        sizeLstNameThings       = lstNameThings.size(); 
        for (i=0 ; i<sizeLstNameThings ; i++)
        {
-               printf("EED wxContourMainFrame::FlipAllContours contour=%d \n", i );
                manualBaseModel * manualModel = kernelManager->getOutlineByKeyName(lstNameThings[i]);
                if (flipDirection==0)
                {
@@ -3500,8 +3395,6 @@ void wxContourMainFrame::FlipLstOfContours(std::vector<std::string> lstNameThing
 
 void wxContourMainFrame::FlipContours(int typeContourFlip, int flipDirection)
 {
-       printf("EED wxContourMainFrame::FlipContours( %d , %d ) \n",typeContourFlip,  flipDirection );
-
        saveState();
        
        if (typeContourFlip==0) // 
@@ -3521,7 +3414,20 @@ void wxContourMainFrame::FlipContours(int typeContourFlip, int flipDirection)
        RefreshInterface();
 }
 
+void wxContourMainFrame::SetXY(int x, int y)
+{
+       _theViewPanel->SetXY(x,y);
+}
 
+wxVtkBaseView* wxContourMainFrame::GetWxVtkBaseView()
+{
+       if (_theViewPanel!=NULL)
+       {
+               return _theViewPanel->getWxVtkBaseView();;
+       } else {
+               return NULL;
+       }
+}