lsst::meas::base::SafeCentroidExtractor Class Reference

#include <InputUtilities.h>

Public Member Functions
	SafeCentroidExtractor (afw::table::Schema &schema, std::string const &name, bool isCentroider=false)
afw::geom::Point2D	operator() (afw::table::SourceRecord &record, FlagHandler const &flags) const

Private Attributes
std::string	_name
bool	_isCentroider

Detailed Description

Utility class for measurement algorithms that extracts a position from the Centroid slot and handles errors in a safe and consistent way.

Definition at line 37 of file InputUtilities.h.

Constructor & Destructor Documentation

lsst::meas::base::SafeCentroidExtractor::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 linking to the slot centroid flag.

Parameters

[out]	schema	Schema to which the alias should be added. The "slot_Centroid" alias must already be present in the Schema's AliasMap.
[in]	name	The name of the algorithm; the flag alias added will be "<name>_flag_badCentroid", or "<name>_flag_badInitialCentroid" if isCentroider=true.
[in]	isCentroider	Indicates whether the calling algorithm is itself a centroid measurement algorithm. If true,, falling back to the Peak because there was no previous centroider or a previous centroider failed will not cause the general failure flag of the current algorithm to be set.

Definition at line 32 of file InputUtilities.cc.

  :
    _name(name),
    _isCentroider(isCentroider)
{
    // Instead of aliasing e.g. MyAlgorithm_flag_badCentroid->slot_Centroid_flag, we actually
    // look up the target of slot_Centroid_flag, and alias that to MyAlgorithm_flag_badCentroid.
    // That way, if someone changes the slots later, after we've already done the measurement,
    // this alias still points to the right thing.
    std::string aliasedFlagName = schema.join("slot", "Centroid", "flag");
    std::string slotFlagName = schema.getAliasMap()->apply(aliasedFlagName);
    if (_isCentroider) {
        if (slotFlagName != schema.join(name, "flag")) {
            // only setup the alias if this isn't the slot algorithm itself (otherwise it'd be circular)
            schema.getAliasMap()->set(schema.join(name, "flag", "badInitialCentroid"), slotFlagName);
        }
    } else {
        if (aliasedFlagName == slotFlagName) {
            throw LSST_EXCEPT(
                pex::exceptions::LogicError,
                (boost::format("Alias for '%s' must be defined before initializing '%s' plugin.")
                % aliasedFlagName % name).str()
            );
        }
        schema.getAliasMap()->set(schema.join(name, "flag", "badCentroid"), slotFlagName);
    }
}

Member Function Documentation

afw::geom::Point2D lsst::meas::base::SafeCentroidExtractor::operator()	(	afw::table::SourceRecord &	record,
		FlagHandler const &	flags
	)		const

Extract a position from the given record.

We use the Centroid slot if it is not NaN, and fall back to the Peak on the Footprint otherwise.

If the Centroid slot is not defined, we throw FatalAlgorithmError, as this indicates a configuration problem.

If the Centroid slot value is NaN and is not flagged (or there is no Centroid slot flag), we throw RuntimeError, which should cause the measurement framework to log a warning and set the current algorithm's general failure flag if it is allowed to propagate out of the algorithm implementation.

If the Centroid slot is NaN and there is no Footprint or Peak (even if the centroid is flagged), we also throw RuntimeError, as this indicates something wrong in another stage of the pipeline that should be addressed before measurement is run.

If the Centroid slot is flagged and we nevertheless obtain a usable position (either from a the Centroid slot itself or from a successful fall-back to the Peak), we set the current algorithm's general failure flag, but return the position as well, allowing it to continue while indicating that the result may not be reliable.

Definition at line 89 of file InputUtilities.cc.

        {
    if (!record.getTable()->getCentroidKey().isValid()) {
        if (_isCentroider) {
            return extractPeak(record, _name);
        } else {
            throw LSST_EXCEPT(
                FatalAlgorithmError,
                (boost::format("%s requires a centroid, but the centroid slot is not defined") % _name).str()
            );
        }
    }
    afw::geom::Point2D result = record.getCentroid();
    if (utils::isnan(result.getX()) || utils::isnan(result.getY())) {
        if (!record.getTable()->getCentroidFlagKey().isValid()) {
            if (_isCentroider) {
                return extractPeak(record, _name);
            } else {
                throw LSST_EXCEPT(
                    pex::exceptions::RuntimeError,
                    (boost::format("%s: Centroid slot value is NaN, but there is no Centroid slot flag "
                                   "(is the executionOrder for %s lower than that of the slot Centroid?)")
                     % _name % _name).str()
                );
            }
        }
        if (!record.getCentroidFlag() && !_isCentroider) {
            throw LSST_EXCEPT(
                pex::exceptions::RuntimeError,
                (boost::format("%s: Centroid slot value is NaN, but the Centroid slot flag is not set "
                               "(is the executionOrder for %s lower than that of the slot Centroid?)")
                 % _name % _name).str()
            );
        }
        result = extractPeak(record, _name);
        if (!_isCentroider) {
            // set the general flag, because using the Peak might affect the current measurement
            flags.setValue(record, FlagHandler::FAILURE, true);
        }
    } else if (!_isCentroider && record.getTable()->getCentroidFlagKey().isValid()
               && record.getCentroidFlag()) {
        // we got a usable value, but the centroid flag is still be set, and that might affect
        // the current measurement
        flags.setValue(record, FlagHandler::FAILURE, true);
    }
    return result;
}

Member Data Documentation

bool lsst::meas::base::SafeCentroidExtractor::_isCentroider

private

Definition at line 85 of file InputUtilities.h.

std::string lsst::meas::base::SafeCentroidExtractor::_name

private

Definition at line 84 of file InputUtilities.h.

---

The documentation for this class was generated from the following files:

• /home/lsstsw/stack/Linux64/meas_base/11.0+1/include/lsst/meas/base/InputUtilities.h
• /home/lsstsw/stack/Linux64/meas_base/11.0+1/src/InputUtilities.cc