RBC.h

 
#pragma once
 
#if defined(RBC_RECURSES)
#error Recursive header files inclusion detected in RBC.h
#else // defined(RBC_RECURSES)
#define RBC_RECURSES
 
#if !defined RBC_h
#define RBC_h
 
#include <iostream>
#include <vector>
#include <map>
#include <chrono>
#include "RBC_vec.h"
 
namespace DGtal {
 
    template<typename TSpace, typename TInputValue = typename TSpace::RealPoint, typename TOutputValue = typename TSpace::Point>
    struct RBC {
        RBC_vec<TSpace,TInputValue,TOutputValue> rot;
        typedef typename RBC_vec<TSpace,TInputValue,TOutputValue>::Circle Circle;
 
        RBC( const RBC_vec<TSpace,TInputValue,TOutputValue>& aRot, const double angle, const TOutputValue center )
          : rot( aRot ),my_angle(angle),my_center(center) {}
 
        template<typename TImage>
        TImage rotateImage( const TImage& img) const
        {
            typedef typename TImage::Domain TDomain;
            typedef DGtal::functors::DomainRigidTransformation2D < typename TImage::Domain, DGtal::RBC_vec<TSpace>> MyDomainTransformer;
            typedef std::pair < typename TSpace::Point, typename TSpace::Point > Bounds;
 
            typename TSpace::Point bottomLeft(-1,-1);
            typename TSpace::Point topRight(1,1);
 
            MyDomainTransformer domainTransformer ( rot );
            Bounds bounds = domainTransformer ( img.domain() );
            TDomain transformedDomain ( bounds.first+bottomLeft, bounds.second+topRight );
            TImage rotatedImage ( transformedDomain );
 
 
 
            // for ( auto r = 1; r < rot.size(); r++ )
            //      rotateCircle( img, rotatedImage, my_center, my_angle, r );
 
            for (typename TDomain::ConstIterator it = img.domain().begin(); it != img.domain().end(); ++it ) {
                rotatedImage.setValue((*this).rot.operator()(*it),img(*it));
            }
            return rotatedImage;
        }
 
        TOutputValue operator()( const TInputValue & aInput ) const
        {
            return (*this).rot.operator()(aInput);
        }
 
        std::string tostring() const {
            return {"RBC"};
        }
 
        void set_angle(const double newAngle) {
            my_angle=newAngle;
            rot.setAngle()= newAngle;
        }
 
        inline TOutputValue  center() const{return my_center;}
 
 
        double my_angle;
        TOutputValue my_center;
 
    };
 
 
}
 
#endif //RBC
 
#undef RBC_RECURSES
#endif // else defined(RBC_RECURSES)