CBDR solver, use a policy to choose the composition of digitized reflections that minimises an error. More...

#include <DGtal/images/bijectiveRotations/CBDRSolver.h>

Public Types
typedef CBDR_naiverotation< TSpace, typename TSpace::RealPoint >	BijectiveReflections

typedef functors::ForwardRigidTransformation2D< TSpace, typename TSpace::RealPoint, typename TSpace::RealPoint, functors::Identity >	RealRotation

typedef std::vector< std::pair< std::vector< GAVector< TSpace > >, GAVector< TSpace > > >	bijectiveReflect

typedef ErrorVectorField< TSpace, TDomain, typename TSpace::RealPoint >	ErrorRealVectors

typedef Reflection< TSpace >	DigitizedReflection

typedef std::vector< std::vector< typename TSpace::RealPoint > >	VectorField

Public Member Functions
	CBDRSolver_GAvec (const size_t km, const double rotAngle, const typename TDomain::Point center)

std::pair< std::vector< GAVector< TSpace > >, double >	outputCompositionReflection (const TDomain &set2d, typename bijectiveReflect::iterator &lowerAngle, typename bijectiveReflect::iterator &upperAngle, const Policy< TSpace, TDomain, BijectiveReflections > &policy)

Data Fields
NBijectiveGenerator< TSpace, typename TSpace::RealPoint >	nBijectiveGenerator

double	my_angle

TDomain::Point	my_center

Detailed Description

template<typename TSpace, typename TDomain>
class DGtal::CBDRSolver_GAvec< TSpace, TDomain >

CBDR solver, use a policy to choose the composition of digitized reflections that minimises an error.

Description of template struct 'CBDR Solver'

Template Parameters

TSpace	a 2 dimensional space.
TDomain	a 2 dimensional domain.

Definition at line 51 of file CBDRSolver.h.

Member Typedef Documentation

◆ bijectiveReflect

template<typename TSpace , typename TDomain >

typedef std::vector<std::pair<std::vector<GAVector<TSpace> >,GAVector<TSpace> > > DGtal::CBDRSolver_GAvec< TSpace, TDomain >::bijectiveReflect

Definition at line 55 of file CBDRSolver.h.

◆ BijectiveReflections

template<typename TSpace , typename TDomain >

typedef CBDR_naiverotation<TSpace,typename TSpace::RealPoint> DGtal::CBDRSolver_GAvec< TSpace, TDomain >::BijectiveReflections

Definition at line 53 of file CBDRSolver.h.

◆ DigitizedReflection

template<typename TSpace , typename TDomain >

typedef Reflection<TSpace> DGtal::CBDRSolver_GAvec< TSpace, TDomain >::DigitizedReflection

Definition at line 57 of file CBDRSolver.h.

◆ ErrorRealVectors

template<typename TSpace , typename TDomain >

typedef ErrorVectorField<TSpace,TDomain,typename TSpace::RealPoint> DGtal::CBDRSolver_GAvec< TSpace, TDomain >::ErrorRealVectors

Definition at line 56 of file CBDRSolver.h.

◆ RealRotation

template<typename TSpace , typename TDomain >

typedef functors::ForwardRigidTransformation2D<TSpace,typename TSpace::RealPoint,typename TSpace::RealPoint,functors::Identity> DGtal::CBDRSolver_GAvec< TSpace, TDomain >::RealRotation

Definition at line 54 of file CBDRSolver.h.

◆ VectorField

template<typename TSpace , typename TDomain >

typedef std::vector<std::vector<typename TSpace::RealPoint> > DGtal::CBDRSolver_GAvec< TSpace, TDomain >::VectorField

Definition at line 58 of file CBDRSolver.h.

Constructor & Destructor Documentation

◆ CBDRSolver_GAvec()

template<typename TSpace , typename TDomain >

DGtal::CBDRSolver_GAvec< TSpace, TDomain >::CBDRSolver_GAvec	(	const size_t	km,
		const double	rotAngle,
		const typename TDomain::Point	center
	)

inline

Definition at line 61 of file CBDRSolver.h.

61 :nBijectiveGenerator(km),my_angle(rotAngle),my_center(center){}

DGtal::CBDRSolver_GAvec::my_center

TDomain::Point my_center

Definition: CBDRSolver.h:109

DGtal::CBDRSolver_GAvec::my_angle

double my_angle

Definition: CBDRSolver.h:108

DGtal::CBDRSolver_GAvec::nBijectiveGenerator

NBijectiveGenerator< TSpace, typename TSpace::RealPoint > nBijectiveGenerator

Definition: CBDRSolver.h:107

Member Function Documentation

◆ outputCompositionReflection()

template<typename TSpace , typename TDomain >

std::pair<std::vector<GAVector<TSpace> >,double> DGtal::CBDRSolver_GAvec< TSpace, TDomain >::outputCompositionReflection	(	const TDomain &	set2d,
		typename bijectiveReflect::iterator &	lowerAngle,
		typename bijectiveReflect::iterator &	upperAngle,
		const Policy< TSpace, TDomain, BijectiveReflections > &	policy
	)

inline

find the composition of bijective reflections that minimize either Linf, L2 or Lcontinuity of error vector field given by the policy

Returns: {vector that minimizes an error given by policy, the error}

Definition at line 66 of file CBDRSolver.h.

                                                                                                                       {
  
             // first error
             auto it = lowerAngle;
  
             typename bijectiveReflect::iterator itMinError;
  
             std::vector<DigitizedReflection> bijectiveReflections;
             std::vector<GAVector<TSpace>> firstReflectionsIndex=(*lowerAngle).first;
  
             for(size_t i = 0 ;i<firstReflectionsIndex.size();++i){
                 bijectiveReflections.push_back(DigitizedReflection(firstReflectionsIndex[i]));//firstReflectionsIndex.size()-1-
             }
  
             CBDR_naiverotation<TSpace> reflections(bijectiveReflections);
  
             // finally get the average error
             double minError = policy.evaluate(set2d, reflections,my_angle, my_center);  //outputError(firsterrors,policy); // should now be policy
  
             itMinError = lowerAngle;
  
             for(it = lowerAngle+1 ; it != upperAngle ; ++it ){
                 std::vector<DigitizedReflection> currentbijectiveReflections;
                 std::vector<GAVector<TSpace>> currentReflectionsIndex=(*it).first;
                 for(size_t i = 0 ;i<currentReflectionsIndex.size();++i){
                     currentbijectiveReflections.push_back(DigitizedReflection(currentReflectionsIndex[i]));//nBijectiveGenerator.BijectiveVectors[currentReflectionsIndex[i]] // currentReflectionsIndex.size()-1-
                 }
                 CBDR_naiverotation<TSpace> currentreflections(currentbijectiveReflections);
  
                 // // finally get the error
                 double error =  policy.evaluate(set2d, currentreflections,my_angle, my_center);
                 if(fabs(error)<fabs(minError)){
                     itMinError = it;
                     minError = error;
                 }
             }
             return {(*itMinError).first,minError};
         }