LSST Applications g1cfbe01172+01aa18f939,g20cdd03214+31e6b93548,g28da252d5a+ea8665a95b,g2bbee38e9b+9ec6cc348d,g2bc492864f+9ec6cc348d,g347aa1857d+9ec6cc348d,g3a166c0a6a+9ec6cc348d,g4322eb9e3a+65eff1e020,g461a3dce89+b86e4b8053,g50ff169b8f+f991eae79d,g52b1c1532d+b86e4b8053,g607f77f49a+31e6b93548,g78056777b3+8ae2798781,g858d7b2824+31e6b93548,g8cd86fa7b1+4851e61ca4,g9ddcbc5298+f24b38b85a,ga1e77700b3+3309dba821,gae0086650b+b86e4b8053,gb0e22166c9+6076c0b52b,gbb886bcc26+dccb771098,gbd462c55f0+dc07f8e65d,gc0c51c7ec2+31e6b93548,gc120e1dc64+a417ce3171,gc28159a63d+9ec6cc348d,gc2a6998b7e+f95f64aeae,gcdd4ae20e8+507450c4cd,gcf0d15dbbd+507450c4cd,gd1535ee943+bcf88ba65f,gd598c5cd71+66126f91fb,gdaeeff99f8+006e14e809,gdbce86181e+39d5515b1a,ge3d4d395c2+b12d4d6a95,ge79ae78c31+9ec6cc348d,gf048a9a2f4+d9c36e6b63,gfbcc870c63+ea41c4420b,w.2024.27
LSST Data Management Base Package
Loading...
Searching...
No Matches
Interp.h
Go to the documentation of this file.
1// -*- LSST-C++ -*-
2
3/*
4 * LSST Data Management System
5 * Copyright 2008-2015 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
25#if !defined(LSST_MEAS_ALGORITHMS_INTERP_H)
26#define LSST_MEAS_ALGORITHMS_INTERP_H
28// Interpolate over defects in a MaskedImage
29//
30#include <limits>
31#include <vector>
32
33#include "lsst/geom/Box.h"
36
37namespace lsst {
38namespace afw {
39namespace detection {
40
41class Psf;
42}
43} // namespace afw
44
45namespace meas {
46namespace algorithms {
47namespace interp {
51double const lpc_1_c1 = 0.7737;
52double const lpc_1_c2 = -0.2737;
57double const lpc_1s2_c1 = 0.7358;
58double const lpc_1s2_c2 = -0.2358;
59/*
60 * Used to debias min(x, y)
61 */
62double const min2GaussianBias = -0.5641895835;
63
64template <typename MaskedImageT>
66 bool horizontal, double minval);
67} // namespace interp
68
73public:
75
88
89 enum { WIDE_DEFECT = 11 };
90
91 explicit Defect(const geom::BoxI &bbox = geom::BoxI()
92 )
93 : afw::image::DefectBase(bbox), _pos(static_cast<DefectPosition>(0)), _type(0) {}
94 virtual ~Defect() {}
95
97 unsigned int type
98 ) {
99 _pos = pos;
100 _type = type;
101 }
102
103 unsigned int getType() const { return _type; }
104 DefectPosition getPos() const { return _pos; }
105private:
106 DefectPosition _pos;
107 unsigned int _type;
108};
109
110template <typename MaskedImageT>
111void interpolateOverDefects(MaskedImageT &image, afw::detection::Psf const &psf,
112 std::vector<Defect::Ptr> &badList, double fallbackValue = 0.0,
113 bool useFallbackValueAtEdge = false);
114
115} // namespace algorithms
116} // namespace meas
117} // namespace lsst
118
119#endif
AmpInfoBoxKey bbox
Definition Amplifier.cc:117
table::Key< int > type
Definition Detector.cc:163
int y
Definition SpanSet.cc:48
A polymorphic base class for representing an image's Point Spread Function.
Definition Psf.h:82
Encapsulate information about a bad portion of a detector.
Definition Defect.h:39
DefectBase(const lsst::geom::Box2I &bbox)
Definition Defect.h:41
An integer coordinate rectangle.
Definition Box.h:55
Encapsulate information about a bad portion of a detector.
Definition Interp.h:72
unsigned int getType() const
Return the defect's interpolation type.
Definition Interp.h:103
void classify(DefectPosition pos, unsigned int type)
Definition Interp.h:96
std::shared_ptr< Defect > Ptr
shared pointer to Defect
Definition Interp.h:74
@ WIDE
defect is in middle, and wide
Definition Interp.h:82
@ NEAR_RIGHT
defect is near right boundary
Definition Interp.h:84
@ NEAR_LEFT
defect is near left boundary
Definition Interp.h:78
@ WIDE_LEFT
defect is wide at left boundary
Definition Interp.h:79
@ WIDE_RIGHT
defect is wide at right boundary
Definition Interp.h:85
@ MIDDLE
defect is in middle of frame
Definition Interp.h:80
@ WIDE_NEAR_LEFT
defect is near left, and wide
Definition Interp.h:81
@ RIGHT
defect is at right boundary
Definition Interp.h:86
@ LEFT
defect is at left boundary
Definition Interp.h:77
@ WIDE_NEAR_RIGHT
defect is near right, and wide
Definition Interp.h:83
Defect(const geom::BoxI &bbox=geom::BoxI())
Definition Interp.h:91
DefectPosition getPos() const
Return the position of the defect.
Definition Interp.h:104
double const min2GaussianBias
Mean value of the minimum of two N(0,1) variates.
Definition Interp.h:62
double const lpc_1_c1
LPC coefficients for sigma = 1, S/N = infty.
Definition Interp.h:51
double const lpc_1s2_c1
LPC coefficients for sigma = 1/sqrt(2), S/N = infty.
Definition Interp.h:57
std::pair< bool, typename MaskedImageT::Image::Pixel > singlePixel(int x, int y, MaskedImageT const &image, bool horizontal, double minval)
Return a boolean status (true: interpolation is OK) and the interpolated value for a pixel,...
Definition Interp.cc:2125
void interpolateOverDefects(MaskedImageT &image, afw::detection::Psf const &psf, std::vector< Defect::Ptr > &badList, double fallbackValue=0.0, bool useFallbackValueAtEdge=false)
Process a set of known bad pixels in an image.
Definition Interp.cc:2047