LSST Applications  21.0.0+04719a4bac,21.0.0-1-ga51b5d4+f5e6047307,21.0.0-11-g2b59f77+a9c1acf22d,21.0.0-11-ga42c5b2+86977b0b17,21.0.0-12-gf4ce030+76814010d2,21.0.0-13-g1721dae+760e7a6536,21.0.0-13-g3a573fe+768d78a30a,21.0.0-15-g5a7caf0+f21cbc5713,21.0.0-16-g0fb55c1+b60e2d390c,21.0.0-19-g4cded4ca+71a93a33c0,21.0.0-2-g103fe59+bb20972958,21.0.0-2-g45278ab+04719a4bac,21.0.0-2-g5242d73+3ad5d60fb1,21.0.0-2-g7f82c8f+8babb168e8,21.0.0-2-g8f08a60+06509c8b61,21.0.0-2-g8faa9b5+616205b9df,21.0.0-2-ga326454+8babb168e8,21.0.0-2-gde069b7+5e4aea9c2f,21.0.0-2-gecfae73+1d3a86e577,21.0.0-2-gfc62afb+3ad5d60fb1,21.0.0-25-g1d57be3cd+e73869a214,21.0.0-3-g357aad2+ed88757d29,21.0.0-3-g4a4ce7f+3ad5d60fb1,21.0.0-3-g4be5c26+3ad5d60fb1,21.0.0-3-g65f322c+e0b24896a3,21.0.0-3-g7d9da8d+616205b9df,21.0.0-3-ge02ed75+a9c1acf22d,21.0.0-4-g591bb35+a9c1acf22d,21.0.0-4-g65b4814+b60e2d390c,21.0.0-4-gccdca77+0de219a2bc,21.0.0-4-ge8a399c+6c55c39e83,21.0.0-5-gd00fb1e+05fce91b99,21.0.0-6-gc675373+3ad5d60fb1,21.0.0-64-g1122c245+4fb2b8f86e,21.0.0-7-g04766d7+cd19d05db2,21.0.0-7-gdf92d54+04719a4bac,21.0.0-8-g5674e7b+d1bd76f71f,master-gac4afde19b+a9c1acf22d,w.2021.13
LSST Data Management Base Package
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"
34 #include "lsst/afw/image/Defect.h"
36 
37 namespace lsst {
38 namespace afw {
39 namespace detection {
40 
41 class Psf;
42 }
43 } // namespace afw
44 
45 namespace meas {
46 namespace algorithms {
47 namespace interp {
51 double const lpc_1_c1 = 0.7737;
52 double const lpc_1_c2 = -0.2737;
57 double const lpc_1s2_c1 = 0.7358;
58 double const lpc_1s2_c2 = -0.2358;
59 /*
60  * Used to debias min(x, y)
61  */
62 double const min2GaussianBias = -0.5641895835;
63 
64 template <typename MaskedImageT>
66  bool horizontal, double minval);
67 } // namespace interp
68 
73 public:
75 
77  LEFT = 1,
82  WIDE,
86  RIGHT
87  };
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; }
105 private:
106  DefectPosition _pos;
107  unsigned int _type;
108 };
109 
110 template <typename MaskedImageT>
111 void 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
double x
table::Key< int > type
Definition: Detector.cc:163
int y
Definition: SpanSet.cc:49
A polymorphic base class for representing an image's Point Spread Function.
Definition: Psf.h:76
Encapsulate information about a bad portion of a detector.
Definition: Defect.h:41
DefectBase(const lsst::geom::Box2I &bbox)
Definition: Defect.h:43
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
Backwards-compatibility support for depersisting the old Calib (FluxMag0/FluxMag0Err) objects.
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
A base class for image defects.
Key< int > psf
Definition: Exposure.cc:65