Loading [MathJax]/extensions/tex2jax.js
LSST Applications g04a91732dc+7fec47d7bc,g07dc498a13+5ab4d22ec3,g0fba68d861+565de8e5d5,g1409bbee79+5ab4d22ec3,g1a7e361dbc+5ab4d22ec3,g1fd858c14a+11200c7927,g20f46db602+25d63fd678,g35bb328faa+fcb1d3bbc8,g4d2262a081+61302e889d,g4d39ba7253+d05e267ece,g4e0f332c67+5d362be553,g53246c7159+fcb1d3bbc8,g60b5630c4e+d05e267ece,g78460c75b0+2f9a1b4bcd,g786e29fd12+cf7ec2a62a,g7b71ed6315+fcb1d3bbc8,g8048e755c2+a1301e4c20,g8852436030+163ceb82d8,g89139ef638+5ab4d22ec3,g89e1512fd8+fbb808ebce,g8d6b6b353c+d05e267ece,g9125e01d80+fcb1d3bbc8,g989de1cb63+5ab4d22ec3,g9f33ca652e+8abe617c77,ga9baa6287d+d05e267ece,gaaedd4e678+5ab4d22ec3,gabe3b4be73+1e0a283bba,gb1101e3267+fefe9ce5b1,gb58c049af0+f03b321e39,gb90eeb9370+824c420ec4,gc741bbaa4f+77ddc33078,gcf25f946ba+163ceb82d8,gd315a588df+0f88d5218e,gd6cbbdb0b4+c8606af20c,gd9a9a58781+fcb1d3bbc8,gde0f65d7ad+e6bd566e97,ge278dab8ac+932305ba37,ge82c20c137+76d20ab76d,w.2025.10
LSST Data Management Base Package
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
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
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 */
72 LSST_DECLARE_CONTROL_FIELD(clsSimpleShapeControl, SimpleShapeControl, sigma);
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
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.
Struct to hold the results of SimpleShape when we don't run it as a plugin.
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