CBDR_naiverotation.h

 
#pragma once
 
#if defined(CBDRNAIVEROTATION_RECURSES)
#error Recursive header files inclusion detected in CBDR_naiverotation.h
#else // defined(CBDRNAIVEROTATION_RECURSES)
#define CBDRNAIVEROTATION_RECURSES
 
#if !defined CBDRNAIVEROTATION_h
#define CBDRNAIVEROTATION_h
 
// Inclusions
#include "DigitizedReflection.h"
#include "GAVector.h"
namespace DGtal {
    template<typename TSpace, typename TInputValue = typename TSpace::RealPoint, typename TOutputValue = typename TSpace::Point>
struct CBDR_naiverotation{
        typedef Reflection<TSpace,TInputValue> DigitizedReflection;
 
        std::vector<DigitizedReflection> bijectiveNormalVectors;
 
        CBDR_naiverotation(const std::vector<DigitizedReflection>& bijectiveReflections={} ):bijectiveNormalVectors(bijectiveReflections){    }
 
        explicit CBDR_naiverotation(const std::vector<GAVector<TSpace>>& bijectiveGAvec ){
             bijectiveNormalVectors.resize(bijectiveGAvec.size());
             std::transform(bijectiveGAvec.begin(), bijectiveGAvec.end(), bijectiveNormalVectors.begin(),
                [](const GAVector<TSpace>& p) { return DigitizedReflection(p); });
        }
 
        TOutputValue operator()( const TInputValue & aInput ) const
        {
            if(bijectiveNormalVectors.size()<=0){
                return (TOutputValue)aInput;
            }
            int i =0;
            DigitizedReflection firstReflection(bijectiveNormalVectors[i]);
            TOutputValue resultat = firstReflection(aInput);
 
            i+=1;
            for(; i < bijectiveNormalVectors.size(); ++i)
            {
                DigitizedReflection currentReflection(bijectiveNormalVectors[i]);
                resultat = currentReflection(resultat);
            }
            return resultat;
        }
 
        template<typename TImage>
        TImage rotateImage(TImage img) const {
            typedef typename TImage::Domain TDomain;
            typedef DGtal::functors::DomainRigidTransformation2D < typename TImage::Domain, DGtal::CBDR_naiverotation<TSpace>> MyDomainTransformer;
            typedef std::pair < typename TSpace::Point, typename TSpace::Point > Bounds;
 
            MyDomainTransformer domainTransformer ( *this );
            Bounds bounds = domainTransformer ( img.domain() );
            TDomain transformedDomain ( bounds.first, bounds.second );
            TImage rotatedImage ( transformedDomain );
 
            for (typename TDomain::ConstIterator it = img.domain().begin(); it != img.domain().end(); ++it )
            {
                rotatedImage.setValue((*this)(*it),img(*it));
            }
            return rotatedImage;
        }
 
    };
}
 
 
#endif //CBDRNAIVEROTATION
#undef CBDRNAIVEROTATION_RECURSES
#endif // else defined(CBDRNAIVEROTATION_RECURSES)