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
StarMatch.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_STAR_MATCH_H
26 #define LSST_JOINTCAL_STAR_MATCH_H
27 
28 #include <iostream>
29 #include <iterator> // for std::ostream_iterator
30 #include <algorithm> // for swap
31 #include <string>
32 #include <list>
33 
34 #include "lsst/jointcal/Point.h"
35 #include "lsst/jointcal/BaseStar.h" // class definition used in inlined functions
36 #include "lsst/jointcal/AstrometryTransform.h" // inlined function calls AstrometryTransform::apply()
37 
38 namespace lsst {
39 namespace jointcal {
40 
52 
54 class StarMatch {
55  friend class StarMatchList;
56 
57 public:
58  /* if one sets that private, then fitting routines will be a nightmare. we could set all the Transform
59  * classes friend... */
63  double distance;
64  double chi2;
65 
67 
69  StarMatch(const FatPoint &point1, const FatPoint &point2, std::shared_ptr<const BaseStar> star1,
71  : point1(point1), point2(point2), s1(std::move(star1)), s2(std::move(star2)), distance(0.){};
72 
73  // the next one would require that StarMatch knows BaseStar which is not mandatory for StarMatch to work
74  // StarMatch(BaseStar *star1, BaseStar *star2) : point1(*star1), point2(*star2), s1(star1), s2(star2) {};
75 
78  return point2.Distance(transform.apply(point1));
79  };
80 
82  double computeChi2(const AstrometryTransform &transform) const;
83 
85  void setDistance(const AstrometryTransform &transform) { distance = computeDistance(transform); };
87  double getDistance() const { return distance; }
88 
89  void swap() {
90  std::swap(point1, point2);
91  std::swap(s1, s2);
92  }
93 
94  /* comparison that ensures that after a sort, duplicates are next one another */
95  explicit StarMatch(){};
96 
97  friend std::ostream &operator<<(std::ostream &stream, const StarMatch &Match);
98 
100 
101 private:
102  // bool operator < (const StarMatch & other) const { return (s1 > other.s1) ? s1 > other.s1 : s2 >
103  // other.s2;}
104 
105  /* for unique to remove duplicates */
106  bool operator==(const StarMatch &other) const { return (s1 == other.s1 && s2 == other.s2); };
107  bool operator!=(const StarMatch &other) const { return (s1 != other.s1 || s2 != other.s2); };
108 
109  friend bool compareStar1(const StarMatch &one, const StarMatch &two);
110  friend bool sameStar1(const StarMatch &one, const StarMatch &two);
111  friend bool compareStar2(const StarMatch &one, const StarMatch &two);
112  friend bool sameStar2(const StarMatch &one, const StarMatch &two);
113 };
114 
115 inline bool compareStar1(const StarMatch &one, const StarMatch &two) {
116  return ((one.s1 == two.s1) ? (one.distance < two.distance) : (&(*one.s1) > &(*two.s1)));
117 }
118 
119 inline bool sameStar1(const StarMatch &one, const StarMatch &two) { return (one.s1 == two.s1); }
120 
121 inline bool compareStar2(const StarMatch &one, const StarMatch &two) {
122  return ((one.s2 == two.s2) ? (one.distance < two.distance) : (&(*one.s2) > &(*two.s2)));
123 }
124 
125 inline bool sameStar2(const StarMatch &one, const StarMatch &two) { return (one.s2 == two.s2); }
126 
127 /* =================================== StarMatchList
128  * ============================================================ */
129 
131 
132 //#ifdef TO_BE_FIXED
133 typedef ::std::list<StarMatch>::iterator StarMatchIterator;
134 typedef ::std::list<StarMatch>::const_iterator StarMatchCIterator;
135 //#endif
136 
138 
147 std::ostream &operator<<(std::ostream &stream, const StarMatchList &starMatchList);
148 
149 class StarMatchList : public std::list<StarMatch> {
150 public:
151  void refineTransform(double nSigmas);
152 
155  void applyTransform(StarMatchList &transformed, const AstrometryTransform *priorTransform,
156  const AstrometryTransform *posteriorTransform = nullptr) const;
157 
158  /* constructor */
159  StarMatchList() : _order(0), _chi2(0){};
160 
162 
165 
167  double getDist2() const { return _dist2; }
168 
170  double getChi2() const { return _chi2; }
171 
173  int getTransformOrder() const { return _order; }
174 
176  void swap();
177 
179  double computeResidual() const;
180 
184  unsigned removeAmbiguities(const AstrometryTransform &transform, int which = 3);
185 
187  void setTransform(const AstrometryTransform *transform) { _transform = transform->clone(); }
189  void setTransform(const AstrometryTransform &transform) { _transform = transform.clone(); }
190  void setTransform(std::shared_ptr<AstrometryTransform> transform) { _transform = std::move(transform); }
191 
193  void setTransformOrder(int order);
194 
197  std::unique_ptr<AstrometryTransform> inverseTransform();
198 
200  void setDistance(const AstrometryTransform &transform);
201 
203  void cutTail(int nKeep);
204 
206  int recoveredNumber(double mindist) const;
207 
209  void dumpTransform(std::ostream &stream = std::cout) const;
210 
211  ~StarMatchList(){/* should delete the transform.... or use counted refs*/};
212 
213 private:
214  StarMatchList(const StarMatchList &); // copies nor properly handled
215  void operator=(const StarMatchList &);
216 
217  int _order;
218  double _chi2;
219  double _dist2;
221 };
222 
225 
227 double computeChi2(const StarMatchList &L, const AstrometryTransform &transform);
228 } // namespace jointcal
229 } // namespace lsst
230 
231 #endif // LSST_JOINTCAL_STAR_MATCH_H
A hanger for star associations.
Definition: StarMatch.h:54
void setTransform(std::shared_ptr< AstrometryTransform > transform)
Definition: StarMatch.h:190
T swap(T... args)
friend class StarMatchList
Definition: StarMatch.h:55
FatPoint point2
2 points
Definition: StarMatch.h:60
double computeChi2(const AstrometryTransform &transform) const
returns the chi2 (using errors in the FatPoint&#39;s)
Definition: StarMatch.cc:44
STL namespace.
void setDistance(const AstrometryTransform &transform)
to be used before sorting on distances.
Definition: StarMatch.h:85
ItemVariant const * other
Definition: Schema.cc:56
double Distance(const Point &other) const
Definition: Point.h:50
A Point with uncertainties.
Definition: FatPoint.h:34
double getDist2() const
access to the sum of squared residuals of the last call to refineTransform.
Definition: StarMatch.h:167
double getChi2() const
access to the chi2 of the last call to refineTransform.
Definition: StarMatch.h:170
A base class for image defects.
::std::list< StarMatch >::const_iterator StarMatchCIterator
Definition: StarMatch.h:134
friend bool sameStar2(const StarMatch &one, const StarMatch &two)
Definition: StarMatch.h:125
StarMatch(const FatPoint &point1, const FatPoint &point2, std::shared_ptr< const BaseStar > star1, std::shared_ptr< const BaseStar > star2)
constructor.
Definition: StarMatch.h:69
friend bool sameStar1(const StarMatch &one, const StarMatch &two)
Definition: StarMatch.h:119
void setTransform(const AstrometryTransform *transform)
sets a transform between the 2 std::lists and deletes the previous or default one. No fit.
Definition: StarMatch.h:187
void setTransform(const AstrometryTransform &transform)
Definition: StarMatch.h:189
::std::list< StarMatch >::iterator StarMatchIterator
Definition: StarMatch.h:133
STL class.
T move(T... args)
friend bool compareStar2(const StarMatch &one, const StarMatch &two)
Definition: StarMatch.h:121
double computeDistance(const AstrometryTransform &transform) const
returns the distance from transform(point1) to point2.
Definition: StarMatch.h:77
std::shared_ptr< const AstrometryTransform > getTransform() const
carries out a fit with outlier rejection
Definition: StarMatch.h:164
STL class.
a virtual (interface) class for geometric transformations.
double computeDist2(const StarMatchList &S, const AstrometryTransform &transform)
sum of distance squared
Definition: StarMatch.cc:239
friend std::ostream & operator<<(std::ostream &stream, const StarMatch &Match)
Definition: StarMatch.cc:56
std::shared_ptr< const BaseStar > s2
Definition: StarMatch.h:62
std::shared_ptr< const BaseStar > s1
Definition: StarMatch.h:62
friend bool compareStar1(const StarMatch &one, const StarMatch &two)
Definition: StarMatch.h:115
STL class.
virtual std::unique_ptr< AstrometryTransform > clone() const =0
returns a copy (allocated by new) of the transformation.
int getTransformOrder() const
returns the order of the used transform
Definition: StarMatch.h:173
double getDistance() const
returns the value computed by the above one.
Definition: StarMatch.h:87
virtual void apply(const double xIn, const double yIn, double &xOut, double &yOut) const =0