GWcsHPX.hpp
5.18 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
/***************************************************************************
* GWcsHPX.hpp - Healpix projection class *
* ----------------------------------------------------------------------- *
* copyright (C) 2010-2011 by Jurgen Knodlseder *
* ----------------------------------------------------------------------- *
* *
* This program is free software: you can redistribute it and/or modify *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation, either version 3 of the License, or *
* (at your option) any later version. *
* *
* This program is distributed in the hope that it will be useful, *
* but WITHOUT ANY WARRANTY; without even the implied warranty of *
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
* GNU General Public License for more details. *
* *
* You should have received a copy of the GNU General Public License *
* along with this program. If not, see <http://www.gnu.org/licenses/>. *
* *
***************************************************************************/
/**
* @file GWcsHPX.hpp
* @brief HealPix projection class definition
* @author J. Knodlseder
*/
#ifndef GWCSHPX_HPP
#define GWCSHPX_HPP
/* __ Includes ___________________________________________________________ */
#include "GWcs.hpp"
#include "GFitsHDU.hpp"
#include "GSkyDir.hpp"
#include "GSkyPixel.hpp"
/***********************************************************************//**
* @class GWcsHPX
*
* @brief HealPix projection class interface defintion
*
* The HealPix projection class has been implemented by adapting code from
* the HealPix library (version 2.1). For more information about HEALPix, see
* http://healpix.jpl.nasa.gov
***************************************************************************/
class GWcsHPX : public GWcs {
public:
// Constructors and destructors
GWcsHPX(void);
explicit GWcsHPX(const int& nside, const std::string& ordering = "NESTED",
const std::string& coordsys = "GAL");
explicit GWcsHPX(const GFitsHDU* hdu);
GWcsHPX(const GWcsHPX& wcs);
virtual ~GWcsHPX(void);
// Operators
GWcsHPX& operator= (const GWcsHPX& wcs);
// Implemented pure virtual base class methods
virtual void clear(void);
virtual GWcsHPX* clone(void) const;
virtual std::string code(void) const { return "HPX"; }
virtual std::string name(void) const { return "HealPix"; }
virtual void read(const GFitsHDU* hdu);
virtual void write(GFitsHDU* hdu) const;
virtual double omega(const int& pix) const;
virtual double omega(const GSkyPixel& pix) const;
virtual GSkyDir pix2dir(const int& pix) const;
virtual int dir2pix(const GSkyDir& dir) const;
virtual GSkyDir xy2dir(const GSkyPixel& pix) const;
virtual GSkyPixel dir2xy(const GSkyDir& dir) const;
virtual std::string print(void) const;
// Additional class specific methods
int npix(void) const;
int nside(void) const;
std::string ordering(void) const;
void ordering(const std::string& ordering);
private:
// Private methods
void init_members(void);
void copy_members(const GWcsHPX& wcs);
void free_members(void);
virtual bool compare(const GWcs& wcs) const;
void std2nat(GVector *coord) const { return; }
void nat2std(GVector *coord) const { return; }
int nside2order(int nside);
void pix2xy(const int& ipix, int* x, int* y) const;
int xy2pix(int x, int y) const;
void pix2ang_ring(int ipix, double* theta, double* phi) const;
void pix2ang_nest(int ipix, double* theta, double* phi) const;
int ang2pix_z_phi_ring(double z, double phi) const;
int ang2pix_z_phi_nest(double z, double phi) const;
unsigned int isqrt(unsigned int arg) const;
// NEW VERSION
void prj_set(void);
int prj_x2s(int nx, int ny, int sxy, int spt,
const double* x, const double* y,
double* phi, double* theta, int* stat) const;
int prj_s2x(int nphi, int ntheta, int spt, int sxy,
const double* phi, const double* theta,
double* x, double* y, int* stat) const;
// Private data area
int m_nside; //!< Number of divisions of each base pixel (1-8192)
int m_npface; //!<
int m_ncap; //!<
int m_order; //!< Order
int m_ordering; //!< Pixel ordering (0=ring, 1=nested, -1=?)
int m_num_pixels; //!< Number of pixels
double m_fact1; //!<
double m_fact2; //!<
double m_omega; //!< Solid angle of pixel
};
#endif /* GWCSHPX_HPP */