LSST Applications 27.0.0,g0265f82a02+469cd937ee,g02d81e74bb+21ad69e7e1,g1470d8bcf6+cbe83ee85a,g2079a07aa2+e67c6346a6,g212a7c68fe+04a9158687,g2305ad1205+94392ce272,g295015adf3+81dd352a9d,g2bbee38e9b+469cd937ee,g337abbeb29+469cd937ee,g3939d97d7f+72a9f7b576,g487adcacf7+71499e7cba,g50ff169b8f+5929b3527e,g52b1c1532d+a6fc98d2e7,g591dd9f2cf+df404f777f,g5a732f18d5+be83d3ecdb,g64a986408d+21ad69e7e1,g858d7b2824+21ad69e7e1,g8a8a8dda67+a6fc98d2e7,g99cad8db69+f62e5b0af5,g9ddcbc5298+d4bad12328,ga1e77700b3+9c366c4306,ga8c6da7877+71e4819109,gb0e22166c9+25ba2f69a1,gb6a65358fc+469cd937ee,gbb8dafda3b+69d3c0e320,gc07e1c2157+a98bf949bb,gc120e1dc64+615ec43309,gc28159a63d+469cd937ee,gcf0d15dbbd+72a9f7b576,gdaeeff99f8+a38ce5ea23,ge6526c86ff+3a7c1ac5f1,ge79ae78c31+469cd937ee,gee10cc3b42+a6fc98d2e7,gf1cff7945b+21ad69e7e1,gfbcc870c63+9a11dc8c8f
LSST Data Management Base Package
Loading...
Searching...
No Matches
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
37#include "lsst/geom/Box.h"
41#include "lsst/jointcal/Frame.h"
42
43namespace lsst {
44namespace jointcal {
45
47
48using VisitIdType = long;
49using CcdIdType = int;
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;
58std::ostream &operator<<(std::ostream &out, CcdImageKey const &key);
59
64class CcdImage {
65public:
69 std::shared_ptr<afw::cameraGeom::Detector> detector, int visit, int ccd,
70 std::string const &fluxField);
71
73 CcdImage(CcdImage const &) = delete;
74 CcdImage(CcdImage &&) = delete;
75 CcdImage &operator=(CcdImage const &) = 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
103
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
157 double getEpoch() const { return _epoch; }
158
161
165 lsst::geom::SpherePoint getBoresightRaDec() const { return _boresightRaDec; }
166
167 double getHourAngle() const { return _hourAngle; }
168
169 double getLstObs() const { return _lstObs; }
170
172 double getSinEta() const { return _sinEta; }
173
175 double getCosEta() const { return _cosEta; }
176
178 double getTanZ() const { return _tanZ; }
179
181 Point getRefractionVector() const { return Point(_tanZ * _cosEta, _tanZ * _sinEta); }
182
184 std::string getFilter() const { return _filter; }
185
187 std::shared_ptr<AstrometryTransform> const getReadWcs() const { return _readWcs; }
188
190 Frame const &getImageFrame() const { return _imageFrame; }
191
192private:
194 std::string const &fluxField);
195
196 jointcal::Frame _imageFrame; // in pixels
197
198 MeasuredStarList _wholeCatalog; // the catalog of measured objets
199 MeasuredStarList _catalogForFit;
200
201 std::shared_ptr<AstrometryTransformSkyWcs> _readWcs; // apply goes from pix to sky
202
203 // The following ones should probably be mostly removed.
204 // go from CommonTangentPlane to this tangent plane.
205 std::shared_ptr<AstrometryTransform> _commonTangentPlaneToTangentPlane;
206 std::shared_ptr<AstrometryTransform> _tangentPlaneToCommonTangentPlane; // reverse one
207 std::shared_ptr<AstrometryTransform> _pixelToCommonTangentPlane; // pixels -> CTP
208 std::shared_ptr<AstrometryTransform> _pixelToTangentPlane;
209
210 std::shared_ptr<AstrometryTransform> _skyToTangentPlane;
211
212 std::string _name;
213 CcdIdType _ccdId;
214 VisitIdType _visit;
215
216 lsst::geom::SpherePoint _boresightRaDec;
217 double _airMass; // airmass value.
218 double _epoch; // julian epoch year (e.g. 2000.0 for J2000)
221 // refraction
222 // eta : parallactic angle, z: zenithal angle (X = 1/cos(z))
223 double _sinEta, _cosEta, _tanZ;
224 // Local Sidereal Time and hour angle of observation
225 double _lstObs, _hourAngle;
226
227 std::string _filter;
228
229 jointcal::Point _commonTangentPoint;
230};
231} // namespace jointcal
232} // namespace lsst
233
234// Add our preferred hash of CcdImageKey to the std:: namespace, so it's always available "for free".
235namespace std {
236template <>
244 size_t operator()(lsst::jointcal::CcdImageKey const &key) const {
245 return hash<size_t>()(static_cast<size_t>(key.visit) | (static_cast<size_t>(key.ccd) << 32));
246 }
247};
248} // namespace std
249
250#endif // LSST_JOINTCAL_CCD_IMAGE_H
AmpInfoBoxKey bbox
Definition Amplifier.cc:117
Implementation of the Photometric Calibration class.
Tag types used to declare specialized field types.
Definition misc.h:31
An integer coordinate rectangle.
Definition Box.h:55
Point in an unspecified spherical coordinate system.
Definition SpherePoint.h:57
Handler of an actual image from a single CCD.
Definition CcdImage.h:64
std::shared_ptr< AstrometryTransform > const getPixelToCommonTangentPlane() const
Definition CcdImage.h:128
CcdIdType getCcdId() const
returns ccd ID
Definition CcdImage.h:145
std::shared_ptr< afw::cameraGeom::Detector > getDetector() const
Definition CcdImage.h:150
VisitIdType getVisit() const
returns visit ID
Definition CcdImage.h:148
void setCommonTangentPoint(Point const &commonTangentPoint)
Sets the common tangent point and computes necessary transforms.
Definition CcdImage.cc:170
MeasuredStarList & getCatalogForFit()
Definition CcdImage.h:95
std::shared_ptr< AstrometryTransform > const getPixelToTangentPlane() const
Definition CcdImage.h:140
double getEpoch() const
Definition CcdImage.h:157
std::pair< int, int > countStars() const
Count the number of valid measured and reference stars that fall within this ccdImage.
Definition CcdImage.cc:155
std::shared_ptr< AstrometryTransform > const getReadWcs() const
the wcs read in the header. NOT updated when fitting.
Definition CcdImage.h:187
Frame const & getImageFrame() const
Frame in pixels.
Definition CcdImage.h:190
MeasuredStarList const & getWholeCatalog() const
Gets the as-read catalog.
Definition CcdImage.h:86
std::shared_ptr< AstrometryTransform > const getTangentPlaneToCommonTangentPlane() const
Definition CcdImage.h:136
std::shared_ptr< AstrometryTransform > const getCommonTangentPlaneToTangentPlane() const
Definition CcdImage.h:132
MeasuredStarList const & getCatalogForFit() const
Gets the catalog to be used for fitting, which may have been cleaned-up.
Definition CcdImage.h:94
CcdImage & operator=(CcdImage &&)=delete
void resetCatalogForFit()
Clear the catalog for fitting and set it to a copy of the whole catalog.
Definition CcdImage.h:99
CcdImage(CcdImage const &)=delete
No move or copy: each CCD image is unique to that ccd+visit, and Associations holds all CcdImages.
double getSinEta() const
Parallactic angle.
Definition CcdImage.h:172
Point getRefractionVector() const
Definition CcdImage.h:181
double getHourAngle() const
Definition CcdImage.h:167
std::shared_ptr< afw::image::PhotoCalib > getPhotoCalib() const
Return the exposure's photometric calibration.
Definition CcdImage.h:160
CcdImage(afw::table::SourceCatalog &record, const std::shared_ptr< lsst::afw::geom::SkyWcs > &wcs, const std::shared_ptr< lsst::afw::image::VisitInfo > &visitInfo, lsst::geom::Box2I const &bbox, std::string filter, std::shared_ptr< afw::image::PhotoCalib > photoCalib, std::shared_ptr< afw::cameraGeom::Detector > detector, int visit, int ccd, std::string const &fluxField)
Definition CcdImage.cc:113
std::string getFilter() const
return the CcdImage filter name
Definition CcdImage.h:184
CcdImage(CcdImage &&)=delete
lsst::geom::SpherePoint getBoresightRaDec() const
Gets the boresight RA/Dec.
Definition CcdImage.h:165
jointcal::Point const & getCommonTangentPoint() const
Gets the common tangent point, shared between all ccdImages.
Definition CcdImage.h:126
std::string getName() const
Return the _name that identifies this ccdImage.
Definition CcdImage.h:79
double getLstObs() const
Definition CcdImage.h:169
double getTanZ() const
Parallactic angle.
Definition CcdImage.h:178
double getCosEta() const
Parallactic angle.
Definition CcdImage.h:175
CcdImage & operator=(CcdImage const &)=delete
CcdImageKey getHashKey() const
Definition CcdImage.h:152
double getAirMass() const
Airmass.
Definition CcdImage.h:155
std::shared_ptr< AstrometryTransform > const getSkyToTangentPlane() const
Definition CcdImage.h:142
rectangle with sides parallel to axes.
Definition Frame.h:38
A list of MeasuredStar. They are usually filled in Associations::createCcdImage.
A point in a plane.
Definition Point.h:37
void copyTo(StarList< Star > &copy) const
clears copy and makes a copy of the std::list to copy
Definition StarList.cc:68
T clear(T... args)
T right(T... args)
long VisitIdType
Definition CcdImage.h:48
std::ostream & operator<<(std::ostream &stream, AstrometryMapping const &mapping)
STL namespace.
For hashing a ccdImage: the pair of (visit, ccd) IDs should be unique to each ccdImage.
Definition CcdImage.h:51
bool operator==(CcdImageKey const &right) const
Definition CcdImage.h:55
bool operator!=(CcdImageKey const &right) const
Definition CcdImage.h:54
size_t operator()(lsst::jointcal::CcdImageKey const &key) const
Definition CcdImage.h:244
Key< int > visitInfo
Definition Exposure.cc:70