sutra_centroider_pyr.h

Go to the documentation of this file.

// -----------------------------------------------------------------------------
//  This file is part of COMPASS <https://anr-compass.github.io/compass/>
//
//  Copyright (C) 2011-2023 COMPASS Team <https://github.com/ANR-COMPASS>
//  All rights reserved.
// -----------------------------------------------------------------------------
 
 
#ifndef _SUTRA_CENTROIDER_PYR_H_
#define _SUTRA_CENTROIDER_PYR_H_
 
#include <sutra_centroider.h>
 
struct Method_CoG {
  bool is_local = false;
  bool is_sinus = true;
 
  Method_CoG(bool isLocal_ = false, bool isSinus_ = true)
      : is_local(isLocal_), is_sinus(isSinus_) {}
 
  Method_CoG(uint8_t method) : is_local(method > 1), is_sinus(!(method % 2)) {}
 
  static const char *str(const struct Method_CoG &method) {
    if (method.is_sinus && !method.is_local) return "sinus global";     // 0
    if (!method.is_sinus && !method.is_local) return "nosinus global";  // 1
    if (method.is_sinus && method.is_local) return "sinus local";       // 2
    if (!method.is_sinus && method.is_local) return "nosinus local";    // 3
    throw "method unknown";
  };
};
 
template <class Tin, class T>
class SutraCentroiderPyr : public SutraCentroider<Tin, T> {
 private:
  string pyr_type;
 
 public:
  SutraCentroiderPyr(CarmaContext *context, SutraWfs *wfs, long nvalid,
                       float offset, float scale, bool filter_TT, int device);
  SutraCentroiderPyr(const SutraCentroiderPyr &centroider);
  ~SutraCentroiderPyr();
 
  string get_type();
  int set_valid_thresh(T valid_thresh);
  T get_valid_thresh();
 
  int set_method(Method_CoG method);
  Method_CoG get_method();
  string get_method_str();
 
  int get_pyr(float *cube, float *intensities, T *centroids, int *subindx,
              int *subindy, int nvalid, int ns, int nim);
  int get_cog(float *cube, float *intensities, T *centroids, int nvalid,
              int npix, int ntot, cudaStream_t stream=0);
  int get_cog(float *intensities, T *slopes, bool noise);
  int get_cog();
 
 private:
  T valid_thresh;
  Method_CoG method;
};
 
template <class T>
void pyr_slopes(T *d_odata, T *d_idata, int *subindx, int *subindy,
                float *intensities, int ns, int nvalid, int nim,
                SlopeOrder slope_order, CarmaDevice *device);
 
template <class T>
void pyr2_slopes(T *d_odata, T *ref, T *d_idata, int *subindx, int *subindy,
                 float *intensities, int ns, int nvalid, float scale,
                 T valid_thresh, int do_sin, SlopeOrder slope_order,
                 CarmaDevice *device);
#endif  // _SUTRA_CENTROIDER_PYR_H_