LSST Applications g0f08755f38+c89d42e150,g1635faa6d4+b6cf076a36,g1653933729+a8ce1bb630,g1a0ca8cf93+4c08b13bf7,g28da252d5a+f33f8200ef,g29321ee8c0+0187be18b1,g2bbee38e9b+9634bc57db,g2bc492864f+9634bc57db,g2cdde0e794+c2c89b37c4,g3156d2b45e+41e33cbcdc,g347aa1857d+9634bc57db,g35bb328faa+a8ce1bb630,g3a166c0a6a+9634bc57db,g3e281a1b8c+9f2c4e2fc3,g414038480c+077ccc18e7,g41af890bb2+e740673f1a,g5fbc88fb19+17cd334064,g7642f7d749+c89d42e150,g781aacb6e4+a8ce1bb630,g80478fca09+f8b2ab54e1,g82479be7b0+e2bd23ab8b,g858d7b2824+c89d42e150,g9125e01d80+a8ce1bb630,g9726552aa6+10f999ec6a,ga5288a1d22+065360aec4,gacf8899fa4+9553554aa7,gae0086650b+a8ce1bb630,gb58c049af0+d64f4d3760,gbd46683f8f+ac57cbb13d,gc28159a63d+9634bc57db,gcf0d15dbbd+e37acf7834,gda3e153d99+c89d42e150,gda6a2b7d83+e37acf7834,gdaeeff99f8+1711a396fd,ge2409df99d+cb1e6652d6,ge79ae78c31+9634bc57db,gf0baf85859+147a0692ba,gf3967379c6+02b11634a5,w.2024.45
LSST Data Management Base Package
Loading...
Searching...
No Matches
simpleShape.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/eigen.h"
24
25#include "ndarray/pybind11.h"
26
29
30namespace py = pybind11;
31using namespace pybind11::literals;
32
33namespace lsst {
34namespace meas {
35namespace extensions {
36namespace simpleShape {
37
38namespace {
39
40template <typename T, typename PyClass>
41static void declareMoments(PyClass cls) {
42 cls.def_static("computeMoments",
43 (SimpleShapeResult(*)(afw::geom::ellipses::Ellipse const &,
44 afw::image::MaskedImage<T> const &, double)) &
46 "weight"_a, "image"_a, "nSigmaRegion"_a = 3.0);
47}
48
49} // namespace
50
51PYBIND11_MODULE(simpleShape, mod) {
52 py::module::import("lsst.afw.geom");
53 py::module::import("lsst.afw.image");
54 py::module::import("lsst.afw.table");
55 py::module::import("lsst.meas.base");
56
57 /* Module level */
58 py::class_<SimpleShape, std::shared_ptr<SimpleShape>, base::SimpleAlgorithm> clsSimpleShape(
59 mod, "SimpleShape");
60 py::class_<SimpleShapeControl> clsSimpleShapeControl(mod, "SimpleShapeControl");
61 py::class_<SimpleShapeResult> clsSimpleShapeResult(mod, "SimpleShapeResult");
62 py::class_<SimpleShapeResultKey> clsSimpleShapeResultKey(mod, "SimpleShapeResultKey");
63
64 /* Constructors */
65 clsSimpleShape.def(py::init<SimpleShape::Control const &, std::string const &, afw::table::Schema &>(),
66 "ctrl"_a, "name"_a, "schema"_a);
67 clsSimpleShapeControl.def(py::init<>());
68
69 clsSimpleShapeResultKey.def(py::init<afw::table::SubSchema const &>(), "s"_a);
70
71 /* Members */
73 LSST_DECLARE_CONTROL_FIELD(clsSimpleShapeControl, SimpleShapeControl, nSigmaRegion);
74
75 declareMoments<float>(clsSimpleShape);
76 declareMoments<double>(clsSimpleShape);
77 clsSimpleShape.def_static("correctWeightedMoments", &SimpleShape::correctWeightedMoments, "weight"_a,
78 "ellipse"_a, "center"_a);
79
80 clsSimpleShapeResult.def_readwrite("ellipse", &SimpleShapeResult::ellipse);
81 clsSimpleShapeResult.def_readwrite("center", &SimpleShapeResult::center);
82 clsSimpleShapeResult.def_readwrite("covariance", &SimpleShapeResult::covariance);
83 clsSimpleShapeResult.def("getFlag", &SimpleShapeResult::getFlag);
84
85 clsSimpleShapeResultKey.def("get", &SimpleShapeResultKey::get, "record"_a);
86}
87
88} // namespace simpleShape
89} // namespace extensions
90} // namespace meas
91} // namespace lsst
afw::table::Key< double > sigma
An abstract base classes for which the same implementation can be used for both SingleFrameAlgorithm ...
Definition Algorithm.h:170
A C++ control class to handle SdssShapeAlgorithm's configuration.
Definition simpleShape.h:44
static Eigen::Matrix< double, 5, 5 > correctWeightedMoments(afw::geom::ellipses::Quadrupole const &weight, afw::geom::ellipses::Quadrupole &ellipse, geom::Point2D &center)
Correct moments measured with a Gaussian weight function by assuming the data was also an elliptical ...
static SimpleShapeResult computeMoments(afw::geom::ellipses::Ellipse const &weight, afw::image::MaskedImage< T > const &image, double nSigmaRegion=3.0)
Compute the Gaussian-weighted moments of an image.
geom::Point2D center
Measured first moments, or the input center if !recentroid.
Eigen::Matrix< double, 5, 5 > covariance
Matrix of uncertainties; ordered Ixx, Iyy, Ixy, Ix, Iy.
afw::geom::ellipses::Quadrupole ellipse
Measured second moments.
virtual SimpleShapeResult get(afw::table::BaseRecord const &record) const
#define LSST_DECLARE_CONTROL_FIELD(WRAPPER, CLASS, NAME)
Macro used to wrap fields declared by LSST_CONTROL_FIELD using Pybind11.
Definition python.h:50