Loading [MathJax]/extensions/tex2jax.js
LSST Applications 28.0.0,g1653933729+a8ce1bb630,g1a997c3884+a8ce1bb630,g28da252d5a+5bd70b7e6d,g2bbee38e9b+638fca75ac,g2bc492864f+638fca75ac,g3156d2b45e+07302053f8,g347aa1857d+638fca75ac,g35bb328faa+a8ce1bb630,g3a166c0a6a+638fca75ac,g3e281a1b8c+7bbb0b2507,g4005a62e65+17cd334064,g414038480c+5b5cd4fff3,g41af890bb2+4ffae9de63,g4e1a3235cc+0f1912dca3,g6249c6f860+3c3976f90c,g80478fca09+46aba80bd6,g82479be7b0+77990446f6,g858d7b2824+78ba4d1ce1,g89c8672015+f667a5183b,g9125e01d80+a8ce1bb630,ga5288a1d22+2a6264e9ca,gae0086650b+a8ce1bb630,gb58c049af0+d64f4d3760,gc22bb204ba+78ba4d1ce1,gc28159a63d+638fca75ac,gcf0d15dbbd+32ddb6096f,gd6b7c0dfd1+3e339405e9,gda3e153d99+78ba4d1ce1,gda6a2b7d83+32ddb6096f,gdaeeff99f8+1711a396fd,gdd5a9049c5+b18c39e5e3,ge2409df99d+a5e4577cdc,ge33fd446bb+78ba4d1ce1,ge79ae78c31+638fca75ac,gf0baf85859+64e8883e75,gf5289d68f6+e1b046a8d7,gfa443fc69c+91d9ed1ecf,gfda6b12a05+8419469a56
LSST Data Management Base Package
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
WcsMap.h
Go to the documentation of this file.
1/*
2 * LSST Data Management System
3 * Copyright 2017 AURA/LSST.
4 *
5 * This product includes software developed by the
6 * LSST Project (http://www.lsst.org/).
7 *
8 * This program is free software: you can redistribute it and/or modify
9 * it under the terms of the GNU General Public License as published by
10 * the Free Software Foundation, either version 3 of the License, or
11 * (at your option) any later version.
12 *
13 * This program is distributed in the hope that it will be useful,
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16 * GNU General Public License for more details.
17 *
18 * You should have received a copy of the LSST License Statement and
19 * the GNU General Public License along with this program. If not,
20 * see <https://www.lsstcorp.org/LegalNotices/>.
21 */
22#ifndef ASTSHIM_WCSMAP_H
23#define ASTSHIM_WCSMAP_H
24
25#include <memory>
26#include <sstream>
27#include <utility>
28
29#include "astshim/base.h"
31#include "astshim/Mapping.h"
32
33namespace ast {
34
48enum class WcsType {
49 AZP = AST__AZP,
50 SZP = AST__SZP,
51 TAN = AST__TAN,
52 STG = AST__STG,
53 SIN = AST__SIN,
54 ARC = AST__ARC,
55 ZPN = AST__ZPN,
56 ZEA = AST__ZEA,
57 AIR = AST__AIR,
58 CYP = AST__CYP,
59 CEA = AST__CEA,
60 CAR = AST__CAR,
61 MER = AST__MER,
62 SFL = AST__SFL,
63 PAR = AST__PAR,
64 MOL = AST__MOL,
65 AIT = AST__AIT,
66 COP = AST__COP,
67 COE = AST__COE,
68 COD = AST__COD,
69 COO = AST__COO,
70 BON = AST__BON,
71 PCO = AST__PCO,
72 TSC = AST__TSC,
73 CSC = AST__CSC,
74 QSC = AST__QSC,
75 NCP = AST__NCP,
76 GLS = AST__GLS,
77 TPN = AST__TPN,
78 HPX = AST__HPX,
79 XPH = AST__XPH,
80 WCSBAD = AST__WCSBAD,
81};
82
146class WcsMap : public Mapping {
147 friend class Object;
148
149public:
209 explicit WcsMap(int ncoord, WcsType type, int lonax, int latax, std::string const &options = "")
210 : Mapping(reinterpret_cast<AstMapping *>(
211 astWcsMap(ncoord, static_cast<int>(type), lonax, latax, "%s", options.c_str()))) {
212 assertOK();
213 }
214
215 virtual ~WcsMap() {}
216
218 WcsMap(WcsMap const &) = default;
219 WcsMap(WcsMap &&) = default;
220 WcsMap &operator=(WcsMap const &) = delete;
221 WcsMap &operator=(WcsMap &&) = default;
222
225
227 double getNatLat() const { return getD("NatLat"); }
228
230 double getNatLon() const { return getD("NatLon"); }
231
238 double getPVi_m(int i, int m) const {
240 os << "PV" << i << "_" << m;
241 return getD(os.str());
242 }
243
245 int getPVMax(int axis) const { return getI(detail::formatAxisAttr("PVMax", axis)); }
246
251 return std::make_pair(getI(detail::formatAxisAttr("WcsAxis", 1)),
252 getI(detail::formatAxisAttr("WcsAxis", 2)));
253 }
254
256 WcsType getWcsType() const { return static_cast<WcsType>(getI("WcsType")); }
257
258protected:
259 virtual std::shared_ptr<Object> copyPolymorphic() const override {
261 }
262
264 explicit WcsMap(AstWcsMap *rawptr) : Mapping(reinterpret_cast<AstMapping *>(rawptr)) {
265 if (!astIsAWcsMap(getRawPtr())) {
267 os << "this is a " << getClassName() << ", which is not a WcsMap";
268 throw std::invalid_argument(os.str());
269 }
270 }
271};
272
273} // namespace ast
274
275#endif
std::ostream * os
Definition Schema.cc:557
int m
Definition SpanSet.cc:48
An abstract base class for objects which transform one set of coordinates to another.
Definition Mapping.h:59
Abstract base class for all AST objects.
Definition Object.h:51
std::shared_ptr< T > copyImpl() const
Implementation of deep copy.
Definition Object.h:320
double getD(std::string const &attrib) const
Get the value of an attribute as a double.
Definition Object.h:374
std::string getClassName() const
Get Class: the name of the class (e.g.
Definition Object.h:139
int getI(std::string const &attrib) const
Get the value of an attribute as an int.
Definition Object.cc:178
AstObject const * getRawPtr() const
Get the raw AST pointer.
Definition Object.h:292
Map from a spherical system to a cartesian system using standard FITS sky coordinate projections.
Definition WcsMap.h:146
WcsType getWcsType() const
Get WcsType: FITS-WCS projection type.
Definition WcsMap.h:256
double getNatLat() const
get NatLat: native latitude of the reference point of a FITS-WCS projection.
Definition WcsMap.h:227
double getNatLon() const
get NatLon: native longitude of the reference point of a FITS-WCS projection.
Definition WcsMap.h:230
WcsMap(WcsMap const &)=default
Copy constructor: make a deep copy.
WcsMap & operator=(WcsMap const &)=delete
WcsMap(int ncoord, WcsType type, int lonax, int latax, std::string const &options="")
Create a WcsMap.
Definition WcsMap.h:209
std::shared_ptr< WcsMap > copy() const
Return a deep copy of this object.
Definition WcsMap.h:224
WcsMap(AstWcsMap *rawptr)
Construct a WcsMap from a raw AST pointer.
Definition WcsMap.h:264
int getPVMax(int axis) const
Get PVMax(axis) for one axis: maximum number of FITS-WCS projection parameters.
Definition WcsMap.h:245
std::pair< int, int > getWcsAxis() const
Get WcsAxis(lonlat) FITS-WCS projection axis for longitude, latitude.
Definition WcsMap.h:250
WcsMap & operator=(WcsMap &&)=default
virtual ~WcsMap()
Definition WcsMap.h:215
double getPVi_m(int i, int m) const
Get PVi_m for one value of i and m: a FITS-WCS projection parameter.
Definition WcsMap.h:238
WcsMap(WcsMap &&)=default
virtual std::shared_ptr< Object > copyPolymorphic() const override
Return a deep copy of this object.
Definition WcsMap.h:259
T make_pair(T... args)
std::string formatAxisAttr(std::string const &name, int axis)
Format an axis-specific attribute by appending the axis index.
Definition utils.h:79
AST wrapper classes and functions.
WcsType
WCS types that give the projection type code (in upper case) as used in the FITS-WCS "CTYPEi" keyword...
Definition WcsMap.h:48
void assertOK(AstObject *rawPtr1=nullptr, AstObject *rawPtr2=nullptr)
Throw std::runtime_error if AST's state is bad.
Definition base.cc:49
T static_pointer_cast(T... args)