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
copyGoodPixels.cc
Go to the documentation of this file.
1// -*- LSST-C++ -*-
2
3/*
4 * LSST Data Management System
5 * Copyright 2008, 2009, 2010 LSST Corporation.
6 *
7 * This product includes software developed by the
8 * LSST Project (http://www.lsst.org/).
9 *
10 * This program is free software: you can redistribute it and/or modify
11 * it under the terms of the GNU General Public License as published by
12 * the Free Software Foundation, either version 3 of the License, or
13 * (at your option) any later version.
14 *
15 * This program is distributed in the hope that it will be useful,
16 * but WITHOUT ANY WARRANTY; without even the implied warranty of
17 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18 * GNU General Public License for more details.
19 *
20 * You should have received a copy of the LSST License Statement and
21 * the GNU General Public License along with this program. If not,
22 * see <http://www.lsstcorp.org/LegalNotices/>.
23 */
24
30#include <cstdint>
31#include <limits>
32
33#include "boost/format.hpp"
34
35#include "lsst/pex/exceptions.h"
36#include "lsst/geom.h"
38
40namespace geom = lsst::geom;
41namespace afwImage = lsst::afw::image;
43
44namespace {
45 /*
46 * A boolean functor to test if a MaskedImage pixel is valid (mask & badPixel == 0)
47 */
48 struct CheckMask {
49 CheckMask(lsst::afw::image::MaskPixel badPixel) : _badPixel(badPixel) {}
50
51 template<typename T>
52 bool operator()(T val) const {
53 return ((val.mask() & _badPixel) == 0) ? true : false;
54 }
55 private:
57 };
58
59 /*
60 * A boolean functor to test if an Image pixel has known value (not NaN)
61 */
62 struct CheckKnownValue {
63 CheckKnownValue(lsst::afw::image::MaskPixel) {}
64
65 template<typename T>
66 bool operator()(T val) const {
67 return !std::isnan(static_cast<float>(*val));
68 }
69 };
70
71 /*
72 * Implementation of copyGoodPixels
73 *
74 * The template parameter isValidPixel is a functor with operator()
75 * which returns true if a given image pixel is valid.
76 * This allows us to support multiple image types including
77 * MaskedImage with a mask bit and Image with a check for NaN.
78 *
79 * @return overlapping bounding box, relative to parent image
80 */
81 template <typename ImageT, typename isValidPixel>
82 int copyGoodPixelsImpl(
83 ImageT &destImage,
84 ImageT const &srcImage,
85 lsst::afw::image::MaskPixel const badPixelMask
86 ) {
87 geom::Box2I overlapBBox = destImage.getBBox();
88 overlapBBox.clip(srcImage.getBBox());
89 if (overlapBBox.isEmpty()) {
90 return 0;
91 }
92
93 ImageT destView(destImage, overlapBBox, afwImage::PARENT, false);
94 ImageT srcView(srcImage, overlapBBox, afwImage::PARENT, false);
95
96 isValidPixel const isValid(badPixelMask); // functor to check if a pixel is good
97 int numGoodPix = 0;
98 for (int y = 0, endY = srcView.getHeight(); y != endY; ++y) {
99 typename ImageT::const_x_iterator srcIter = srcView.row_begin(y);
100 typename ImageT::const_x_iterator const srcEndIter = srcView.row_end(y);
101 typename ImageT::x_iterator destIter = destView.row_begin(y);
102 for (; srcIter != srcEndIter; ++srcIter, ++destIter) {
103 if (isValid(srcIter)) {
104 *destIter = *srcIter;
105// typename ImageT::SinglePixel pix = *srcIter;
106// *destIter = pix;
107 ++numGoodPix;
108 }
109 }
110 }
111 return numGoodPix;
112 }
113} // anonymous namespace
114
115template <typename ImagePixelT>
117 // spell out lsst:afw::image to make Doxygen happy
120) {
122 return copyGoodPixelsImpl<Image, CheckKnownValue>(destImage, srcImage, 0x0);
123}
124
125template <typename ImagePixelT>
127 // spell out lsst:afw::image to make Doxygen happy
131 lsst::afw::image::VariancePixel> const &srcImage,
132 lsst::afw::image::MaskPixel const badPixelMask
133) {
135 return copyGoodPixelsImpl<Image, CheckMask>(destImage, srcImage, badPixelMask);
136}
137
138// Explicit instantiations
139
141#define MASKEDIMAGE(IMAGEPIXEL) afwImage::MaskedImage<IMAGEPIXEL, \
142 afwImage::MaskPixel, afwImage::VariancePixel>
143#define INSTANTIATE(IMAGEPIXEL) \
144 template int coaddUtils::copyGoodPixels<IMAGEPIXEL>( \
145 afwImage::Image<IMAGEPIXEL> &destImage, \
146 afwImage::Image<IMAGEPIXEL> const &srcImage \
147 ); \
148 \
149 template int coaddUtils::copyGoodPixels<IMAGEPIXEL>( \
150 MASKEDIMAGE(IMAGEPIXEL) &destImage, \
151 MASKEDIMAGE(IMAGEPIXEL) const &srcImage, \
152 afwImage::MaskPixel const badPixelMask \
153 );
154
155INSTANTIATE(double);
156INSTANTIATE(float);
157INSTANTIATE(int);
#define INSTANTIATE(FROMSYS, TOSYS)
Definition Detector.cc:509
int y
Definition SpanSet.cc:48
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
An integer coordinate rectangle.
Definition Box.h:55
void clip(Box2I const &other) noexcept
Shrink this to ensure that other.contains(*this).
Definition Box.cc:189
bool isEmpty() const noexcept
Return true if the box contains no points.
Definition Box.h:213
bool isValid
Definition fits.cc:404
T isnan(T... args)
float VariancePixel
default type for MaskedImage variance images
int copyGoodPixels(lsst::afw::image::Image< ImagePixelT > &destImage, lsst::afw::image::Image< ImagePixelT > const &srcImage)
copy good pixels from one image to another
ImageT val
Definition CR.cc:146