vv/vvToolBinarize.cxx

   1 /*=========================================================================
   2
   3   Program:   vv
   4   Module:    $RCSfile: vvToolBinarize.cxx,v $
   5   Language:  C++
   6   Date:      $Date: 2010/03/05 10:32:33 $
   7   Version:   $Revision: 1.11 $
   8   Author :   David Sarrut (david.sarrut@creatis.insa-lyon.fr)
   9
  10   Copyright (C) 2008
  11   Léon Bérard cancer center http://oncora1.lyon.fnclcc.fr
  12   CREATIS-LRMN http://www.creatis.insa-lyon.fr
  13
  14   This program is free software: you can redistribute it and/or modify
  15   it under the terms of the GNU General Public License as published by
  16   the Free Software Foundation, version 3 of the License.
  17
  18   This program is distributed in the hope that it will be useful,
  19   but WITHOUT ANY WARRANTY; without even the implied warranty of
  20   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  21   GNU General Public License for more details.
  22
  23   You should have received a copy of the GNU General Public License
  24   along with this program.  If not, see <http://www.gnu.org/licenses/>.
  25
  26   =========================================================================*/
  27
  28 #include "vvToolBinarize.h"
  29 #include "vvSlicerManager.h"
  30 #include "vvSlicer.h"
  31 #include "vvToolInputSelectorWidget.h"
  32
  33 #include "clitkBinarizeImageGenericFilter.h"
  34
  35 #include <vtkImageActor.h>
  36 #include <vtkCamera.h>
  37 #include <vtkImageClip.h>
  38
  39 //------------------------------------------------------------------------------
  40 // Create the tool and automagically (I like this word) insert it in
  41 // the main window menu.
  42 ADD_TOOL(vvToolBinarize);
  43 //------------------------------------------------------------------------------
  44
  45
  46 //------------------------------------------------------------------------------
  47 void vvToolBinarize::Initialize() {
  48   SetToolName("Binarize");
  49   SetToolMenuName("Binarize");
  50   SetToolIconFilename(":/new/prefix1/icons/binarize.png");
  51   SetToolTip("Image interactive binarization with thresholds.");
  52 }
  53 //------------------------------------------------------------------------------
  54
  55
  56 //------------------------------------------------------------------------------
  57 vvToolBinarize::vvToolBinarize(vvMainWindowBase * parent, Qt::WindowFlags f)
  58   :vvToolWidgetBase(parent,f),
  59    vvToolBase<vvToolBinarize>(parent),
  60    Ui::vvToolBinarize()
  61 {
  62   // GUI Initialization
  63   Ui_vvToolBinarize::setupUi(mToolWidget);
  64   mInteractiveDisplayIsEnabled = mCheckBoxInteractiveDisplay->isChecked();
  65
  66   // Connect signals & slots
  67   connect(mRadioButtonLowerThan, SIGNAL(toggled(bool)), this, SLOT(enableLowerThan(bool)));
  68   connect(mCheckBoxUseFG, SIGNAL(toggled(bool)), this, SLOT(useFGBGtoggled(bool)));
  69   connect(mCheckBoxUseBG, SIGNAL(toggled(bool)), this, SLOT(useFGBGtoggled(bool)));
  70   connect(mCheckBoxInteractiveDisplay, SIGNAL(toggled(bool)), this, SLOT(InteractiveDisplayToggled(bool)));
  71
  72   // Initialize some widget
  73   mThresholdSlider1->SetText("");
  74   mThresholdSlider2->SetText("");
  75   mFGSlider->SetText("Foreground value");
  76   mBGSlider->SetText("Background value");
  77
  78   // Main filter
  79   mFilter = new clitk::BinarizeImageGenericFilter<args_info_clitkBinarizeImage>; //DS PUT IN BASECLASS ?
  80 }
  81 //------------------------------------------------------------------------------
  82
  83
  84 //------------------------------------------------------------------------------
  85 vvToolBinarize::~vvToolBinarize() {
  86 }
  87 //------------------------------------------------------------------------------
  88
  89
  90 //------------------------------------------------------------------------------
  91 void vvToolBinarize::InteractiveDisplayToggled(bool b) {
  92   mInteractiveDisplayIsEnabled = b;
  93   if (!mInteractiveDisplayIsEnabled) {
  94     RemoveVTKObjects();
  95   }
  96   else {
  97     for(unsigned int i=0; i<mImageContour.size(); i++) {
  98       mImageContour[i]->showActors();
  99       mImageContourLower[i]->showActors();
 100     }
 101     if (mCurrentSlicerManager)
 102       mCurrentSlicerManager->Render();
 103   }
 104 }
 105 //------------------------------------------------------------------------------
 106
 107
 108 //------------------------------------------------------------------------------
 109 void vvToolBinarize::RemoveVTKObjects() {
 110   for(unsigned int i=0; i<mImageContour.size(); i++) {
 111     mImageContour[i]->hideActors();
 112     mImageContourLower[i]->hideActors();
 113   }
 114   if (mCurrentSlicerManager)
 115     mCurrentSlicerManager->Render();
 116 }
 117 //------------------------------------------------------------------------------
 118
 119
 120 //------------------------------------------------------------------------------
 121 bool vvToolBinarize::close() {
 122   RemoveVTKObjects();
 123   return vvToolWidgetBase::close();
 124 }
 125 //------------------------------------------------------------------------------
 126
 127
 128 //------------------------------------------------------------------------------
 129 void vvToolBinarize::reject() {
 130   RemoveVTKObjects();
 131   return vvToolWidgetBase::reject();
 132 }
 133 //------------------------------------------------------------------------------
 134
 135
 136 //------------------------------------------------------------------------------
 137 void vvToolBinarize::enableLowerThan(bool b) {
 138   if (!b) {
 139     mThresholdSlider1->resetMaximum();
 140     for(unsigned int i=0; i<mImageContour.size(); i++) {
 141       mImageContourLower[i]->hideActors();
 142     }
 143     mCurrentSlicerManager->Render();
 144   }
 145   else {
 146     valueChangedT1(mThresholdSlider1->GetValue());
 147     valueChangedT2(mThresholdSlider2->GetValue());
 148     for(unsigned int i=0; i<mImageContour.size(); i++) {
 149       mImageContourLower[i]->showActors();
 150     }
 151     mCurrentSlicerManager->Render();
 152   }
 153 }
 154 //------------------------------------------------------------------------------
 155
 156
 157 //------------------------------------------------------------------------------
 158 void vvToolBinarize::useFGBGtoggled(bool) {
 159   if (!mCheckBoxUseBG->isChecked() && !mCheckBoxUseFG->isChecked())
 160     mCheckBoxUseBG->toggle();
 161 }
 162 //------------------------------------------------------------------------------
 163
 164
 165 //------------------------------------------------------------------------------
 166 void vvToolBinarize::InputIsSelected(vvSlicerManager * m) {
 167   mCurrentSlicerManager = m;
 168
 169   // Specific for this gui
 170   mThresholdSlider1->SetValue(0);
 171   mThresholdSlider2->SetValue(0);
 172   mThresholdSlider1->SetImage(mCurrentImage);
 173   mThresholdSlider2->SetImage(mCurrentImage);
 174   mFGSlider->SetImage(mCurrentImage);
 175   mBGSlider->SetImage(mCurrentImage);
 176   //  DD(mCurrentSlicerManager->GetFileName().c_str());
 177   mFGSlider->SetMaximum(mCurrentImage->GetFirstVTKImageData()->GetScalarTypeMax());
 178   mFGSlider->SetMinimum(mCurrentImage->GetFirstVTKImageData()->GetScalarTypeMin());
 179   mBGSlider->SetMaximum(mCurrentImage->GetFirstVTKImageData()->GetScalarTypeMax());
 180   mBGSlider->SetMinimum(mCurrentImage->GetFirstVTKImageData()->GetScalarTypeMin());
 181   mFGSlider->SetValue(1);
 182   mBGSlider->SetValue(0);
 183   mFGSlider->SetSingleStep(1);
 184   mBGSlider->SetSingleStep(1);
 185
 186   // VTK objects for interactive display
 187   for(int i=0;i<mCurrentSlicerManager->NumberOfSlicers(); i++) {
 188     mImageContour.push_back(new vvImageContour);
 189     mImageContour[i]->setSlicer(mCurrentSlicerManager->GetSlicer(i));
 190     mImageContour[i]->setColor(1.0, 0.0, 0.0);
 191     mImageContourLower.push_back(new vvImageContour);
 192     mImageContourLower[i]->setSlicer(mCurrentSlicerManager->GetSlicer(i));
 193     mImageContourLower[i]->setColor(0.0, 0.0, 1.0);
 194   }
 195   valueChangedT1(mThresholdSlider1->GetValue());
 196
 197   connect(mThresholdSlider1, SIGNAL(valueChanged(double)), this, SLOT(valueChangedT1(double)));
 198   connect(mThresholdSlider2, SIGNAL(valueChanged(double)), this, SLOT(valueChangedT2(double)));
 199
 200   connect(mCurrentSlicerManager,SIGNAL(UpdateSlice(int,int)),this,SLOT(UpdateSlice(int, int)));
 201   connect(mCurrentSlicerManager,SIGNAL(UpdateTSlice(int,int)),this,SLOT(UpdateSlice(int, int)));
 202
 203 }
 204 //------------------------------------------------------------------------------
 205
 206
 207 //------------------------------------------------------------------------------
 208 void vvToolBinarize::UpdateSlice(int slicer,int slices) {
 209   if (!mInteractiveDisplayIsEnabled) return;
 210   if (!mCurrentSlicerManager) close();
 211   for(int i=0;i<mCurrentSlicerManager->NumberOfSlicers(); i++) {
 212     mImageContour[i]->update(mThresholdSlider1->GetValue());
 213     if (mRadioButtonLowerThan->isChecked())
 214       mImageContourLower[i]->update(mThresholdSlider2->GetValue());
 215
 216   }
 217   mCurrentSlicerManager->Render();
 218 }
 219 //------------------------------------------------------------------------------
 220
 221 //------------------------------------------------------------------------------
 222 void vvToolBinarize::GetArgsInfoFromGUI() {
 223
 224   /* //KEEP THIS FOR READING GGO FROM FILE
 225     int argc=1;
 226     std::string a = "toto";
 227     char * const* argv = new char*;
 228     //a.c_str();
 229     struct cmdline_parser_params p;
 230     p.check_required = 0;
 231     int good = cmdline_parser_ext(argc, argv, &args_info, &p);
 232   DD(good);
 233   */
 234
 235   mArgsInfo.imagetypes_flag = 0;
 236   mArgsInfo.upper_given = 0;
 237   mArgsInfo.lower_given = 0;
 238   bool inverseBGandFG = false;
 239
 240   mArgsInfo.lower_given = 1;
 241   mArgsInfo.lower_arg = mThresholdSlider1->GetValue();
 242   // DD(mArgsInfo.lower_arg);
 243   if (mRadioButtonLowerThan->isChecked()) {
 244     mArgsInfo.upper_given = 1;
 245     mArgsInfo.upper_arg = mThresholdSlider2->GetValue();
 246     if (mArgsInfo.upper_arg<mArgsInfo.lower_arg) {
 247       mArgsInfo.upper_given = 0;
 248       DD("TODO : lower thres greater than greater thres ! Ignoring ");
 249     }
 250   }
 251
 252   mArgsInfo.fg_arg = mFGSlider->GetValue();
 253   mArgsInfo.bg_arg = mBGSlider->GetValue();
 254
 255   if (inverseBGandFG) {
 256     mArgsInfo.fg_arg = mFGSlider->GetValue();
 257     mArgsInfo.bg_arg = mBGSlider->GetValue();
 258   }
 259   mArgsInfo.fg_given = 1;
 260   mArgsInfo.bg_given = 1;
 261
 262   if (mCheckBoxUseBG->isChecked()) {
 263     if (mCheckBoxUseFG->isChecked()) mArgsInfo.mode_arg = (char*)"both";
 264     else mArgsInfo.mode_arg = (char*)"BG";
 265   }
 266   else mArgsInfo.mode_arg = (char*)"FG";
 267
 268   mArgsInfo.verbose_flag = false;
 269
 270   // Required (even if not used)
 271   mArgsInfo.input_given = 0;
 272   mArgsInfo.output_given = 0;
 273
 274   mArgsInfo.input_arg = new char;
 275   mArgsInfo.output_arg = new char;
 276 }
 277 //------------------------------------------------------------------------------
 278
 279
 280 //------------------------------------------------------------------------------
 281 void vvToolBinarize::apply() {
 282   if (!mCurrentSlicerManager) close();
 283   QApplication::setOverrideCursor(QCursor(Qt::WaitCursor));
 284   GetArgsInfoFromGUI();
 285
 286   // Main filter
 287   clitk::BinarizeImageGenericFilter<args_info_clitkBinarizeImage>::Pointer filter =
 288      clitk::BinarizeImageGenericFilter<args_info_clitkBinarizeImage>::New();
 289   filter->SetArgsInfo(mArgsInfo);
 290   filter->SetInputVVImage(mCurrentImage);
 291   filter->Update();
 292
 293   // Output
 294   vvImage::Pointer output = filter->GetOutputVVImage();
 295   std::ostringstream osstream;
 296   osstream << "Binarized_" << mCurrentSlicerManager->GetSlicer(0)->GetFileName() << ".mhd";
 297   AddImage(output,osstream.str());
 298   QApplication::restoreOverrideCursor();
 299   close();
 300 }
 301 //------------------------------------------------------------------------------
 302
 303
 304 //------------------------------------------------------------------------------
 305 void vvToolBinarize::valueChangedT2(double v) {
 306   //  DD("valueChangedT2");
 307   if (mRadioButtonLowerThan->isChecked()) {
 308     mThresholdSlider1->SetMaximum(v);
 309     if (!mInteractiveDisplayIsEnabled) return;
 310     for(int i=0;i<mCurrentSlicerManager->NumberOfSlicers(); i++) {
 311       mImageContourLower[i]->update(v);
 312     }
 313     mCurrentSlicerManager->Render();
 314   }
 315 }
 316 //------------------------------------------------------------------------------
 317
 318
 319 //------------------------------------------------------------------------------
 320 void vvToolBinarize::valueChangedT1(double v) {
 321   //  DD("valueChangedT1");
 322   if (!mCurrentSlicerManager) close();
 323   mThresholdSlider2->SetMinimum(v);
 324   //  int m1 = (int)lrint(v);
 325   if (!mInteractiveDisplayIsEnabled) return;
 326   for(int i=0;i<mCurrentSlicerManager->NumberOfSlicers(); i++) {
 327     mImageContour[i]->update(v);
 328   }
 329   mCurrentSlicerManager->Render();
 330 }
 331 //------------------------------------------------------------------------------