--- /dev/null
+// =========================================================================
+// @author Leonardo Florez Valencia
+// @email florez-l@javeriana.edu.co
+// =========================================================================
+#ifndef __fpa__Filters__Mori__h__
+#define __fpa__Filters__Mori__h__
+
+#include <deque>
+#include <set>
+#include <itkConceptChecking.h>
+#include <fpa/Common/PeakDetector.h>
+#include <fpa/Filters/Algorithm.h>
+#include <fpa/Functors/RegionGrow/BinaryThreshold.h>
+
+namespace fpa
+{
+ namespace Filters
+ {
+ /**
+ */
+ template< class _TTraits >
+ class Mori
+ : public fpa::Filters::Algorithm< _TTraits >
+ {
+ public:
+ typedef _TTraits TTraits;
+ fpaTraitsMacro( typename, TTraits );
+
+ typedef fpa::Filters::Algorithm< TTraits > Superclass;
+ typedef Mori Self;
+ typedef itk::SmartPointer< Self > Pointer;
+ typedef itk::SmartPointer< const Self > ConstPointer;
+
+ typedef std::deque< TNode > TQueue;
+ typedef std::set< TInputValue > TThresholds;
+ typedef fpa::Common::PeakDetector TPeakDetector;
+ typedef TPeakDetector::TPeak TPeak;
+ typedef fpa::Functors::RegionGrow::BinaryThreshold< TInputValue > TPredicate;
+
+ protected:
+ itkConceptMacro(
+ Check_TOutputValue,
+ ( itk::Concept::IsUnsignedInteger< TOutputValue > )
+ );
+
+ public:
+ itkTypeMacro( fpa::Filters::Image::Mori, fpa::Filters::Algorithm );
+
+ itkGetConstMacro( InsideValue, TOutputValue );
+ itkSetMacro( InsideValue, TOutputValue );
+
+ itkGetConstMacro( MinimumThreshold, TInputValue );
+ itkSetMacro( MinimumThreshold, TInputValue );
+
+ public:
+ TOutputValue GetOutsideValue( ) const;
+ void SetOutsideValue( const TOutputValue& v );
+
+ unsigned long GetSignalKernelSize( ) const;
+ void SetSignalKernelSize( unsigned long k );
+
+ double GetSignalThreshold( ) const;
+ void SetSignalThreshold( double t );
+
+ double GetSignalInfluence( ) const;
+ void SetSignalInfluence( double i );
+
+ void ClearThresholds( );
+ void AddThreshold( const TInputValue& thr );
+ void SetThresholds(
+ const TInputValue& lower,
+ const TInputValue& upper,
+ const TInputValue& delta
+ );
+ const TThresholds& GetThresholds( ) const;
+ unsigned long GetNumberOfEvaluatedThresholds( ) const;
+ TInputValue GetOptimumThreshold( ) const;
+ void GetSignalValues(
+ unsigned long i, double& x, double& y, TPeak& p
+ ) const;
+
+ protected:
+ Mori( );
+ virtual ~Mori( );
+
+ virtual void _BeforeGenerateData( );
+ virtual void _Reinitialize( );
+ virtual void _UpdateOutputValue( TNode& n ) override;
+ virtual void _QueueClear( ) override;
+ virtual TNode _QueuePop( ) override;
+ virtual void _QueuePush( const TNode& n ) override;
+ virtual unsigned long _QueueSize( ) const override;
+ virtual void _ComputeOutputValue( TNode& n ) override;
+
+ private:
+ // Purposely not implemented.
+ Mori( const Self& other );
+ Self& operator=( const Self& other );
+
+ protected:
+ typename TPredicate::Pointer m_Predicate;
+
+ TOutputValue m_InsideValue;
+ TInputValue m_MinimumThreshold;
+ TThresholds m_Thresholds;
+ typename TThresholds::const_iterator m_CurrThr;
+
+ TQueue m_Queues[ 2 ];
+ unsigned int m_CurrQueue;
+ double m_CurrCount;
+
+ TPeakDetector m_PeakDetector;
+ };
+
+ } // ecapseman
+
+} // ecapseman
+
+#ifndef ITK_MANUAL_INSTANTIATION
+# include <fpa/Filters/Mori.hxx>
+#endif // ITK_MANUAL_INSTANTIATION
+
+#endif // __fpa__Filters__Mori__h__
+// eof - $RCSfile$