ErrorVectorField.h

 
#pragma once
 
#if defined(ERRORVECTORFIELD_RECURSES)
#error Recursive header files inclusion detected in ErrorVectorField.h
#else // defined(ERRORVECTORFIELD_RECURSES)
#define ERRORVECTORFIELD_RECURSES
 
#if !defined ERRORVECTORFIELD_h
#define ERRORVECTORFIELD_h
 
// Inclusions
#include <DGtal/images/RigidTransformation2D.h>
 
namespace DGtal{
    template<typename TSpace, typename TDomain, typename TBijectiveReflections, typename TInputValue = typename TSpace::RealPoint, typename TOutputValue = typename TSpace::Point>
    struct ErrorVectorField{
 
        typedef functors::ForwardRigidTransformation2D<TSpace,TInputValue,typename TSpace::RealPoint,functors::Identity> RealRotation;
        typedef std::vector<std::vector<typename TSpace::RealPoint>> VectorField;
 
    protected:
        TBijectiveReflections normalVectors;
        typename TSpace::Point my_center;
        RealRotation targetRotation;
        RealRotation originCenteredRotation;
 
    public:
        ErrorVectorField(const TBijectiveReflections& reflections,const double theta,const typename TSpace::Point center):normalVectors(reflections),targetRotation(center,theta,{0,0}),my_center(center),originCenteredRotation({0,0},theta,{0,0})
        {}
 
        VectorField getOutputVectorFieldFromContour(const TDomain& set2dContour, bool continuityVecField =false ){
            VectorField outVecField;
            // std::cout << "DSL output vector error : normalVectors "<<normalVectors.my_startingAngle<<std::endl;
            // std::cout << "DSL output vector error : rotation angle "<<normalVectors.my_angle<<std::endl;
            // std::cout << "input test point = "<<*(set2dContour.begin()+30)<<std::endl;
            // std::cout << "output double reflections = "<<normalVectors(*(set2dContour.begin()+30)-my_center)+my_center<<std::endl;
            // std::cout << "output rotation = "<<targetRotation(*(set2dContour.begin()+30))<<std::endl;
 
            for (typename TDomain::ConstIterator it = set2dContour.begin(); it != set2dContour.end(); ++it ) {
                std::vector<typename TSpace::RealPoint> pixelError;
                typename TSpace::Point p = *it;
 
                typename TSpace::Point preflections = normalVectors(p-my_center)+my_center;
                typename TSpace::RealPoint protation = targetRotation(p);
 
                // compute the error field
                typename TSpace::RealPoint error = protation - preflections ;
                pixelError.push_back(error);
 
                // compute the eucliean rotation of the neighbors of p
                if(continuityVecField) {
                    // for the 8-Neighbor, compute the rotation
                    for(int veci = -1 ; veci <2 ; ++veci) {
                        for(int vecj = -1 ; vecj<2 ; ++vecj) {
                            if(veci!=0 || vecj!=0){
                                typename TSpace::RealPoint vecij_rot = originCenteredRotation({static_cast<double>(veci),static_cast<double>(vecj)});
                                typename TSpace::RealPoint neigh_rot = protation+vecij_rot;
                                pixelError.push_back(neigh_rot-preflections);
                            }
 
                        }
                    }
 
                }
 
 
                outVecField.push_back(pixelError);
            }
 
            return outVecField;
        }
 
    };
}
 
 
#endif //ERRORVECTORFIELD
#undef ERRORVECTORFIELD_RECURSES
#endif // else defined(ERRORVECTORFIELD_RECURSES)