LSSTApplications  17.0+124,17.0+14,17.0+73,18.0.0+37,18.0.0+80,18.0.0-4-g68ffd23+4,18.1.0-1-g0001055+12,18.1.0-1-g03d53ef+5,18.1.0-1-g1349e88+55,18.1.0-1-g2505f39+44,18.1.0-1-g5315e5e+4,18.1.0-1-g5e4b7ea+14,18.1.0-1-g7e8fceb+4,18.1.0-1-g85f8cd4+48,18.1.0-1-g8ff0b9f+4,18.1.0-1-ga2c679d+1,18.1.0-1-gd55f500+35,18.1.0-10-gb58edde+2,18.1.0-11-g0997b02+4,18.1.0-13-gfe4edf0b+12,18.1.0-14-g259bd21+21,18.1.0-19-gdb69f3f+2,18.1.0-2-g5f9922c+24,18.1.0-2-gd3b74e5+11,18.1.0-2-gfbf3545+32,18.1.0-26-g728bddb4+5,18.1.0-27-g6ff7ca9+2,18.1.0-3-g52aa583+25,18.1.0-3-g8ea57af+9,18.1.0-3-gb69f684+42,18.1.0-3-gfcaddf3+6,18.1.0-32-gd8786685a,18.1.0-4-gf3f9b77+6,18.1.0-5-g1dd662b+2,18.1.0-5-g6dbcb01+41,18.1.0-6-gae77429+3,18.1.0-7-g9d75d83+9,18.1.0-7-gae09a6d+30,18.1.0-9-gc381ef5+4,w.2019.45
LSSTDataManagementBasePackage
CcdImage.h
Go to the documentation of this file.
1 // -*- LSST-C++ -*-
2 /*
3  * This file is part of jointcal.
4  *
5  * Developed for the LSST Data Management System.
6  * This product includes software developed by the LSST Project
7  * (https://www.lsst.org).
8  * See the COPYRIGHT file at the top-level directory of this distribution
9  * for details of code ownership.
10  *
11  * This program is free software: you can redistribute it and/or modify
12  * it under the terms of the GNU General Public License as published by
13  * the Free Software Foundation, either version 3 of the License, or
14  * (at your option) any later version.
15  *
16  * This program is distributed in the hope that it will be useful,
17  * but WITHOUT ANY WARRANTY; without even the implied warranty of
18  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
19  * GNU General Public License for more details.
20  *
21  * You should have received a copy of the GNU General Public License
22  * along with this program. If not, see <https://www.gnu.org/licenses/>.
23  */
24 
25 #ifndef LSST_JOINTCAL_CCD_IMAGE_H
26 #define LSST_JOINTCAL_CCD_IMAGE_H
27 
28 #include <list>
29 #include <string>
30 
32 #include "lsst/afw/table/Source.h"
33 #include "lsst/afw/geom/SkyWcs.h"
37 #include "lsst/geom/Box.h"
38 #include "lsst/geom/SpherePoint.h"
41 #include "lsst/jointcal/Frame.h"
42 
43 namespace lsst {
44 namespace jointcal {
45 
47 
48 typedef int VisitIdType;
49 typedef int CcdIdType;
51 struct CcdImageKey {
52  VisitIdType visit;
53  CcdIdType ccd;
54  bool operator!=(CcdImageKey const &right) const { return !(*this == right); }
55  bool operator==(CcdImageKey const &right) const { return (visit == right.visit) && (ccd == right.ccd); }
56 };
57 // typedef std::pair<VisitIdType, CcdIdType> CcdImageKey;
59 
64 class CcdImage {
65 public:
70  std::string const &fluxField);
71 
73  CcdImage(CcdImage const &) = delete;
74  CcdImage(CcdImage &&) = delete;
75  CcdImage &operator=(CcdImage const &) = delete;
76  CcdImage &operator=(CcdImage &&) = delete;
77 
79  std::string getName() const { return _name; }
80 
86  MeasuredStarList const &getWholeCatalog() const { return _wholeCatalog; }
87 
89 
94  MeasuredStarList const &getCatalogForFit() const { return _catalogForFit; }
95  MeasuredStarList &getCatalogForFit() { return _catalogForFit; }
97 
100  getCatalogForFit().clear();
101  getWholeCatalog().copyTo(getCatalogForFit());
102  }
103 
112  std::pair<int, int> countStars() const;
113 
119  void setCommonTangentPoint(Point const &commonTangentPoint);
120 
126  jointcal::Point const &getCommonTangentPoint() const { return _commonTangentPoint; }
127 
129  return _pixelToCommonTangentPlane;
130  }
131 
133  return _commonTangentPlaneToTangentPlane;
134  }
135 
137  return _tangentPlaneToCommonTangentPlane;
138  }
139 
140  std::shared_ptr<AstrometryTransform> const getPixelToTangentPlane() const { return _pixelToTangentPlane; }
141 
142  std::shared_ptr<AstrometryTransform> const getSkyToTangentPlane() const { return _skyToTangentPlane; }
143 
145  CcdIdType getCcdId() const { return _ccdId; }
146 
148  VisitIdType getVisit() const { return _visit; }
149 
151 
152  CcdImageKey getHashKey() const { return CcdImageKey{_visit, _ccdId}; }
153 
155  double getAirMass() const { return _airMass; }
156 
158  double getMjd() const { return _mjd; }
159 
162 
166  lsst::geom::SpherePoint getBoresightRaDec() const { return _boresightRaDec; }
167 
168  double getHourAngle() const { return _hourAngle; }
169 
170  double getLstObs() const { return _lstObs; }
171 
173  double getSinEta() const { return _sinEta; }
174 
176  double getCosEta() const { return _cosEta; }
177 
179  double getTanZ() const { return _tanZ; }
180 
182  Point getRefractionVector() const { return Point(_tanZ * _cosEta, _tanZ * _sinEta); }
183 
185  std::string getFilter() const { return _filter; }
186 
188  std::shared_ptr<AstrometryTransform> const getReadWcs() const { return _readWcs; }
189 
191  Frame const &getImageFrame() const { return _imageFrame; }
192 
193 private:
195  std::string const &fluxField);
196 
197  jointcal::Frame _imageFrame; // in pixels
198 
199  MeasuredStarList _wholeCatalog; // the catalog of measured objets
200  MeasuredStarList _catalogForFit;
201 
202  std::shared_ptr<AstrometryTransformSkyWcs> _readWcs; // apply goes from pix to sky
203 
204  // The following ones should probably be mostly removed.
205  // go from CommonTangentPlane to this tangent plane.
206  std::shared_ptr<AstrometryTransform> _commonTangentPlaneToTangentPlane;
207  std::shared_ptr<AstrometryTransform> _tangentPlaneToCommonTangentPlane; // reverse one
208  std::shared_ptr<AstrometryTransform> _pixelToCommonTangentPlane; // pixels -> CTP
209  std::shared_ptr<AstrometryTransform> _pixelToTangentPlane;
210 
211  std::shared_ptr<AstrometryTransform> _skyToTangentPlane;
212 
213  std::string _name;
214  CcdIdType _ccdId;
215  VisitIdType _visit;
216 
217  lsst::geom::SpherePoint _boresightRaDec;
218  double _airMass; // airmass value.
219  double _mjd; // modified julian date
222  // refraction
223  // eta : parallactic angle, z: zenithal angle (X = 1/cos(z))
224  double _sinEta, _cosEta, _tanZ;
225  // Local Sidereal Time and hour angle of observation
226  double _lstObs, _hourAngle;
227 
228  std::string _filter;
229 
230  jointcal::Point _commonTangentPoint;
231 };
232 } // namespace jointcal
233 } // namespace lsst
234 
235 // Add our preferred hash of CcdImageKey to the std:: namespace, so it's always available "for free".
236 namespace std {
237 template <>
246  return hash<size_t>()(static_cast<size_t>(key.visit) | (static_cast<size_t>(key.ccd) << 32));
247  }
248 };
249 } // namespace std
250 
251 #endif // LSST_JOINTCAL_CCD_IMAGE_H
AmpInfoBoxKey bbox
Definition: Amplifier.cc:117
VisitIdType getVisit() const
returns visit ID
Definition: CcdImage.h:148
A point in a plane.
Definition: Point.h:36
bool operator==(CcdImageKey const &right) const
Definition: CcdImage.h:55
double getCosEta() const
Parallactic angle.
Definition: CcdImage.h:176
std::ostream & operator<<(std::ostream &stream, AstrometryTransform const &transform)
Delegates to transform.dump()
std::shared_ptr< AstrometryTransform > const getPixelToTangentPlane() const
Definition: CcdImage.h:140
std::string getName() const
Return the _name that identifies this ccdImage.
Definition: CcdImage.h:79
std::shared_ptr< AstrometryTransform > const getReadWcs() const
the wcs read in the header. NOT updated when fitting.
Definition: CcdImage.h:188
For hashing a ccdImage: the pair of (visit, ccd) IDs should be unique to each ccdImage.
Definition: CcdImage.h:51
double getMjd() const
Julian Date.
Definition: CcdImage.h:158
STL namespace.
std::shared_ptr< AstrometryTransform > const getTangentPlaneToCommonTangentPlane() const
Definition: CcdImage.h:136
CcdIdType getCcdId() const
returns ccd ID
Definition: CcdImage.h:145
double getSinEta() const
Parallactic angle.
Definition: CcdImage.h:173
A list of MeasuredStar. They are usually filled in Associations::createCcdImage.
Definition: MeasuredStar.h:146
table::Key< table::Array< std::uint8_t > > wcs
Definition: SkyWcs.cc:71
jointcal::Point const & getCommonTangentPoint() const
Gets the common tangent point, shared between all ccdImages.
Definition: CcdImage.h:126
STL class.
void resetCatalogForFit()
Clear the catalog for fitting and set it to a copy of the whole catalog.
Definition: CcdImage.h:99
std::shared_ptr< afw::image::PhotoCalib > getPhotoCalib() const
Return the exposure&#39;s photometric calibration.
Definition: CcdImage.h:161
rectangle with sides parallel to axes.
Definition: Frame.h:38
MeasuredStarList & getCatalogForFit()
Gets the catalog to be used for fitting, which may have been cleaned-up.
Definition: CcdImage.h:95
A base class for image defects.
Custom catalog class for record/table subclasses that are guaranteed to have an ID, and should generally be sorted by that ID.
Definition: fwd.h:63
size_t operator()(lsst::jointcal::CcdImageKey const &key) const
Definition: CcdImage.h:245
Key< U > key
Definition: Schema.cc:281
double getLstObs() const
Definition: CcdImage.h:170
table::Key< int > detector
double getHourAngle() const
Definition: CcdImage.h:168
STL class.
lsst::geom::SpherePoint getBoresightRaDec() const
Gets the boresight RA/Dec.
Definition: CcdImage.h:166
double getAirMass() const
Airmass.
Definition: CcdImage.h:155
Point in an unspecified spherical coordinate system.
Definition: SpherePoint.h:57
double getTanZ() const
Parallactic angle.
Definition: CcdImage.h:179
Point getRefractionVector() const
Definition: CcdImage.h:182
MeasuredStarList const & getCatalogForFit() const
Gets the catalog to be used for fitting, which may have been cleaned-up.
Definition: CcdImage.h:94
std::string getFilter() const
return the CcdImage filter name
Definition: CcdImage.h:185
bool operator!=(CcdImageKey const &right) const
Definition: CcdImage.h:54
std::shared_ptr< AstrometryTransform > const getPixelToCommonTangentPlane() const
Definition: CcdImage.h:128
std::shared_ptr< AstrometryTransform > const getCommonTangentPlaneToTangentPlane() const
Definition: CcdImage.h:132
std::shared_ptr< AstrometryTransform > const getSkyToTangentPlane() const
Definition: CcdImage.h:142
std::list< std::shared_ptr< CcdImage > > CcdImageList
Definition: CcdImage.h:46
int VisitIdType
Definition: CcdImage.h:48
CcdImageKey getHashKey() const
Definition: CcdImage.h:152
std::shared_ptr< afw::cameraGeom::Detector > getDetector() const
Definition: CcdImage.h:150
Handler of an actual image from a single CCD.
Definition: CcdImage.h:64
An integer coordinate rectangle.
Definition: Box.h:55
STL class.
Implementation of the Photometric Calibration class.
MeasuredStarList const & getWholeCatalog() const
Gets the as-read catalog.
Definition: CcdImage.h:86
Frame const & getImageFrame() const
Frame in pixels.
Definition: CcdImage.h:191