LSST Applications 27.0.0,g0265f82a02+469cd937ee,g02d81e74bb+21ad69e7e1,g1470d8bcf6+cbe83ee85a,g2079a07aa2+e67c6346a6,g212a7c68fe+04a9158687,g2305ad1205+94392ce272,g295015adf3+81dd352a9d,g2bbee38e9b+469cd937ee,g337abbeb29+469cd937ee,g3939d97d7f+72a9f7b576,g487adcacf7+71499e7cba,g50ff169b8f+5929b3527e,g52b1c1532d+a6fc98d2e7,g591dd9f2cf+df404f777f,g5a732f18d5+be83d3ecdb,g64a986408d+21ad69e7e1,g858d7b2824+21ad69e7e1,g8a8a8dda67+a6fc98d2e7,g99cad8db69+f62e5b0af5,g9ddcbc5298+d4bad12328,ga1e77700b3+9c366c4306,ga8c6da7877+71e4819109,gb0e22166c9+25ba2f69a1,gb6a65358fc+469cd937ee,gbb8dafda3b+69d3c0e320,gc07e1c2157+a98bf949bb,gc120e1dc64+615ec43309,gc28159a63d+469cd937ee,gcf0d15dbbd+72a9f7b576,gdaeeff99f8+a38ce5ea23,ge6526c86ff+3a7c1ac5f1,ge79ae78c31+469cd937ee,gee10cc3b42+a6fc98d2e7,gf1cff7945b+21ad69e7e1,gfbcc870c63+9a11dc8c8f
LSST Data Management Base Package
Loading...
Searching...
No Matches
_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>
24#include <lsst/utils/python.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::utils::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 &,
56 convolveWithBruteForce<OutImageT, InImageT>);
57 });
58}
59template <typename PixelType1, typename PixelType2>
60void declareAll(lsst::utils::python::WrapperCollection &wrappers) {
63
64 declareByType<image::Image<PixelType1>, image::Image<PixelType2>>(wrappers);
65 declareByType<M1, M2>(wrappers);
66}
67} // namespace
68
69void declareConvolve(lsst::utils::python::WrapperCollection &wrappers) {
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}
118void wrapConvolve(lsst::utils::python::WrapperCollection &wrappers) {
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 class to represent a 2-dimensional array of pixels.
Definition Image.h:51
A class to manipulate images, masks, and variance as a single object.
Definition MaskedImage.h:74
Parameters to control convolution.
A kernel that has only one non-zero pixel (of value 1)
Definition Kernel.h:643
Kernels are used for convolution with MaskedImages and (eventually) Images.
Definition Kernel.h:110
A kernel that is a linear combination of fixed basis kernels.
Definition Kernel.h:704
A kernel described by a pair of functions: func(x, y) = colFunc(x) * rowFunc(y)
Definition Kernel.h:860
void wrapConvolve(lsst::utils::python::WrapperCollection &wrappers)
Definition _convolve.cc:118
void declareConvolve(lsst::utils::python::WrapperCollection &wrappers)
Definition _convolve.cc:69
py::class_< PixelAreaBoundedField, std::shared_ptr< PixelAreaBoundedField >, BoundedField > PyClass