Loading [MathJax]/extensions/tex2jax.js
LSST Applications g0485b4d2cb+c8d56b10d4,g0fba68d861+c513f10703,g1ec0fe41b4+3e153770da,g1fd858c14a+fca1f36da6,g2440f9efcc+8c5ae1fdc5,g35bb328faa+8c5ae1fdc5,g4d2262a081+73db4aee98,g53246c7159+8c5ae1fdc5,g56a49b3a55+8699aedcf1,g60b5630c4e+0e13f5d8e4,g65f11fe441+d5838b1573,g67b6fd64d1+035c836e50,g6ba5f4ee81+5218eae6c4,g78460c75b0+7e33a9eb6d,g786e29fd12+668abc6043,g8352419a5c+8c5ae1fdc5,g8852436030+5e76f8a839,g89139ef638+035c836e50,g94187f82dc+0e13f5d8e4,g989de1cb63+035c836e50,g9d31334357+0e13f5d8e4,g9f33ca652e+190f04d267,gabe3b4be73+8856018cbb,gabf8522325+21619da9f3,gb1101e3267+e14fd59f59,gb89ab40317+035c836e50,gc91f06edcd+a47578e525,gcef618a4dd+0e13f5d8e4,gcf25f946ba+5e76f8a839,gd6cbbdb0b4+958adf5c1f,gdb8242c116+18fb55f1cc,gde0f65d7ad+fe97e93be4,ge278dab8ac+83c63f4893,ge410e46f29+035c836e50,gf35d7ec915+97dd712d81,gf5e32f922b+8c5ae1fdc5,gf67bdafdda+035c836e50,gf6800124b1+a9633b82fa,w.2025.19
LSST Data Management Base Package
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
spatialModelPsf.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"
24#include "pybind11/stl.h"
25
27
28namespace py = pybind11;
29using namespace pybind11::literals;
30
31namespace lsst {
32namespace meas {
33namespace algorithms {
34namespace {
35
36template <typename PixelT>
37void declareFunctions(lsst::cpputils::python::WrapperCollection &wrappers) {
38 using MaskedImageT = afw::image::MaskedImage<PixelT, afw::image::MaskPixel, afw::image::VariancePixel>;
39 auto &mod = wrappers.module;
40 mod.def("createKernelFromPsfCandidates", createKernelFromPsfCandidates<PixelT>, "psfCells"_a, "dims"_a,
41 "xy0"_a, "nEigenComponents"_a, "spatialOrder"_a, "ksize"_a, "nStarPerCell"_a = -1,
42 "constantWeight"_a = true, "border"_a = 3);
43 mod.def("countPsfCandidates", countPsfCandidates<PixelT>, "psfCells"_a, "nStarPerCell"_a = -1);
44 mod.def("fitSpatialKernelFromPsfCandidates",
45 (std::pair<bool, double>(*)(afw::math::Kernel *, afw::math::SpatialCellSet const &, int const,
46 double const, double const))fitSpatialKernelFromPsfCandidates<PixelT>,
47 "kernel"_a, "psfCells"_a, "nStarPerCell"_a = -1, "tolerance"_a = 1e-5, "lambda"_a = 0.0);
48 mod.def("fitSpatialKernelFromPsfCandidates",
49 (std::pair<bool, double>(*)(afw::math::Kernel *, afw::math::SpatialCellSet const &, bool const,
50 int const, double const,
52 "kernel"_a, "psfCells"_a, "doNonLinearFit"_a, "nStarPerCell"_a = -1, "tolerance"_a = 1e-5,
53 "lambda"_a = 0.0);
54 mod.def("subtractPsf", subtractPsf<MaskedImageT>, "psf"_a, "data"_a, "x"_a, "y"_a,
56 mod.def("fitKernelParamsToImage", fitKernelParamsToImage<MaskedImageT>, "kernel"_a, "image"_a, "pos"_a);
57 mod.def("fitKernelToImage", fitKernelToImage<MaskedImageT>, "kernel"_a, "image"_a, "pos"_a);
58}
59} // namespace
60
62 declareFunctions<float>(wrappers);
63}
64
65} // namespace algorithms
66} // namespace meas
67} // namespace lsst
Class used by SpatialCell for spatial PSF fittig.
A helper class for subdividing pybind11 module across multiple translation units (i....
Definition python.h:242
pybind11::module module
The module object passed to the PYBIND11_MODULE block that contains this WrapperCollection.
Definition python.h:448
void wrapSpatialModelPsf(WrapperCollection &wrappers)
std::pair< std::shared_ptr< afw::math::LinearCombinationKernel >, std::vector< double > > createKernelFromPsfCandidates(afw::math::SpatialCellSet const &psfCells, geom::Extent2I const &dims, geom::Point2I const &xy0, int const nEigenComponents, int const spatialOrder, int const ksize, int const nStarPerCell=-1, bool const constantWeight=true, int const border=3)
Return a Kernel pointer and a list of eigenvalues resulting from analysing the provided SpatialCellSe...
int countPsfCandidates(afw::math::SpatialCellSet const &psfCells, int const nStarPerCell=-1)
Count the number of candidates in use.
std::pair< std::vector< double >, afw::math::KernelList > fitKernelParamsToImage(afw::math::LinearCombinationKernel const &kernel, Image const &image, geom::Point2D const &pos)
Fit a LinearCombinationKernel to an Image, allowing the coefficients of the components to vary.
std::pair< std::shared_ptr< afw::math::Kernel >, std::pair< double, double > > fitKernelToImage(afw::math::LinearCombinationKernel const &kernel, Image const &image, geom::Point2D const &pos)
Fit a LinearCombinationKernel to an Image, allowing the coefficients of the components to vary.
std::pair< bool, double > fitSpatialKernelFromPsfCandidates(afw::math::Kernel *kernel, afw::math::SpatialCellSet const &psfCells, int const nStarPerCell=-1, double const tolerance=1e-5, double const lambda=0.0)
Fit spatial kernel using full-nonlinear optimization to estimate candidate amplitudes.
double subtractPsf(afw::detection::Psf const &psf, ImageT *data, double x, double y, double psfFlux=std::numeric_limits< double >::quiet_NaN())
T quiet_NaN(T... args)