better manage resizing roi tool

[clitk.git] / vv / vvToolResample.cxx

diff --git a/vv/vvToolResample.cxx b/vv/vvToolResample.cxx
index 531a03b879098dcd7566f68008ed2ea135b98c89..03ed9cb32f5a8da11483442a441b7e425302c88d 100644 (file)
--- a/vv/vvToolResample.cxx
+++ b/vv/vvToolResample.cxx
@@ -1,7 +1,7 @@
 /*=========================================================================
   Program:   vv                     http://www.creatis.insa-lyon.fr/rio/vv
 
-  Authors belong to: 
+  Authors belong to:
   - University of LYON              http://www.universite-lyon.fr/
   - Léon Bérard cancer center       http://oncora1.lyon.fnclcc.fr
   - CREATIS CNRS laboratory         http://www.creatis.insa-lyon.fr
@@ -16,10 +16,12 @@
   - CeCILL-B   http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
   ======================================================================-====*/
 
+// vv
 #include "vvToolResample.h"
 #include "vvSlicer.h"
 #include "vvToolCreator.h"
 
+// qt
 #include <QFileInfo>
 #include <QMessageBox>
 
@@ -31,9 +33,11 @@ ADD_TOOL(vvToolResample);
 vvToolResample::vvToolResample(vvMainWindowBase * parent, Qt::WindowFlags f):
   vvToolWidgetBase(parent,f),
   vvToolBase<vvToolResample>(parent),
-  Ui::vvToolResample() {
+  Ui::vvToolResample()
+{
   Ui_vvToolResample::setupUi(mToolWidget);
 
+  //  mFilter = clitk::ResampleImageGenericFilter::New();
   mFilter = clitk::ImageResampleGenericFilter::New();
   mLastError ="";
 
@@ -115,7 +119,8 @@ vvToolResample::vvToolResample(vvMainWindowBase * parent, Qt::WindowFlags f):
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::Initialize() {
+void vvToolResample::Initialize()
+{
   SetToolName("Resample Image");
   SetToolMenuName("Resample");
   SetToolIconFilename(":/common/icons/resample.png");
@@ -125,7 +130,8 @@ void vvToolResample::Initialize() {
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::InputIsSelected(vvSlicerManager* m) {
+void vvToolResample::InputIsSelected(vvSlicerManager* m)
+{
 
   mCurrentSlicerManager = m;
   mCurrentImage = mCurrentSlicerManager->GetSlicer(0)->GetImage();
@@ -133,7 +139,7 @@ void vvToolResample::InputIsSelected(vvSlicerManager* m) {
   mInputFileName = mCurrentSlicerManager->GetFileName().c_str();
 
   // Set current information
-  mPixelType = mCurrentImage->GetScalarTypeAsString().c_str();
+  mPixelType = mCurrentImage->GetScalarTypeAsITKString().c_str();
   //ds      ComponentType = mCurrentImageGetNumberOfScalarComponents();
   mDimension = mCurrentImage->GetNumberOfDimensions();
 
@@ -141,7 +147,7 @@ void vvToolResample::InputIsSelected(vvSlicerManager* m) {
   mInputOrigin.resize(mDimension);
   mInputSpacing.resize(mDimension);
   mInputSize.resize(mDimension);
-  for (int i = 0; i < mDimension;i++) {
+  for (int i = 0; i < mDimension; i++) {
     mInputOrigin[i] = mCurrentImage->GetOrigin()[i];
     mInputSpacing[i] = mCurrentImage->GetSpacing()[i];
     mInputSize[i] = mCurrentImage->GetSize()[i];
@@ -171,7 +177,8 @@ void vvToolResample::InputIsSelected(vvSlicerManager* m) {
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::UpdateOutputInfo() {
+void vvToolResample::UpdateOutputInfo()
+{
   mOutputSizeLabel->setText(GetVectorIntAsString(mOutputSize));
   mOutputSpacingLabel->setText(GetVectorDoubleAsString(mOutputSpacing));
   mOutputMemoryLabel->setText(GetSizeInBytes(mOutputSize));
@@ -180,7 +187,8 @@ void vvToolResample::UpdateOutputInfo() {
 
 
 //------------------------------------------------------------------------------
-QString vvToolResample::GetSizeInBytes(std::vector<int> & size) {
+QString vvToolResample::GetSizeInBytes(std::vector<int> & size)
+{
   int t = 1;
   for (unsigned int i=0; i<size.size(); i++) t *= size[i];
   t *= mCurrentImage->GetScalarSize()*mCurrentImage->GetNumberOfScalarComponents();
@@ -190,25 +198,23 @@ QString vvToolResample::GetSizeInBytes(std::vector<int> & size) {
     t /= 1000000000;
     result += QString::number(t);
     result += " GB)";
-  }
-  else if (t > 1000000) {
+  } else if (t > 1000000) {
     t /= 1000000;
     result += QString::number(t);
     result += " MB)";
-  }
-  else if (t > 1000) {
+  } else if (t > 1000) {
     t /= 1000;
     result += QString::number(t);
     result += " KB)";
-  }
-  else result += ")";
+  } else result += ")";
   return result;
 }
 //------------------------------------------------------------------------------
 
 
 //------------------------------------------------------------------------------
-QString vvToolResample::GetVectorDoubleAsString(std::vector<double> vectorDouble) {
+QString vvToolResample::GetVectorDoubleAsString(std::vector<double> vectorDouble)
+{
   QString result;
   for (unsigned int i= 0; i<vectorDouble.size(); i++) {
     if (i != 0)
@@ -221,7 +227,8 @@ QString vvToolResample::GetVectorDoubleAsString(std::vector<double> vectorDouble
 
 
 //------------------------------------------------------------------------------
-QString vvToolResample::GetVectorIntAsString(std::vector<int> vectorInt) {
+QString vvToolResample::GetVectorIntAsString(std::vector<int> vectorInt)
+{
   QString result;
   for (unsigned int i= 0; i<vectorInt.size(); i++) {
     if (i != 0)
@@ -234,7 +241,8 @@ QString vvToolResample::GetVectorIntAsString(std::vector<int> vectorInt) {
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::FillSizeEdit(std::vector<int> size) {
+void vvToolResample::FillSizeEdit(std::vector<int> size)
+{
   xSizeLineEdit->setText(QString::number(size[0]));
   ySizeLineEdit->setText(QString::number(size[1]));
   if (size.size() > 2)
@@ -244,7 +252,8 @@ void vvToolResample::FillSizeEdit(std::vector<int> size) {
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::FillSpacingEdit(std::vector<double> spacing) {
+void vvToolResample::FillSpacingEdit(std::vector<double> spacing)
+{
   xSpacingLineEdit->setText(QString::number(spacing[0]));
   ySpacingLineEdit->setText(QString::number(spacing[1]));
   if (spacing.size() > 2)
@@ -254,7 +263,8 @@ void vvToolResample::FillSpacingEdit(std::vector<double> spacing) {
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::UpdateOutputSizeAndSpacing() {
+void vvToolResample::UpdateOutputSizeAndSpacing()
+{
   mOutputSize.resize(mDimension);
   mOutputSize = mInputSize;
   mOutputSpacing.resize(mDimension);
@@ -275,7 +285,8 @@ void vvToolResample::UpdateOutputSizeAndSpacing() {
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::UpdateControlSizeAndSpacing() {
+void vvToolResample::UpdateControlSizeAndSpacing()
+{
   scaleSizeLineEdit->setText("");
   scaleSpacingLineEdit->setText("");
   isoSizeLineEdit->setText("");
@@ -298,15 +309,13 @@ void vvToolResample::UpdateControlSizeAndSpacing() {
     ySizeLineEdit->setReadOnly(0);
     if (mDimension > 2)
       zSizeLineEdit->setReadOnly(0);
-  }
-  else {
+  } else {
     if (spacingRadioButton->isChecked()) {
       xSpacingLineEdit->setReadOnly(0);
       ySpacingLineEdit->setReadOnly(0);
       if (mDimension > 2)
         zSpacingLineEdit->setReadOnly(0);
-    }
-    else if (scaleSizeRadioButton->isChecked())
+    } else if (scaleSizeRadioButton->isChecked())
       scaleSizeLineEdit->setReadOnly(0);
     else if (scaleSpacingRadioButton->isChecked())
       scaleSpacingLineEdit->setReadOnly(0);
@@ -320,39 +329,40 @@ void vvToolResample::UpdateControlSizeAndSpacing() {
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::ComputeNewSpacingFromSize() {
+void vvToolResample::ComputeNewSpacingFromSize()
+{
   double newSpacing = mInputSpacing[0]*mInputSize[0];
   xSpacingLineEdit->setText(QString::number(newSpacing/xSizeLineEdit->text().toDouble()));
   newSpacing = mInputSpacing[1]*mInputSize[1];
   ySpacingLineEdit->setText(QString::number(newSpacing/ySizeLineEdit->text().toDouble()));
-  if (mDimension > 2)
-    {
-      newSpacing = mInputSpacing[2]*mInputSize[2];
-      zSpacingLineEdit->setText(QString::number(newSpacing/zSizeLineEdit->text().toDouble()));
-    }
+  if (mDimension > 2) {
+    newSpacing = mInputSpacing[2]*mInputSize[2];
+    zSpacingLineEdit->setText(QString::number(newSpacing/zSizeLineEdit->text().toDouble()));
+  }
   UpdateOutputSizeAndSpacing();
 }
 //------------------------------------------------------------------------------
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::ComputeNewSizeFromSpacing() {
+void vvToolResample::ComputeNewSizeFromSpacing()
+{
   double newSize = mInputSpacing[0]*mInputSize[0];
   xSizeLineEdit->setText(QString::number(newSize/xSpacingLineEdit->text().toDouble()));
   newSize = mInputSpacing[1]*mInputSize[1];
   ySizeLineEdit->setText(QString::number(newSize/ySpacingLineEdit->text().toDouble()));
-  if (mDimension > 2)
-    {
-      newSize = mInputSpacing[2]*mInputSize[2];
-      zSizeLineEdit->setText(QString::number(newSize/zSpacingLineEdit->text().toDouble()));
-    }
+  if (mDimension > 2) {
+    newSize = mInputSpacing[2]*mInputSize[2];
+    zSizeLineEdit->setText(QString::number(newSize/zSpacingLineEdit->text().toDouble()));
+  }
   UpdateOutputSizeAndSpacing();
 }
 //------------------------------------------------------------------------------
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::ComputeNewSpacingFromScale() {
+void vvToolResample::ComputeNewSpacingFromScale()
+{
   xSpacingLineEdit->setText(QString::number(mInputSpacing[0]*scaleSpacingLineEdit->text().toDouble()/100));
   ySpacingLineEdit->setText(QString::number(mInputSpacing[1]*scaleSpacingLineEdit->text().toDouble()/100));
   if (mDimension > 2)
@@ -363,7 +373,8 @@ void vvToolResample::ComputeNewSpacingFromScale() {
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::ComputeNewSizeFromScale() {
+void vvToolResample::ComputeNewSizeFromScale()
+{
   xSizeLineEdit->setText(QString::number(mInputSize[0]*scaleSizeLineEdit->text().toDouble()/100));
   ySizeLineEdit->setText(QString::number(mInputSize[1]*scaleSizeLineEdit->text().toDouble()/100));
   if (mDimension > 2)
@@ -374,7 +385,8 @@ void vvToolResample::ComputeNewSizeFromScale() {
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::ComputeNewSpacingFromIso() {
+void vvToolResample::ComputeNewSpacingFromIso()
+{
   xSpacingLineEdit->setText(QString::number(isoSpacingLineEdit->text().toDouble()));
   ySpacingLineEdit->setText(QString::number(isoSpacingLineEdit->text().toDouble()));
   if (mDimension > 2)
@@ -385,7 +397,8 @@ void vvToolResample::ComputeNewSpacingFromIso() {
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::ComputeNewSizeFromIso() {
+void vvToolResample::ComputeNewSizeFromIso()
+{
   xSizeLineEdit->setText(QString::number(isoSizeLineEdit->text().toDouble()));
   ySizeLineEdit->setText(QString::number(isoSizeLineEdit->text().toDouble()));
   if (mDimension > 2)
@@ -396,20 +409,19 @@ void vvToolResample::ComputeNewSizeFromIso() {
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::UpdateInterpolation() {
+void vvToolResample::UpdateInterpolation()
+{
   if (interpolationComboBox->currentText() == "BSpline") {
     bSplineLabel->show();
     bSplineOrderSpinBox->show();
     bLUTFactorLabel->hide();
     bLUTSpinBox->hide();
-  }
-  else if (interpolationComboBox->currentText() == "B-LUT (faster BSpline)")   {
+  } else if (interpolationComboBox->currentText() == "Blut (faster BSpline)")   {
     bSplineLabel->show();
     bSplineOrderSpinBox->show();
     bLUTFactorLabel->show();
     bLUTSpinBox->show();
-  }
-  else {
+  } else {
     bSplineLabel->hide();
     bSplineOrderSpinBox->hide();
     bLUTFactorLabel->hide();
@@ -420,15 +432,15 @@ void vvToolResample::UpdateInterpolation() {
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::UpdateGaussianFilter() {
+void vvToolResample::UpdateGaussianFilter()
+{
   if (gaussianFilterCheckBox->isChecked()) {
     gaussianFilterLabel->show();
     xGaussianLineEdit->show();
     yGaussianLineEdit->show();
     if (mDimension > 2)
       zGaussianLineEdit->show();
-  }
-  else {
+  } else {
     gaussianFilterLabel->hide();
     xGaussianLineEdit->hide();
     yGaussianLineEdit->hide();
@@ -439,7 +451,8 @@ void vvToolResample::UpdateGaussianFilter() {
 
 
 //------------------------------------------------------------------------------
-void vvToolResample::apply() {
+void vvToolResample::apply()
+{
 
   // Get resampler options
   std::vector<double> sigma;
@@ -448,13 +461,22 @@ void vvToolResample::apply() {
   if (mDimension > 2) sigma.push_back(zGaussianLineEdit->text().toDouble());
   if (mDimension == 4) sigma.push_back(0.01); //FIXME Don't filter along the temporal direction
 
+  /*
+  // Build ArgsInfo
+  typename clitk::ResampleImageGenericFilter::ArgsInfoType mArgsInfo;
+  mArgsInfo.size_given = mDimension;
+  mArgsInfo.size_arg = mDimension;
+  for(int i=0; i<mDimension; i++) {
+    mArgsInfo.size = mOutputSize;
+  */
+
   mFilter->SetOutputSize(mOutputSize);
   mFilter->SetOutputSpacing(mOutputSpacing);
   mFilter->SetInterpolationName(interpolationComboBox->currentText().toLower().toStdString());
 
   if (interpolationComboBox->currentText() == "BSpline")
     mFilter->SetBSplineOrder(bSplineOrderSpinBox->value());
-  else if (interpolationComboBox->currentText() == "B-LUT (faster BSpline)") {
+  else if (interpolationComboBox->currentText() == "Blut (faster BSpline)") {
     mFilter->SetInterpolationName("blut");
     mFilter->SetBSplineOrder(bSplineOrderSpinBox->value());
     mFilter->SetBLUTSampling(bLUTSpinBox->value());
@@ -475,7 +497,8 @@ void vvToolResample::apply() {
 
 
 //------------------------------------------------------------------------------
-std::string vvToolResample::GetOutputFileName() {
+std::string vvToolResample::GetOutputFileName()
+{
   QFileInfo info(QString(mCurrentSlicerManager->GetFileName().c_str()));
   return (info.path().toStdString() + "/resampled_" + info.fileName().toStdString());
 }