LSSTApplications  19.0.0-14-gb0260a2+72efe9b372,20.0.0+7927753e06,20.0.0+8829bf0056,20.0.0+995114c5d2,20.0.0+b6f4b2abd1,20.0.0+bddc4f4cbe,20.0.0-1-g253301a+8829bf0056,20.0.0-1-g2b7511a+0d71a2d77f,20.0.0-1-g5b95a8c+7461dd0434,20.0.0-12-g321c96ea+23efe4bbff,20.0.0-16-gfab17e72e+fdf35455f6,20.0.0-2-g0070d88+ba3ffc8f0b,20.0.0-2-g4dae9ad+ee58a624b3,20.0.0-2-g61b8584+5d3db074ba,20.0.0-2-gb780d76+d529cf1a41,20.0.0-2-ged6426c+226a441f5f,20.0.0-2-gf072044+8829bf0056,20.0.0-2-gf1f7952+ee58a624b3,20.0.0-20-geae50cf+e37fec0aee,20.0.0-25-g3dcad98+544a109665,20.0.0-25-g5eafb0f+ee58a624b3,20.0.0-27-g64178ef+f1f297b00a,20.0.0-3-g4cc78c6+e0676b0dc8,20.0.0-3-g8f21e14+4fd2c12c9a,20.0.0-3-gbd60e8c+187b78b4b8,20.0.0-3-gbecbe05+48431fa087,20.0.0-38-ge4adf513+a12e1f8e37,20.0.0-4-g97dc21a+544a109665,20.0.0-4-gb4befbc+087873070b,20.0.0-4-gf910f65+5d3db074ba,20.0.0-5-gdfe0fee+199202a608,20.0.0-5-gfbfe500+d529cf1a41,20.0.0-6-g64f541c+d529cf1a41,20.0.0-6-g9a5b7a1+a1cd37312e,20.0.0-68-ga3f3dda+5fca18c6a4,20.0.0-9-g4aef684+e18322736b,w.2020.45
LSSTDataManagementBasePackage
Blendedness.h
Go to the documentation of this file.
1 // -*- LSST-C++ -*-
2 /*
3  * This file is part of meas_base.
4  *
5  * Developed for the LSST Data Management System.
6  * This product includes software developed by the LSST Project
7  * (https://www.lsst.org).
8  * See the COPYRIGHT file at the top-level directory of this distribution
9  * for details of code ownership.
10  *
11  * This program is free software: you can redistribute it and/or modify
12  * it under the terms of the GNU General Public License as published by
13  * the Free Software Foundation, either version 3 of the License, or
14  * (at your option) any later version.
15  *
16  * This program is distributed in the hope that it will be useful,
17  * but WITHOUT ANY WARRANTY; without even the implied warranty of
18  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
19  * GNU General Public License for more details.
20  *
21  * You should have received a copy of the GNU General Public License
22  * along with this program. If not, see <https://www.gnu.org/licenses/>.
23  */
24 
25 #ifndef LSST_MEAS_BASE_Blendedness_h_INCLUDED
26 #define LSST_MEAS_BASE_Blendedness_h_INCLUDED
27 
28 #include "lsst/pex/config.h"
31 #include "lsst/afw/table/Source.h"
32 #include "lsst/afw/image/Image.h"
35 
36 namespace lsst {
37 namespace meas {
38 namespace base {
39 
41 public:
43  doOld, bool,
44  "Whether to compute HeavyFootprint dot products (the old deblend.blendedness parameter)");
45 
46  LSST_CONTROL_FIELD(doFlux, bool, "Whether to compute quantities related to the Gaussian-weighted flux");
47 
48  LSST_CONTROL_FIELD(doShape, bool, "Whether to compute quantities related to the Gaussian-weighted shape");
49 
51  "Radius factor that sets the maximum extent of the weight function (and hence the "
52  "flux measurements)");
53 
54  BlendednessControl() : doOld(true), doFlux(true), doShape(true), nSigmaWeightMax(3.0) {}
55 };
56 
67 public:
68  // Structures and routines to manage flaghandler
69  static FlagDefinitionList const& getFlagDefinitions();
70  static FlagDefinition const FAILURE;
72  static FlagDefinition const NO_SHAPE;
73 
75 
77 
94  static float computeAbsExpectation(float data, float variance);
95 
112  static float computeAbsBias(float mu, float variance);
113 
115  afw::table::SourceRecord& child) const;
116 
118  afw::table::SourceRecord& child) const;
119 
120  virtual void measure(afw::table::SourceRecord& measRecord,
121  afw::image::Exposure<float> const& exposure) const {}
122 
123  virtual void fail(afw::table::SourceRecord& measRecord, MeasurementError* error = nullptr) const {}
124 
125 private:
126  void _measureMoments(afw::image::MaskedImage<float> const& image, afw::table::SourceRecord& child,
127  afw::table::Key<double> const& instFluxRawKey,
128  afw::table::Key<double> const& instFluxAbsKey, ShapeResultKey const& _shapeRawKey,
129  ShapeResultKey const& _shapeAbsKey) const;
130 
131  Control const _ctrl;
134  afw::table::Key<double> _instFluxChildRaw;
135  afw::table::Key<double> _instFluxParentRaw;
137  afw::table::Key<double> _instFluxChildAbs;
138  afw::table::Key<double> _instFluxParentAbs;
139  ShapeResultKey _shapeChildRaw;
140  ShapeResultKey _shapeParentRaw;
141  ShapeResultKey _shapeChildAbs;
142  ShapeResultKey _shapeParentAbs;
143  FlagHandler _flagHandler;
144 };
145 
146 } // namespace base
147 } // namespace meas
148 } // namespace lsst
149 
150 #endif // LSST_MEAS_BASE_Blendedness_h_INCLUDED
schema
table::Schema schema
Definition: Amplifier.cc:115
lsst::afw::image
Backwards-compatibility support for depersisting the old Calib (FluxMag0/FluxMag0Err) objects.
Definition: imageAlgorithm.dox:1
lsst::meas::base::SimpleAlgorithm
An abstract base classes for which the same implementation can be used for both SingleFrameAlgorithm ...
Definition: Algorithm.h:170
lsst::meas::base::BlendednessAlgorithm::FAILURE
static FlagDefinition const FAILURE
Definition: Blendedness.h:70
lsst::meas::base::BlendednessAlgorithm
Compute metrics that measure how blended objects are.
Definition: Blendedness.h:66
lsst::meas::base::BlendednessAlgorithm::measure
virtual void measure(afw::table::SourceRecord &measRecord, afw::image::Exposure< float > const &exposure) const
Called to measure a single child source in an image.
Definition: Blendedness.h:120
lsst::meas::base::BlendednessControl
Definition: Blendedness.h:40
std::string
STL class.
lsst::log.log.logContinued.error
def error(fmt, *args)
Definition: logContinued.py:213
lsst::afw::table::SourceRecord
Record class that contains measurements made on a single exposure.
Definition: Source.h:80
lsst::meas::base::BlendednessAlgorithm::measureChildPixels
void measureChildPixels(afw::image::MaskedImage< float > const &image, afw::table::SourceRecord &child) const
Definition: Blendedness.cc:346
lsst::afw::image::Exposure< float >
lsst::meas::base::MeasurementError
Exception to be thrown when a measurement algorithm experiences a known failure mode.
Definition: exceptions.h:48
lsst::afw::table::Schema
Defines the fields and offsets for a table.
Definition: Schema.h:50
lsst::meas::base::ShapeResultKey
A FunctorKey for ShapeResult.
Definition: ShapeUtilities.h:113
lsst::afw::geom.transform.transformContinued.name
string name
Definition: transformContinued.py:32
lsst::meas::base::BlendednessAlgorithm::BlendednessAlgorithm
BlendednessAlgorithm(Control const &ctrl, std::string const &name, afw::table::Schema &schema)
Definition: Blendedness.cc:195
lsst::meas::base::BlendednessControl::doOld
bool doOld
"Whether to compute HeavyFootprint dot products (the old deblend.blendedness parameter)" ;
Definition: Blendedness.h:44
Image.h
lsst::meas::base::BlendednessAlgorithm::Control
BlendednessControl Control
Definition: Blendedness.h:74
lsst::meas::base::FlagDefinitionList
vector-type utility class to build a collection of FlagDefinitions
Definition: FlagHandler.h:60
lsst::meas::base::BlendednessAlgorithm::measureParentPixels
void measureParentPixels(afw::image::MaskedImage< float > const &image, afw::table::SourceRecord &child) const
Definition: Blendedness.cc:351
lsst::meas::base::BlendednessAlgorithm::computeAbsBias
static float computeAbsBias(float mu, float variance)
Compute the bias induced by using the absolute value of a pixel instead of its value.
Definition: Blendedness.cc:270
lsst::meas::base::BlendednessAlgorithm::getFlagDefinitions
static FlagDefinitionList const & getFlagDefinitions()
Definition: Blendedness.cc:49
data
char * data
Definition: BaseRecord.cc:62
lsst::afw::image::MaskedImage
A class to manipulate images, masks, and variance as a single object.
Definition: MaskedImage.h:73
lsst::meas::base::BlendednessControl::doFlux
bool doFlux
"Whether to compute quantities related to the Gaussian-weighted flux" ;
Definition: Blendedness.h:46
lsst::meas::base::BlendednessControl::BlendednessControl
BlendednessControl()
Definition: Blendedness.h:54
lsst::meas::base::FlagHandler
Utility class for handling flag fields that indicate the failure modes of an algorithm.
Definition: FlagHandler.h:148
lsst::afw::table::Key< double >
LSST_CONTROL_FIELD
#define LSST_CONTROL_FIELD(NAME, TYPE, DOC)
A preprocessor macro used to define fields in C++ "control object" structs.
Definition: config.h:43
Algorithm.h
lsst::meas::base::BlendednessAlgorithm::fail
virtual void fail(afw::table::SourceRecord &measRecord, MeasurementError *error=nullptr) const
Handle an exception thrown by the current algorithm by setting flags in the given record.
Definition: Blendedness.h:123
Source.h
lsst
A base class for image defects.
Definition: imageAlgorithm.dox:1
ShapeUtilities.h
lsst::meas::base::FlagDefinition
Simple class used to define and document flags The name and doc constitute the identity of the FlagDe...
Definition: FlagHandler.h:40
lsst::meas::base::BlendednessAlgorithm::NO_CENTROID
static FlagDefinition const NO_CENTROID
Definition: Blendedness.h:71
variance
afw::table::Key< afw::table::Array< VariancePixelT > > variance
Definition: HeavyFootprint.cc:218
Transform.h
FlagHandler.h
config.h
lsst::meas::base::BlendednessControl::nSigmaWeightMax
double nSigmaWeightMax
"Radius factor that sets the maximum extent of the weight function (and hence the " "flux measurement...
Definition: Blendedness.h:52
lsst::meas::base::BlendednessControl::doShape
bool doShape
"Whether to compute quantities related to the Gaussian-weighted shape" ;
Definition: Blendedness.h:48
lsst::meas::base::BlendednessAlgorithm::NO_SHAPE
static FlagDefinition const NO_SHAPE
Definition: Blendedness.h:72
lsst::meas::base::BlendednessAlgorithm::computeAbsExpectation
static float computeAbsExpectation(float data, float variance)
Compute the posterior expectation value of the true instrumental flux in a pixel from its (Gaussian) ...
Definition: Blendedness.cc:260