LSST Applications  21.0.0-147-g0e635eb1+1acddb5be5,22.0.0+052faf71bd,22.0.0+1ea9a8b2b2,22.0.0+6312710a6c,22.0.0+729191ecac,22.0.0+7589c3a021,22.0.0+9f079a9461,22.0.1-1-g7d6de66+b8044ec9de,22.0.1-1-g87000a6+536b1ee016,22.0.1-1-g8e32f31+6312710a6c,22.0.1-10-gd060f87+016f7cdc03,22.0.1-12-g9c3108e+df145f6f68,22.0.1-16-g314fa6d+c825727ab8,22.0.1-19-g93a5c75+d23f2fb6d8,22.0.1-19-gb93eaa13+aab3ef7709,22.0.1-2-g8ef0a89+b8044ec9de,22.0.1-2-g92698f7+9f079a9461,22.0.1-2-ga9b0f51+052faf71bd,22.0.1-2-gac51dbf+052faf71bd,22.0.1-2-gb66926d+6312710a6c,22.0.1-2-gcb770ba+09e3807989,22.0.1-20-g32debb5+b8044ec9de,22.0.1-23-gc2439a9a+fb0756638e,22.0.1-3-g496fd5d+09117f784f,22.0.1-3-g59f966b+1e6ba2c031,22.0.1-3-g849a1b8+f8b568069f,22.0.1-3-gaaec9c0+c5c846a8b1,22.0.1-32-g5ddfab5d3+60ce4897b0,22.0.1-4-g037fbe1+64e601228d,22.0.1-4-g8623105+b8044ec9de,22.0.1-5-g096abc9+d18c45d440,22.0.1-5-g15c806e+57f5c03693,22.0.1-7-gba73697+57f5c03693,master-g6e05de7fdc+c1283a92b8,master-g72cdda8301+729191ecac,w.2021.39
LSST Data Management Base Package
matchOptimisticB.h
Go to the documentation of this file.
1 // -*- lsst-c++ -*-
2 #if !defined(LSST_MEAS_ASTROM_MATCHOPTIMISTICB_H)
3 #define LSST_MEAS_ASTROM_MATCHOPTIMISTICB_H
4 
5 #include <cmath>
6 #include <string>
7 #include <vector>
8 
9 #include "lsst/pex/config.h"
10 #include "lsst/geom/Point.h"
11 #include "lsst/afw/table/Source.h"
12 #include "lsst/afw/table/Match.h"
13 
14 namespace lsst {
15 namespace meas {
16 namespace astrom {
17 
22 struct RecordProxy {
25  mutable bool used; // set true if this star object has already been used for a match
26  // mutable to allow freezing the fundamental object data using const
27  // while keeping track of which objects have been used in each search
28 
29  double getX() const { return position.getX(); }
30  double getY() const { return position.getY(); }
31 
36 
37  bool operator==(RecordProxy const& other) const { return record == other.record; }
38  bool operator!=(RecordProxy const& other) const { return record != other.record; }
39 
48 
49  explicit RecordProxy() {} // default constructor needed so we can call ProxyVector::resize()
50 };
51 
53 
54 ProxyVector makeProxies(afw::table::SourceCatalog const& sourceCat, afw::geom::SkyWcs const& distortedWcs,
55  afw::geom::SkyWcs const& tanWcs);
56 
58 
59 struct ProxyPair {
62  double distance;
63  double pa;
64 
65  ProxyPair(RecordProxy const& s1, RecordProxy const& s2) : first(s1), second(s2) {
66  double x1 = first.position.getX();
67  double y1 = first.position.getY();
68  double x2 = second.position.getX();
69  double y2 = second.position.getY();
70  distance = std::hypot(x2 - x1, y2 - y1);
71  pa = std::atan2(y2 - y1, x2 - x1);
72  }
73 };
74 
76  LSST_CONTROL_FIELD(refFluxField, std::string, "name of flux field in reference catalog");
77  LSST_CONTROL_FIELD(sourceFluxField, std::string, "name of flux field in source catalog");
78  LSST_CONTROL_FIELD(numBrightStars, int, "maximum number of bright reference stars to use");
79  LSST_CONTROL_FIELD(minMatchedPairs, int, "minimum number of matches");
81  "maximum allowed distance between reference objects and sources (pixels)");
82  LSST_CONTROL_FIELD(maxOffsetPix, double, "maximum allowed frame translation (pixels)");
83  LSST_CONTROL_FIELD(maxRotationDeg, double, "maximum allowed frame rotation (deg)");
84  LSST_CONTROL_FIELD(allowedNonperpDeg, double, "allowed non-perpendicularity of x and y axes (deg)");
85  LSST_CONTROL_FIELD(numPointsForShape, int, "number of points in a matching shape");
87 
89  : refFluxField("r_flux"),
90  sourceFluxField("slot_ApFlux_instFlux"),
91  numBrightStars(100),
92  minMatchedPairs(50),
94  maxOffsetPix(300),
95  maxRotationDeg(1.0),
96  allowedNonperpDeg(3.0),
98  maxDeterminant(0.02) {
99  validate();
100  }
101 
102  void validate() const;
103 
105 };
106 
127  afw::table::SourceCatalog const& sourceCat,
128  MatchOptimisticBControl const& control,
129  afw::geom::SkyWcs const& wcs, int posRefBegInd = 0,
130  bool verbose = false);
131 
132 } // namespace astrom
133 } // namespace meas
134 } // namespace lsst
135 
136 #endif
table::Key< table::Array< std::uint8_t > > wcs
Definition: SkyWcs.cc:66
T atan2(T... args)
A 2-dimensional celestial WCS that transform pixels to ICRS RA/Dec, using the LSST standard for pixel...
Definition: SkyWcs.h:117
Custom catalog class for record/table subclasses that are guaranteed to have an ID,...
Definition: SortedCatalog.h:42
#define LSST_CONTROL_FIELD(NAME, TYPE, DOC)
A preprocessor macro used to define fields in C++ "control object" structs.
Definition: config.h:43
T hypot(T... args)
std::vector< RecordProxy > ProxyVector
afw::table::ReferenceMatchVector matchOptimisticB(afw::table::SimpleCatalog const &posRefCat, afw::table::SourceCatalog const &sourceCat, MatchOptimisticBControl const &control, afw::geom::SkyWcs const &wcs, int posRefBegInd=0, bool verbose=false)
Match sources to stars in a position reference catalog using optimistic pattern matching B.
ProxyVector makeProxies(afw::table::SourceCatalog const &sourceCat, afw::geom::SkyWcs const &distortedWcs, afw::geom::SkyWcs const &tanWcs)
A base class for image defects.
std::string refFluxField
"name of flux field in reference catalog" ;
double maxRotationDeg
"maximum allowed frame rotation (deg)" ;
double matchingAllowancePix
"maximum allowed distance between reference objects and sources (pixels)" ;
int numPointsForShape
"number of points in a matching shape" ;
int minMatchedPairs
"minimum number of matches" ;
int numBrightStars
"maximum number of bright reference stars to use" ;
std::string sourceFluxField
"name of flux field in source catalog" ;
double maxOffsetPix
"maximum allowed frame translation (pixels)" ;
double allowedNonperpDeg
"allowed non-perpendicularity of x and y axes (deg)" ;
ProxyPair(RecordProxy const &s1, RecordProxy const &s2)
A wrapper around a SimpleRecord or SourceRecord that allows us to record a pixel position in a way th...
bool operator==(RecordProxy const &other) const
bool operator!=(RecordProxy const &other) const
RecordProxy(std::shared_ptr< afw::table::SimpleRecord > record, geom::Point2D const &position)
Construct a RecordProxy.
std::shared_ptr< afw::table::SimpleRecord > record