yet another refactoring

[cpPlugins.git] / lib / cpExtensions / Algorithms / ParallelImageMean.hxx

diff --git a/lib/cpExtensions/Algorithms/ParallelImageMean.hxx b/lib/cpExtensions/Algorithms/ParallelImageMean.hxx
new file mode 100644 (file)
index 0000000..b3e4716
--- /dev/null
+++ b/lib/cpExtensions/Algorithms/ParallelImageMean.hxx
@@ -0,0 +1,65 @@
+// -------------------------------------------------------------------------
+// @author Leonardo Florez-Valencia (florez-l@javeriana.edu.co)
+// -------------------------------------------------------------------------
+
+#ifndef __CPEXTENSIONS__ALGORITHMS__PARALLELIMAGEMEAN__HXX__
+#define __CPEXTENSIONS__ALGORITHMS__PARALLELIMAGEMEAN__HXX__
+
+#include <itkImageRegionConstIterator.h>
+
+// -------------------------------------------------------------------------
+template< class I >
+cpExtensions::Algorithms::ParallelImageMean< I >::
+ParallelImageMean( )
+  : Superclass( )
+{
+}
+
+// -------------------------------------------------------------------------
+template< class I >
+cpExtensions::Algorithms::ParallelImageMean< I >::
+~ParallelImageMean( )
+{
+}
+
+// -------------------------------------------------------------------------
+template< class I >
+void cpExtensions::Algorithms::ParallelImageMean< I >::
+BeforeThreadedExecution( )
+{
+  this->m_Mean.SetSize( _TPixelTraits::GetLength( ) );
+  this->m_Mean.Fill( double( 0 ) );
+  this->m_N = 0;
+}
+
+// -------------------------------------------------------------------------
+template< class I >
+void cpExtensions::Algorithms::ParallelImageMean< I >::
+ThreadedExecution( const DomainType& region, const itk::ThreadIdType id )
+{
+  itk::ImageRegionConstIterator< I > i( this->m_Associate, region );
+  for( i.GoToBegin( ); !i.IsAtEnd( ); ++i )
+  {
+    for( unsigned int d = 0; d < _TPixelTraits::GetLength( ); ++d )
+    {
+      this->m_Mean.SetElement(
+        d, this->m_Mean.GetElement( d ) + double( i.Get( )[ d ] )
+        );
+      this->m_N++;
+
+    } // rof
+
+  } // rof
+}
+
+// -------------------------------------------------------------------------
+template< class I >
+void cpExtensions::Algorithms::ParallelImageMean< I >::
+AfterThreadedExecution( )
+{
+  this->m_Mean /= double( this->m_N );
+}
+
+#endif // __CPEXTENSIONS__ALGORITHMS__PARALLELIMAGEMEAN__HXX__
+
+// eof - $RCSfile$