LSST Applications  21.0.0-147-g0e635eb1+1acddb5be5,22.0.0+052faf71bd,22.0.0+1ea9a8b2b2,22.0.0+6312710a6c,22.0.0+729191ecac,22.0.0+7589c3a021,22.0.0+9f079a9461,22.0.1-1-g7d6de66+b8044ec9de,22.0.1-1-g87000a6+536b1ee016,22.0.1-1-g8e32f31+6312710a6c,22.0.1-10-gd060f87+016f7cdc03,22.0.1-12-g9c3108e+df145f6f68,22.0.1-16-g314fa6d+c825727ab8,22.0.1-19-g93a5c75+d23f2fb6d8,22.0.1-19-gb93eaa13+aab3ef7709,22.0.1-2-g8ef0a89+b8044ec9de,22.0.1-2-g92698f7+9f079a9461,22.0.1-2-ga9b0f51+052faf71bd,22.0.1-2-gac51dbf+052faf71bd,22.0.1-2-gb66926d+6312710a6c,22.0.1-2-gcb770ba+09e3807989,22.0.1-20-g32debb5+b8044ec9de,22.0.1-23-gc2439a9a+fb0756638e,22.0.1-3-g496fd5d+09117f784f,22.0.1-3-g59f966b+1e6ba2c031,22.0.1-3-g849a1b8+f8b568069f,22.0.1-3-gaaec9c0+c5c846a8b1,22.0.1-32-g5ddfab5d3+60ce4897b0,22.0.1-4-g037fbe1+64e601228d,22.0.1-4-g8623105+b8044ec9de,22.0.1-5-g096abc9+d18c45d440,22.0.1-5-g15c806e+57f5c03693,22.0.1-7-gba73697+57f5c03693,master-g6e05de7fdc+c1283a92b8,master-g72cdda8301+729191ecac,w.2021.39
LSST Data Management Base Package
InputUtilities.cc
Go to the documentation of this file.
1 // -*- lsst-c++ -*-
2 /*
3  * LSST Data Management System
4  * Copyright 2008-2014 LSST Corporation.
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 #include <cmath>
25 
26 #include "lsst/afw/table/Source.h"
30 
31 namespace lsst {
32 namespace meas {
33 namespace base {
34 
36  bool isCentroider)
37  : _name(name), _isCentroider(isCentroider) {
38  // Instead of aliasing e.g. MyAlgorithm_flag_badCentroid->slot_Centroid_flag, we actually
39  // look up the target of slot_Centroid_flag, and alias that to MyAlgorithm_flag_badCentroid.
40  // That way, if someone changes the slots later, after we've already done the measurement,
41  // this alias still points to the right thing.
42  std::string aliasedFlagName = schema.join("slot", "Centroid", "flag");
43  std::string slotFlagName = schema.getAliasMap()->apply(aliasedFlagName);
44  if (_isCentroider) {
45  if (slotFlagName != schema.join(name, "flag")) {
46  // only setup the alias if this isn't the slot algorithm itself (otherwise it'd be circular)
47  schema.getAliasMap()->set(schema.join(name, "flag", "badInitialCentroid"), slotFlagName);
48  }
49  } else {
50  if (aliasedFlagName == slotFlagName) {
51  throw LSST_EXCEPT(
53  (boost::format("Alias for '%s' must be defined before initializing '%s' plugin.") %
54  aliasedFlagName % name)
55  .str());
56  }
57  schema.getAliasMap()->set(schema.join(name, "flag", "badCentroid"), slotFlagName);
58  }
59 }
60 
61 namespace {
62 
63 geom::Point2D extractPeak(afw::table::SourceRecord const& record, std::string const& name) {
66  if (!footprint) {
67  throw LSST_EXCEPT(
69  (boost::format("%s: Centroid slot value is NaN, but no Footprint attached to record") % name)
70  .str());
71  }
72  if (footprint->getPeaks().empty()) {
73  throw LSST_EXCEPT(
75  (boost::format("%s: Centroid slot value is NaN, but Footprint has no Peaks") % name).str());
76  }
77  result.setX(footprint->getPeaks().front().getFx());
78  result.setY(footprint->getPeaks().front().getFy());
79  return result;
80 }
81 
82 } // namespace
83 
85  FlagHandler const& flags) const {
86  if (!record.getTable()->getCentroidSlot().getMeasKey().isValid()) {
87  if (_isCentroider) {
88  return extractPeak(record, _name);
89  } else {
90  throw LSST_EXCEPT(
92  (boost::format("%s requires a centroid, but the centroid slot is not defined") % _name)
93  .str());
94  }
95  }
96  geom::Point2D result = record.getCentroid();
97  if (std::isnan(result.getX()) || std::isnan(result.getY())) {
98  if (!record.getTable()->getCentroidSlot().getFlagKey().isValid()) {
99  if (_isCentroider) {
100  return extractPeak(record, _name);
101  } else {
102  throw LSST_EXCEPT(
104  (boost::format(
105  "%s: Centroid slot value is NaN, but there is no Centroid slot flag "
106  "(is the executionOrder for %s lower than that of the slot Centroid?)") %
107  _name % _name)
108  .str());
109  }
110  }
111  if (!record.getCentroidFlag() && !_isCentroider) {
112  throw LSST_EXCEPT(
114  (boost::format("%s: Centroid slot value is NaN, but the Centroid slot flag is not set "
115  "(is the executionOrder for %s lower than that of the slot Centroid?)") %
116  _name % _name)
117  .str());
118  }
119  result = extractPeak(record, _name);
120  if (!_isCentroider) {
121  // set the general flag, because using the Peak might affect the current measurement
122  flags.setValue(record, flags.getFailureFlagNumber(), true);
123  }
124  } else if (!_isCentroider && record.getTable()->getCentroidSlot().getFlagKey().isValid() &&
125  record.getCentroidFlag()) {
126  // we got a usable value, but the centroid flag is still be set, and that might affect
127  // the current measurement
128  flags.setValue(record, flags.getFailureFlagNumber(), true);
129  }
130  return result;
131 }
132 
134  // Instead of aliasing e.g. MyAlgorithm_flag_badShape->slot_Shape_flag, we actually
135  // look up the target of slot_Shape_flag, and alias that to MyAlgorithm_flag_badCentroid.
136  // That way, if someone changes the slots later, after we've already done the measurement,
137  // this alias still points to the right thing.
138  std::string aliasedFlagName = schema.join("slot", "Shape", "flag");
139  std::string slotFlagName = schema.getAliasMap()->apply(aliasedFlagName);
140  if (aliasedFlagName == slotFlagName) {
142  (boost::format("Alias for '%s' must be defined before initializing '%s' plugin.") %
143  aliasedFlagName % name)
144  .str());
145  }
146  schema.getAliasMap()->set(schema.join(name, "flag", "badShape"), slotFlagName);
147 }
148 
150  FlagHandler const& flags) const {
151  if (!record.getTable()->getCentroidSlot().getMeasKey().isValid()) {
152  throw LSST_EXCEPT(
154  (boost::format("%s requires a shape, but the shape slot is not defined") % _name).str());
155  }
157  if (std::isnan(result.getIxx()) || std::isnan(result.getIyy()) || std::isnan(result.getIxy()) ||
158  result.getIxx() * result.getIyy() < (1.0 + 1.0e-6) * result.getIxy() * result.getIxy()
159  // We are checking that Ixx*Iyy > (1 + epsilon)*Ixy*Ixy where epsilon is suitably small. The
160  // value of epsilon used here is a magic number. DM-5801 is supposed to figure out if we are
161  // to keep this value.
162  ) {
163  if (!record.getTable()->getShapeSlot().getFlagKey().isValid()) {
164  throw LSST_EXCEPT(
166  (boost::format("%s: Shape slot value is NaN, but there is no Shape slot flag "
167  "(is the executionOrder for %s lower than that of the slot Shape?)") %
168  _name % _name)
169  .str());
170  }
171  if (!record.getShapeFlag()) {
172  throw LSST_EXCEPT(
174  (boost::format("%s: Shape slot value is NaN, but the Shape slot flag is not set "
175  "(is the executionOrder for %s lower than that of the slot Shape?)") %
176  _name % _name)
177  .str());
178  }
179  throw LSST_EXCEPT(
181  (boost::format("%s: Shape needed, and Shape slot measurement failed.") % _name).str(),
182  flags.getFailureFlagNumber());
183  } else if (record.getTable()->getShapeSlot().getFlagKey().isValid() && record.getShapeFlag()) {
184  // we got a usable value, but the shape flag might still be set, and that might affect
185  // the current measurement
186  flags.setValue(record, flags.getFailureFlagNumber(), true);
187  }
188  return result;
189 }
190 
191 } // namespace base
192 } // namespace meas
193 } // namespace lsst
py::object result
Definition: _schema.cc:429
table::Key< std::string > name
Definition: Amplifier.cc:116
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
table::Schema schema
Definition: python.h:134
An ellipse core with quadrupole moments as parameters.
Definition: Quadrupole.h:47
Defines the fields and offsets for a table.
Definition: Schema.h:51
Record class that contains measurements made on a single exposure.
Definition: Source.h:78
std::shared_ptr< SourceTable const > getTable() const
Definition: Source.h:102
bool getShapeFlag() const
Return true if the measurement in the Shape slot failed.
Definition: Source.h:590
CentroidSlotDefinition::MeasValue getCentroid() const
Get the value of the Centroid slot measurement.
Definition: Source.h:570
bool getCentroidFlag() const
Return true if the measurement in the Centroid slot failed.
Definition: Source.h:578
ShapeSlotDefinition::MeasValue getShape() const
Get the value of the Shape slot measurement.
Definition: Source.h:582
std::shared_ptr< Footprint > getFootprint() const
Definition: Source.h:98
Exception to be thrown when a measurement algorithm experiences a fatal error.
Definition: exceptions.h:76
Utility class for handling flag fields that indicate the failure modes of an algorithm.
Definition: FlagHandler.h:148
std::size_t getFailureFlagNumber() const
Get the index of the General Failure flag, if one is defined.
Definition: FlagHandler.h:286
void setValue(afw::table::BaseRecord &record, std::size_t i, bool value) const
Set the flag field corresponding to the given flag index.
Definition: FlagHandler.h:262
Exception to be thrown when a measurement algorithm experiences a known failure mode.
Definition: exceptions.h:48
SafeCentroidExtractor(afw::table::Schema &schema, std::string const &name, bool isCentroider=false)
Construct the extractor, creating a flag alias that indicates failure in the input centroid by linkin...
geom::Point2D operator()(afw::table::SourceRecord &record, FlagHandler const &flags) const
Extract a position from the given record.
afw::geom::ellipses::Quadrupole operator()(afw::table::SourceRecord &record, FlagHandler const &flags) const
Extract a shape from the given record.
SafeShapeExtractor(afw::table::Schema &schema, std::string const &name)
Construct the extractor, creating a flag alias that indicates failure in the input centroid by linkin...
Reports errors in the logical structure of the program.
Definition: Runtime.h:46
Reports errors that are due to events beyond the control of the program.
Definition: Runtime.h:104
T isnan(T... args)
def format(config, name=None, writeSourceLine=True, prefix="", verbose=False)
Definition: history.py:174
A base class for image defects.