threshold_algorithm.h

Go to the documentation of this file.

//--------------------------------------------------
//--------------------------------------------------

#ifndef THRESHOLD_ALGORITHM_H__
#define THRESHOLD_ALGORITHM_H__

#include "../algorithm.h"
#include "threshold_configuration.h"

#include <itkImageRegionIterator.h>


namespace imedgine
{

//--------------------------------------------------
//--------------------------------------------------

  class ThresholdAlgorithm: public Algorithm
  {
  public:

    //--------------------------------------------------

    ThresholdAlgorithm();


    //--------------------------------------------------

    virtual ~ThresholdAlgorithm();

    //--------------------------------------------------

    virtual void run();
    
  protected:

    //--------------------------------------------------

    ThresholdAlgorithm(ThresholdAlgorithm const& src);
    
    //--------------------------------------------------

    template <class DatasetPointerType>
    void applyAlgorithm(DatasetPointerType dataset); 
    
    //--------------------------------------------------
    
    imedgine::dataset_key_type dataset_key_;

  private:
    
  //----------------------------------------

    ThresholdAlgorithm& operator = (ThresholdAlgorithm const& src);

  }; // class ThresholdAlgorithm
  
  //--------------------------------------------------
  
  template <class DatasetPointerType>
      void ThresholdAlgorithm::applyAlgorithm(DatasetPointerType dataset)
  {
    ThresholdConfiguration* config = 
        dynamic_cast<ThresholdConfiguration*>(configuration_.get());
    
    if (config == 0)
    {
      return;
    }
    
    typedef typename DatasetPointerType::ObjectType DatasetType;
    typedef typename DatasetType::PixelType pixel_value_type;
    
    pixel_value_type lower_bound = 
        static_cast<pixel_value_type>(config->getLowerBound());
    pixel_value_type upper_bound = 
        static_cast<pixel_value_type>(config->getUpperBound());
    pixel_value_type default_pixel_value = 
        static_cast<pixel_value_type>(config->getDefaultPixelValue());
    
    std::cout << "lower bound: " << (unsigned int)lower_bound << std::endl
        << "upper bound: " << (unsigned int)upper_bound << std::endl;
    
    typedef itk::ImageRegionIterator<DatasetType> iterator_type;
    
    iterator_type dataset_iter(
        dataset, dataset->GetLargestPossibleRegion());
    
    for (dataset_iter.GoToBegin(); !dataset_iter.IsAtEnd(); ++dataset_iter)
    {
      if ((dataset_iter.Value() > upper_bound) ||
            (dataset_iter.Value() < lower_bound))
      {
        dataset_iter.Value() = default_pixel_value;
      }  
    }
    
    emit datasetChangedSignal(dataset_key_);
  }

} // namespace imedgine
  
#endif // THRESHOLD_ALGORITHM_H__

---