LSSTApplications  18.1.0
LSSTDataManagementBasePackage
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 
35  operator PTR(afw::table::SimpleRecord)() const { return record; }
36 
37  bool operator==(RecordProxy const& other) const { return record == other.record; }
38  bool operator!=(RecordProxy const& other) const { return record != other.record; }
39 
47  : record(record), position(position) {}
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 
57 ProxyVector makeProxies(afw::table::SimpleCatalog const& posRefCat, afw::geom::SkyWcs const& tanWcs);
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");
80  LSST_CONTROL_FIELD(matchingAllowancePix, double,
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");
86  LSST_CONTROL_FIELD(maxDeterminant, double, "?");
87 
89  : refFluxField("r_flux"),
90  sourceFluxField("slot_ApFlux_instFlux"),
91  numBrightStars(100),
92  minMatchedPairs(50),
93  matchingAllowancePix(10.0),
94  maxOffsetPix(300),
95  maxRotationDeg(1.0),
96  allowedNonperpDeg(3.0),
97  numPointsForShape(6),
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
bool operator==(RecordProxy const &other) const
A 2-dimensional celestial WCS that transform pixels to ICRS RA/Dec, using the LSST standard for pixel...
Definition: SkyWcs.h:117
#define PTR(...)
Definition: base.h:41
bool operator!=(RecordProxy const &other) const
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...
#define LSST_CONTROL_FIELD(NAME, TYPE, DOC)
A preprocessor macro used to define fields in C++ "control object" structs.
Definition: config.h:36
table::Key< table::Array< std::uint8_t > > wcs
Definition: SkyWcs.cc:71
std::vector< RecordProxy > ProxyVector
ProxyPair(RecordProxy const &s1, RecordProxy const &s2)
T atan2(T... args)
STL class.
A base class for image defects.
RecordProxy(boost::shared_ptr< afw::table::SimpleRecord > record, geom::Point2D const &position)
Construct a RecordProxy.
Custom catalog class for record/table subclasses that are guaranteed to have an ID, and should generally be sorted by that ID.
Definition: fwd.h:63
ProxyVector makeProxies(afw::table::SourceCatalog const &sourceCat, afw::geom::SkyWcs const &distortedWcs, afw::geom::SkyWcs const &tanWcs)
T hypot(T... args)
STL class.
A wrapper around a SimpleRecord or SourceRecord that allows us to record a pixel position in a way th...
ItemVariant const * other
Definition: Schema.cc:56
Record class that must contain a unique ID field and a celestial coordinate field.
Definition: Simple.h:48
boost::shared_ptr< afw::table::SimpleRecord > record