LSST Applications g063fba187b+fee0456c91,g0f08755f38+ea96e5a5a3,g1653933729+a8ce1bb630,g168dd56ebc+a8ce1bb630,g1a2382251a+90257ff92a,g20f6ffc8e0+ea96e5a5a3,g217e2c1bcf+937a289c59,g28da252d5a+daa7da44eb,g2bbee38e9b+253935c60e,g2bc492864f+253935c60e,g3156d2b45e+6e55a43351,g32e5bea42b+31359a2a7a,g347aa1857d+253935c60e,g35bb328faa+a8ce1bb630,g3a166c0a6a+253935c60e,g3b1af351f3+a8ce1bb630,g3e281a1b8c+c5dd892a6c,g414038480c+416496e02f,g41af890bb2+afe91b1188,g599934f4f4+0db33f7991,g7af13505b9+e36de7bce6,g80478fca09+da231ba887,g82479be7b0+a4516e59e3,g858d7b2824+ea96e5a5a3,g89c8672015+f4add4ffd5,g9125e01d80+a8ce1bb630,ga5288a1d22+bc6ab8dfbd,gb58c049af0+d64f4d3760,gc28159a63d+253935c60e,gcab2d0539d+3f2b72788c,gcf0d15dbbd+4ea9c45075,gda6a2b7d83+4ea9c45075,gdaeeff99f8+1711a396fd,ge79ae78c31+253935c60e,gef2f8181fd+3031e3cf99,gf0baf85859+c1f95f4921,gfa517265be+ea96e5a5a3,gfa999e8aa5+17cd334064,w.2024.50
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