...

[FrontAlgorithms.git] / lib / fpa / Image / Dijkstra.h

diff --git a/lib/fpa/Image/Dijkstra.h b/lib/fpa/Image/Dijkstra.h
index 37c28e808458eeed85c4713f5b8bf1a123fc36a7..7992842295bb400e78461cb18f2e1255047d3d5f 100644 (file)
--- a/lib/fpa/Image/Dijkstra.h
+++ b/lib/fpa/Image/Dijkstra.h
@@ -1,10 +1,18 @@
+// =========================================================================
+// @author Leonardo Florez Valencia
+// @email florez-l@javeriana.edu.co
+// =========================================================================
+
 #ifndef __fpa__Image__Dijkstra__h__
 #define __fpa__Image__Dijkstra__h__
 
 #include <fpa/Base/Dijkstra.h>
-#include <fpa/Image/Algorithm.h>
+#include <fpa/Base/SeedsInterface.h>
+#include <fpa/Image/MarksInterface.h>
+#include <fpa/Image/Filter.h>
 #include <fpa/Image/MinimumSpanningTree.h>
-#include <fpa/Image/Functors/Base.h>
+#include <fpa/Image/Functors/VertexParentBase.h>
+#include <fpa/Image/Functors/VertexIdentity.h>
 
 namespace fpa
 {
@@ -14,36 +22,70 @@ namespace fpa
      */
     template< class _TInputImage, class _TOutputImage >
     class Dijkstra
-      : public fpa::Base::Dijkstra< fpa::Image::Algorithm< _TInputImage, _TOutputImage >, fpa::Image::MinimumSpanningTree< _TInputImage::ImageDimension > >
+      : public fpa::Base::Dijkstra< fpa::Image::Filter< _TInputImage, _TOutputImage >, fpa::Image::MarksInterface< _TInputImage::ImageDimension >, fpa::Base::SeedsInterface< typename _TInputImage::IndexType, typename _TInputImage::IndexType::LexicographicCompare >, fpa::Image::MinimumSpanningTree< _TInputImage::ImageDimension > >
     {
     public:
-      typedef Dijkstra                                             Self;
-      typedef fpa::Image::Algorithm< _TInputImage, _TOutputImage > TAlgorithm;
+      // Interfaces
+      typedef fpa::Image::Filter< _TInputImage, _TOutputImage > TFilter;
+      typedef fpa::Image::MarksInterface< _TInputImage::ImageDimension > TMarksInterface;
+      typedef fpa::Base::SeedsInterface< typename _TInputImage::IndexType, typename _TInputImage::IndexType::LexicographicCompare > TSeedsInterface;
       typedef fpa::Image::MinimumSpanningTree< _TInputImage::ImageDimension > TMST;
-      typedef fpa::Base::Dijkstra< TAlgorithm, TMST > Superclass;
-      typedef itk::SmartPointer< Self >               Pointer;
-      typedef itk::SmartPointer< const Self >         ConstPointer;
-
-      typedef typename Superclass::TOutput TOutput;
-      typedef typename Superclass::TVertex TVertex;
 
-      typedef fpa::Image::Functors::Base< _TInputImage, fpa::Base::DijkstraCostFunctionBase< TVertex, TOutput > > TCostFunction;
+      // Smart pointers
+      typedef Dijkstra Self;
+      typedef fpa::Base::Dijkstra< TFilter, TMarksInterface, TSeedsInterface, TMST > Superclass;
+      typedef itk::SmartPointer< Self > Pointer;
+      typedef itk::SmartPointer< const Self > ConstPointer;
 
-    protected:
-      typedef typename Superclass::_TQueueNode _TQueueNode;
+      typedef typename TFilter::TInputImage TInputImage;
+      typedef typename TFilter::TOutputValue TOutputValue;
+      typedef typename TFilter::TVertex TVertex;
 
     public:
       itkNewMacro( Self );
       itkTypeMacro( fpa::Image::Dijkstra, fpa::Base::Dijkstra );
 
     protected:
-      Dijkstra( );
-      virtual ~Dijkstra( );
+      Dijkstra( )
+        : Superclass( )
+        {
+          this->SetFunctor(
+            fpa::Image::Functors::
+            VertexIdentity< _TInputImage, typename TFilter::TOutputValue >::
+            New( )
+            );
+        }
+      virtual ~Dijkstra( )        { }
+
+      virtual void _ConfigureOutputs( const TOutputValue& init_value ) override
+        {
+          this->Superclass::_ConfigureOutputs( init_value );
 
-      virtual void _BeforeGenerateData( ) fpa_OVERRIDE;
+          typename TVertex::OffsetType o;
+          o.Fill( 0 );
+          const TInputImage* input = this->GetInput( );
+          TMST* mst = this->GetMinimumSpanningTree( );
+          mst->CopyInformation( input );
+          mst->SetBufferedRegion( input->GetRequestedRegion( ) );
+          mst->Allocate( );
+          mst->FillBuffer( o );
+        }
+
+      virtual void GenerateData( ) override
+        {
+          // Configure functors with input image
+          typedef typename TFilter::TOutputValue _TOutputValue;
+          typedef fpa::Image::Functors::VertexParentBase< _TInputImage, _TOutputValue > _TVFunc;
+          _TVFunc* vfunc =
+            dynamic_cast< _TVFunc* >( this->m_VertexFunctor.GetPointer( ) );
+          if( vfunc != NULL )
+            vfunc->SetImage( this->GetInput( ) );
+
+          // Ok, continue
+          this->Superclass::GenerateData( );
+        }
 
     private:
-      // Purposely not defined
       Dijkstra( const Self& other );
       Self& operator=( const Self& other );
     };
@@ -52,10 +94,6 @@ namespace fpa
 
 } // ecapseman
 
-#ifndef ITK_MANUAL_INSTANTIATION
-#  include <fpa/Image/Dijkstra.hxx>
-#endif // ITK_MANUAL_INSTANTIATION
-
 #endif // __fpa__Image__Dijkstra__h__
 
 // eof - $RCSfile$