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
ApertureFlux.h
Go to the documentation of this file.
1 // -*- lsst-c++ -*-
2 /*
3  * LSST Data Management System
4  * Copyright 2008-2016 AURA/LSST.
5  *
6  * This product includes software developed by the
7  * LSST Project (http://www.lsst.org/).
8  *
9  * This program is free software: you can redistribute it and/or modify
10  * it under the terms of the GNU General Public License as published by
11  * the Free Software Foundation, either version 3 of the License, or
12  * (at your option) any later version.
13  *
14  * This program is distributed in the hope that it will be useful,
15  * but WITHOUT ANY WARRANTY; without even the implied warranty of
16  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17  * GNU General Public License for more details.
18  *
19  * You should have received a copy of the LSST License Statement and
20  * the GNU General Public License along with this program. If not,
21  * see <http://www.lsstcorp.org/LegalNotices/>.
22  */
23 
24 #ifndef LSST_MEAS_BASE_ApertureFlux_h_INCLUDED
25 #define LSST_MEAS_BASE_ApertureFlux_h_INCLUDED
26 
27 #include <array>
28 #include <string>
29 #include <vector>
30 #include <bitset>
31 
32 #include "lsst/pex/config.h"
34 #include "lsst/afw/table/arrays.h"
41 
42 namespace lsst {
43 namespace meas {
44 namespace base {
45 
50 public:
52 
53  LSST_CONTROL_FIELD(radii, std::vector<double>, "Radius (in pixels) of apertures.");
54 
56  maxSincRadius, double,
57  "Maximum radius (in pixels) for which the sinc algorithm should be used instead of the "
58  "faster naive algorithm. For elliptical apertures, this is the minor axis radius.");
59 
62  "Warping kernel used to shift Sinc photometry coefficients to different center positions");
63 };
64 
65 struct ApertureFluxResult;
66 
81 public:
82  // Structures and routines to manage flaghandler
83  static FlagDefinitionList const& getFlagDefinitions();
84  static unsigned int const N_FLAGS = 3;
85  static FlagDefinition const FAILURE;
88 
90 
93 
95 
106  template <typename T>
108  afw::geom::ellipses::Ellipse const& ellipse,
109  Control const& ctrl = Control());
110  template <typename T>
112  afw::geom::ellipses::Ellipse const& ellipse,
113  Control const& ctrl = Control());
115 
117 
126  template <typename T>
128  afw::geom::ellipses::Ellipse const& ellipse,
129  Control const& ctrl = Control());
130  template <typename T>
132  afw::geom::ellipses::Ellipse const& ellipse,
133  Control const& ctrl = Control());
135 
137 
148  template <typename T>
150  Control const& ctrl = Control());
151 
152  template <typename T>
154  afw::geom::ellipses::Ellipse const& ellipse, Control const& ctrl = Control());
156 
161  daf::base::PropertySet& metadata);
162 
164 
173  virtual void measure(afw::table::SourceRecord& record,
174  afw::image::Exposure<float> const& exposure) const = 0;
175 
177  virtual void fail(afw::table::SourceRecord& measRecord, MeasurementError* error = nullptr) const;
178 
191  static std::string makeFieldPrefix(std::string const& name, double radius);
192 
193 protected:
194  void copyResultToRecord(Result const& result, afw::table::SourceRecord& record, int index) const;
195 
196  FlagHandler const& getFlagHandler(int index) const { return _keys[index].flags; }
197 
198  Control const _ctrl;
200 
201 private:
202  struct Keys {
203  FluxResultKey instFluxKey;
204  FlagHandler flags;
205 
206  Keys(afw::table::Schema& schema, std::string const& prefix, std::string const& doc, bool isSinc);
207  };
208 
209  std::vector<Keys> _keys;
210 };
211 
219  bool getFlag(unsigned int index) const { return _flags[index]; }
220 
222  bool getFlag(std::string const& name) const {
224  }
225 
227  void setFlag(unsigned int index, bool value = true) { _flags[index] = value; }
228 
230  void unsetFlag(unsigned int index) { _flags[index] = false; }
231 
232 private:
234 };
235 
243 public:
246 
247  /*
248  * @brief Perform transformation from inputCatalog to outputCatalog.
249  *
250  * @param[in] inputCatalog Source of data to be transformed
251  * @param[in,out] outputCatalog Container for transformed results
252  * @param[in] wcs World coordinate system under which transformation will take place
253  * @param[in] photoCalib Photometric calibration under which transformation will take place
254  * @throws LengthError Catalog sizes do not match
255  */
256  virtual void operator()(afw::table::SourceCatalog const& inputCatalog,
257  afw::table::BaseCatalog& outputCatalog, afw::geom::SkyWcs const& wcs,
258  afw::image::PhotoCalib const& photoCalib) const;
259 
260 private:
261  std::vector<MagResultKey> _magKeys;
262  Control _ctrl;
263 };
264 
265 } // namespace base
266 } // namespace meas
267 } // namespace lsst
268 
269 #endif // !LSST_MEAS_BASE_ApertureFlux_h_INCLUDED
py::object result
Definition: _schema.cc:430
table::Key< std::string > name
Definition: Amplifier.cc:116
std::string prefix
Definition: SchemaMapper.cc:79
SchemaMapper * mapper
Definition: SchemaMapper.cc:78
table::Key< table::Array< std::uint8_t > > wcs
Definition: SkyWcs.cc:71
table::Schema schema
Definition: python.h:134
A 2-dimensional celestial WCS that transform pixels to ICRS RA/Dec, using the LSST standard for pixel...
Definition: SkyWcs.h:117
An ellipse defined by an arbitrary BaseCore and a center point.
Definition: Ellipse.h:51
A class to represent a 2-dimensional array of pixels.
Definition: Image.h:58
A class to manipulate images, masks, and variance as a single object.
Definition: MaskedImage.h:73
The photometric calibration of an exposure.
Definition: PhotoCalib.h:114
Defines the fields and offsets for a table.
Definition: Schema.h:50
A mapping between the keys of two Schemas, used to copy data between them.
Definition: SchemaMapper.h:21
Record class that contains measurements made on a single exposure.
Definition: Source.h:80
Class for storing generic metadata.
Definition: PropertySet.h:67
Base class for multiple-aperture photometry algorithms.
Definition: ApertureFlux.h:80
SafeCentroidExtractor _centroidExtractor
Definition: ApertureFlux.h:199
ApertureFluxAlgorithm(Control const &ctrl, std::string const &name, afw::table::Schema &schema, daf::base::PropertySet &metadata)
Construct the algorithm and add its fields to the given Schema.
Definition: ApertureFlux.cc:71
void copyResultToRecord(Result const &result, afw::table::SourceRecord &record, int index) const
Definition: ApertureFlux.cc:96
static FlagDefinitionList const & getFlagDefinitions()
Definition: ApertureFlux.cc:49
static FlagDefinition const FAILURE
Definition: ApertureFlux.h:85
static unsigned int const N_FLAGS
Definition: ApertureFlux.h:84
static FlagDefinition const SINC_COEFFS_TRUNCATED
Definition: ApertureFlux.h:87
static std::string makeFieldPrefix(std::string const &name, double radius)
Construct an appropriate prefix for table fields.
Definition: ApertureFlux.cc:57
static Result computeSincFlux(afw::image::Image< T > const &image, afw::geom::ellipses::Ellipse const &ellipse, Control const &ctrl=Control())
Compute the instFlux (and optionally, uncertanties) within an aperture using Sinc photometry.
static Result computeFlux(afw::image::Image< T > const &image, afw::geom::ellipses::Ellipse const &ellipse, Control const &ctrl=Control())
Compute the instFlux (and optionally, uncertanties) within an aperture using the algorithm determined...
static FlagDefinition const APERTURE_TRUNCATED
Definition: ApertureFlux.h:86
ApertureFluxResult Result
Result object returned by static methods.
Definition: ApertureFlux.h:92
FlagHandler const & getFlagHandler(int index) const
Definition: ApertureFlux.h:196
static Result computeNaiveFlux(afw::image::Image< T > const &image, afw::geom::ellipses::Ellipse const &ellipse, Control const &ctrl=Control())
Compute the instFlux (and optionally, uncertanties) within an aperture using naive photometry.
virtual void measure(afw::table::SourceRecord &record, afw::image::Exposure< float > const &exposure) const =0
Measure the configured apertures on the given image.
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: ApertureFlux.cc:87
Configuration object for multiple-aperture flux algorithms.
Definition: ApertureFlux.h:49
std::string shiftKernel
"Warping kernel used to shift Sinc photometry coefficients to different center positions" ;
Definition: ApertureFlux.h:62
double maxSincRadius
"Maximum radius (in pixels) for which the sinc algorithm should be used instead of the " "faster naiv...
Definition: ApertureFlux.h:58
std::vector< double > radii
"Radius (in pixels) of apertures." ;
Definition: ApertureFlux.h:53
Measurement transformation for aperture fluxes.
Definition: ApertureFlux.h:242
ApertureFluxTransform(Control const &ctrl, std::string const &name, afw::table::SchemaMapper &mapper)
virtual void operator()(afw::table::SourceCatalog const &inputCatalog, afw::table::BaseCatalog &outputCatalog, afw::geom::SkyWcs const &wcs, afw::image::PhotoCalib const &photoCalib) const
Abstract base class for all C++ measurement transformations.
Definition: Transform.h:86
vector-type utility class to build a collection of FlagDefinitions
Definition: FlagHandler.h:60
FlagDefinition getDefinition(std::size_t index) const
get a reference to the FlagDefinition with specified index.
Definition: FlagHandler.h:83
Utility class for handling flag fields that indicate the failure modes of an algorithm.
Definition: FlagHandler.h:148
A FunctorKey for FluxResult.
Definition: FluxUtilities.h:59
Exception to be thrown when a measurement algorithm experiences a known failure mode.
Definition: exceptions.h:48
Utility class for measurement algorithms that extracts a position from the Centroid slot and handles ...
An abstract base classes for which the same implementation can be used for both SingleFrameAlgorithm ...
Definition: Algorithm.h:170
#define LSST_CONTROL_FIELD(NAME, TYPE, DOC)
A preprocessor macro used to define fields in C++ "control object" structs.
Definition: config.h:43
This defines the base of measurement transformations.
Backwards-compatibility support for depersisting the old Calib (FluxMag0/FluxMag0Err) objects.
A base class for image defects.
A Result struct for running an aperture flux algorithm with a single radius.
Definition: ApertureFlux.h:217
bool getFlag(std::string const &name) const
Return the flag value associated with the given flag name.
Definition: ApertureFlux.h:222
bool getFlag(unsigned int index) const
Return the flag value associated with the given bit.
Definition: ApertureFlux.h:219
void setFlag(unsigned int index, bool value=true)
Set the flag value associated with the given bit.
Definition: ApertureFlux.h:227
void unsetFlag(unsigned int index)
Clear (i.e. set to false) the flag associated with the given bit.
Definition: ApertureFlux.h:230
Simple class used to define and document flags The name and doc constitute the identity of the FlagDe...
Definition: FlagHandler.h:40
A reusable result struct for instFlux measurements.
Definition: FluxUtilities.h:41
Key< int > photoCalib
Definition: Exposure.cc:67