LSST Applications g0f08755f38+82efc23009,g12f32b3c4e+e7bdf1200e,g1653933729+a8ce1bb630,g1a0ca8cf93+50eff2b06f,g28da252d5a+52db39f6a5,g2bbee38e9b+37c5a29d61,g2bc492864f+37c5a29d61,g2cdde0e794+c05ff076ad,g3156d2b45e+41e33cbcdc,g347aa1857d+37c5a29d61,g35bb328faa+a8ce1bb630,g3a166c0a6a+37c5a29d61,g3e281a1b8c+fb992f5633,g414038480c+7f03dfc1b0,g41af890bb2+11b950c980,g5fbc88fb19+17cd334064,g6b1c1869cb+12dd639c9a,g781aacb6e4+a8ce1bb630,g80478fca09+72e9651da0,g82479be7b0+04c31367b4,g858d7b2824+82efc23009,g9125e01d80+a8ce1bb630,g9726552aa6+8047e3811d,ga5288a1d22+e532dc0a0b,gae0086650b+a8ce1bb630,gb58c049af0+d64f4d3760,gc28159a63d+37c5a29d61,gcf0d15dbbd+2acd6d4d48,gd7358e8bfb+778a810b6e,gda3e153d99+82efc23009,gda6a2b7d83+2acd6d4d48,gdaeeff99f8+1711a396fd,ge2409df99d+6b12de1076,ge79ae78c31+37c5a29d61,gf0baf85859+d0a5978c5a,gf3967379c6+4954f8c433,gfb92a5be7c+82efc23009,gfec2e1e490+2aaed99252,w.2024.46
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