LSSTApplications  20.0.0
LSSTDataManagementBasePackage
SimpleAstrometryModel.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_SIMPLE_ASTROMETRY_MODEL_H
26 #define LSST_JOINTCAL_SIMPLE_ASTROMETRY_MODEL_H
27 
28 #include "memory"
29 
31 
36 #include <map>
37 
38 namespace lsst {
39 namespace jointcal {
40 
41 class CcdImage;
42 
43 /* We deal here with coordinate transforms which are fitted
44  and/or necessary to AstrometryFit. The classes SimpleAstrometryModel
45 and SimplePolyMapping implement a model where there is one
46 separate transfrom per CcdImage. One could chose other setups.
47 
48 */
49 
63 public:
64  SimpleAstrometryModel(CcdImageList const &ccdImageList,
65  const std::shared_ptr<ProjectionHandler const> projectionHandler, bool initFromWCS,
66  unsigned nNotFit = 0, unsigned order = 3);
67 
73 
74  // The following routines are the interface to AstrometryFit
76  const AstrometryMapping *getMapping(CcdImage const &) const override;
77 
79  Eigen::Index assignIndices(std::string const &whatToFit, Eigen::Index firstIndex) override;
80 
81  // dispaches the offsets after a fit step into the actual locations of parameters
82  void offsetParams(Eigen::VectorXd const &delta) override;
83 
88  CcdImage const &ccdImage) const override {
89  return _skyToTangentPlane->getSkyToTangentPlane(ccdImage);
90  }
91 
93  void freezeErrorTransform() override;
94 
96  std::size_t getTotalParameters() const override;
97 
98  void print(std::ostream &out) const override;
99 
101  AstrometryTransform const &getTransform(CcdImage const &ccdImage) const;
102 
104  std::shared_ptr<afw::geom::SkyWcs> makeSkyWcs(CcdImage const &ccdImage) const override;
105 
107 
108 private:
110  const std::shared_ptr<ProjectionHandler const> _skyToTangentPlane;
111 
113  AstrometryMapping *findMapping(CcdImage const &ccdImage) const override;
114 };
115 } // namespace jointcal
116 } // namespace lsst
117 
118 #endif // LSST_JOINTCAL_SIMPLE_ASTROMETRY_MODEL_H
lsst::jointcal::SimpleAstrometryModel::operator=
SimpleAstrometryModel & operator=(SimpleAstrometryModel const &)=delete
std::string
STL class.
std::shared_ptr
STL class.
std::list< std::shared_ptr< CcdImage > >
lsst::jointcal::SimpleAstrometryModel::makeSkyWcs
std::shared_ptr< afw::geom::SkyWcs > makeSkyWcs(CcdImage const &ccdImage) const override
Make a SkyWcs that contains this model.
Definition: SimpleAstrometryModel.cc:150
lsst::jointcal::SimpleAstrometryModel::freezeErrorTransform
void freezeErrorTransform() override
Definition: SimpleAstrometryModel.cc:123
lsst::jointcal::SimpleAstrometryModel::offsetParams
void offsetParams(Eigen::VectorXd const &delta) override
Offset the parameters by the provided amounts (by -delta).
Definition: SimpleAstrometryModel.cc:116
AstrometryTransform.h
Eigenstuff.h
lsst::jointcal::SimpleAstrometryModel::operator=
SimpleAstrometryModel & operator=(SimpleAstrometryModel &&)=delete
lsst::jointcal::SimpleAstrometryModel::SimpleAstrometryModel
SimpleAstrometryModel(CcdImageList const &ccdImageList, const std::shared_ptr< ProjectionHandler const > projectionHandler, bool initFromWCS, unsigned nNotFit=0, unsigned order=3)
Definition: SimpleAstrometryModel.cc:41
lsst::jointcal::SimpleAstrometryModel::SimpleAstrometryModel
SimpleAstrometryModel(SimpleAstrometryModel const &)=delete
No copy or move: there is only ever one instance of a given model (i.e.. per ccd+visit)
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::SimpleAstrometryModel
A model where there is one independent transform per CcdImage.
Definition: SimpleAstrometryModel.h:62
lsst::jointcal::SimpleAstrometryModel::getTotalParameters
std::size_t getTotalParameters() const override
Return the total number of parameters in this model.
Definition: SimpleAstrometryModel.cc:127
lsst::jointcal
Definition: Associations.h:49
lsst
A base class for image defects.
Definition: imageAlgorithm.dox:1
AstrometryModel.h
lsst::jointcal::SimpleAstrometryModel::getMapping
const AstrometryMapping * getMapping(CcdImage const &) const override
Mapping associated to a given CcdImage.
Definition: SimpleAstrometryModel.cc:97
lsst::jointcal::SimpleAstrometryModel::getSkyToTangentPlane
const std::shared_ptr< AstrometryTransform const > getSkyToTangentPlane(CcdImage const &ccdImage) const override
the mapping of sky coordinates (i.e.
Definition: SimpleAstrometryModel.h:87
lsst::jointcal::AstrometryMapping
virtual class needed in the abstraction of the distortion model
Definition: AstrometryMapping.h:39
std::size_t
lsst::jointcal::AstrometryTransform
a virtual (interface) class for geometric transformations.
Definition: AstrometryTransform.h:65
lsst::jointcal::SimpleAstrometryModel::print
void print(std::ostream &out) const override
Print a string representation of the contents of this mapping, for debugging.
Definition: SimpleAstrometryModel.cc:135
std::unordered_map
STL class.
lsst::jointcal::SimpleAstrometryModel::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: SimpleAstrometryModel.cc:101
lsst::jointcal::SimpleAstrometryModel::getTransform
AstrometryTransform const & getTransform(CcdImage const &ccdImage) const
Access to mappings.
Definition: SimpleAstrometryModel.cc:146
lsst::jointcal::SimpleAstrometryModel::~SimpleAstrometryModel
~SimpleAstrometryModel()
Definition: SimpleAstrometryModel.h:106
lsst::jointcal::SimpleAstrometryModel::SimpleAstrometryModel
SimpleAstrometryModel(SimpleAstrometryModel &&)=delete
SimpleAstrometryMapping.h