transformation_d.h

Go to the documentation of this file.

/* Copyright (C) 2000-2013 CEA
 *
 * This software and supporting documentation were developed by
 *     bioPICSEL
 *     CEA/DSV/I²BM/MIRCen/LMN, Batiment 61,
 *     18, route du Panorama
 *     92265 Fontenay-aux-Roses
 *     France
 */


#ifndef AIMS_REGISTRATION_TRANSFORMATION_D_H
#define AIMS_REGISTRATION_TRANSFORMATION_D_H


#include <aims/data/data_g.h>
#include <aims/resampling/motion.h>
#include <aims/registration/transformation.h>
#include <iostream>
#include <math.h>


template <class T>
AimsData<T> Transformation::application(const AimsData<T>& test,Point3df p)
{
  
  float angle_rad = p[2] /180.0 * M_PI;
  float c = cos( angle_rad );
  float s = sin( angle_rad );
  AimsData<T>    testtrans(test.dimX(), test.dimY());
  
  
  int x,y;
  ForEach2d(testtrans, x, y)
    {
      testtrans(x,y)      = T( 195 );

      if(((_cx+c*(x-_cx)+s*(y-_cy)-p[1])<test.dimX())&&((_cx+c*(x-_cx)+s*(y-_cy)-p[1])>=0)&&((_cy-s*(x-_cx)+c*(y-_cy)-p[0])<test.dimY())&&((_cy-s*(x-_cx)+c*(y-_cy)-p[0])>=0))
        testtrans(x,y)  = test( unsigned(_cx+c*(x-_cx)+s*(y-_cy)-p[1]), unsigned(_cy-s*(x-_cx)+c*(y-_cy)-p[0]) ) ;
    }

  return (testtrans);
} 


#endif