45 template <
typename KeyType,
typename ValueType>
49 for (
auto const &i : map) {
63 ConstrainedAstrometryModel::ConstrainedAstrometryModel(
65 int chipOrder,
int visitOrder)
67 _skyToTangentPlane(projectionHandler) {
73 for (
auto &
ccdImage : ccdImageList) {
77 auto visitp = _visitMap.find(visit);
78 if (visitp == _visitMap.end()) {
82 auto chipp = _chipMap.find(chip);
83 if (chipp == _chipMap.end()) {
86 if (radius2 < minRadius2) {
88 constrainedChip = chip;
98 _chipMap[chip] = std::make_shared<SimplePolyMapping>(shiftAndNormalize,
99 pol * shiftAndNormalize.inverted());
104 _chipMap.at(constrainedChip)->setToBeFit(
false);
107 for (
auto &
ccdImage : ccdImageList) {
114 if (_chipMap.find(chip) == _chipMap.end()) {
115 LOGLS_WARN(
_log,
"Chip " << chip <<
" is missing in the reference exposure, expect troubles.");
121 std::make_unique<ChipVisitAstrometryMapping>(_chipMap[chip], _visitMap[visit]);
123 LOGLS_INFO(
_log,
"Got " << _chipMap.size() <<
" chip mappings and " << _visitMap.size()
124 <<
" visit mappings; holding chip " << constrainedChip <<
" fixed (" 126 LOGLS_DEBUG(
_log,
"CcdImage map has " << _mappings.size() <<
" mappings, with " 127 << _mappings.bucket_count() <<
" buckets and a load factor of " 128 << _mappings.load_factor());
132 return findMapping(ccdImage);
141 Eigen::Index firstIndex
143 Eigen::Index index = firstIndex;
144 if (whatToFit.
find(
"Distortions") == std::string::npos) {
145 LOGLS_ERROR(
_log,
"assignIndices was called and Distortions is *not* in whatToFit");
149 _fittingChips = (whatToFit.
find(
"DistortionsChip") != std::string::npos);
150 _fittingVisits = (whatToFit.
find(
"DistortionsVisit") != std::string::npos);
152 if ((!_fittingChips) && (!_fittingVisits)) {
153 _fittingChips = _fittingVisits =
true;
156 for (
auto &i : _chipMap) {
157 i.second->setIndex(index);
158 index += i.second->getNpar();
161 for (
auto &i : _visitMap) {
162 i.second->setIndex(index);
163 index += i.second->getNpar();
166 for (
auto &i : _mappings) {
167 i.second->setWhatToFit(_fittingChips, _fittingVisits);
174 for (
auto &i : _chipMap) {
175 auto mapping = i.second.get();
176 mapping->offsetParams(delta.segment(mapping->getIndex(), mapping->getNpar()));
179 for (
auto &i : _visitMap) {
180 auto mapping = i.second.get();
181 mapping->offsetParams(delta.segment(mapping->getIndex(), mapping->getNpar()));
186 for (
auto i = _visitMap.begin(); i != _visitMap.end(); ++i) i->second->freezeErrorTransform();
187 for (
auto i = _chipMap.begin(); i != _chipMap.end(); ++i) i->second->freezeErrorTransform();
191 auto chipp = _chipMap.find(chip);
192 if (chipp == _chipMap.end()) {
194 errMsg <<
"No such chipId: " << chip <<
" among ";
195 outputMapKeys(_chipMap, errMsg);
199 return chipp->second->getTransform();
206 for (
auto i = _visitMap.begin(); i != _visitMap.end(); ++i) res.
push_back(i->first);
211 auto visitp = _visitMap.find(visit);
212 if (visitp == _visitMap.end()) {
214 errMsg <<
"No such visitId: " << visit <<
" among ";
215 outputMapKeys(_visitMap, errMsg);
219 return visitp->second->getTransform();
224 for (
auto &i : _chipMap) {
225 total += i.second->getNpar();
227 for (
auto &i : _visitMap) {
228 total += i.second->getNpar();
248 afw::geom::Point2D(0, 0),
249 afw::geom::SpherePoint(tangentPoint.x, tangentPoint.y, afw::geom::degrees),
251 auto iwcToSkyMap = iwcToSkyWcs->getFrameDict()->getMapping(
"PIXELS",
"SKY");
252 auto skyFrame = iwcToSkyWcs->getFrameDict()->getFrame(
"SKY");
255 frameDict.
addFrame(
"PIXELS", *pixelsToFocal, focalFrame);
256 frameDict.
addFrame(
"FOCAL", *focalToIwc, iwcFrame);
257 frameDict.
addFrame(
"IWC", *iwcToSkyMap, *skyFrame);
258 return std::make_shared<afw::geom::SkyWcs>(frameDict);
262 auto i = _mappings.find(ccdImage.
getHashKey());
263 if (i == _mappings.end())
265 "ConstrainedAstrometryModel cannot find CcdImage " + ccdImage.
getName());
266 return i->second.get();
AstrometryTransform const & getVisitTransform(VisitIdType const &visit) const
Access to mappings.
VisitIdType getVisit() const
returns visit ID
Eigen::Index assignIndices(std::string const &whatToFit, Eigen::Index firstIndex) override
Positions the various parameter sets into the parameter vector, starting at firstIndex.
LOG_LOGGER _log
lsst.logging instance, to be created by a subclass so that messages have consistent name...
std::string getName() const
Return the _name that identifies this ccdImage.
CameraSysPrefix const PIXELS
const std::shared_ptr< AstrometryTransform const > getSkyToTangentPlane(CcdImage const &ccdImage) const override
The mapping of sky coordinates (i.e.
Interface between AstrometryFit and the combinations of Mappings from pixels to some tangent plane (a...
std::vector< VisitIdType > getVisits() const
Access to array of visits involved in the solution.
std::size_t getTotalParameters() const override
Return the total number of parameters in this model.
CcdIdType getCcdId() const
returns ccd ID
AstrometryMapping const * getMapping(CcdImage const &) const override
Mapping associated to a given CcdImage.
#define LOGLS_INFO(logger, message)
void addFrame(int iframe, Mapping const &map, Frame const &frame) override
rectangle with sides parallel to axes.
AstrometryTransformLinear normalizeCoordinatesTransform(const Frame &frame)
Returns the transformation that maps the input frame along both axes to [-1,1].
Class for a simple mapping implementing a generic AstrometryTransform.
std::shared_ptr< SkyWcs > makeSkyWcs(daf::base::PropertySet &metadata, bool strip=false)
std::shared_ptr< afw::geom::SkyWcs > makeSkyWcs(CcdImage const &ccdImage) const override
Make a SkyWcs that contains this model.
void freezeErrorTransform() override
From there on, measurement errors are propagated using the current transforms (and no longer evolve)...
This one is the Tangent Plane (called gnomonic) projection (from celestial sphere to tangent plane) ...
T dynamic_pointer_cast(T... args)
#define LSST_EXCEPT(type,...)
void offsetParams(Eigen::VectorXd const &Delta) override
Dispaches the offsets after a fit step into the actual locations of parameters.
AstrometryTransform const & getChipTransform(CcdIdType const chip) const
Access to mappings.
#define LOGLS_DEBUG(logger, message)
Eigen::Matrix2d makeCdMatrix(lsst::geom::Angle const &scale, lsst::geom::Angle const &orientation=0 *lsst::geom::degrees, bool flipX=false)
virtual class needed in the abstraction of the distortion model
CcdImageKey getHashKey() const
std::shared_ptr< afw::cameraGeom::Detector > getDetector() const
Handler of an actual image from a single CCD.
CameraSys const FOCAL_PLANE
Frame const & getImageFrame() const
Frame in pixels.
#define LOGLS_ERROR(logger, message)
#define LOGLS_WARN(logger, message)