images/exampleRigidtransformation2d.cpp

Example of 2D rigid transformation using forward and backward model.

See also

Geometric transformations

Result for backward model

 

#include <iostream>
#include <cmath>
#include "DGtal/images/ImageSelector.h"
#include "DGtal/images/ImageContainerBySTLVector.h"
#include "DGtal/images/ConstImageAdapter.h"
#include "ConfigExamples.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/base/Common.h"
#include "DGtal/io/readers/PGMReader.h"
#include "DGtal/io/writers/GenericWriter.h"
#include "DGtal/images/RigidTransformation2D.h"
 
using namespace std;
using namespace DGtal;
using namespace functors;
using namespace Z2i;

 
int main( int , char** )
{
  typedef ImageSelector<Domain, unsigned char >::Type Image;
  typedef ForwardRigidTransformation2D < Space > ForwardTrans;
  typedef BackwardRigidTransformation2D < Space > BackwardTrans;
  typedef ConstImageAdapter<Image, Domain, BackwardTrans, Image::Value, Identity > MyImageBackwardAdapter;
  typedef DomainRigidTransformation2D < Domain, ForwardTrans > MyDomainTransformer;
  typedef MyDomainTransformer::Bounds Bounds;
  trace.beginBlock ( "Example rigidtransformation2d" );
    ForwardTrans forwardTrans( RealPoint ( 5, 5 ), M_PI_4, RealVector( 3, -3 ) );
    BackwardTrans backwardTrans( RealPoint ( 5, 5 ), M_PI_4, RealVector( 3, -3 ) );
    MyDomainTransformer domainTransformer ( forwardTrans );
    Identity idD;
 
    Image image = PGMReader<Image>::importPGM ( examplesPath + "samples/church.pgm" );
    Bounds bounds = domainTransformer ( image.domain() );
    Domain transformedDomain ( bounds.first, bounds.second );
    
    trace.beginBlock ( "Backward - Eulerian model" );
      MyImageBackwardAdapter backwardImageAdapter ( image, transformedDomain , backwardTrans, idD );
      backwardImageAdapter >> "backward_transform.pgm";
    trace.endBlock();
    
    trace.beginBlock( "Forward - Lagrangian model" );
      Image forwardTransformedImage ( transformedDomain );
      for ( Domain::ConstIterator it = image.domain().begin(); it != image.domain().end(); ++it )
      {
        forwardTransformedImage.setValue ( forwardTrans ( *it ), image ( *it ) );
      }
      forwardTransformedImage >> "forward_transform.pgm";
    trace.endBlock();
  trace.endBlock();
  return 0;
}
//                                                                           //