GloriaMedianMap.h

// -*- mode: c++; mode: auto-fill; mode: flyspell-prog; -*-


#ifndef GloriaMedianMap_h
#define GloriaMedianMap_h 1

#include "GloriaValueMap.h"
#include "gloria/GloriaFitsImage.h"

/*
 *   This source code is part of the LUIZA software package for GLORIA.
 *   You are free to use this source files for your own development as
 *   long as it stays in a public research context. You are not
 *   allowed to use it for commercial purpose. You must put this
 *   header with author names in all development based on this file.
 *
 */

namespace gloria {


  class GloriaMedianMap : public GloriaValueMap {

  public:

    GloriaMedianMap() : GloriaValueMap() { };


    GloriaMedianMap(int Nbx, int Nby, int Nvar = 1);



    GloriaMedianMap(const GloriaFitsImage* fitsImage, int Nbx, int Nby, int Nvar=1);


    GloriaMedianMap(GloriaFitsTable* mapTable) : GloriaValueMap(mapTable) { };


    void AddMedian();


    void AddMedian(string name);


    void SetMedianName(int ivar, string name);
 

    inline string GetMedianName(int ivar) {return _valueNames.at(2*ivar); } ;



    long CalculateMap(int iVar, int nVal, double* Xpos, double* Ypos, double* Val);


    long CalculateMap(int iVar, vector<double>* Xpos, vector<double>* Ypos, vector<double>* Val);




    void CalculateBgMap(const GloriaFitsImage* fitsImage, int Nbx, int Nby, int sampling=1);

    void CalculateBgMap(const GloriaFitsImage* fitsImage, int sampling=1);



    inline double GetMedian(int ivar, int ix, int iy, int ipol=1)
                         { return GetMapValue(2*ivar,ix,iy,ipol); } ;

    inline double GetMedian(int ivar, double x, double y, int ipol=1)
                         { return GetMapValue(2*ivar,x,y,ipol); } ;

    inline double GetMedian(int ivar, long ipix, int ipol=1)
                        { return GetMapValue(2*ivar,ipix,ipol); } ;

    inline double GetSpread(int ivar, int ix, int iy, int ipol=1)
                         { return GetMapValue(2*ivar+1,ix,iy,ipol); } ;

    inline double GetSpread(int ivar, double x, double y, int ipol=1)
                         { return GetMapValue(2*ivar+1,x,y,ipol); } ;

    inline double GetSpread(int ivar, long ipix, int ipol=1)
                        { return GetMapValue(2*ivar+1,ipix,ipol); } ;



    inline double GetBgLevel(int ix, int iy, int ipol=1)
                         { return GetMapValue(0,ix,iy,ipol); } ;

    inline double GetBgLevel(double x, double y, int ipol=1)
                         { return GetMapValue(0,x,y,ipol); } ;

    inline double GetBgLevel(long ipix, int ipol=1)
                        { return GetMapValue(0,ipix,ipol); } ;

    inline double GetNoise(int ix, int iy, int ipol=1)
                         { return GetMapValue(1,ix,iy,ipol); } ;

    inline double GetNoise(double x, double y, int ipol=1)
                         { return GetMapValue(1,x,y,ipol); } ;

    inline double GetNoise(long ipix, int ipol=1)
                        { return GetMapValue(1,ipix,ipol); } ;



    inline double* GetGloriaMedianMap(int ivar=0) {return _valueMaps.at(2*ivar); };

    inline double* GetSpreadMap(int ivar=0) {return _valueMaps.at(2*ivar+1); };

    inline double* GetBgMap() {return _valueMaps.at(0); };

    inline double* GetNoiseMap() {return _valueMaps.at(1); };



    inline  double GetMedianMean(int ivar=0){ return GetMapMean(2*ivar);} ;

    inline double GetSpreadMean(int ivar=0){ return GetMapMean(2*ivar+1);} ;

    inline double GetBgMean() { return GetMapMean(0); };

    inline double GetNoiseMean() { return GetMapMean(1); } ;



    inline GloriaFitsImage* GetMedianImage(int ivar=0, int ipol=1)
                    { return GetMapImage(2*ivar, ipol);}  ;

    inline GloriaFitsImage* GetSpreadImage(int ivar=0, int ipol=1)
                    { return GetMapImage(2*ivar+1, ipol);}  ;


    inline GloriaFitsImage* GetBgImage(int ipol=1)
                    { return GetMapImage(0, ipol);}  ;

    inline GloriaFitsImage* GetNoiseImage(int ipol=1)
                    { return GetMapImage(1, ipol);}  ;

    void InitializeMaps(int Nbx, int Nby, int Nvar);

  } ;


}

#endif