GloriaValueMap.h

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


#ifndef GloriaValueMap_h
#define GloriaValueMap_h 1

#include "gloria/GloriaFitsTable.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 GloriaValueMap {

  public:


    GloriaValueMap();


    GloriaValueMap(int Nbx, int Nby, int Nmap=1);



     GloriaValueMap(const GloriaFitsImage* fitsImage, int Nbx, int Nby, int Nmap=1);



     GloriaValueMap(GloriaFitsTable* mapTable);


    virtual ~GloriaValueMap();

    void SetImageSize(const GloriaFitsImage* fitsImage);

    void SetImageSize(int SizeX, int SizeY);



    void AddMap();


    void AddMap(string name);


    inline void SetMapName(int imap, string name) 
                                    { _valueNames.at(imap) = name ; } ;


    double GetMapValue(int imap, int ix, int iy, int ipol=1);

    double GetMapValue(int imap, long ipix, int ipol=1);

    double GetMapValue(int imap, double x, double y, int ipol=1);

    inline int Nmap() { return _valueMaps.size(); };

    inline int NMapX() { return _nMapX; };

    inline int NMapY() { return _nMapY; };

    inline int ImgSizeX() { return _sizeX; };

    inline int ImgSizeY() { return _sizeY; };

    inline double* GetMap(int imap) {return _valueMaps.at(imap); };

    inline string GetMapName(int imap) {return _valueNames.at(imap); } ;

    double GetMapMean(int imap);

    GloriaFitsImage* GetMapImage(int imap, int ipol=1);

    GloriaFitsTable* GetTable();

    bool LoadFromTable(GloriaFitsTable* mapTable);

    void ClearMaps();

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

protected:

    // Different interpolation orders
    void FillMapImage0(double* map, GloriaFitsImage* mapImage);
    void FillMapImage1(double* map, GloriaFitsImage* mapImage);
    void FillMapImage2(double* map, GloriaFitsImage* mapImage);

    // Number of map bins in X,Y
    int _nMapX, _nMapY, _nMapXY;

    // Number of bins in the input image
    int _sizeX, _sizeY;
    long _sizeXY;

    // Size of map bin
    int _binX, _binY, _binXY;

    // vector of map arrays
    std::vector<double*> _valueMaps;

    // vector of map names
    std::vector<string> _valueNames;

  } ;


}

#endif