FluxMeasurement.h

Go to the documentation of this file.

/*
 * FluxMeasurement.h
 *
 *  Created on: Oct 19, 2018
 *      Author: Alejandro Alvarez
 */
 
#ifndef _SEFRAMEWORK_SEFRAMEWORK_APERTURE_MEASUREFLUX_H
#define _SEFRAMEWORK_SEFRAMEWORK_APERTURE_MEASUREFLUX_H
 
#include "Aperture.h"
#include "SEFramework/Image/WriteableImage.h"
#include "SEFramework/Source/SourceFlags.h"
 
namespace SourceXtractor {
 
struct FluxMeasurement {
  SeFloat m_flux;
  SeFloat m_variance;
  SeFloat m_total_area, m_bad_area;
  Flags m_flags;
 
  FluxMeasurement() : m_flux(0.), m_variance(0.), m_total_area(0.), m_bad_area(0.), m_flags(Flags::NONE) {}
};

FluxMeasurement measureFlux(const std::shared_ptr<Aperture> &aperture, SeFloat centroid_x, SeFloat centroid_y,
                            const std::shared_ptr<Image<SeFloat>> &img,
                            const std::shared_ptr<Image<SeFloat>> &variance_map, SeFloat variance_threshold,
                            bool use_symmetry);

template <typename T>
void fillAperture(const std::shared_ptr<Aperture> &aperture, SeFloat centroid_x, SeFloat centroid_y,
                  const std::shared_ptr<WriteableImage<T>> &img, T value) {
  auto min_pixel = aperture->getMinPixel(centroid_x, centroid_y);
  auto max_pixel = aperture->getMaxPixel(centroid_x, centroid_y);
 
  for (int y = min_pixel.m_y; y <= max_pixel.m_y; ++y) {
    for (int x = min_pixel.m_x; x <= max_pixel.m_x; ++x) {
      if (aperture->getArea(centroid_x, centroid_y, x, y) > 0) {
        if (x >= 0 && y >= 0 && x < img->getWidth() && y < img->getHeight()) {
          img->setValue(x, y, value);
        }
      }
    }
  }
}

template <typename T>
void drawAperture(const std::shared_ptr<Aperture> &aperture, SeFloat centroid_x, SeFloat centroid_y,
                  const std::shared_ptr<WriteableImage<T>> &img, T value) {
  auto min_pixel = aperture->getMinPixel(centroid_x, centroid_y);
  auto max_pixel = aperture->getMaxPixel(centroid_x, centroid_y);
 
  for (int y = min_pixel.m_y; y <= max_pixel.m_y; ++y) {
    for (int x = min_pixel.m_x; x <= max_pixel.m_x; ++x) {
      if (aperture->drawArea(centroid_x, centroid_y, x, y) > 0) {
        if (x >= 0 && y >= 0 && x < img->getWidth() && y < img->getHeight()) {
          img->setValue(x, y, value);
        }
      }
    }
  }
}
 
} // end SourceXtractor
 
#endif // _SEFRAMEWORK_SEFRAMEWORK_APERTURE_MEASUREFLUX_H