VariablePsfStack.h

Go to the documentation of this file.

/*
 * VariablePsf.h
 *
 *  Created on: Jun 25, 2018
 *      Author:
 */

#ifndef _SEIMPLEMENTATION_PSF_VARIABLEPSFSTACK_H_
#define _SEIMPLEMENTATION_PSF_VARIABLEPSFSTACK_H_

#include <CCfits/CCfits>
#include <SEFramework/Image/VectorImage.h>
#include <SEFramework/Psf/Psf.h>

namespace SourceXtractor {

class VariablePsfStack final : public Psf {
public:
  explicit VariablePsfStack(std::shared_ptr<CCfits::FITS> pFits) : m_pFits(pFits), m_psf_size(0), mm_pixel_sampling(0.0) {
    setup(pFits);
    selfTest();
  };

  virtual ~VariablePsfStack() = default;

  int getWidth() const {
    return m_psf_size;
  };

  int getHeight() const {
    return m_psf_size;
  };

  double getPixelSampling() const {
    return mm_pixel_sampling;
  };

  const std::vector<std::string>& getComponents() const {
    return m_components;
  };

  virtual std::shared_ptr<VectorImage<SeFloat>> getPsf(const std::vector<double>& values) const;

private:
  mutable std::mutex m_mutex;
  std::shared_ptr<CCfits::FITS> m_pFits;

  int m_psf_size;
  int m_grid_offset;

  double mm_pixel_sampling;

  long m_nrows;

  std::vector<SeFloat> m_ra_values;
  std::vector<SeFloat> m_dec_values;
  std::vector<SeFloat> m_x_values;
  std::vector<SeFloat> m_y_values;
  std::vector<int>     m_gridx_values;
  std::vector<int>     m_gridy_values;

  std::vector<std::string> m_components = {"X_IMAGE", "Y_IMAGE"};

  /*
   * Check the file, load the positions and so on
   */
  void setup(std::shared_ptr<CCfits::FITS> pFits);

  /*
   * consistency of the stackedPSF
   */
  void selfTest();
};

}  // namespace SourceXtractor

#endif  //_SEIMPLEMENTATION_PSF_VARIABLEPSFSTACK_H_