LSST Applications g063fba187b+fee0456c91,g0f08755f38+ea96e5a5a3,g1653933729+a8ce1bb630,g168dd56ebc+a8ce1bb630,g1a2382251a+90257ff92a,g20f6ffc8e0+ea96e5a5a3,g217e2c1bcf+937a289c59,g28da252d5a+daa7da44eb,g2bbee38e9b+253935c60e,g2bc492864f+253935c60e,g3156d2b45e+6e55a43351,g32e5bea42b+31359a2a7a,g347aa1857d+253935c60e,g35bb328faa+a8ce1bb630,g3a166c0a6a+253935c60e,g3b1af351f3+a8ce1bb630,g3e281a1b8c+c5dd892a6c,g414038480c+416496e02f,g41af890bb2+afe91b1188,g599934f4f4+0db33f7991,g7af13505b9+e36de7bce6,g80478fca09+da231ba887,g82479be7b0+a4516e59e3,g858d7b2824+ea96e5a5a3,g89c8672015+f4add4ffd5,g9125e01d80+a8ce1bb630,ga5288a1d22+bc6ab8dfbd,gb58c049af0+d64f4d3760,gc28159a63d+253935c60e,gcab2d0539d+3f2b72788c,gcf0d15dbbd+4ea9c45075,gda6a2b7d83+4ea9c45075,gdaeeff99f8+1711a396fd,ge79ae78c31+253935c60e,gef2f8181fd+3031e3cf99,gf0baf85859+c1f95f4921,gfa517265be+ea96e5a5a3,gfa999e8aa5+17cd334064,w.2024.50
LSST Data Management Base Package
Loading...
Searching...
No Matches
binImage.cc
Go to the documentation of this file.
1/*
2 * LSST Data Management System
3 * Copyright 2008, 2009, 2010 LSST Corporation.
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 <http://www.lsstcorp.org/LegalNotices/>.
21 */
22
23/*
24 * Bin an Image or MaskedImage by an integral factor (the same in x and y)
25 */
26#include <memory>
27#include <cstdint>
28
29#include "lsst/pex/exceptions.h"
31
33
34namespace lsst {
35namespace afw {
36namespace math {
37
38template <typename ImageT>
39std::shared_ptr<ImageT> binImage(ImageT const& in, int const binsize, lsst::afw::math::Property const flags) {
40 return binImage(in, binsize, binsize, flags);
41}
42
43template <typename ImageT>
44std::shared_ptr<ImageT> binImage(ImageT const& in, int const binX, int const binY,
45 lsst::afw::math::Property const flags) {
46 if (flags != lsst::afw::math::MEAN) {
48 (boost::format("Only afwMath::MEAN is supported, saw 0x%x") % flags).str());
49 }
50 if (binX <= 0 || binY <= 0) {
52 (boost::format("Binning must be >= 0, saw %dx%d") % binX % binY).str());
53 }
54
55 int const outWidth = in.getWidth() / binX;
56 int const outHeight = in.getHeight() / binY;
57
59 std::shared_ptr<ImageT>(new ImageT(lsst::geom::Extent2I(outWidth, outHeight)));
60 out->setXY0(in.getXY0());
61 *out = typename ImageT::SinglePixel(0);
62
63 for (int oy = 0, iy = 0; oy < out->getHeight(); ++oy) {
64 for (int i = 0; i != binY; ++i, ++iy) {
65 typename ImageT::x_iterator optr = out->row_begin(oy);
66 for (typename ImageT::x_iterator iptr = in.row_begin(iy), iend = iptr + binX * outWidth;
67 iptr < iend;) {
68 typename ImageT::SinglePixel val = *iptr;
69 ++iptr;
70 for (int j = 1; j != binX; ++j, ++iptr) {
71 val += *iptr;
72 }
73 *optr += val;
74 ++optr;
75 }
76 }
77 for (typename ImageT::x_iterator ptr = out->row_begin(oy), end = out->row_end(oy); ptr != end;
78 ++ptr) {
79 *ptr /= binX * binY;
80 }
81 }
82
83 return out;
84}
85
86//
87// Explicit instantiations
88//
90#define INSTANTIATE(TYPE) \
91 template std::shared_ptr<image::Image<TYPE>> binImage(image::Image<TYPE> const&, int, \
92 lsst::afw::math::Property const); \
93 template std::shared_ptr<image::Image<TYPE>> binImage(image::Image<TYPE> const&, int, int, \
94 lsst::afw::math::Property const); \
95 template std::shared_ptr<image::MaskedImage<TYPE>> binImage(image::MaskedImage<TYPE> const&, int, \
96 lsst::afw::math::Property const); \
97 template std::shared_ptr<image::MaskedImage<TYPE>> binImage(image::MaskedImage<TYPE> const&, int, int, \
98 lsst::afw::math::Property const);
99
101INSTANTIATE(int)
102INSTANTIATE(float)
103INSTANTIATE(double)
105} // namespace math
106} // namespace afw
107} // namespace lsst
int end
#define INSTANTIATE(FROMSYS, TOSYS)
Definition Detector.cc:509
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition Exception.h:48
std::uint64_t * ptr
Definition RangeSet.cc:95
Reports arguments outside the domain of an operation.
Definition Runtime.h:57
Reports invalid arguments.
Definition Runtime.h:66
Property
control what is calculated
Definition Statistics.h:53
@ MEAN
estimate sample mean
Definition Statistics.h:57
std::shared_ptr< ImageT > binImage(ImageT const &inImage, int const binX, int const binY, lsst::afw::math::Property const flags=lsst::afw::math::MEAN)
Definition binImage.cc:44
ImageT val
Definition CR.cc:146