]> Creatis software - creaContours.git/blobdiff - lib/Interface_ManagerContour_NDimensions/wxContourMainFrame.cxx
*** empty log message ***
[creaContours.git] / lib / Interface_ManagerContour_NDimensions / wxContourMainFrame.cxx
index 3070d426e1bf0325c8ae1fa9655246512d548af6..4ca7174761dd6c4e30e7a090ef083ab15a51120e 100644 (file)
@@ -58,6 +58,8 @@ char wxContourMainFrame::COPY = 'C';
                _refLineControl                         = NULL;
                _refLineModel                           = NULL;
                _refLineView                            = NULL;
+               _thresholdActor                         = NULL;
+               _actorPresent                           = false;
 
                // set up default notebook style
 //             m_notebook_style =wxAUI_NB_TAB_SPLIT | wxAUI_NB_TAB_EXTERNAL_MOVE | wxNO_BORDER;
@@ -77,6 +79,8 @@ char wxContourMainFrame::COPY = 'C';
                _refLineControl                         = NULL;
                _refLineModel                           = NULL;
                _refLineView                            = NULL;
+               _thresholdActor                         = NULL;
+               _actorPresent                           = false;
 
                // set up default notebook style
                m_notebook_style =wxAUI_NB_TAB_SPLIT | wxAUI_NB_TAB_EXTERNAL_MOVE | wxNO_BORDER;
@@ -592,11 +596,12 @@ void wxContourMainFrame::changeInstant()
                updateInstantOutlines();
                updateInstantImageData();
                updateInstantAxes();
-
-               
+               updateActor();
 }
-void wxContourMainFrame :: updateInstantOutlines(){
-       
+
+
+void wxContourMainFrame :: updateInstantOutlines()
+{
        _theViewPanel->removeSceneContours();
        _theViewPanel->addNameWrapperToScene();
 
@@ -613,7 +618,8 @@ std::string wxContourMainFrame::getNameWrapping(int i){
        return kernelManager->getNameWrapping(i);
 }
 
-void wxContourMainFrame::updateInstantImageData(){
+void wxContourMainFrame::updateInstantImageData()
+{
        std::vector<int> inst;
        _instantPanel->getInstant(inst);
 
@@ -627,7 +633,84 @@ void wxContourMainFrame::updateInstantImageData(){
        }
 }
 
-void wxContourMainFrame::updateInstantAxes(){
+void wxContourMainFrame::updateInstantAxes()
+{
+}
+
+void wxContourMainFrame :: updateActor()
+{
+               if (_actorPresent)
+               {
+                       int minMax[2];
+                       _pannew->onThresholdInstantChange(minMax);
+
+                       int     z = _theViewPanel->GetZ();
+
+                       vtkImageData * img = getImageData();
+                       vtkImageReslice *imageReslice = vtkImageReslice::New();
+
+                       imageReslice->SetInput( img );
+                       imageReslice->SetInformationInput(img);
+                       imageReslice->SetResliceAxesDirectionCosines(1,0,0, 0,1,0 ,0,0,1);
+                       imageReslice->SetResliceAxesOrigin(0,0,z);
+                       imageReslice->SetOutputDimensionality(2);
+                       imageReslice->SetInterpolationModeToLinear();
+
+                       img = imageReslice->GetOutput();
+                       img->Update();
+                       img->UpdateInformation();
+
+                       wxVtkBaseView * baseView = _theViewPanel->getWxVtkBaseView();
+                       
+                       double range[2];
+                       img->GetScalarRange(range);
+
+                       int min = floor (range[0]);
+                       int max = ceil (range[1]);
+
+                       //Lookup Table
+                       vtkLookupTable *lookup = vtkLookupTable::New();
+                       lookup->SetNumberOfTableValues(max+1);
+                       lookup->SetTableRange(range); 
+                       lookup->SetAlphaRange(0, 1);
+                       lookup->SetValueRange(0, 1);
+                       lookup->SetSaturationRange(0, 0); 
+                       lookup->SetRampToLinear( );
+
+                       //Assign a fake color for the upper image, and set the white as transparent
+                       int i;
+                       int minVal = minMax[0];
+                       int maxVal = minMax[1];
+
+                       for(i = min; i <= max; i++)
+                       {
+                               if( i >= minVal && i <= maxVal )
+                               {
+                                       lookup->SetTableValue(i, 1.0, 0.0, 0.0, 1);
+                               }
+                               else if( i >= min && i < minVal )
+                               {
+                                       lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                               }
+                               else if( i > maxVal && i < max )
+                               {
+                                       lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                               }
+                               else
+                               {
+                                       lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                               }
+                       }
+                               
+                       lookup->Build( );
+
+                       vtkImageMapToColors *mapperImage = vtkImageMapToColors::New( );
+                       mapperImage->SetLookupTable( lookup );
+                       mapperImage->SetInput( img );
+                       //mapperImage->SetOutputFormatToRGBA( );
+
+                       _thresholdActor->SetInput( mapperImage->GetOutput() );
+               }
 }
 
 void wxContourMainFrame::onChangeDeep(int val){
@@ -1576,49 +1659,52 @@ void wxContourMainFrame::onSegmentationAllSlice(int minZ,int maxZ,int isovalue,i
 //Creates a reference line to be used as the axis of the mirroring
 void wxContourMainFrame::referenceLine()
 {
-    wxBusyCursor wait;
-       double spc[3];
+       if (_refLineControl==NULL)
+       {
+               wxBusyCursor wait;
+               double spc[3];
 
-       vtkImageData *vtkimagedata = _theViewPanel->getImageData();
-       vtkimagedata->GetSpacing(spc);
+               vtkImageData *vtkimagedata = _theViewPanel->getImageData();
+               vtkimagedata->GetSpacing(spc);
 
-       wxVtkBaseView * viewer2D = _theViewPanel->getWxVtkBaseView();
+               wxVtkBaseView * viewer2D = _theViewPanel->getWxVtkBaseView();
 
-       _refLineControl = new manualLineControler();
-       _refLineModel   = new manualContourModelLine();
-       _refLineView = new manualViewLine();
-       _refLineView->SetModel( _refLineModel );
-       _refLineView->SetWxVtkBaseView( viewer2D );
-       _refLineView->SetRange( 2 );
-       _refLineView->SetZ( 1200 );
+               _refLineControl = new manualLineControler();
+               _refLineModel   = new manualContourModelLine();
+               _refLineView = new manualViewLine();
+               _refLineView->SetModel( _refLineModel );
+               _refLineView->SetWxVtkBaseView( viewer2D );
+               _refLineView->SetRange( 2 );
+               _refLineView->SetZ( 1200 );
 
-       _refLineView->SetSpacing(spc);
+               _refLineView->SetSpacing(spc);
 
-       _refLineView->SetColorNormalContour(0, 0, 1);
-       _refLineView->SetColorEditContour(0, 0.5, 0.5);
-       _refLineView->SetWidthLine(4);
+               _refLineView->SetColorNormalContour(0, 0, 1);
+               _refLineView->SetColorEditContour(0, 0.5, 0.5);
+               _refLineView->SetWidthLine(4);
 
-       _refLineControl->SetModelView( _refLineModel , _refLineView );
-       viewer2D->GetInteractorStyleBaseView()->AddInteractorStyleMaracas( _refLineControl );
-       _refLineModel->SetCloseContour(false);
-       _refLineControl->CreateNewManualContour();
+               _refLineControl->SetModelView( _refLineModel , _refLineView );
+               viewer2D->GetInteractorStyleBaseView()->AddInteractorStyleMaracas( _refLineControl );
+               _refLineModel->SetCloseContour(false);
+               _refLineControl->CreateNewManualContour();
 
-       double z = _refLineControl->GetZ();
-/*
-       _refLineControl->AddPoint(50,100,z);
-       _refLineControl->AddPoint(80,40,z);
-       _refLineControl->SetCompleteCreation(true);
+               double z = _refLineControl->GetZ();
+       /*
+               _refLineControl->AddPoint(50,100,z);
+               _refLineControl->AddPoint(80,40,z);
+               _refLineControl->SetCompleteCreation(true);
 
-       _refLineModel->AddPoint(50,100,z);
-       _refLineModel->AddPoint(80,40,z);
-       _refLineModel->UpdateSpline();
+               _refLineModel->AddPoint(50,100,z);
+               _refLineModel->AddPoint(80,40,z);
+               _refLineModel->UpdateSpline();
 
-       _refLineView->UpdateViewPoint(0);
-       _refLineView->UpdateViewPoint(1);
-*/
+               _refLineView->UpdateViewPoint(0);
+               _refLineView->UpdateViewPoint(1);
+       */
 
-       _refLineControl->SetActive(true);
-       _refLineView->RefreshContour();
+               _refLineControl->SetActive(true);
+               _refLineView->RefreshContour();
+       }
 }
 
 //Hides the referenceLine
@@ -1626,9 +1712,14 @@ void wxContourMainFrame::refLineHide()
 {
        if (_refLineControl!=NULL)
        {
-               _refLineView->RemoveCompleteContourActor();
-               _refLineControl->SetEditable(false);
-               _refLineControl->SetActive(false);
+               if (_refLineControl->GetActive())
+               {
+                       _refLineView->RemoveCompleteContourActor();
+                       _refLineControl->SetEditable(false);
+                       _refLineControl->SetActive(false);
+
+                       RefreshInterface();
+               }
        }
 }
 
@@ -1637,18 +1728,23 @@ void wxContourMainFrame::refLineShow()
 {
        if (_refLineControl!=NULL)
        {
-               wxVtkBaseView * viewer2D = _theViewPanel->getWxVtkBaseView();
+               if (_refLineControl->GetActive()==false)
+               {
+                       wxVtkBaseView * viewer2D = _theViewPanel->getWxVtkBaseView();
 
-               _refLineView->SetModel( _refLineModel );
-               _refLineView->SetWxVtkBaseView( viewer2D );
-               _refLineView->SetRange( 2 );
-               _refLineView->SetZ( 1200 );
+                       _refLineView->SetModel( _refLineModel );
+                       _refLineView->SetWxVtkBaseView( viewer2D );
+                       _refLineView->SetRange( 2 );
+                       _refLineView->SetZ( 1200 );
 
-               _refLineControl->SetModelView( _refLineModel , _refLineView );
-               viewer2D->GetInteractorStyleBaseView()->AddInteractorStyleMaracas( _refLineControl );
-               _refLineControl->CreateNewManualContour();
-               _refLineControl->SetActive(true);
-               _refLineView->RefreshContour();
+                       _refLineControl->SetModelView( _refLineModel , _refLineView );
+                       viewer2D->GetInteractorStyleBaseView()->AddInteractorStyleMaracas( _refLineControl );
+                       _refLineControl->CreateNewManualContour();
+                       _refLineControl->SetActive(true);
+                       _refLineView->RefreshContour();
+
+                       RefreshInterface();
+               }
        }
 }
 
@@ -1697,42 +1793,157 @@ void wxContourMainFrame::onMirror()
 
 void wxContourMainFrame::onThreshold(int minVal, int maxVal)
 {
-       
-       vtkImageData * img = _theViewPanel->getImageData();
-       
-       double range[2];
-       img->GetScalarRange(range);
-
-       //vtkImageData * fakeImage = new vtkImageData();
+       if (!_actorPresent)
+       {
+               int     z = _theViewPanel->GetZ();
 
-       //Lookup Table
-       vtkLookupTable *lookup = vtkLookupTable::New();
-       lookup->SetNumberOfTableValues(range[1]+1);
-       lookup->SetTableRange(range); 
-       lookup->SetAlphaRange(0, 1);
-       lookup->SetValueRange(0, 1);
-       lookup->SetSaturationRange(0, 0); 
-       lookup->SetRampToLinear( );
+               vtkImageData * img = getImageData();
+               vtkImageReslice *imageReslice = vtkImageReslice::New();
 
+               imageReslice->SetInput( img );
+               imageReslice->SetInformationInput(img);
+               imageReslice->SetResliceAxesDirectionCosines(1,0,0, 0,1,0 ,0,0,1);
+               imageReslice->SetResliceAxesOrigin(0,0,z);
+               imageReslice->SetOutputDimensionality(2);
+               imageReslice->SetInterpolationModeToLinear();
 
+               img = imageReslice->GetOutput();
+               img->Update();
+               img->UpdateInformation();
 
-       //Assign a fake color for the upper image, and set the white as transparent
-       for(int i = range[0]; i < range[1]; i++)
-       {
-               if( i >= minVal && i <= maxVal )
+               wxVtkBaseView * baseView = _theViewPanel->getWxVtkBaseView();
+               
+               double range[2];
+               img->GetScalarRange(range);
+
+               int min = floor (range[0]);
+               int max = ceil (range[1]);
+
+               //Lookup Table
+               vtkLookupTable *lookup = vtkLookupTable::New();
+               lookup->SetNumberOfTableValues(max+1);
+               lookup->SetTableRange(range); 
+               lookup->SetAlphaRange(0, 1);
+               lookup->SetValueRange(0, 1);
+               lookup->SetSaturationRange(0, 0); 
+               lookup->SetRampToLinear( );
+
+               //Assign a fake color for the upper image, and set the white as transparent
+               int i;
+               for(i = min; i <= max; i++)
                {
-                       lookup->SetTableValue(i, 1.0, 0.0, 0.0, 1);
+                       if( i >= minVal && i <= maxVal )
+                       {
+                               lookup->SetTableValue(i, 1.0, 0.0, 0.0, 1);
+                       }
+                       else if( i >= min && i < minVal )
+                       {
+                               lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                       }
+                       else if( i > maxVal && i < max )
+                       {
+                               lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                       }
+                       else
+                       {
+                               lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                       }
                }
-               else if( i >= range[0] && i < minVal )
+                       
+               lookup->Build( );
+
+               vtkImageMapToColors *mapperImage = vtkImageMapToColors::New( );
+               mapperImage->SetLookupTable( lookup );
+               mapperImage->SetInput( img );
+               //mapperImage->SetOutputFormatToRGBA( );
+
+               if (_thresholdActor==NULL)
                {
-                       lookup->SetTableValue(i, 0.0, 0.0, 0.0, 0); //transparent
+                       _thresholdActor = vtkImageActor::New( );
+                       _thresholdActor->SetInput( mapperImage->GetOutput() );
+                       _thresholdActor->SetOpacity( 0.6 );
+                       _thresholdActor->InterpolateOn( ); 
                }
-               else if( i > maxVal && i < range[1] )
+               else
                {
-                       lookup->SetTableValue(i, 0.0, 0.0, 0.0, 0); //transparent
+                       _thresholdActor->SetInput( mapperImage->GetOutput() );
                }
+
+               baseView->GetRenderer()->AddActor( _thresholdActor );
+               _actorPresent = true;
        }
+
+       RefreshInterface();
+/*
+       vtkImageViewer2* viewer = vtkImageViewer2::New();        
+       viewer->SetInput( upperImageActor->GetInput() ); 
+       viewer->SetColorLevel((range[1]-range[0])/2);
+       viewer->SetColorWindow(range[1]);
+       viewer->GetRenderer()->AddActor( upperImageActor );
+       viewer->Render();
+*/
+
+}
+
+void wxContourMainFrame::onThresholdChange(int minVal, int maxVal)
+{
+       if (_actorPresent)
+       {
+               int     z = _theViewPanel->GetZ();
+
+               vtkImageData * img = getImageData();
+               vtkImageReslice *imageReslice = vtkImageReslice::New();
+
+               imageReslice->SetInput( img );
+               imageReslice->SetInformationInput(img);
+               imageReslice->SetResliceAxesDirectionCosines(1,0,0, 0,1,0 ,0,0,1);
+               imageReslice->SetResliceAxesOrigin(0,0,z);
+               imageReslice->SetOutputDimensionality(2);
+               imageReslice->SetInterpolationModeToLinear();
+
+               img = imageReslice->GetOutput();
+               img->Update();
+               img->UpdateInformation();
+
+               wxVtkBaseView * baseView = _theViewPanel->getWxVtkBaseView();
                
+               double range[2];
+               img->GetScalarRange(range);
+
+               int min = floor (range[0]);
+               int max = ceil (range[1]);
+
+               //Lookup Table
+               vtkLookupTable *lookup = vtkLookupTable::New();
+               lookup->SetNumberOfTableValues(max+1);
+               lookup->SetTableRange(range); 
+               lookup->SetAlphaRange(0, 1);
+               lookup->SetValueRange(0, 1);
+               lookup->SetSaturationRange(0, 0); 
+               lookup->SetRampToLinear( );
+
+               //Assign a fake color for the upper image, and set the white as transparent
+               int i;
+               for(i = min; i <= max; i++)
+               {
+                       if( i >= minVal && i <= maxVal )
+                       {
+                               lookup->SetTableValue(i, 1.0, 0.0, 0.0, 1);
+                       }
+                       else if( i >= min && i < minVal )
+                       {
+                               lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                       }
+                       else if( i > maxVal && i < max )
+                       {
+                               lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                       }
+                       else
+                       {
+                               lookup->SetTableValue(i, 1.0, 1.0, 1.0, 0); //transparent
+                       }
+               }
+                       
                lookup->Build( );
 
                vtkImageMapToColors *mapperImage = vtkImageMapToColors::New( );
@@ -1740,20 +1951,51 @@ void wxContourMainFrame::onThreshold(int minVal, int maxVal)
                mapperImage->SetInput( img );
                //mapperImage->SetOutputFormatToRGBA( );
 
-               vtkImageActor * upperImageActor = vtkImageActor::New( );
-               upperImageActor->SetInput( mapperImage->GetOutput() );
-               upperImageActor->SetOpacity( 1 );
+               _thresholdActor->SetInput( mapperImage->GetOutput() );
+       }
 
-               vtkImageViewer2* viewer = vtkImageViewer2::New();        
-               viewer->SetInput( upperImageActor->GetInput() ); 
-               viewer->SetColorLevel((range[1]-range[0])/2);
-               viewer->SetColorWindow(range[1]);
-               viewer->GetRenderer()->AddActor( upperImageActor );
-               viewer->Render();
+       RefreshInterface();
+}
+
+void wxContourMainFrame::onThresholdInterpolation(bool interpolate)
+{
+       if (_thresholdActor!=NULL)
+       {
+               if (interpolate)
+               {
+                       _thresholdActor->InterpolateOn( );
+               }
+
+               else
+               {
+                       _thresholdActor->InterpolateOff( );
+               }
 
                RefreshInterface();
+       }
 }
 
+void wxContourMainFrame::onThresholdChangeOpacity (int opacity)
+{
+       if (_actorPresent)
+       {
+               _thresholdActor->SetOpacity(opacity*0.1);
+       }
+
+       RefreshInterface();
+}
+
+void wxContourMainFrame::onThresholdRemove()
+{
+       if (_actorPresent)
+       {
+               wxVtkBaseView * baseView = _theViewPanel->getWxVtkBaseView();
+               baseView->GetRenderer()->RemoveActor( _thresholdActor );
+               _actorPresent = false;
+       }
+
+       RefreshInterface();
+}
 
 void wxContourMainFrame::showAxis(bool show)
 {
@@ -1802,7 +2044,6 @@ int wxContourMainFrame::getWindowLevel(){
 void wxContourMainFrame::onBrigthnessColorWindowLevel(int colorwindow,int windowlevel){
        _theViewPanel->onBrigthnessColorWindowLevel(colorwindow, windowlevel);
 
-       
        RefreshInterface();
 }