imageUtils_test.cpp

Go to the documentation of this file.

/** \file imageUtils_test.cpp
 */
#include "../../catch2/catch.hpp"
 
#include <vector>
#include <Eigen/Dense>
 
#define MX_NO_ERROR_REPORTS
 
#include "../../../include/math/func/gaussian.hpp"
#include "../../../include/improc/imageUtils.hpp"
#include "../../../include/improc/eigenCube.hpp"
 
/** Scenario: centroiding Gaussians with center of light
  * 
  * Verify center of light calculation
  * 
  * \anchor tests_improc_imageUtils_imageCenterOfLight
  */
SCENARIO( "Verify center of light calculation", "[improc::imageCenterOfLight]" ) 
{
   GIVEN("a Gaussian")
   {
      WHEN("geometric center")
      {
         mx::improc::eigenImage<double> im;
         im.resize(64,64);
         
         mx::math::func::gaussian2D<double>(im.data(), im.rows(), im.cols(), 0., 1.0, 31.5, 31.5, 2);
         
         double x, y;
         mx::improc::imageCenterOfLight( x, y, im);
         
         REQUIRE(fabs(x-31.5) < 1e-8 );
         REQUIRE(fabs(y-31.5) < 1e-8 );
      }
      WHEN("geometric quarter")
      {
         mx::improc::eigenImage<double> im;
         im.resize(64,64);
         
         mx::math::func::gaussian2D<double>(im.data(), im.rows(), im.cols(), 0., 1.0, 15.5, 15.5, 2);
         
         double x, y;
         mx::improc::imageCenterOfLight( x, y, im);
         
         REQUIRE(fabs(x-15.5) < 1e-8 );
         REQUIRE(fabs(y-15.5) < 1e-8 );
      }
   }
}