LSST Applications 27.0.0,g0265f82a02+469cd937ee,g02d81e74bb+21ad69e7e1,g1470d8bcf6+cbe83ee85a,g2079a07aa2+e67c6346a6,g212a7c68fe+04a9158687,g2305ad1205+94392ce272,g295015adf3+81dd352a9d,g2bbee38e9b+469cd937ee,g337abbeb29+469cd937ee,g3939d97d7f+72a9f7b576,g487adcacf7+71499e7cba,g50ff169b8f+5929b3527e,g52b1c1532d+a6fc98d2e7,g591dd9f2cf+df404f777f,g5a732f18d5+be83d3ecdb,g64a986408d+21ad69e7e1,g858d7b2824+21ad69e7e1,g8a8a8dda67+a6fc98d2e7,g99cad8db69+f62e5b0af5,g9ddcbc5298+d4bad12328,ga1e77700b3+9c366c4306,ga8c6da7877+71e4819109,gb0e22166c9+25ba2f69a1,gb6a65358fc+469cd937ee,gbb8dafda3b+69d3c0e320,gc07e1c2157+a98bf949bb,gc120e1dc64+615ec43309,gc28159a63d+469cd937ee,gcf0d15dbbd+72a9f7b576,gdaeeff99f8+a38ce5ea23,ge6526c86ff+3a7c1ac5f1,ge79ae78c31+469cd937ee,gee10cc3b42+a6fc98d2e7,gf1cff7945b+21ad69e7e1,gfbcc870c63+9a11dc8c8f
LSST Data Management Base Package
Loading...
Searching...
No Matches
NaiveCentroid.cc
Go to the documentation of this file.
1// -*- lsst-c++ -*-
2/*
3 * LSST Data Management System
4 * Copyright 2008-2013 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#include "ndarray/eigen.h"
26
27namespace lsst {
28namespace meas {
29namespace base {
30namespace {
31FlagDefinitionList flagDefinitions;
32} // namespace
33
34FlagDefinition const NaiveCentroidAlgorithm::FAILURE = flagDefinitions.addFailureFlag();
35FlagDefinition const NaiveCentroidAlgorithm::NO_COUNTS =
36 flagDefinitions.add("flag_noCounts", "Object to be centroided has no counts");
37FlagDefinition const NaiveCentroidAlgorithm::EDGE =
38 flagDefinitions.add("flag_edge", "Object too close to edge");
39
41
43 afw::table::Schema& schema)
44 : _ctrl(ctrl),
45 _centroidKey(CentroidResultKey::addFields(schema, name, "centroid from Naive Centroid algorithm",
47 _flagHandler(FlagHandler::addFields(schema, name, getFlagDefinitions())),
48 _centroidExtractor(schema, name, true),
49 _centroidChecker(schema, name, ctrl.doFootprintCheck, ctrl.maxDistToPeak) {}
50
52 afw::image::Exposure<float> const& exposure) const {
53 geom::Point2D center = _centroidExtractor(measRecord, _flagHandler);
55 result.x = center.getX();
56 result.y = center.getY();
57 measRecord.set(_centroidKey, result); // better than NaN
58
59 typedef afw::image::Image<float> ImageT;
60 ImageT const& image = *exposure.getMaskedImage().getImage();
61
62 int x = center.getX(); // FIXME: this is different from GaussianCentroid and SdssCentroid here,
63 int y = center.getY(); // and probably shouldn't be.
64
65 x -= image.getX0(); // work in image Pixel coordinates
66 y -= image.getY0();
67
68 if (x < 1 || x >= image.getWidth() - 1 || y < 1 || y >= image.getHeight() - 1) {
69 _flagHandler.setValue(measRecord, EDGE.number, true);
70 _flagHandler.setValue(measRecord, NaiveCentroidAlgorithm::FAILURE.number, true);
71 return;
72 }
73
74 ImageT::xy_locator im = image.xy_at(x, y);
75
76 double const sum = (im(-1, 1) + im(0, 1) + im(1, 1) + im(-1, 0) + im(0, 0) + im(1, 0) + im(-1, -1) +
77 im(0, -1) + im(1, -1)) -
78 9 * _ctrl.background;
79
80 if (sum == 0.0) {
81 _flagHandler.setValue(measRecord, NO_COUNTS.number, true);
82 _flagHandler.setValue(measRecord, NaiveCentroidAlgorithm::FAILURE.number, true);
83 return;
84 }
85
86 double const sum_x = -im(-1, 1) + im(1, 1) + -im(-1, 0) + im(1, 0) + -im(-1, -1) + im(1, -1);
87 double const sum_y = (im(-1, 1) + im(0, 1) + im(1, 1)) - (im(-1, -1) + im(0, -1) + im(1, -1));
88
89 result.x = afw::image::indexToPosition(x + image.getX0()) + sum_x / sum;
90 result.y = afw::image::indexToPosition(y + image.getY0()) + sum_y / sum;
91 measRecord.set(_centroidKey, result);
92 _centroidChecker(measRecord);
93}
94
96 _flagHandler.handleFailure(measRecord, error);
97}
98
101 : CentroidTransform{name, mapper} {
104 if (flag == NaiveCentroidAlgorithm::FAILURE) continue;
105 if (mapper.getInputSchema().getNames().count(mapper.getInputSchema().join(name, flag.name)) == 0)
106 continue;
108 mapper.getInputSchema().find<afw::table::Flag>(name + "_" + flag.name).key;
109 mapper.addMapping(key);
110 }
111}
112
113} // namespace base
114} // namespace meas
115} // namespace lsst
SchemaMapper * mapper
int y
Definition SpanSet.cc:48
A class to contain the data, WCS, and other information needed to describe an image of the sky.
Definition Exposure.h:72
A class to represent a 2-dimensional array of pixels.
Definition Image.h:51
Tag types used to declare specialized field types.
Definition misc.h:31
void set(Key< T > const &key, U const &value)
Set value of a field for the given key.
Definition BaseRecord.h:164
Defines the fields and offsets for a table.
Definition Schema.h:51
A mapping between the keys of two Schemas, used to copy data between them.
Record class that contains measurements made on a single exposure.
Definition Source.h:78
A FunctorKey for CentroidResult.
Base for centroid measurement transformations.
vector-type utility class to build a collection of FlagDefinitions
Definition FlagHandler.h:60
std::size_t size() const
return the current size (number of defined elements) of the collection
Utility class for handling flag fields that indicate the failure modes of an algorithm.
void handleFailure(afw::table::BaseRecord &record, MeasurementError const *error=nullptr) const
Handle an expected or unexpected Exception thrown by a measurement algorithm.
void setValue(afw::table::BaseRecord &record, std::size_t i, bool value) const
Set the flag field corresponding to the given flag index.
Exception to be thrown when a measurement algorithm experiences a known failure mode.
Definition exceptions.h:48
static FlagDefinition const NO_COUNTS
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.
virtual void measure(afw::table::SourceRecord &measRecord, afw::image::Exposure< float > const &exposure) const
Called to measure a single child source in an image.
static FlagDefinition const FAILURE
NaiveCentroidAlgorithm(Control const &ctrl, std::string const &name, afw::table::Schema &schema)
static FlagDefinition const EDGE
static FlagDefinitionList const & getFlagDefinitions()
A C++ control class to handle NaiveCentroidAlgorithm's configuration.
double background
"Value to subtract from the image pixel values" ;
NaiveCentroidTransform(Control const &ctrl, std::string const &name, afw::table::SchemaMapper &mapper)
double indexToPosition(double ind)
Convert image index to image position.
Definition ImageUtils.h:55
@ NO_UNCERTAINTY
Algorithm provides no uncertainy information at all.
Definition constants.h:44
A reusable struct for centroid measurements.
CentroidElement x
x (column) coordinate of the measured position
Simple class used to define and document flags The name and doc constitute the identity of the FlagDe...
Definition FlagHandler.h:40