LSSTApplications  20.0.0
LSSTDataManagementBasePackage
ConstrainedAstrometryModel.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_CONSTRAINED_ASTROMETRY_MODEL_H
26 #define LSST_JOINTCAL_CONSTRAINED_ASTROMETRY_MODEL_H
27 
28 #include "memory" // for std::*_ptr
29 
31 
32 class CcdImage;
33 
36 #include "lsst/jointcal/Frame.h"
40 #include "lsst/jointcal/CcdImage.h"
41 
42 #include <map>
43 
44 namespace lsst {
45 namespace jointcal {
46 
63 public:
64  ConstrainedAstrometryModel(CcdImageList const &ccdImageList,
65  std::shared_ptr<ProjectionHandler const> projectionHandler, int chipOrder,
66  int visitOrder);
67 
73 
74  // The following routines are the interface to AstrometryFit
76  AstrometryMapping const *getMapping(CcdImage const &) const override;
77 
82  Eigen::Index assignIndices(std::string const &whatToFit, Eigen::Index firstIndex) override;
83 
88  void offsetParams(Eigen::VectorXd const &Delta) override;
89 
94  void freezeErrorTransform() override;
95 
97  std::size_t getTotalParameters() const override;
98 
100  AstrometryTransform const &getChipTransform(CcdIdType const chip) const;
101 
104 
107 
114  CcdImage const &ccdImage) const override {
115  return _skyToTangentPlane->getSkyToTangentPlane(ccdImage);
116  }
117 
119  std::shared_ptr<afw::geom::SkyWcs> makeSkyWcs(CcdImage const &ccdImage) const override;
120 
121  void print(std::ostream &out) const override;
122 
123 private:
127  const std::shared_ptr<ProjectionHandler const> _skyToTangentPlane;
128  bool _fittingChips, _fittingVisits;
129 
131  AstrometryMapping *findMapping(CcdImage const &ccdImage) const override;
132 };
133 } // namespace jointcal
134 } // namespace lsst
135 
136 #endif // LSST_JOINTCAL_CONSTRAINED_ASTROMETRY_MODEL_H
std::string
STL class.
std::shared_ptr
STL class.
lsst::jointcal::ConstrainedAstrometryModel::ConstrainedAstrometryModel
ConstrainedAstrometryModel(CcdImageList const &ccdImageList, std::shared_ptr< ProjectionHandler const > projectionHandler, int chipOrder, int visitOrder)
Definition: ConstrainedAstrometryModel.cc:64
lsst::jointcal::ConstrainedAstrometryModel::getVisits
std::vector< VisitIdType > getVisits() const
Access to array of visits involved in the solution.
Definition: ConstrainedAstrometryModel.cc:202
std::list< std::shared_ptr< CcdImage > >
CcdImage.h
lsst::jointcal::ConstrainedAstrometryModel::getTotalParameters
std::size_t getTotalParameters() const override
Return the total number of parameters in this model.
Definition: ConstrainedAstrometryModel.cc:221
lsst::jointcal::CcdIdType
int CcdIdType
Definition: CcdImage.h:49
AstrometryTransform.h
Eigenstuff.h
Frame.h
std::vector
STL class.
lsst::jointcal::ConstrainedAstrometryModel::freezeErrorTransform
void freezeErrorTransform() override
From there on, measurement errors are propagated using the current transforms (and no longer evolve).
Definition: ConstrainedAstrometryModel.cc:184
lsst::jointcal::ConstrainedAstrometryModel::ConstrainedAstrometryModel
ConstrainedAstrometryModel(ConstrainedAstrometryModel const &)=delete
No copy or move: there is only ever one instance of a given model (i.e. per ccd+visit)
lsst::jointcal::ConstrainedAstrometryModel::assignIndices
Eigen::Index assignIndices(std::string const &whatToFit, Eigen::Index firstIndex) override
Positions the various parameter sets into the parameter vector, starting at firstIndex.
Definition: ConstrainedAstrometryModel.cc:140
lsst::jointcal::AstrometryModel
Interface between AstrometryFit and the combinations of Mappings from pixels to some tangent plane (a...
Definition: AstrometryModel.h:50
ProjectionHandler.h
std::ostream
STL class.
lsst::jointcal::CcdImage
Handler of an actual image from a single CCD.
Definition: CcdImage.h:64
lsst::jointcal::ConstrainedAstrometryModel::getMapping
AstrometryMapping const * getMapping(CcdImage const &) const override
Mapping associated to a given CcdImage.
Definition: ConstrainedAstrometryModel.cc:132
lsst::jointcal::ConstrainedAstrometryModel::getVisitTransform
AstrometryTransform const & getVisitTransform(VisitIdType const &visit) const
Access to mappings.
Definition: ConstrainedAstrometryModel.cc:209
std::map
STL class.
lsst::jointcal
Definition: Associations.h:49
lsst::jointcal::ConstrainedAstrometryModel::ConstrainedAstrometryModel
ConstrainedAstrometryModel(ConstrainedAstrometryModel &&)=delete
lsst
A base class for image defects.
Definition: imageAlgorithm.dox:1
lsst::jointcal::ConstrainedAstrometryModel::operator=
ConstrainedAstrometryModel & operator=(ConstrainedAstrometryModel &&)=delete
AstrometryModel.h
lsst::jointcal::ConstrainedAstrometryModel::makeSkyWcs
std::shared_ptr< afw::geom::SkyWcs > makeSkyWcs(CcdImage const &ccdImage) const override
Make a SkyWcs that contains this model.
Definition: ConstrainedAstrometryModel.cc:232
lsst::jointcal::ConstrainedAstrometryModel::getSkyToTangentPlane
const std::shared_ptr< AstrometryTransform const > getSkyToTangentPlane(CcdImage const &ccdImage) const override
The mapping of sky coordinates (i.e.
Definition: ConstrainedAstrometryModel.h:113
lsst::jointcal::AstrometryMapping
virtual class needed in the abstraction of the distortion model
Definition: AstrometryMapping.h:39
lsst::jointcal::ConstrainedAstrometryModel
A multi-component model, fitting mappings for sensors and visits simultaneously.
Definition: ConstrainedAstrometryModel.h:62
std::size_t
lsst::jointcal::AstrometryTransform
a virtual (interface) class for geometric transformations.
Definition: AstrometryTransform.h:65
lsst.pipe.tasks.assembleCoadd.visit
visit
Definition: assembleCoadd.py:336
lsst::jointcal::VisitIdType
int VisitIdType
Definition: CcdImage.h:48
lsst::jointcal::ConstrainedAstrometryModel::offsetParams
void offsetParams(Eigen::VectorXd const &Delta) override
Dispaches the offsets after a fit step into the actual locations of parameters.
Definition: ConstrainedAstrometryModel.cc:171
std::unordered_map
STL class.
lsst::jointcal::ConstrainedAstrometryModel::print
void print(std::ostream &out) const override
Print a string representation of the contents of this mapping, for debugging.
Definition: ConstrainedAstrometryModel.cc:259
ChipVisitAstrometryMapping.h
lsst::jointcal::ConstrainedAstrometryModel::getChipTransform
AstrometryTransform const & getChipTransform(CcdIdType const chip) const
Access to mappings.
Definition: ConstrainedAstrometryModel.cc:189
lsst::jointcal::ConstrainedAstrometryModel::operator=
ConstrainedAstrometryModel & operator=(ConstrainedAstrometryModel const &)=delete
SimpleAstrometryMapping.h