LSST Applications g0b6bd0c080+a72a5dd7e6,g1182afd7b4+2a019aa3bb,g17e5ecfddb+2b8207f7de,g1d67935e3f+06cf436103,g38293774b4+ac198e9f13,g396055baef+6a2097e274,g3b44f30a73+6611e0205b,g480783c3b1+98f8679e14,g48ccf36440+89c08d0516,g4b93dc025c+98f8679e14,g5c4744a4d9+a302e8c7f0,g613e996a0d+e1c447f2e0,g6c8d09e9e7+25247a063c,g7271f0639c+98f8679e14,g7a9cd813b8+124095ede6,g9d27549199+a302e8c7f0,ga1cf026fa3+ac198e9f13,ga32aa97882+7403ac30ac,ga786bb30fb+7a139211af,gaa63f70f4e+9994eb9896,gabf319e997+ade567573c,gba47b54d5d+94dc90c3ea,gbec6a3398f+06cf436103,gc6308e37c7+07dd123edb,gc655b1545f+ade567573c,gcc9029db3c+ab229f5caf,gd01420fc67+06cf436103,gd877ba84e5+06cf436103,gdb4cecd868+6f279b5b48,ge2d134c3d5+cc4dbb2e3f,ge448b5faa6+86d1ceac1d,gecc7e12556+98f8679e14,gf3ee170dca+25247a063c,gf4ac96e456+ade567573c,gf9f5ea5b4d+ac198e9f13,gff490e6085+8c2580be5c,w.2022.27
LSST Data Management Base Package
multiShapeletBasis.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
24#include "ndarray/pybind11.h"
25
27
28namespace py = pybind11;
29using namespace pybind11::literals;
30
31namespace lsst {
32namespace shapelet {
33
34PYBIND11_MODULE(multiShapeletBasis, mod) {
35 py::class_<MultiShapeletBasisComponent, std::shared_ptr<MultiShapeletBasisComponent>>
36 clsMultiShapeletBasisComponent(mod, "MultiShapeletBasisComponent");
37
38 clsMultiShapeletBasisComponent.def(py::init<double, int, ndarray::Array<double const, 2, 2> const &>(),
39 "radius"_a, "order"_a, "matrix"_a);
40
41 clsMultiShapeletBasisComponent.def("getRadius", &MultiShapeletBasisComponent::getRadius);
42 clsMultiShapeletBasisComponent.def("getOrder", &MultiShapeletBasisComponent::getOrder);
43 clsMultiShapeletBasisComponent.def("getMatrix", &MultiShapeletBasisComponent::getMatrix);
44
45 py::class_<MultiShapeletBasis, std::shared_ptr<MultiShapeletBasis>> clsMultiShapeletBasis(
46 mod, "MultiShapeletBasis");
47
48 clsMultiShapeletBasis.attr("Component") = clsMultiShapeletBasisComponent;
49
50 clsMultiShapeletBasis.def(py::init<int>());
51 clsMultiShapeletBasis.def(py::init<MultiShapeletBasis const &>());
52
53 clsMultiShapeletBasis.def("getSize", &MultiShapeletBasis::getSize);
54 clsMultiShapeletBasis.def("getComponentCount", &MultiShapeletBasis::getComponentCount);
55 clsMultiShapeletBasis.def("addComponent", &MultiShapeletBasis::addComponent);
56 clsMultiShapeletBasis.def("scale", &MultiShapeletBasis::scale);
57 clsMultiShapeletBasis.def("normalize", &MultiShapeletBasis::normalize);
58 clsMultiShapeletBasis.def("merge", &MultiShapeletBasis::merge);
59 clsMultiShapeletBasis.def("makeFunction", &MultiShapeletBasis::makeFunction);
60}
61
62} // shapelet
63} // lsst
int getOrder() const
Order of this shapelet expansion.
double getRadius() const
Return the radius of this shapelet expansion.
ndarray::Array< double const, 2, 2 > getMatrix() const
Matrix whose elements [i,j] map MultiShapeletBasis elements j to shapelet terms i.
void merge(MultiShapeletBasis const &other)
Combine the given basis with this (in place), by appending its elements.
void normalize()
Rescale all matrices so each element has unit flux.
int getComponentCount() const
Return the number of components (distinct shapelet bases) in the MultiShapeletBasis.
int getSize() const
Return the number of elements (i.e. free amplitudes) in the MultiShapeletBasis.
void scale(double factor)
Multiply the radius of all basis elements by the given factor.
MultiShapeletFunction makeFunction(afw::geom::ellipses::Ellipse const &ellipse, ndarray::Array< double const, 1, 1 > const &coefficients) const
Create a MultiShapeletFunction from the basis.
void addComponent(double radius, int order, ndarray::Array< double const, 2, 2 > const &matrix)
Add a new component (shapelet basis) to the MultiShapeletBasis.
PYBIND11_MODULE(basisEvaluator, mod)
A base class for image defects.