...

[FrontAlgorithms.git] / lib / fpa / Filters / Image / Mori.hxx

diff --git a/lib/fpa/Filters/Image/Mori.hxx b/lib/fpa/Filters/Image/Mori.hxx
new file mode 100644 (file)
index 0000000..8b44364
--- /dev/null
+++ b/lib/fpa/Filters/Image/Mori.hxx
@@ -0,0 +1,46 @@
+// =========================================================================
+// @author Leonardo Florez Valencia
+// @email florez-l@javeriana.edu.co
+// =========================================================================
+#ifndef __fpa__Filters__Image__Mori__hxx__
+#define __fpa__Filters__Image__Mori__hxx__
+
+#include <itkBinaryThresholdImageFilter.h>
+
+// -------------------------------------------------------------------------
+template< class _TInputImage, class _TOutputImage, class _TTraits >
+fpa::Filters::Image::Mori< _TInputImage, _TOutputImage, _TTraits >::
+Mori( )
+  : Superclass( )
+{
+}
+
+// -------------------------------------------------------------------------
+template< class _TInputImage, class _TOutputImage, class _TTraits >
+fpa::Filters::Image::Mori< _TInputImage, _TOutputImage, _TTraits >::
+~Mori( )
+{
+}
+
+// -------------------------------------------------------------------------
+template< class _TInputImage, class _TOutputImage, class _TTraits >
+void fpa::Filters::Image::Mori< _TInputImage, _TOutputImage, _TTraits >::
+_AfterGenerateData( )
+{
+  this->Superclass::_AfterGenerateData( );
+
+  // Extract correct mark
+  typedef typename Superclass::TMarksImage _TMarks;
+  typedef itk::BinaryThresholdImageFilter< _TMarks, TOutputImage > _TFilter;
+  typename _TFilter::Pointer filter = _TFilter::New( );
+  filter->SetInput( this->GetMarks( ) );
+  filter->SetLowerThreshold( 1 );
+  filter->SetUpperThreshold( this->m_PeakDetector.GetNumberOfSamples( ) - 2 );
+  filter->SetInsideValue( this->GetInsideValue( ) );
+  filter->SetOutsideValue( this->GetOutsideValue( ) );
+  filter->Update( );
+  this->GraftOutput( filter->GetOutput( ) );
+}
+
+#endif // __fpa__Filters__Image__Mori__hxx__
+// eof - $RCSfile$