Loading [MathJax]/extensions/tex2jax.js
LSST Applications g04a91732dc+cc8870d3f5,g07dc498a13+5aa0b8792f,g0fba68d861+488cddfaa9,g1409bbee79+5aa0b8792f,g1a7e361dbc+5aa0b8792f,g1fd858c14a+f64bc332a9,g35bb328faa+fcb1d3bbc8,g4d2262a081+b1c1982739,g4d39ba7253+9633a327c1,g4e0f332c67+5d362be553,g53246c7159+fcb1d3bbc8,g60b5630c4e+9633a327c1,g668ecb457e+25d63fd678,g78460c75b0+2f9a1b4bcd,g786e29fd12+cf7ec2a62a,g7b71ed6315+fcb1d3bbc8,g8852436030+8b64ca622a,g89139ef638+5aa0b8792f,g89e1512fd8+37f975783e,g8d6b6b353c+9633a327c1,g9125e01d80+fcb1d3bbc8,g989de1cb63+5aa0b8792f,g9f33ca652e+b196626af7,ga9baa6287d+9633a327c1,gaaedd4e678+5aa0b8792f,gabe3b4be73+1e0a283bba,gb1101e3267+71e32094df,gb58c049af0+f03b321e39,gb90eeb9370+2807b1ad02,gc741bbaa4f+1ae86710ed,gcf25f946ba+8b64ca622a,gd315a588df+a39986a76f,gd6cbbdb0b4+c8606af20c,gd9a9a58781+fcb1d3bbc8,gde0f65d7ad+94dfc458f4,ge278dab8ac+932305ba37,ge82c20c137+76d20ab76d,gfe73954cf8+a1301e4c20,w.2025.11
LSST Data Management Base Package
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
_convolve.cc
Go to the documentation of this file.
1/*
2 * LSST Data Management System
3 * Copyright 2008-2016 AURA/LSST.
4 *
5 * This product includes software developed by the
6 * LSST Project (http://www.lsst.org/).
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
23#include <pybind11/pybind11.h>
25
27
28namespace py = pybind11;
29
30using namespace py::literals;
31
32namespace lsst {
33namespace afw {
34namespace math {
35namespace detail {
36
37namespace {
38template <typename OutImageT, typename InImageT>
39void declareByType(lsst::cpputils::python::WrapperCollection &wrappers) {
40 wrappers.wrap([](auto &mod) {
41 mod.def("basicConvolve",
42 (void (*)(OutImageT &, InImageT const &, lsst::afw::math::Kernel const &,
43 lsst::afw::math::ConvolutionControl const &))basicConvolve<OutImageT, InImageT>);
44 mod.def("basicConvolve",
45 (void (*)(OutImageT &, InImageT const &, lsst::afw::math::DeltaFunctionKernel const &,
46 lsst::afw::math::ConvolutionControl const &))basicConvolve<OutImageT, InImageT>);
47 mod.def("basicConvolve",
48 (void (*)(OutImageT &, InImageT const &, lsst::afw::math::LinearCombinationKernel const &,
49 lsst::afw::math::ConvolutionControl const &))basicConvolve<OutImageT, InImageT>);
50 mod.def("basicConvolve",
51 (void (*)(OutImageT &, InImageT const &, lsst::afw::math::SeparableKernel const &,
52 lsst::afw::math::ConvolutionControl const &))basicConvolve<OutImageT, InImageT>);
53 mod.def("convolveWithBruteForce",
54 (void (*)(OutImageT &, InImageT const &, lsst::afw::math::Kernel const &,
55 lsst::afw::math::ConvolutionControl const &))
57 });
58}
59template <typename PixelType1, typename PixelType2>
60void declareAll(lsst::cpputils::python::WrapperCollection &wrappers) {
61 using M1 = image::MaskedImage<PixelType1, image::MaskPixel, image::VariancePixel>;
62 using M2 = image::MaskedImage<PixelType2, image::MaskPixel, image::VariancePixel>;
63
64 declareByType<image::Image<PixelType1>, image::Image<PixelType2>>(wrappers);
65 declareByType<M1, M2>(wrappers);
66}
67} // namespace
68
70 using PyClass = py::class_<KernelImagesForRegion, std::shared_ptr<KernelImagesForRegion>>;
71 auto clsKernelImagesForRegion =
72 wrappers.wrapType(PyClass(wrappers.module, "KernelImagesForRegion"), [](auto &mod, auto &cls) {
73 cls.def(py::init<KernelImagesForRegion::KernelConstPtr, lsst::geom::Box2I const &,
74 lsst::geom::Point2I const &, bool>(),
75 "kernelPtr"_a, "bbox"_a, "xy0"_a, "doNormalize"_a);
76 cls.def(py::init<KernelImagesForRegion::KernelConstPtr, lsst::geom::Box2I const &,
77 lsst::geom::Point2I const &, bool, KernelImagesForRegion::ImagePtr,
78 KernelImagesForRegion::ImagePtr, KernelImagesForRegion::ImagePtr,
79 KernelImagesForRegion::ImagePtr>(),
80 "kernelPtr"_a, "bbox"_a, "xy0"_a, "doNormalize"_a, "bottomLeftImagePtr"_a,
81 "bottomRightImagePtr"_a, "topLeftImagePtr"_a, "topRightImagePtr"_a);
82
83 cls.def("getBBox", &KernelImagesForRegion::getBBox);
84 cls.def("getXY0", &KernelImagesForRegion::getXY0);
85 cls.def("getDoNormalize", &KernelImagesForRegion::getDoNormalize);
86 cls.def("getImage", &KernelImagesForRegion::getImage);
87 cls.def("getKernel", &KernelImagesForRegion::getKernel);
88 cls.def("getPixelIndex", &KernelImagesForRegion::getPixelIndex);
89 cls.def("computeNextRow", &KernelImagesForRegion::computeNextRow);
90 cls.def_static("getMinInterpolationSize", KernelImagesForRegion::getMinInterpolationSize);
91 });
92
93 wrappers.wrapType(py::enum_<KernelImagesForRegion::Location>(clsKernelImagesForRegion, "Location"),
94 [](auto &mod, auto &enm) {
95 enm.value("BOTTOM_LEFT", KernelImagesForRegion::Location::BOTTOM_LEFT);
96 enm.value("BOTTOM_RIGHT", KernelImagesForRegion::Location::BOTTOM_RIGHT);
97 enm.value("TOP_LEFT", KernelImagesForRegion::Location::TOP_LEFT);
98 enm.value("TOP_RIGHT", KernelImagesForRegion::Location::TOP_RIGHT);
99 enm.export_values();
100 });
101
103 wrappers.module, "RowOfKernelImagesForRegion"),
104 [](auto &mod, auto &cls) {
105 cls.def(py::init<int, int>(), "nx"_a, "ny"_a);
106
107 cls.def("front", &RowOfKernelImagesForRegion::front);
108 cls.def("back", &RowOfKernelImagesForRegion::back);
109 cls.def("getNX", &RowOfKernelImagesForRegion::getNX);
110 cls.def("getNY", &RowOfKernelImagesForRegion::getNY);
111 cls.def("getYInd", &RowOfKernelImagesForRegion::getYInd);
112 cls.def("getRegion", &RowOfKernelImagesForRegion::getRegion);
113 cls.def("hasData", &RowOfKernelImagesForRegion::hasData);
114 cls.def("isLastRow", &RowOfKernelImagesForRegion::isLastRow);
115 cls.def("incrYInd", &RowOfKernelImagesForRegion::incrYInd);
116 });
117}
119 wrappers.addSignatureDependency("lsst.afw.image");
120 wrappers.addSignatureDependency("lsst.afw.math");
121 declareConvolve(wrappers);
122 declareAll<double, double>(wrappers);
123 declareAll<double, float>(wrappers);
124 declareAll<double, int>(wrappers);
125 declareAll<double, std::uint16_t>(wrappers);
126 declareAll<float, float>(wrappers);
127 declareAll<float, int>(wrappers);
128 declareAll<float, std::uint16_t>(wrappers);
129 declareAll<int, int>(wrappers);
130 declareAll<std::uint16_t, std::uint16_t>(wrappers);
131}
132} // namespace detail
133} // namespace math
134} // namespace afw
135} // namespace lsst
A helper class for subdividing pybind11 module across multiple translation units (i....
Definition python.h:242
void addSignatureDependency(std::string const &name)
Indicate an external module that provides a type used in function/method signatures.
Definition python.h:357
void wrap(WrapperCallback function)
Add a set of wrappers without defining a class.
Definition python.h:369
PyType wrapType(PyType cls, ClassWrapperCallback function, bool setModuleName=true)
Add a type (class or enum) wrapper, deferring method and other attribute definitions until finish() i...
Definition python.h:391
pybind11::module module
The module object passed to the PYBIND11_MODULE block that contains this WrapperCollection.
Definition python.h:448
void declareConvolve(lsst::cpputils::python::WrapperCollection &wrappers)
Definition _convolve.cc:69
void wrapConvolve(lsst::cpputils::python::WrapperCollection &wrappers)
Definition _convolve.cc:118
void basicConvolve(OutImageT &convolvedImage, InImageT const &inImage, lsst::afw::math::Kernel const &kernel, lsst::afw::math::ConvolutionControl const &convolutionControl)
Low-level convolution function that does not set edge pixels.
void convolveWithBruteForce(OutImageT &convolvedImage, InImageT const &inImage, lsst::afw::math::Kernel const &kernel, lsst::afw::math::ConvolutionControl const &convolutionControl)
Convolve an Image or MaskedImage with a Kernel by computing the kernel image at every point.
py::class_< PixelAreaBoundedField, std::shared_ptr< PixelAreaBoundedField >, BoundedField > PyClass