LSST Applications  21.0.0+04719a4bac,21.0.0-1-ga51b5d4+f5e6047307,21.0.0-11-g2b59f77+a9c1acf22d,21.0.0-11-ga42c5b2+86977b0b17,21.0.0-12-gf4ce030+76814010d2,21.0.0-13-g1721dae+760e7a6536,21.0.0-13-g3a573fe+768d78a30a,21.0.0-15-g5a7caf0+f21cbc5713,21.0.0-16-g0fb55c1+b60e2d390c,21.0.0-19-g4cded4ca+71a93a33c0,21.0.0-2-g103fe59+bb20972958,21.0.0-2-g45278ab+04719a4bac,21.0.0-2-g5242d73+3ad5d60fb1,21.0.0-2-g7f82c8f+8babb168e8,21.0.0-2-g8f08a60+06509c8b61,21.0.0-2-g8faa9b5+616205b9df,21.0.0-2-ga326454+8babb168e8,21.0.0-2-gde069b7+5e4aea9c2f,21.0.0-2-gecfae73+1d3a86e577,21.0.0-2-gfc62afb+3ad5d60fb1,21.0.0-25-g1d57be3cd+e73869a214,21.0.0-3-g357aad2+ed88757d29,21.0.0-3-g4a4ce7f+3ad5d60fb1,21.0.0-3-g4be5c26+3ad5d60fb1,21.0.0-3-g65f322c+e0b24896a3,21.0.0-3-g7d9da8d+616205b9df,21.0.0-3-ge02ed75+a9c1acf22d,21.0.0-4-g591bb35+a9c1acf22d,21.0.0-4-g65b4814+b60e2d390c,21.0.0-4-gccdca77+0de219a2bc,21.0.0-4-ge8a399c+6c55c39e83,21.0.0-5-gd00fb1e+05fce91b99,21.0.0-6-gc675373+3ad5d60fb1,21.0.0-64-g1122c245+4fb2b8f86e,21.0.0-7-g04766d7+cd19d05db2,21.0.0-7-gdf92d54+04719a4bac,21.0.0-8-g5674e7b+d1bd76f71f,master-gac4afde19b+a9c1acf22d,w.2021.13
LSST Data Management Base Package
SkyFrame.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_SKYFRAME_H
23 #define ASTSHIM_SKYFRAME_H
24 
25 #include <memory>
26 #include <vector>
27 
28 #include "astshim/Frame.h"
29 
30 namespace ast {
31 
66 class SkyFrame : public Frame {
67  friend class Object;
68 
69 public:
124  explicit SkyFrame(std::string const &options = "")
125  : Frame(reinterpret_cast<AstFrame *>(astSkyFrame("%s", options.c_str()))) {
126  assertOK();
127  }
128 
129  virtual ~SkyFrame() {}
130 
132  SkyFrame(SkyFrame const &) = default;
133  SkyFrame(SkyFrame &&) = default;
134  SkyFrame &operator=(SkyFrame const &) = delete;
135  SkyFrame &operator=(SkyFrame &&) = default;
136 
138  std::shared_ptr<SkyFrame> copy() const { return std::static_pointer_cast<SkyFrame>(copyPolymorphic()); }
139 
141  bool getAlignOffset() const { return getB("AlignOffset"); }
142 
144  bool getAsTime(int axis) const { return getB(detail::formatAxisAttr("AsTime", axis)); }
145 
147  double getEquinox() const { return getD("Equinox"); }
148 
150  bool getIsLatAxis(int axis) const { return getB(detail::formatAxisAttr("IsLatAxis", axis)); }
151 
153  bool getIsLonAxis(int axis) const { return getB(detail::formatAxisAttr("IsLonAxis", axis)); }
154 
156  int getLatAxis() const { return getI("LatAxis"); }
157 
159  int getLonAxis() const { return getI("LonAxis"); }
160 
162  bool getNegLon() const { return getB("NegLon"); }
163 
165  std::string getProjection() const { return getC("Projection"); }
166 
171  for (int axis = 1; axis < 3; ++axis) {
172  ret.push_back(getD(detail::formatAxisAttr("SkyRef", axis)));
173  }
174  return ret;
175  }
176 
178  std::string getSkyRefIs() const { return getC("SkyRefIs"); }
179 
184  for (int axis = 1; axis < 3; ++axis) {
185  ret.push_back(getD(detail::formatAxisAttr("SkyRefP", axis)));
186  }
187  return ret;
188  }
189 
191  double getSkyTol() const { return getD("SkyTol"); }
192 
194  void setAlignOffset(bool alignOffset) { setB("AlignOffset", alignOffset); }
195 
197  void setAsTime(int axis, bool asTime) { setB(detail::formatAxisAttr("AsTime", axis), asTime); }
198 
200  void setEquinox(double equinox) { setD("Equinox", equinox); }
201 
203  void setNegLon(bool negLon) { setB("NegLon", negLon); }
204 
206  void setProjection(std::string const &projection) { setC("Projection", projection); }
207 
209  void setSkyRef(std::vector<double> const &skyRef) {
210  detail::assertEqual(skyRef.size(), "skyRef length", 2u, "number of axes");
211  for (int i = 0; i < 2; ++i) {
212  setD(detail::formatAxisAttr("SkyRef", i + 1), skyRef[i]);
213  }
214  }
215 
217  void setSkyRefIs(std::string const &skyRefIs) { setC("SkyRefIs", skyRefIs); }
218 
220  void setSkyRefP(std::vector<double> const &skyRefP) {
221  detail::assertEqual(skyRefP.size(), "skyRefP length", 2u, "number of axes");
222  for (int i = 0; i < 2; ++i) {
223  setD(detail::formatAxisAttr("SkyRefP", i + 1), skyRefP[i]);
224  }
225  }
226 
228  void setSkyTol(double skyTol) { setD("SkyTol", skyTol); }
229 
239  auto *rawMap = reinterpret_cast<AstObject *>(astSkyOffsetMap(getRawPtr()));
240  assertOK();
241  return Object::fromAstObject<Mapping>(rawMap, false);
242  }
243 
244 protected:
245  virtual std::shared_ptr<Object> copyPolymorphic() const override {
246  return copyImpl<SkyFrame, AstSkyFrame>();
247  }
248 
250  explicit SkyFrame(AstSkyFrame *rawptr) : Frame(reinterpret_cast<AstFrame *>(rawptr)) {
251  if (!astIsASkyFrame(getRawPtr())) {
253  os << "this is a " << getClassName() << ", which is not a SkyFrame";
254  throw std::invalid_argument(os.str());
255  }
256  }
257 };
258 
259 } // namespace ast
260 
261 #endif
table::Key< double > equinox
Definition: OldWcs.cc:135
std::ostream * os
Definition: Schema.cc:746
Frame is used to represent a coordinate system.
Definition: Frame.h:157
Abstract base class for all AST objects.
Definition: Object.h:51
void setD(std::string const &attrib, double value)
Set the value of an attribute as a double.
Definition: Object.h:476
double getD(std::string const &attrib) const
Get the value of an attribute as a double.
Definition: Object.h:374
std::string const getC(std::string const &attrib) const
Get the value of an attribute as a string.
Definition: Object.h:361
std::string getClassName() const
Get Class: the name of the class (e.g.
Definition: Object.h:139
bool getB(std::string const &attrib) const
Get the value of an attribute as a bool.
Definition: Object.h:348
AstObject const * getRawPtr() const
Get the raw AST pointer.
Definition: Object.h:292
int getI(std::string const &attrib) const
Get the value of an attribute as an int.
Definition: Object.h:400
void setC(std::string const &attrib, std::string const &value)
Set the value of an attribute as a string.
Definition: Object.h:464
void setB(std::string const &attrib, bool value)
Set the value of an attribute as a bool.
Definition: Object.h:452
SkyFrame is a specialised form of Frame which describes celestial longitude/latitude coordinate syste...
Definition: SkyFrame.h:66
std::shared_ptr< Mapping > skyOffsetMap()
Get a sky offset map.
Definition: SkyFrame.h:238
SkyFrame & operator=(SkyFrame const &)=delete
bool getNegLon() const
Get NegLon: display longitude values in the range [-pi,pi]?
Definition: SkyFrame.h:162
void setNegLon(bool negLon)
Set NegLon: display longitude values in the range [-pi,pi]?
Definition: SkyFrame.h:203
bool getAsTime(int axis) const
Get AsTime(axis) for one axis: format celestial coordinates as times?
Definition: SkyFrame.h:144
SkyFrame(std::string const &options="")
Construct a SkyFrame.
Definition: SkyFrame.h:124
SkyFrame(SkyFrame &&)=default
std::string getProjection() const
Get Projection: sky projection description.
Definition: SkyFrame.h:165
virtual std::shared_ptr< Object > copyPolymorphic() const override
Return a deep copy of this object.
Definition: SkyFrame.h:245
void setSkyRef(std::vector< double > const &skyRef)
Set SkyRef: position defining location of the offset coordinate system.
Definition: SkyFrame.h:209
void setProjection(std::string const &projection)
Set Projection: sky projection description.
Definition: SkyFrame.h:206
std::vector< double > getSkyRefP() const
Get SkyRefP for both axes: position defining orientation of the offset coordinate system.
Definition: SkyFrame.h:182
double getSkyTol() const
Get SkyTol: smallest significant shift in sky coordinates.
Definition: SkyFrame.h:191
virtual ~SkyFrame()
Definition: SkyFrame.h:129
void setEquinox(double equinox)
Set Equinox: epoch of the mean equinox.
Definition: SkyFrame.h:200
SkyFrame & operator=(SkyFrame &&)=default
int getLatAxis() const
Get LatAxis: index of the latitude axis.
Definition: SkyFrame.h:156
bool getIsLonAxis(int axis) const
Get IsLonAxis(axis) for one axis: is the specified axis the longitude axis?
Definition: SkyFrame.h:153
bool getIsLatAxis(int axis) const
Get IsLatAxis(axis) for one axis: is the specified axis the latitude axis?
Definition: SkyFrame.h:150
SkyFrame(AstSkyFrame *rawptr)
Construct a SkyFrame from a raw AST pointer.
Definition: SkyFrame.h:250
void setSkyRefP(std::vector< double > const &skyRefP)
Set SkyRefP: position defining orientation of the offset coordinate system.
Definition: SkyFrame.h:220
void setSkyRefIs(std::string const &skyRefIs)
Set SkyRefIs: selects the nature of the offset coordinate system.
Definition: SkyFrame.h:217
SkyFrame(SkyFrame const &)=default
Copy constructor: make a deep copy.
void setAsTime(int axis, bool asTime)
Set AsTime(axis) for one axis: format celestial coordinates as times?
Definition: SkyFrame.h:197
void setSkyTol(double skyTol)
Set SkyTol: smallest significant shift in sky coordinates.
Definition: SkyFrame.h:228
int getLonAxis() const
Get LonAxis: index of the longitude axis.
Definition: SkyFrame.h:159
std::shared_ptr< SkyFrame > copy() const
Return a deep copy of this object.
Definition: SkyFrame.h:138
bool getAlignOffset() const
Get AlignOffset: align SkyFrames using the offset coordinate system?
Definition: SkyFrame.h:141
void setAlignOffset(bool alignOffset)
Set AlignOffset: align SkyFrames using the offset coordinate system?
Definition: SkyFrame.h:194
double getEquinox() const
Get Equinox: epoch of the mean equinox.
Definition: SkyFrame.h:147
std::string getSkyRefIs() const
Get SkyRefIs: selects the nature of the offset coordinate system.
Definition: SkyFrame.h:178
std::vector< double > getSkyRef() const
Get SkyRef for both axes: position defining location of the offset coordinate system.
Definition: SkyFrame.h:169
std::string formatAxisAttr(std::string const &name, int axis)
Format an axis-specific attribute by appending the axis index.
Definition: utils.h:79
void assertEqual(T1 val1, std::string const &descr1, T2 val2, std::string const &descr2)
Definition: utils.h:48
AST wrapper classes and functions.
void assertOK(AstObject *rawPtr1=nullptr, AstObject *rawPtr2=nullptr)
Throw std::runtime_error if AST's state is bad.
Definition: base.cc:49
T push_back(T... args)
T size(T... args)