LSSTApplications  19.0.0-14-gb0260a2+72efe9b372,20.0.0+7927753e06,20.0.0+8829bf0056,20.0.0+995114c5d2,20.0.0+b6f4b2abd1,20.0.0+bddc4f4cbe,20.0.0-1-g253301a+8829bf0056,20.0.0-1-g2b7511a+0d71a2d77f,20.0.0-1-g5b95a8c+7461dd0434,20.0.0-12-g321c96ea+23efe4bbff,20.0.0-16-gfab17e72e+fdf35455f6,20.0.0-2-g0070d88+ba3ffc8f0b,20.0.0-2-g4dae9ad+ee58a624b3,20.0.0-2-g61b8584+5d3db074ba,20.0.0-2-gb780d76+d529cf1a41,20.0.0-2-ged6426c+226a441f5f,20.0.0-2-gf072044+8829bf0056,20.0.0-2-gf1f7952+ee58a624b3,20.0.0-20-geae50cf+e37fec0aee,20.0.0-25-g3dcad98+544a109665,20.0.0-25-g5eafb0f+ee58a624b3,20.0.0-27-g64178ef+f1f297b00a,20.0.0-3-g4cc78c6+e0676b0dc8,20.0.0-3-g8f21e14+4fd2c12c9a,20.0.0-3-gbd60e8c+187b78b4b8,20.0.0-3-gbecbe05+48431fa087,20.0.0-38-ge4adf513+a12e1f8e37,20.0.0-4-g97dc21a+544a109665,20.0.0-4-gb4befbc+087873070b,20.0.0-4-gf910f65+5d3db074ba,20.0.0-5-gdfe0fee+199202a608,20.0.0-5-gfbfe500+d529cf1a41,20.0.0-6-g64f541c+d529cf1a41,20.0.0-6-g9a5b7a1+a1cd37312e,20.0.0-68-ga3f3dda+5fca18c6a4,20.0.0-9-g4aef684+e18322736b,w.2020.45
LSSTDataManagementBasePackage
FitsChan.cc
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 
23 #include <complex>
24 #include <memory>
25 #include <string>
26 #include <vector>
27 
28 #include "astshim/base.h"
29 #include "astshim/Object.h"
30 #include "astshim/Stream.h"
31 #include "astshim/FitsChan.h"
32 
33 namespace ast {
34 
35 namespace {
36 
40 char const *cstrOrNull(std::string const &str) { return str.empty() ? nullptr : str.c_str(); }
41 
42 } // namespace
43 
44 FitsChan::FitsChan(Stream &stream, std::string const &options)
45  : Channel(reinterpret_cast<AstChannel *>(
46  astFitsChan(detail::source, detail::sink, "%s", options.c_str())),
47  stream, true) {
48  assertOK();
49 }
50 
52  // when an astFitsChan is destroyed it first writes out any cards, but if I let astFitsChan
53  // do this automatically then it occurs while the Channel and its Source are being destroyed,
54  // which is too late
55  astWriteFits(getRawPtr());
56  // No 'assertOK' here: can't throw in a Dtor.
57 }
58 
60  std::complex<double> defval) const {
61  std::complex<double> val = defval;
62  // this use of reinterpret_cast is explicitly permitted, for C compatibility
63  double *rawval = reinterpret_cast<double(&)[2]>(val);
64  bool found = astGetFitsCF(getRawPtr(), cstrOrNull(name), rawval);
65  assertOK();
66  return FoundValue<std::complex<double>>(found, val);
67 }
68 
70  char *rawval; // astGetFitsCN has its own static buffer for the value
71  bool found = astGetFitsCN(getRawPtr(), cstrOrNull(name), &rawval);
72  assertOK();
73  std::string val = found ? rawval : defval;
74  return FoundValue<std::string>(found, val);
75 }
76 
77 FoundValue<double> FitsChan::getFitsF(std::string const &name, double defval) const {
78  double val = defval;
79  bool found = astGetFitsF(getRawPtr(), cstrOrNull(name), &val);
80  assertOK();
81  return FoundValue<double>(found, val);
82 }
83 
85  int val = defval;
86  bool found = astGetFitsI(getRawPtr(), cstrOrNull(name), &val);
87  assertOK();
88  return FoundValue<int>(found, val);
89 }
90 
92  int val = static_cast<int>(defval);
93  bool found = astGetFitsL(getRawPtr(), cstrOrNull(name), &val);
94  assertOK();
95  return FoundValue<bool>(found, static_cast<bool>(val));
96 }
97 
99  char *rawval; // astGetFitsS has its own static buffer for the value
100  bool found = astGetFitsS(getRawPtr(), cstrOrNull(name), &rawval);
101  assertOK();
102  std::string val = found ? rawval : defval;
103  return FoundValue<std::string>(found, val);
104 }
105 
107  int const initialIndex = getCard();
108  int const numCards = getNCard();
109  std::vector<std::string> nameList;
110  nameList.reserve(numCards);
111  try {
112  for (auto i = 1; i <= numCards; ++i) {
113  setCard(i);
114  nameList.emplace_back(getCardName());
115  }
116  } catch (...) {
117  setCard(initialIndex);
118  throw;
119  }
120  setCard(initialIndex);
121  return nameList;
122 }
123 
125  std::unique_ptr<char[]> fitsbuf(new char[detail::FITSLEN + 1]);
126  fitsbuf[0] = '\0'; // in case nothing is found
127  bool success = static_cast<bool>(astFindFits(getRawPtr(), name.c_str(), fitsbuf.get(), inc));
128  assertOK();
129  return FoundValue<std::string>(success, std::string(fitsbuf.get()));
130 }
131 
133  int there;
134  int hasvalue = astTestFits(getRawPtr(), cstrOrNull(name), &there);
135  assertOK();
136  if (hasvalue) {
137  return FitsKeyState::PRESENT;
138  }
140 }
141 
142 } // namespace ast
ast::FitsChan::getFitsCN
FoundValue< std::string > getFitsCN(std::string const &name="", std::string defval="") const
Get the value of a CONTINUE card.
Definition: FitsChan.cc:69
ast::FitsKeyState::NOVALUE
@ NOVALUE
keyword is present, but has no value
std::string
STL class.
ast::assertOK
void assertOK(AstObject *rawPtr1=nullptr, AstObject *rawPtr2=nullptr)
Throw std::runtime_error if AST's state is bad.
Definition: base.cc:49
std::vector::reserve
T reserve(T... args)
ast::FoundValue
A value and associated validity flag.
Definition: FitsChan.h:69
std::vector< std::string >
ast::FitsChan::findFits
FoundValue< std::string > findFits(std::string const &name, bool inc)
Search for a card in a FitsChan by keyword.
Definition: FitsChan.cc:124
std::unique_ptr::get
T get(T... args)
ast::FitsChan::getCardName
std::string getCardName() const
Get CardName: the keyword name of the current card.
Definition: FitsChan.h:496
val
ImageT val
Definition: CR.cc:146
ast::FitsChan::getAllCardNames
std::vector< std::string > getAllCardNames()
Get the name of all cards, in order, starting from the first card.
Definition: FitsChan.cc:106
ast::FitsChan::getNCard
int getNCard() const
Get NCard: the number of cards.
Definition: FitsChan.h:555
lsst::afw::geom.transform.transformContinued.name
string name
Definition: transformContinued.py:32
ast::FitsChan::FitsChan
FitsChan(Stream &stream, std::string const &options="")
Construct a channel that uses a provided Stream.
Definition: FitsChan.cc:44
ast::FitsChan::setCard
void setCard(int ind)
Set Card: the index of the current card, where 1 is the first card.
Definition: FitsChan.h:1036
ast::Channel
Channel provides input/output of AST objects.
Definition: Channel.h:60
ast::Stream
A stream for ast::Channel.
Definition: Stream.h:41
ast::FitsChan::getFitsCF
FoundValue< std::complex< double > > getFitsCF(std::string const &name="", std::complex< double > defval={0, 0}) const
Get the value of a complex double card.
Definition: FitsChan.cc:59
FitsChan.h
ast::detail::source
const char * source()
Source function that allows astChannel to source from a Stream.
Definition: Stream.h:224
ast::FitsChan::getFitsI
FoundValue< int > getFitsI(std::string const &name="", int defval=0) const
Get the value of a int card.
Definition: FitsChan.cc:84
ast::FitsChan::testFits
FitsKeyState testFits(std::string const &name="") const
Determine if a card is present, and if so, whether it has a value.
Definition: FitsChan.cc:132
std::string::c_str
T c_str(T... args)
ast::FitsChan::getFitsS
FoundValue< std::string > getFitsS(std::string const &name="", std::string defval="") const
Get the value of a string card.
Definition: FitsChan.cc:98
base.h
ast
AST wrapper classes and functions.
Definition: attributes_channel.dox:1
std::vector::emplace_back
T emplace_back(T... args)
ast::FitsChan::getFitsF
FoundValue< double > getFitsF(std::string const &name="", double defval=0) const
Get the value of a double card.
Definition: FitsChan.cc:77
ast::FitsChan::~FitsChan
virtual ~FitsChan()
Definition: FitsChan.cc:51
ast::FitsChan::getCard
int getCard() const
Get Card: the index of the current card, where 1 is the first card.
Definition: FitsChan.h:486
std::string::empty
T empty(T... args)
std::complex< double >
ast::Object::getRawPtr
AstObject const * getRawPtr() const
Get the raw AST pointer.
Definition: Object.h:292
ast::FitsKeyState
FitsKeyState
Enums describing the presence or absence of a FITS keyword.
Definition: FitsChan.h:40
Object.h
std::unique_ptr
STL class.
ast::FitsKeyState::PRESENT
@ PRESENT
keyword is present and has a value
ast::FitsChan::getFitsL
FoundValue< bool > getFitsL(std::string const &name="", bool defval=false) const
Get the value of a bool card.
Definition: FitsChan.cc:91
ast::FitsKeyState::ABSENT
@ ABSENT
keyword is not present
ast::detail::sink
void sink(const char *cstr)
Sink function that allows astChannel to sink to a Stream.
Definition: Stream.h:239
Stream.h