lsst.meas.base g62aa8f1a4b+bf10eb883e
Loading...
Searching...
No Matches
GaussianFlux.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
24#include "ndarray/eigen.h"
25
27#include "lsst/geom/Box.h"
32
33namespace lsst {
34namespace meas {
35namespace base {
36namespace {
37FlagDefinitionList flagDefinitions;
38} // namespace
39
40FlagDefinition const GaussianFluxAlgorithm::FAILURE = flagDefinitions.addFailureFlag();
41
43
45 afw::table::Schema& schema)
46 : _ctrl(ctrl),
47 _instFluxResultKey(FluxResultKey::addFields(schema, name, "instFlux from Gaussian Flux algorithm")),
48 _centroidExtractor(schema, name),
49 _shapeExtractor(schema, name) {
50 _flagHandler = FlagHandler::addFields(schema, name, getFlagDefinitions());
51}
52
54 afw::image::Exposure<float> const& exposure) const {
55 geom::Point2D centroid = _centroidExtractor(measRecord, _flagHandler);
56 afw::geom::ellipses::Quadrupole shape = _shapeExtractor(measRecord, _flagHandler);
57
59 SdssShapeAlgorithm::computeFixedMomentsFlux(exposure.getMaskedImage(), shape, centroid);
60
61 measRecord.set(_instFluxResultKey, result);
62 _flagHandler.setValue(measRecord, FAILURE.number, false);
63}
64
68
69} // namespace base
70} // namespace meas
71} // namespace lsst
vector-type utility class to build a collection of FlagDefinitions
Definition FlagHandler.h:60
void handleFailure(afw::table::BaseRecord &record, MeasurementError const *error=nullptr) const
Handle an expected or unexpected Exception thrown by a measurement algorithm.
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.
void setValue(afw::table::BaseRecord &record, std::size_t i, bool value) const
Set the flag field corresponding to the given flag index.
A FunctorKey for FluxResult.
static FlagDefinition const FAILURE
static FlagDefinitionList const & getFlagDefinitions()
GaussianFluxAlgorithm(Control const &ctrl, std::string const &name, afw::table::Schema &schema)
virtual void measure(afw::table::SourceRecord &measRecord, afw::image::Exposure< float > const &exposure) const
Called to measure a single child source in an 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.
A C++ control class to handle GaussianFluxAlgorithm's configuration.
Exception to be thrown when a measurement algorithm experiences a known failure mode.
Definition exceptions.h:48
static FluxResult computeFixedMomentsFlux(ImageT const &image, afw::geom::ellipses::Quadrupole const &shape, geom::Point2D const &position)
Compute the instFlux within a fixed Gaussian aperture.
Definition SdssShape.cc:834
A reusable result struct for instFlux measurements.