ObjectInfo.cpp

Go to the documentation of this file.

/* 
 * @file ObjectInfo.cpp
 * @author Nikolaos Apostolakos <nikoapos@gmail.com>
 */

#include <SEImplementation/PythonConfig/ObjectInfo.h>
#include <SEImplementation/Plugin/PixelCentroid/PixelCentroid.h>
#include <SEImplementation/Plugin/IsophotalFlux/IsophotalFlux.h>
#include <SEImplementation/Plugin/ShapeParameters/ShapeParameters.h>
#include <SEImplementation/Plugin/AssocMode/AssocMode.h>

namespace SourceXtractor {

ObjectInfo::ObjectInfo() {
  emplace(std::make_pair("centroid_x", 0.));
  emplace(std::make_pair("centroid_y", 0.));
  emplace(std::make_pair("isophotal_flux", 0.));
  emplace(std::make_pair("radius", 0.));
  emplace(std::make_pair("angle", 0.));
  emplace(std::make_pair("aspect_ratio", 0.));
  emplace(std::make_pair("assoc_match", false));
  emplace(std::make_pair("assoc_size", 0.));
  for (int i=0; i<100; i++) {
    std::stringstream label;
    label << "assoc_value_" << i;
    emplace(std::make_pair(label.str(), 0.));
  }
}

ObjectInfo::ObjectInfo(const SourceInterface& source) {
  auto centroid = source.getProperty<PixelCentroid>();
  auto iso_flux = source.getProperty<IsophotalFlux>();
  auto shape = source.getProperty<ShapeParameters>();
  auto assoc = source.getProperty<AssocMode>();

  double aspect_guess = std::max<double>(shape.getEllipseB() / shape.getEllipseA(), 0.01);

  emplace(std::make_pair("centroid_x", centroid.getCentroidX() + 1.0));
  emplace(std::make_pair("centroid_y", centroid.getCentroidY() + 1.0));
  emplace(std::make_pair("isophotal_flux", std::max<double>(iso_flux.getFlux(), 0.0001)));
  emplace(std::make_pair("radius", std::max<double>(shape.getEllipseA() / 2.0, 0.01)));
  emplace(std::make_pair("angle", shape.getEllipseTheta()));
  emplace(std::make_pair("aspect_ratio", aspect_guess));

  emplace(std::make_pair("assoc_match", assoc.getMatch()));
  if (assoc.getMatch()) {
    emplace(std::make_pair("assoc_size", (double) assoc.getAssocValues().shape()[0]));
    int i=0;
    for (auto assoc_value : assoc.getAssocValues()) {
      std::stringstream label;
      label << "assoc_value_" << i;
      emplace(std::make_pair(label.str(), assoc_value));
      i++;
    }

  }
}

}