CTK/Libs/Visualization/VTK/Widgets/ctkVTKThresholdWidget.cpp

/*=========================================================================

  Library:   CTK

  Copyright (c) Kitware Inc.

  Licensed under the Apache License, Version 2.0 (the "License");
  you may not use this file except in compliance with the License.
  You may obtain a copy of the License at

      http://www.commontk.org/LICENSE

  Unless required by applicable law or agreed to in writing, software
  distributed under the License is distributed on an "AS IS" BASIS,
  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
  See the License for the specific language governing permissions and
  limitations under the License.

=========================================================================*/

// Qt includes
#include <QDebug>

// CTK includes
#include "ctkLogger.h"
#include "ctkVTKThresholdWidget.h"
#include "ctkUtils.h"
#include "ui_ctkVTKThresholdWidget.h"

// VTK includes
#include <vtkPiecewiseFunction.h>

//----------------------------------------------------------------------------
static ctkLogger logger("org.commontk.visualization.vtk.widgets.ctkVTKThresholdWidget");
//----------------------------------------------------------------------------

class ctkVTKThresholdWidgetPrivate:
  public Ui_ctkVTKThresholdWidget
{
   Q_DECLARE_PUBLIC(ctkVTKThresholdWidget);
protected:
  ctkVTKThresholdWidget* const q_ptr;
public:
  ctkVTKThresholdWidgetPrivate(ctkVTKThresholdWidget& object);
  void setupUi(QWidget* widget);
  void setThreshold(double min, double max, double opacity);
  void guessThreshold(double& min, double& max, double& opacity);
  
  vtkPiecewiseFunction* PiecewiseFunction;
};

// ----------------------------------------------------------------------------
// ctkVTKThresholdWidgetPrivate methods

// ----------------------------------------------------------------------------
ctkVTKThresholdWidgetPrivate::ctkVTKThresholdWidgetPrivate(
  ctkVTKThresholdWidget& object)
  : q_ptr(&object)
{
  this->PiecewiseFunction = 0;
}

// ----------------------------------------------------------------------------
void ctkVTKThresholdWidgetPrivate::setupUi(QWidget* widget)
{
  Q_Q(ctkVTKThresholdWidget);
  Q_ASSERT(q == widget);
  this->Ui_ctkVTKThresholdWidget::setupUi(widget);

  QObject::connect(this->ThresholdSliderWidget, SIGNAL(valuesChanged(double, double)),
                   q, SLOT(setThresholdValues(double, double)));
  QObject::connect(this->OpacitySliderWidget, SIGNAL(valueChanged(double)),
                   q, SLOT(setOpacity(double)));
}

// ----------------------------------------------------------------------------
void ctkVTKThresholdWidgetPrivate::guessThreshold(double& min, double& max, double& opacity)
{
  min = this->ThresholdSliderWidget->minimum();
  max = this->ThresholdSliderWidget->maximum();
  opacity = 1.;
  if (!this->PiecewiseFunction || this->PiecewiseFunction->GetSize() == 0)
    {
    return;
    }
  int minIndex; 
  for (minIndex=0; minIndex < this->PiecewiseFunction->GetSize(); ++minIndex)
    {
    double node[4];
    this->PiecewiseFunction->GetNodeValue(minIndex, node);
    if (node[1] > 0.)
      {
      min = node[0];
      opacity = node[1];
      break;
      }
    }
  if (minIndex == this->PiecewiseFunction->GetSize())
    {
    return;
    }
  int maxIndex;
  for (maxIndex = minIndex + 1;
       maxIndex < this->PiecewiseFunction->GetSize();
       ++maxIndex)
    {
    double node[4];
    this->PiecewiseFunction->GetNodeValue(maxIndex, node);
    // average the opacities
    opacity += node[1];
    if (node[1] == 0.)
      {
      max = node[0];
      opacity /= maxIndex - minIndex + 1;
      break;
      }
    }
  // couldn't find the upper threshold value, use the upper range
  if (maxIndex == this->PiecewiseFunction->GetSize())
    {
    opacity /= maxIndex - minIndex;
    }
}

// ----------------------------------------------------------------------------
void ctkVTKThresholdWidgetPrivate::setThreshold(double min, double max, double opacity)
{
  if (!this->PiecewiseFunction)
    {
    return;
    }
  double range[2];
  this->ThresholdSliderWidget->range(range);
 
  double node[4];
  node[0] = range[0];
  node[1] = 0.;
  node[2] = 0.5;
  node[3] = 0.;
  if (this->PiecewiseFunction->GetSize() < 1)
    {
    this->PiecewiseFunction->AddPoint(node[0], node[1], node[2], node[3]);
    }
  else
    {
    this->PiecewiseFunction->SetNodeValue(0, node);
    }

  node[0] = min;
  node[1] = 0.;
  node[2] = 0.;
  node[3] = 1.;
  if (this->PiecewiseFunction->GetSize() < 2)
    {
    this->PiecewiseFunction->AddPoint(node[0], node[1], node[2], node[3]);
    }
  else
    {
    this->PiecewiseFunction->SetNodeValue(1, node);
    }
  node[0] = max;
  node[1] = opacity;
  if (this->PiecewiseFunction->GetSize() < 3)
    {
    this->PiecewiseFunction->AddPoint(node[0], node[1], node[2], node[3]);
    }
  else
    {
    this->PiecewiseFunction->SetNodeValue(2, node);
    }
  node[0] = (max == range[1]) ? range[1] + 0.00000000000001 : range[1];
  node[1] = 0.;
  node[2] = 0.5;
  node[3] = 0.;
  if (this->PiecewiseFunction->GetSize() < 4)
    {
    this->PiecewiseFunction->AddPoint(node[0], node[1], node[2], node[3]);
    }
  else
    {
    this->PiecewiseFunction->SetNodeValue(3, node);
    }

}

// ----------------------------------------------------------------------------
// ctkVTKThresholdWidget methods

// ----------------------------------------------------------------------------
ctkVTKThresholdWidget::ctkVTKThresholdWidget(QWidget* parentWidget)
  :QWidget(parentWidget)
   , d_ptr(new ctkVTKThresholdWidgetPrivate(*this))
{
  Q_D(ctkVTKThresholdWidget);
  d->setupUi(this);
}

// ----------------------------------------------------------------------------
ctkVTKThresholdWidget::~ctkVTKThresholdWidget()
{
}

// ----------------------------------------------------------------------------
vtkPiecewiseFunction* ctkVTKThresholdWidget::piecewiseFunction()const
{
  Q_D(const ctkVTKThresholdWidget);
  return d->PiecewiseFunction;
}

// ----------------------------------------------------------------------------
void ctkVTKThresholdWidget
::setPiecewiseFunction(vtkPiecewiseFunction* newFunction)
{
  Q_D(ctkVTKThresholdWidget);
  this->qvtkReconnect(d->PiecewiseFunction, newFunction, vtkCommand::ModifiedEvent,
                      this, SLOT(updateFromPiecewiseFunction()));
  d->PiecewiseFunction = newFunction;
  
  double range[2] = {0., 1.};
  if (d->PiecewiseFunction)
    {
    d->PiecewiseFunction->GetRange(range);
    }
  this->setRange(range[0], range[1]);
  if (d->PiecewiseFunction)
    {
    double min, max, value;
    d->guessThreshold(min, max, value);
    d->setThreshold(min, max, value);
    }
}

// ----------------------------------------------------------------------------
void ctkVTKThresholdWidget::updateFromPiecewiseFunction()
{
  Q_D(ctkVTKThresholdWidget);
  
  if (!d->PiecewiseFunction)
    {
    return;
    }
  double range[2];
  d->ThresholdSliderWidget->range(range);
  double minThreshold = range[0];
  double maxThreshold = range[1];
  double opacity = 1.;
  double node[4];
  if (d->PiecewiseFunction->GetSize() > 1)
    {
    d->PiecewiseFunction->GetNodeValue(1, node);
    minThreshold = node[0];
    }
  if (d->PiecewiseFunction->GetSize() > 2)
    {
    d->PiecewiseFunction->GetNodeValue(2, node);
    maxThreshold = node[0];
    opacity = node[1];
    }
  bool wasBlocking = d->ThresholdSliderWidget->blockSignals(true);
  d->ThresholdSliderWidget->setValues(minThreshold, maxThreshold);
  d->ThresholdSliderWidget->blockSignals(wasBlocking);
  d->OpacitySliderWidget->blockSignals(true);
  d->OpacitySliderWidget->setValue(opacity);
  d->OpacitySliderWidget->blockSignals(wasBlocking);
}

// ----------------------------------------------------------------------------
void ctkVTKThresholdWidget::setThresholdValues(double min, double max)
{
  Q_D(ctkVTKThresholdWidget);
  d->setThreshold(min, max, d->OpacitySliderWidget->value());
}

// ----------------------------------------------------------------------------
void ctkVTKThresholdWidget::setOpacity(double opacity)
{
  Q_D(ctkVTKThresholdWidget);
  d->setThreshold(d->ThresholdSliderWidget->minimumValue(),
                  d->ThresholdSliderWidget->maximumValue(), opacity);
}

// ----------------------------------------------------------------------------
void ctkVTKThresholdWidget::range(double* range)const
{
  Q_D(const ctkVTKThresholdWidget);
  d->ThresholdSliderWidget->range(range);
}

// ----------------------------------------------------------------------------
void ctkVTKThresholdWidget::setRange(double min, double max)
{
  Q_D(ctkVTKThresholdWidget);
  bool wasBlocking = d->ThresholdSliderWidget->blockSignals(true);
  int decimals = qMax(0, -ctk::orderOfMagnitude(max - min) + 2);
  d->ThresholdSliderWidget->setDecimals(decimals);
  d->ThresholdSliderWidget->setSingleStep(pow(10., -decimals));
  d->ThresholdSliderWidget->setRange(min, max);
  d->ThresholdSliderWidget->blockSignals(wasBlocking);
}