LSSTApplications  21.0.0+1b62c9342b,21.0.0+45a059f35e,21.0.0-1-ga51b5d4+ceb9cf20a3,21.0.0-2-g103fe59+4d13aff7ba,21.0.0-2-g1367e85+571a348718,21.0.0-2-g2909d54+45a059f35e,21.0.0-2-g45278ab+1b62c9342b,21.0.0-2-g4bc9b9f+ebfe466dad,21.0.0-2-g5242d73+571a348718,21.0.0-2-g54e2caa+07cebfb09d,21.0.0-2-g66bcc37+0b2c5d3971,21.0.0-2-g7f82c8f+08f1f55c36,21.0.0-2-g8dde007+5d1b9cb3f5,21.0.0-2-g8f08a60+73884b2cf5,21.0.0-2-ga326454+08f1f55c36,21.0.0-2-ga63a54e+458e82fbcd,21.0.0-2-gc738bc1+8c4731df06,21.0.0-2-gde069b7+5a8f2956b8,21.0.0-2-ge17e5af+571a348718,21.0.0-2-ge712728+cfa36ee5f9,21.0.0-2-gecfae73+e597808034,21.0.0-2-gfc62afb+571a348718,21.0.0-20-g4449a12+6d1341e0f3,21.0.0-22-gf0532904+1cd928f0c5,21.0.0-3-g4c5b185+c3794955c6,21.0.0-3-g6d51c4a+0b2c5d3971,21.0.0-3-g8076721+5adeb471db,21.0.0-3-gaa929c8+01f4b7cfca,21.0.0-3-gd222c45+afc8332dbe,21.0.0-4-g1383c07+0b2c5d3971,21.0.0-4-g3300ddd+1b62c9342b,21.0.0-4-g5873dc9+9a92674037,21.0.0-4-g8a80011+bd904b6426,21.0.0-5-gcff38f6+844b7f7b93,21.0.0-6-g463d161+18af5fb57b,21.0.0-6-gd3283ba+01f4b7cfca,21.0.0-8-g19111d86+8234efb485,21.0.0-9-g7bed000b9+c7d3cce47e,w.2021.04
LSSTDataManagementBasePackage
matchOptimisticB.cc
Go to the documentation of this file.
1 /*
2  * LSST Data Management System
3  *
4  * This product includes software developed by the
5  * LSST Project (http://www.lsst.org/).
6  * See the COPYRIGHT file
7  *
8  * This program is free software: you can redistribute it and/or modify
9  * it under the terms of the GNU General Public License as published by
10  * the Free Software Foundation, either version 3 of the License, or
11  * (at your option) any later version.
12  *
13  * This program is distributed in the hope that it will be useful,
14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16  * GNU General Public License for more details.
17  *
18  * You should have received a copy of the LSST License Statement and
19  * the GNU General Public License along with this program. If not,
20  * see <https://www.lsstcorp.org/LegalNotices/>.
21  */
22 #include "pybind11/pybind11.h"
23 #include "pybind11/stl.h"
24 
25 #include <memory>
26 
27 #include "lsst/geom/Point.h"
28 #include "lsst/afw/geom/SkyWcs.h"
29 #include "lsst/pex/config/python.h" // defines LSST_DECLARE_CONTROL_FIELD
31 
32 namespace py = pybind11;
33 using namespace pybind11::literals;
34 
35 namespace lsst {
36 namespace meas {
37 namespace astrom {
38 namespace {
39 
40 static void declareRecordProxy(py::module &mod) {
41  py::class_<RecordProxy, std::shared_ptr<RecordProxy>> cls(mod, "RecordProxy");
42 
43  cls.def_readwrite("record", &RecordProxy::record);
44  cls.def_readwrite("position", &RecordProxy::position);
45  cls.def_readwrite("used", &RecordProxy::used);
46 
48  "position"_a);
49 
50  // TO DO: decide if we need to wrap operator PTR(lsst::afw::table::SimpleRecord)()
51 
52  cls.def("__eq__", &RecordProxy::operator==, py::is_operator());
53  cls.def("__ne__", &RecordProxy::operator!=, py::is_operator());
54 
55  cls.def("getX", &RecordProxy::getX);
56  cls.def("getY", &RecordProxy::getY);
57 }
58 
59 static void declareProxyPair(py::module &mod) {
60  py::class_<ProxyPair, std::shared_ptr<ProxyPair>> cls(mod, "ProxyPair");
61 
62  cls.def_readwrite("first", &ProxyPair::first);
63  cls.def_readwrite("second", &ProxyPair::second);
64  cls.def_readwrite("distance", &ProxyPair::distance);
65  cls.def_readwrite("pa", &ProxyPair::pa);
66 
67  cls.def(py::init<RecordProxy const &, RecordProxy const &>(), "s1"_a, "s2"_a);
68 }
69 
70 static void declareMatchOptimisticBControl(py::module &mod) {
71  py::class_<MatchOptimisticBControl> cls(mod, "MatchOptimisticBControl");
72 
73  cls.def(py::init<>());
74 
75  LSST_DECLARE_CONTROL_FIELD(cls, MatchOptimisticBControl, refFluxField);
76  LSST_DECLARE_CONTROL_FIELD(cls, MatchOptimisticBControl, sourceFluxField);
77  LSST_DECLARE_CONTROL_FIELD(cls, MatchOptimisticBControl, numBrightStars);
78  LSST_DECLARE_CONTROL_FIELD(cls, MatchOptimisticBControl, minMatchedPairs);
79  LSST_DECLARE_CONTROL_FIELD(cls, MatchOptimisticBControl, matchingAllowancePix);
80  LSST_DECLARE_CONTROL_FIELD(cls, MatchOptimisticBControl, maxOffsetPix);
81  LSST_DECLARE_CONTROL_FIELD(cls, MatchOptimisticBControl, maxRotationDeg);
82  LSST_DECLARE_CONTROL_FIELD(cls, MatchOptimisticBControl, allowedNonperpDeg);
83  LSST_DECLARE_CONTROL_FIELD(cls, MatchOptimisticBControl, numPointsForShape);
84  LSST_DECLARE_CONTROL_FIELD(cls, MatchOptimisticBControl, maxDeterminant);
85 
86  cls.def("validate", &MatchOptimisticBControl::validate);
87 }
88 
89 } // namespace
90 
92  declareRecordProxy(mod);
93  declareProxyPair(mod);
94  declareMatchOptimisticBControl(mod);
95 
96  mod.def("makeProxies",
98  afw::geom::SkyWcs const &)) &
100  "sourceCat"_a, "distortedWcs"_a, "tanWcs"_a);
101  mod.def("makeProxies",
103  "posRefCat"_a, "tanWcs"_a);
104 
105  mod.def("matchOptimisticB", &matchOptimisticB, "posRefCat"_a, "sourceCat"_a, "control"_a, "wcs"_a,
106  "posRefBegInd"_a = 0, "verbose"_a = false);
107 }
108 
109 } // namespace astrom
110 } // namespace meas
111 } // namespace lsst
lsst::meas::astrom::PYBIND11_MODULE
PYBIND11_MODULE(matchOptimisticB, mod)
Definition: matchOptimisticB.cc:91
std::shared_ptr
STL class.
lsst::afw::table._match.second
second
Definition: _match.py:76
lsst::meas::astrom::makeProxies
ProxyVector makeProxies(afw::table::SourceCatalog const &sourceCat, afw::geom::SkyWcs const &distortedWcs, afw::geom::SkyWcs const &tanWcs)
Definition: matchOptimisticB.cc:436
lsst::afw::table._match.first
first
Definition: _match.py:74
SkyWcs.h
lsst::afw::geom::SkyWcs
A 2-dimensional celestial WCS that transform pixels to ICRS RA/Dec, using the LSST standard for pixel...
Definition: SkyWcs.h:117
LSST_DECLARE_CONTROL_FIELD
#define LSST_DECLARE_CONTROL_FIELD(WRAPPER, CLASS, NAME)
Macro used to wrap fields declared by LSST_CONTROL_FIELD using Pybind11.
Definition: python.h:50
lsst::afw::geom.transform.transformContinued.cls
cls
Definition: transformContinued.py:33
lsst::afw::table._source.SourceCatalog
Definition: _source.py:32
lsst::meas::astrom::matchOptimisticB
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.
Definition: matchOptimisticB.cc:459
lsst::afw::table._match.distance
distance
Definition: _match.py:78
lsst
A base class for image defects.
Definition: imageAlgorithm.dox:1
matchOptimisticB.h
lsst::geom::Point< double, 2 >
Point.h
pybind11
Definition: _GenericMap.cc:40
lsst::utils.tests.init
def init()
Definition: tests.py:59
lsst.pex.config.wrap.validate
validate
Definition: wrap.py:295
python.h
lsst::meas::modelfit.psf.psfContinued.module
module
Definition: psfContinued.py:42
lsst::afw::table::SortedCatalogT
Custom catalog class for record/table subclasses that are guaranteed to have an ID,...
Definition: SortedCatalog.h:42
lsst::meas::astrom::ProxyVector
std::vector< RecordProxy > ProxyVector
Definition: matchOptimisticB.h:52