example_ImageGaussianModelEstimator
example_ReadQuadEdgeMeshWithoutPlugins
example_RenderQuadEdgeMeshWithoutPlugins
+ example_ParallelImageMean
+ example_LightCompensation
)
FOREACH(prog ${NOPLUGINS_EXAMPLES_PROGRAMS})
ADD_EXECUTABLE(
--- /dev/null
+#include <cstdlib>
+#include <iostream>
+#include <string>
+
+#include <itkImage.h>
+#include <itkImageFileReader.h>
+#include <itkRGBPixel.h>
+
+#include <cpPlugins/Extensions/Algorithms/LightCompensationFilter.h>
+
+// -------------------------------------------------------------------------
+const unsigned int Dim = 2;
+typedef itk::RGBPixel< unsigned char > TPixel;
+typedef itk::Image< TPixel, Dim > TImage;
+
+// -------------------------------------------------------------------------
+int main( int argc, char* argv[] )
+{
+ if( argc < 3 )
+ {
+ std::cerr
+ << "Usage: " << argv[ 0 ]
+ << " input_image output_image"
+ << std::endl;
+ return( 1 );
+
+ } // fi
+ std::string input_image_fn = argv[ 1 ];
+ std::string output_image_fn = argv[ 2 ];
+
+ // Read image
+ itk::ImageFileReader< TImage >::Pointer input_image_reader =
+ itk::ImageFileReader< TImage >::New( );
+ input_image_reader->SetFileName( input_image_fn );
+ try
+ {
+ input_image_reader->Update( );
+ }
+ catch( itk::ExceptionObject& err )
+ {
+ std::cerr << "Error caught: " << err << std::endl;
+ return( 1 );
+
+ } // yrt
+ TImage::Pointer input_image = input_image_reader->GetOutput( );
+
+ cpPlugins::Extensions::Algorithms::
+ LightCompensationFilter< TImage >::Pointer filter =
+ cpPlugins::Extensions::Algorithms::
+ LightCompensationFilter< TImage >::New( );
+ filter->SetInput( input_image );
+ filter->Update( );
+
+ return( 0 );
+}
+
+// eof - $RCSfile$
--- /dev/null
+#include <cstdlib>
+#include <iostream>
+#include <string>
+
+#include <itkImage.h>
+#include <itkImageFileReader.h>
+#include <itkRGBPixel.h>
+
+#include <cpPlugins/Extensions/Algorithms/ParallelImageMean.h>
+
+// -------------------------------------------------------------------------
+const unsigned int Dim = 2;
+typedef itk::RGBPixel< unsigned char > TPixel;
+typedef itk::Image< TPixel, Dim > TImage;
+
+// -------------------------------------------------------------------------
+int main( int argc, char* argv[] )
+{
+ if( argc < 2 )
+ {
+ std::cerr << "Usage: " << argv[ 0 ] << " image" << std::endl;
+ return( 1 );
+
+ } // fi
+ std::string image_fn = argv[ 1 ];
+
+ // Read image
+ itk::ImageFileReader< TImage >::Pointer image_reader =
+ itk::ImageFileReader< TImage >::New( );
+ image_reader->SetFileName( image_fn );
+ try
+ {
+ image_reader->Update( );
+ }
+ catch( itk::ExceptionObject& err )
+ {
+ std::cerr << "Error caught: " << err << std::endl;
+ return( 1 );
+
+ } // yrt
+ TImage::Pointer image = image_reader->GetOutput( );
+
+ // Compute mean
+ cpPlugins::Extensions::Algorithms::ParallelImageMean< TImage >::Pointer mean =
+ cpPlugins::Extensions::Algorithms::ParallelImageMean< TImage >::New( );
+ mean->Execute( image, image->GetRequestedRegion( ) );
+ std::cout << "Mean: " << mean->GetMean( ) << std::endl;
+
+ return( 0 );
+}
+
+// eof - $RCSfile$
--- /dev/null
+// -------------------------------------------------------------------------
+// @author Leonardo Florez-Valencia (florez-l@javeriana.edu.co)
+// -------------------------------------------------------------------------
+
+#ifndef __CPPLUGINS__EXTENSIONS__ALGORITHMS__LIGHTCOMPENSATIONFILTER__H__
+#define __CPPLUGINS__EXTENSIONS__ALGORITHMS__LIGHTCOMPENSATIONFILTER__H__
+
+#include <itkInPlaceImageFilter.h>
+#include <cpPlugins/Extensions/Algorithms/ParallelImageMean.h>
+
+namespace cpPlugins
+{
+ namespace Extensions
+ {
+ namespace Algorithms
+ {
+ /**
+ */
+ template< class I >
+ class LightCompensationFilter
+ : public itk::InPlaceImageFilter< I, I >
+ {
+ public:
+ typedef LightCompensationFilter Self;
+ typedef itk::InPlaceImageFilter< I, I > Superclass;
+ typedef itk::SmartPointer< Self > Pointer;
+ typedef itk::SmartPointer< const Self > ConstPointer;
+
+ typedef I TImage;
+ typedef typename I::RegionType TRegion;
+
+ typedef cpPlugins::Extensions::Algorithms::ParallelImageMean< I > TMeanCalculator;
+ typedef typename TMeanCalculator::TMean TMean;
+
+ public:
+ itkNewMacro( Self );
+ itkTypeMacro( LightCompensationFilter, itkInPlaceImageFilter );
+
+ itkGetConstMacro( Mean, TMean );
+
+ protected:
+ LightCompensationFilter( );
+ virtual ~LightCompensationFilter( );
+
+ private:
+ virtual void BeforeThreadedGenerateData( );
+ virtual void ThreadedGenerateData(
+ const TRegion& region, itk::ThreadIdType id
+ );
+ virtual void AfterThreadedGenerateData( );
+
+ protected:
+ TMean m_Mean;
+ TMean m_Coefficient;
+ };
+
+ } // ecapseman
+
+ } // ecapseman
+
+} // ecapseman
+
+#include <cpPlugins/Extensions/Algorithms/LightCompensationFilter.hxx>
+
+#endif // __CPPLUGINS__EXTENSIONS__ALGORITHMS__LIGHTCOMPENSATIONFILTER__H__
+
+// eof - $RCSfile$
--- /dev/null
+// -------------------------------------------------------------------------
+// @author Leonardo Florez-Valencia (florez-l@javeriana.edu.co)
+// -------------------------------------------------------------------------
+
+#ifndef __CPPLUGINS__EXTENSIONS__ALGORITHMS__LIGHTCOMPENSATIONFILTER__HXX__
+#define __CPPLUGINS__EXTENSIONS__ALGORITHMS__LIGHTCOMPENSATIONFILTER__HXX__
+
+#include <itkImageRegionIterator.h>
+#include <itkImageRegionConstIterator.h>
+#include <itkNumericTraits.h>
+
+// -------------------------------------------------------------------------
+template< class I >
+cpPlugins::Extensions::Algorithms::LightCompensationFilter< I >::
+LightCompensationFilter( )
+ : Superclass( )
+{
+ this->SetNumberOfRequiredInputs( 1 );
+ this->InPlaceOff( );
+}
+
+// -------------------------------------------------------------------------
+template< class I >
+cpPlugins::Extensions::Algorithms::LightCompensationFilter< I >::
+~LightCompensationFilter( )
+{
+}
+
+// -------------------------------------------------------------------------
+template< class I >
+void cpPlugins::Extensions::Algorithms::LightCompensationFilter< I >::
+BeforeThreadedGenerateData( )
+{
+ I* input = const_cast< I* >( this->GetInput( ) );
+
+ typename TMeanCalculator::Pointer mc = TMeanCalculator::New( );
+ mc->Execute( input, input->GetRequestedRegion( ) );
+ this->m_Mean = mc->GetMean( );
+
+ double gm = this->m_Mean[ 0 ] + this->m_Mean[ 1 ] + this->m_Mean[ 2 ];
+ gm /= double( 3 );
+ this->m_Coefficient = this->m_Mean;
+ this->m_Coefficient.Fill( gm );
+ this->m_Coefficient[ 0 ] /= this->m_Mean[ 0 ];
+ this->m_Coefficient[ 1 ] /= this->m_Mean[ 1 ];
+ this->m_Coefficient[ 2 ] /= this->m_Mean[ 2 ];
+}
+
+// -------------------------------------------------------------------------
+template< class I >
+void cpPlugins::Extensions::Algorithms::LightCompensationFilter< I >::
+ThreadedGenerateData( const TRegion& region, itk::ThreadIdType id )
+{
+ typedef itk::NumericTraits< typename I::PixelType > _TPixelTraits;
+ typedef typename _TPixelTraits::ValueType _TPixelChannel;
+ typedef itk::NumericTraits< _TPixelChannel > _TChannelTraits;
+ static const double max_value = double( _TChannelTraits::max( ) );
+
+ itk::ImageRegionConstIterator< I > iIt( this->GetInput( ), region );
+ itk::ImageRegionIterator< I > oIt( this->GetOutput( ), region );
+
+ iIt.GoToBegin( );
+ oIt.GoToBegin( );
+ for( ; !iIt.IsAtEnd( ); ++iIt, ++oIt )
+ {
+ double r = double( iIt.Get( )[ 0 ] ) * this->m_Coefficient[ 0 ];
+ double g = double( iIt.Get( )[ 1 ] ) * this->m_Coefficient[ 1 ];
+ double b = double( iIt.Get( )[ 2 ] ) * this->m_Coefficient[ 2 ];
+ if( r > max_value ) r = max_value;
+ if( g > max_value ) g = max_value;
+ if( b > max_value ) b = max_value;
+ oIt.Get( ).Set(
+ _TPixelChannel( r ), _TPixelChannel( g ), _TPixelChannel( b )
+ );
+
+ } // rof
+}
+
+// -------------------------------------------------------------------------
+template< class I >
+void cpPlugins::Extensions::Algorithms::LightCompensationFilter< I >::
+AfterThreadedGenerateData( )
+{
+}
+
+#endif // __CPPLUGINS__EXTENSIONS__ALGORITHMS__LIGHTCOMPENSATIONFILTER__HXX__
+
+// eof - $RCSfile$
--- /dev/null
+// -------------------------------------------------------------------------
+// @author Leonardo Florez-Valencia (florez-l@javeriana.edu.co)
+// -------------------------------------------------------------------------
+
+#ifndef __CPPLUGINS__EXTENSIONS__ALGORITHMS__PARALLELIMAGEMEAN__H__
+#define __CPPLUGINS__EXTENSIONS__ALGORITHMS__PARALLELIMAGEMEAN__H__
+
+#include <itkDomainThreader.h>
+#include <itkThreadedImageRegionPartitioner.h>
+#include <itkArray.h>
+#include <itkNumericTraits.h>
+
+namespace cpPlugins
+{
+ namespace Extensions
+ {
+ namespace Algorithms
+ {
+ /**
+ */
+ template< class I >
+ class ParallelImageMean
+ : public itk::DomainThreader< itk::ThreadedImageRegionPartitioner< I::ImageDimension >, I >
+ {
+ public:
+ // Standard ITK typedefs.
+ typedef itk::DomainThreader< itk::ThreadedImageRegionPartitioner< I::ImageDimension >, I > Superclass;
+ typedef ParallelImageMean Self;
+ typedef itk::SmartPointer< Self > Pointer;
+ typedef itk::SmartPointer< const Self > ConstPointer;
+
+ typedef typename Superclass::DomainType DomainType;
+
+ typedef itk::Array< double > TMean;
+
+ protected:
+ typedef itk::NumericTraits< typename I::PixelType > _TPixelTraits;
+
+
+ public:
+ itkNewMacro( Self );
+ itkTypeMacro( ParallelImageMean, itkDomainThreader );
+
+ itkGetConstMacro( Mean, TMean );
+
+ protected:
+ ParallelImageMean( );
+ virtual ~ParallelImageMean( );
+
+ private:
+ virtual void BeforeThreadedExecution( );
+ virtual void ThreadedExecution(
+ const DomainType& region, const itk::ThreadIdType id
+ );
+ virtual void AfterThreadedExecution( );
+
+ protected:
+ itk::Array< double > m_Mean;
+ unsigned long m_N;
+ };
+
+ } // ecapseman
+
+ } // ecapseman
+
+} // ecapseman
+
+#include <cpPlugins/Extensions/Algorithms/ParallelImageMean.hxx>
+
+#endif // __CPPLUGINS__EXTENSIONS__ALGORITHMS__PARALLELIMAGEMEAN__H__
+
+// eof - $RCSfile$
--- /dev/null
+// -------------------------------------------------------------------------
+// @author Leonardo Florez-Valencia (florez-l@javeriana.edu.co)
+// -------------------------------------------------------------------------
+
+#ifndef __CPPLUGINS__EXTENSIONS__ALGORITHMS__PARALLELIMAGEMEAN__HXX__
+#define __CPPLUGINS__EXTENSIONS__ALGORITHMS__PARALLELIMAGEMEAN__HXX__
+
+#include <itkImageRegionConstIterator.h>
+
+// -------------------------------------------------------------------------
+template< class I >
+cpPlugins::Extensions::Algorithms::ParallelImageMean< I >::
+ParallelImageMean( )
+ : Superclass( )
+{
+}
+
+// -------------------------------------------------------------------------
+template< class I >
+cpPlugins::Extensions::Algorithms::ParallelImageMean< I >::
+~ParallelImageMean( )
+{
+}
+
+// -------------------------------------------------------------------------
+template< class I >
+void cpPlugins::Extensions::Algorithms::ParallelImageMean< I >::
+BeforeThreadedExecution( )
+{
+ this->m_Mean.SetSize( _TPixelTraits::GetLength( ) );
+ this->m_Mean.Fill( double( 0 ) );
+ this->m_N = 0;
+}
+
+// -------------------------------------------------------------------------
+template< class I >
+void cpPlugins::Extensions::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 cpPlugins::Extensions::Algorithms::ParallelImageMean< I >::
+AfterThreadedExecution( )
+{
+ this->m_Mean /= double( this->m_N );
+}
+
+#endif // __CPPLUGINS__EXTENSIONS__ALGORITHMS__PARALLELIMAGEMEAN__HXX__
+
+// eof - $RCSfile$
git://git.creatis.insa-lyon.fr / cpPlugins.git / commitdiff