LSSTApplications  21.0.0+75b29a8a7f,21.0.0+e70536a077,21.0.0-1-ga51b5d4+62c747d40b,21.0.0-11-ga6ea59e8e+47cba9fc36,21.0.0-2-g103fe59+914993bf7c,21.0.0-2-g1367e85+e2614ded12,21.0.0-2-g45278ab+e70536a077,21.0.0-2-g4bc9b9f+7b2b5f8678,21.0.0-2-g5242d73+e2614ded12,21.0.0-2-g54e2caa+6403186824,21.0.0-2-g7f82c8f+3ac4acbffc,21.0.0-2-g8dde007+04a6aea1af,21.0.0-2-g8f08a60+9402881886,21.0.0-2-ga326454+3ac4acbffc,21.0.0-2-ga63a54e+81dd751046,21.0.0-2-gc738bc1+5f65c6e7a9,21.0.0-2-gde069b7+26c92b3210,21.0.0-2-gecfae73+0993ddc9bd,21.0.0-2-gfc62afb+e2614ded12,21.0.0-21-gba890a8+5a4f502a26,21.0.0-23-g9966ff26+03098d1af8,21.0.0-3-g357aad2+8ad216c477,21.0.0-3-g4be5c26+e2614ded12,21.0.0-3-g6d51c4a+4d2fe0280d,21.0.0-3-g7d9da8d+75b29a8a7f,21.0.0-3-gaa929c8+522e0f12c2,21.0.0-3-ge02ed75+4d2fe0280d,21.0.0-4-g3300ddd+e70536a077,21.0.0-4-gc004bbf+eac6615e82,21.0.0-4-gccdca77+f94adcd104,21.0.0-4-gd1c1571+18b81799f9,21.0.0-5-g7b47fff+4d2fe0280d,21.0.0-5-gb155db7+d2632f662b,21.0.0-5-gdf36809+637e4641ee,21.0.0-6-g722ad07+28c848f42a,21.0.0-7-g959bb79+522e0f12c2,21.0.0-7-gfd72ab2+cf01990774,21.0.0-9-g87fb7b8d+e2ab11cdd6,w.2021.04
LSSTDataManagementBasePackage
Blendedness.cc
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 #include <cmath>
26 
27 #include "boost/math/constants/constants.hpp"
28 
35 
36 namespace lsst {
37 namespace meas {
38 namespace base {
39 namespace {
40 FlagDefinitionList flagDefinitions;
41 } // namespace
42 
43 FlagDefinition const BlendednessAlgorithm::FAILURE = flagDefinitions.addFailureFlag();
44 FlagDefinition const BlendednessAlgorithm::NO_CENTROID =
45  flagDefinitions.add("flag_noCentroid", "Object has no centroid");
46 FlagDefinition const BlendednessAlgorithm::NO_SHAPE =
47  flagDefinitions.add("flag_noShape", "Object has no shape");
48 
50 
51 namespace {
52 
53 double computeOldBlendedness(PTR(afw::detection::Footprint const) childFootprint,
54  afw::image::Image<float> const& parentImage) {
55  if (!childFootprint) {
56  throw LSST_EXCEPT(pex::exceptions::LogicError, "blendedness_old requires a Footprint.");
57  }
58 
60  childHeavy = std::dynamic_pointer_cast<afw::detection::HeavyFootprint<float> const>(childFootprint);
61 
62  if (!childHeavy) {
63  return 0.0; // if it's not a HeavyFootprint, it's not blended.
64  }
65 
66  if (!parentImage.getBBox(afw::image::PARENT).contains(childHeavy->getBBox())) {
67  throw LSST_EXCEPT(pex::exceptions::LogicError, "Child footprint extends beyond image.");
68  }
69 
70  // Iterate over all the spans in the child HeavyFootprint,
71  // along with iterators for the child pixels (from the HeavyFootprint)
72  // and parent pixels (from the Exposure).
73  typedef afw::image::Image<float>::const_x_iterator ParentPixIter;
75  auto spanIter = childHeavy->getSpans()->begin();
76  auto const spanEnd = childHeavy->getSpans()->end();
77  ChildPixIter childPixIter = childHeavy->getImageArray().begin();
78  double cp = 0.0; // child.dot(parent)
79  double cc = 0.0; // child.dot(child)
80  while (spanIter != spanEnd) {
81  afw::geom::Span const& span = *spanIter;
82  ParentPixIter parentPixIter =
83  parentImage.x_at(span.getBeginX() - parentImage.getX0(), span.getY() - parentImage.getY0());
84  int const width = span.getWidth();
85  // Iterate over the pixels within the span, updating the dot products.
86  for (int x = 0; x < width; ++parentPixIter, ++childPixIter, ++x) {
87  cp += (*childPixIter) * ((*parentPixIter) - (*childPixIter));
88  cc += (*childPixIter) * (*childPixIter);
89  }
90  ++spanIter;
91  }
92  if (cc > 0.0) {
93  return cp / cc;
94  }
95  return 0.0;
96 }
97 
98 class FluxAccumulator {
99 public:
100  FluxAccumulator() : _w(0.0), _ww(0.0), _wd(0.0) {}
101 
102  void operator()(double, double, float weight, float data) {
103  _w += weight;
104  _ww += weight * weight;
105  _wd += weight * data;
106  }
107 
108  double getFlux() const { return _w * _wd / _ww; }
109 
110 protected:
111  double _w;
112  double _ww;
113  double _wd;
114 };
115 
116 class ShapeAccumulator : public FluxAccumulator {
117 public:
118  ShapeAccumulator() : FluxAccumulator(), _wdxx(0.0), _wdyy(0.0), _wdxy(0.0) {}
119 
120  void operator()(double x, double y, float weight, float data) {
121  FluxAccumulator::operator()(x, y, weight, data);
122  _wdxx += x * x * weight * data;
123  _wdyy += y * y * weight * data;
124  _wdxy += x * y * weight * data;
125  }
126 
127  ShapeResult getShape() const {
128  // Factor of 2 corrects for bias from weight function (correct is exact for an object
129  // with a Gaussian profile.)
130  ShapeResult result;
131  result.xx = 2.0 * _wdxx / _wd;
132  result.yy = 2.0 * _wdyy / _wd;
133  result.xy = 2.0 * _wdxy / _wd;
134  return result;
135  }
136 
137 private:
138  double _wdxx;
139  double _wdyy;
140  double _wdxy;
141 };
142 
143 template <typename Accumulator>
144 void computeMoments(afw::image::MaskedImage<float> const& image, geom::Point2D const& centroid,
145  afw::geom::ellipses::Quadrupole const& shape, double nSigmaWeightMax,
146  Accumulator& accumulatorRaw, Accumulator& accumulatorAbs) {
148 
149  afw::geom::ellipses::Ellipse ellipse(shape, centroid);
150  ellipse.getCore().scale(nSigmaWeightMax);
151 
152  // To evaluate an elliptically-symmetric function, we transform points
153  // by the following transform, then evaluate a circularly-symmetric function
154  // at the transformed positions.
155  geom::LinearTransform transform = shape.getGridTransform();
156 
157  typedef afw::geom::ellipses::PixelRegion::Iterator SpanIter; // yields Spans
158  typedef afw::geom::Span::Iterator PointIter; // yields Point2I positions
159  typedef afw::image::MaskedImage<float>::const_x_iterator PixelIter; // yields pixel values
160 
161  afw::geom::ellipses::PixelRegion region(ellipse);
162  bool isContained = bbox.contains(region.getBBox());
163  SpanIter const spanEnd = region.end();
164  for (SpanIter spanIter = region.begin(); spanIter != spanEnd; ++spanIter) {
165  afw::geom::Span span = *spanIter;
166  if (!isContained) {
167  if (span.getY() < bbox.getMinY() || span.getY() > bbox.getMaxY()) {
168  continue;
169  }
170  span = afw::geom::Span(span.getY(), std::max(span.getMinX(), bbox.getMinX()),
171  std::min(span.getMaxX(), bbox.getMaxX()));
172  if (span.getMinX() > span.getMaxX()) {
173  continue;
174  }
175  }
176  PixelIter pixelIter = image.x_at(span.getBeginX() - image.getX0(), span.getY() - image.getY0());
177  PointIter const pointEnd = span.end();
178  for (PointIter pointIter = span.begin(); pointIter != pointEnd; ++pointIter, ++pixelIter) {
179  geom::Extent2D d = geom::Point2D(*pointIter) - centroid;
180  geom::Extent2D td = transform(d);
181  // use single precision for faster exp, erf
182  float weight = std::exp(static_cast<float>(-0.5 * td.computeSquaredNorm()));
183  float data = pixelIter.image();
184  accumulatorRaw(d.getX(), d.getY(), weight, data);
185  float variance = pixelIter.variance();
188  accumulatorAbs(d.getX(), d.getY(), weight, std::abs(data) - bias);
189  }
190  }
191 }
192 
193 } // namespace
194 
197  : _ctrl(ctrl) {
198  if (_ctrl.doOld) {
199  _old = schema.addField<double>(
200  schema.join(name, "old"),
201  "Blendedness from dot products: (child.dot(parent)/child.dot(child) - 1)");
202  }
203  if (_ctrl.doFlux) {
204  _raw = schema.addField<double>(
205  schema.join(name, "raw"),
206  "Measure of how much the flux is affected by neighbors: "
207  "(1 - child_instFlux/parent_instFlux). Operates on the \"raw\" pixel values.");
208  _instFluxChildRaw = schema.addField<double>(
209  schema.join(name, "raw_child_instFlux"),
210  "Instrumental flux of the child, measured with a Gaussian weight matched to the child. "
211  "Operates on the \"raw\" pixel values.", "count");
212  _instFluxParentRaw = schema.addField<double>(
213  schema.join(name, "raw_parent_instFlux"),
214  "Instrumental flux of the parent, measured with a Gaussian weight matched to the child. "
215  "Operates on the \"raw\" pixel values.", "count");
216  _abs = schema.addField<double>(
217  schema.join(name, "abs"),
218  "Measure of how much the flux is affected by neighbors: "
219  "(1 - child_instFlux/parent_instFlux). "
220  "Operates on the absolute value of the pixels to try to obtain a \"de-noised\" value. "
221  "See section 4.9.11 of Bosch et al. 2018, PASJ, 70, S5 for details.");
222  _instFluxChildAbs = schema.addField<double>(
223  schema.join(name, "abs_child_instFlux"),
224  "Instrumental flux of the child, measured with a Gaussian weight matched to the child. "
225  "Operates on the absolute value of the pixels to try to obtain a \"de-noised\" value. "
226  "See section 4.9.11 of Bosch et al. 2018, PASJ, 70, S5 for details.", "count");
227  _instFluxParentAbs = schema.addField<double>(
228  schema.join(name, "abs_parent_instFlux"),
229  "Instrumental flux of the parent, measured with a Gaussian weight matched to the child. "
230  "Operates on the absolute value of the pixels to try to obtain a \"de-noised\" value. "
231  "See section 4.9.11 of Bosch et al. 2018, PASJ, 70, S5 for details.", "count");
232  }
233  if (_ctrl.doShape) {
234  _shapeChildRaw = ShapeResultKey::addFields(
235  schema, schema.join(name, "raw_child"),
236  "Shape of the child, measured with a Gaussian weight matched to the child. "
237  "Operates on the \"raw\" pixel values.", NO_UNCERTAINTY);
238  _shapeParentRaw = ShapeResultKey::addFields(
239  schema, schema.join(name, "raw_parent"),
240  "Shape of the parent, measured with a Gaussian weight matched to the child. "
241  "Operates on the \"raw\" pixel values.", NO_UNCERTAINTY);
242  _shapeChildAbs = ShapeResultKey::addFields(
243  schema, schema.join(name, "abs_child"),
244  "Shape of the child, measured with a Gaussian weight matched to the child. "
245  "Operates on the absolute value of the pixels to try to obtain a \"de-noised\" value. "
246  "See section 4.9.11 of Bosch et al. 2018, PASJ, 70, S5 for details.",
248  _shapeParentAbs = ShapeResultKey::addFields(
249  schema, schema.join(name, "abs_parent"),
250  "Shape of the parent, measured with a Gaussian weight matched to the child. "
251  "Operates on the absolute value of the pixels to try to obtain a \"de-noised\" value. "
252  "See section 4.9.11 of Bosch et al. 2018, PASJ, 70, S5 for details.",
254  }
255  if (_ctrl.doShape || _ctrl.doFlux) {
257  }
258 }
259 
261  float normalization = 0.5f * std::erfc(-data / std::sqrt(2.0f * variance));
262  if (!(normalization > 0)) {
263  // avoid division by zero; we know the limit at data << -sigma -> 0.
264  return 0.0;
265  }
266  return data + (std::sqrt(0.5f * variance / boost::math::constants::pi<float>()) *
267  std::exp(-0.5f * (data * data) / variance) / normalization);
268 }
269 
271  return (std::sqrt(2.0f * variance / boost::math::constants::pi<float>()) *
272  std::exp(-0.5f * (mu * mu) / variance)) -
273  mu * std::erfc(mu / std::sqrt(2.0f * variance));
274 }
275 
276 void BlendednessAlgorithm::_measureMoments(afw::image::MaskedImage<float> const& image,
278  afw::table::Key<double> const& instFluxRawKey,
279  afw::table::Key<double> const& instFluxAbsKey,
280  ShapeResultKey const& _shapeRawKey,
281  ShapeResultKey const& _shapeAbsKey) const {
282  if (_ctrl.doFlux || _ctrl.doShape) {
283  if (!child.getTable()->getCentroidSlot().getMeasKey().isValid()) {
285  "Centroid Key must be defined to measure the blendedness instFlux");
286  }
287  }
288  if (_ctrl.doShape) {
289  if (!child.getTable()->getCentroidSlot().getMeasKey().isValid()) {
291  "Shape Key must be defined to measure the blendedness shape");
292  }
293  }
294  if (_ctrl.doShape || _ctrl.doFlux) {
295  bool fatal = false;
296  if (child.getTable()->getCentroidSlot().getFlagKey().isValid()) {
297  if (child.getCentroidFlag()) {
298  // don't set general flag, because even a failed centroid should
299  // just fall back to the peak, and that should be fine for this
300  // measurement.
301  _flagHandler.setValue(child, NO_CENTROID.number, true);
302  }
303  }
304  if (child.getTable()->getShapeSlot().getFlagKey().isValid()) {
305  if (child.getShapeFlag()) {
306  _flagHandler.setValue(child, NO_SHAPE.number, true);
307  _flagHandler.setValue(child, FAILURE.number, true);
308  }
309  }
310  if (!(child.getShape().getDeterminant() >= 0.0)) {
311  // shape flag should have been set already, but we're paranoid
312  _flagHandler.setValue(child, NO_SHAPE.number, true);
313  _flagHandler.setValue(child, FAILURE.number, true);
314  fatal = true;
315  }
316  if (!(std::isfinite(child.getX()) && std::isfinite(child.getY()))) {
317  // shape flag should have been set already, but we're paranoid
318  _flagHandler.setValue(child, NO_CENTROID.number, true);
319  _flagHandler.setValue(child, FAILURE.number, true);
320  fatal = true;
321  }
322  if (fatal) return;
323  }
324 
325  if (_ctrl.doShape) {
326  ShapeAccumulator accumulatorRaw;
327  ShapeAccumulator accumulatorAbs;
328  computeMoments(image, child.getCentroid(), child.getShape(), _ctrl.nSigmaWeightMax, accumulatorRaw,
329  accumulatorAbs);
330  if (_ctrl.doFlux) {
331  child.set(instFluxRawKey, accumulatorRaw.getFlux());
332  child.set(instFluxAbsKey, std::max(accumulatorAbs.getFlux(), 0.0));
333  }
334  _shapeRawKey.set(child, accumulatorRaw.getShape());
335  _shapeAbsKey.set(child, accumulatorAbs.getShape());
336  } else if (_ctrl.doFlux) {
337  FluxAccumulator accumulatorRaw;
338  FluxAccumulator accumulatorAbs;
339  computeMoments(image, child.getCentroid(), child.getShape(), _ctrl.nSigmaWeightMax, accumulatorRaw,
340  accumulatorAbs);
341  child.set(instFluxRawKey, accumulatorRaw.getFlux());
342  child.set(instFluxAbsKey, std::max(accumulatorAbs.getFlux(), 0.0));
343  }
344 }
345 
347  afw::table::SourceRecord& child) const {
348  _measureMoments(image, child, _instFluxChildRaw, _instFluxChildAbs, _shapeChildRaw, _shapeChildAbs);
349 }
350 
352  afw::table::SourceRecord& child) const {
353  if (_ctrl.doOld) {
354  child.set(_old, computeOldBlendedness(child.getFootprint(), *image.getImage()));
355  }
356  _measureMoments(image, child, _instFluxParentRaw, _instFluxParentAbs, _shapeParentRaw, _shapeParentAbs);
357  if (_ctrl.doFlux) {
358  child.set(_raw, 1.0 - child.get(_instFluxChildRaw) / child.get(_instFluxParentRaw));
359  child.set(_abs, 1.0 - child.get(_instFluxChildAbs) / child.get(_instFluxParentAbs));
360  if (child.get(_instFluxParentAbs) == 0.0) {
361  // We can get NaNs in the absolute measure if both parent and child have only negative
362  // biased-corrected instFluxes (which we clip to zero). We can't really recover from this,
363  // so we should set the flag.
364  _flagHandler.setValue(child, FAILURE.number, true);
365  }
366  }
367 }
368 
369 } // namespace base
370 } // namespace meas
371 } // namespace lsst
y
int y
Definition: SpanSet.cc:49
lsst::afw::geom::Span::Iterator
SpanPixelIterator Iterator
An iterator over points in the Span.
Definition: Span.h:50
lsst::afw::image
Backwards-compatibility support for depersisting the old Calib (FluxMag0/FluxMag0Err) objects.
Definition: imageAlgorithm.dox:1
lsst::meas::base::BlendednessAlgorithm::FAILURE
static FlagDefinition const FAILURE
Definition: Blendedness.h:70
lsst::meas::base::BlendednessControl
Definition: Blendedness.h:40
lsst::meas::base::FlagDefinition::number
std::size_t number
Definition: FlagHandler.h:54
std::string
STL class.
_w
double _w
Definition: Blendedness.cc:111
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::table::BaseRecord::get
Field< T >::Value get(Key< T > const &key) const
Return the value of a field for the given key.
Definition: BaseRecord.h:151
lsst::meas::base::NO_UNCERTAINTY
@ NO_UNCERTAINTY
Algorithm provides no uncertainy information at all.
Definition: constants.h:44
lsst::afw::geom::ellipses::PixelRegion::Iterator
std::vector< Span >::const_iterator Iterator
Iterator type used by begin() and end().
Definition: PixelRegion.h:50
lsst::afw::table::SourceRecord::getCentroid
CentroidSlotDefinition::MeasValue getCentroid() const
Get the value of the Centroid slot measurement.
Definition: Source.h:572
lsst::sphgeom::abs
Angle abs(Angle const &a)
Definition: Angle.h:106
lsst::geom::LinearTransform
A 2D linear coordinate transformation.
Definition: LinearTransform.h:69
GridTransform.h
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::geom::Point2D
Point< double, 2 > Point2D
Definition: Point.h:324
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
std::sqrt
T sqrt(T... args)
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::afw::table::SourceRecord::getFootprint
std::shared_ptr< Footprint > getFootprint() const
Definition: Source.h:100
lsst::afw::image::ImageBase::x_at
x_iterator x_at(int x, int y) const
Return an x_iterator to the point (x, y) in the image.
Definition: ImageBase.h:405
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
schema
table::Schema schema
Definition: python.h:134
lsst::meas::base::BlendednessAlgorithm::getFlagDefinitions
static FlagDefinitionList const & getFlagDefinitions()
Definition: Blendedness.cc:49
lsst::afw::image::MaskedImage::const_x_iterator
const_MaskedImageIterator< typename Image::x_iterator, typename Mask::x_iterator, typename Variance::x_iterator > const_x_iterator
A const_iterator to a row of a MaskedImage.
Definition: MaskedImage.h:560
lsst::meas::base::FlagHandler::setValue
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
lsst::afw::table::SourceRecord::getShape
ShapeSlotDefinition::MeasValue getShape() const
Get the value of the Shape slot measurement.
Definition: Source.h:584
lsst::geom::Extent::computeSquaredNorm
T computeSquaredNorm() const
Return the squared L2 norm of the Extent (x^2 + y^2 + ...).
Definition: Extent.h:240
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
std::isfinite
T isfinite(T... args)
lsst::sphgeom::detail::centroid
UnitVector3d centroid(VertexIterator const begin, VertexIterator const end)
Definition: ConvexPolygonImpl.h:48
lsst::log.log.logContinued.fatal
def fatal(fmt, *args)
Definition: logContinued.py:217
lsst::meas::base::BlendednessControl::doFlux
bool doFlux
"Whether to compute quantities related to the Gaussian-weighted flux" ;
Definition: Blendedness.h:46
lsst::afw::detection::HeavyFootprint
A set of pixels in an Image, including those pixels' actual values.
Definition: HeavyFootprint.h:49
x
double x
Definition: ChebyshevBoundedField.cc:277
lsst::meas::base::ShapeResultKey::addFields
static ShapeResultKey addFields(afw::table::Schema &schema, std::string const &name, std::string const &doc, UncertaintyEnum uncertainty, afw::table::CoordinateType coordType=afw::table::CoordinateType::PIXEL)
Add the appropriate fields to a Schema, and return a ShapeResultKey that manages them.
Definition: ShapeUtilities.cc:74
_wd
double _wd
Definition: Blendedness.cc:113
lsst::afw::table::SourceRecord::getX
double getX() const
Return the centroid slot x coordinate.
Definition: Source.h:596
lsst::geom::Box2I::contains
bool contains(Point2I const &point) const noexcept
Return true if the box contains the point.
Definition: Box.cc:114
lsst::jointcal::Iterator
FastFinder::Iterator Iterator
Definition: FastFinder.cc:178
lsst.pex::exceptions::LogicError
Reports errors in the logical structure of the program.
Definition: Runtime.h:46
lsst::afw::table::Key< double >
PTR
#define PTR(...)
Definition: base.h:41
result
py::object result
Definition: _schema.cc:429
weight
afw::table::Key< double > weight
Definition: CoaddBoundedField.cc:166
lsst::afw::table::SourceRecord::getCentroidFlag
bool getCentroidFlag() const
Return true if the measurement in the Centroid slot failed.
Definition: Source.h:580
exceptions.h
lsst
A base class for image defects.
Definition: imageAlgorithm.dox:1
LSST_EXCEPT
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
lsst::afw::table::SourceRecord::getY
double getY() const
Return the centroid slot y coordinate.
Definition: Source.h:597
std::min
T min(T... args)
lsst::afw::table::SourceRecord::getTable
std::shared_ptr< SourceTable const > getTable() const
Definition: Source.h:104
PixelRegion.h
lsst::afw::image::ImageBase::getX0
int getX0() const
Return the image's column-origin.
Definition: ImageBase.h:304
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
std::exp
T exp(T... args)
std::erfc
T erfc(T... args)
lsst::geom::Point< double, 2 >
_ww
double _ww
Definition: Blendedness.cc:112
lsst::afw::geom::ellipses::Quadrupole::getDeterminant
double getDeterminant() const
Return the determinant of the matrix representation.
Definition: Quadrupole.h:83
lsst::geom::Box2I
An integer coordinate rectangle.
Definition: Box.h:55
lsst::afw::image::PARENT
@ PARENT
Definition: ImageBase.h:94
transform
table::Key< int > transform
Definition: TransformMap.cc:299
lsst::afw::image::ImageBase::getY0
int getY0() const
Return the image's row-origin.
Definition: ImageBase.h:312
lsst::afw::table::SourceRecord::getShapeFlag
bool getShapeFlag() const
Return true if the measurement in the Shape slot failed.
Definition: Source.h:592
Blendedness.h
lsst::afw::table::BaseRecord::set
void set(Key< T > const &key, U const &value)
Set value of a field for the given key.
Definition: BaseRecord.h:164
lsst::afw::image::ImageBase::getBBox
lsst::geom::Box2I getBBox(ImageOrigin origin=PARENT) const
Definition: ImageBase.h:443
HeavyFootprint.h
std::max
T max(T... args)
lsst::meas::base::ShapeResultKey::set
virtual void set(afw::table::BaseRecord &record, ShapeResult const &value) const
Set a ShapeResult in the given record.
Definition: ShapeUtilities.cc:136
lsst::afw::detection::Footprint
Class to describe the properties of a detected object from an image.
Definition: Footprint.h:63
lsst::afw::image::Image< float >
lsst::afw::image::ImageBase< float >::const_x_iterator
_const_view_t::x_iterator const_x_iterator
A const iterator for traversing the pixels in a row.
Definition: ImageBase.h:141
lsst::geom::Extent< double, 2 >
Ellipse.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::FlagHandler::addFields
static FlagHandler addFields(afw::table::Schema &schema, std::string const &prefix, FlagDefinitionList const &flagDefs, FlagDefinitionList const &exclDefs=FlagDefinitionList::getEmptyList())
Add Flag fields to a schema, creating a FlagHandler object to manage them.
Definition: FlagHandler.cc:37
lsst::meas::base::BlendednessAlgorithm::NO_SHAPE
static FlagDefinition const NO_SHAPE
Definition: Blendedness.h:72
bbox
AmpInfoBoxKey bbox
Definition: Amplifier.cc:117
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