28#include "ndarray/eigen.h"
40FlagDefinitionList flagDefinitions;
45 flagDefinitions.add(
"flag_noGoodPixels",
"not enough non-rejected pixels in data to attempt the fit");
47 "flag_edge",
"object was too close to the edge of the image to use the full PSF model");
57 schema, name,
"instFlux derived from linear least-squares fit of PSF model")),
58 _areaKey(schema.addField<
float>(name +
"_area",
"effective area of PSF",
"pixel")),
59 _chi2Key(schema.addField<
float>(name +
"_chi2",
"chi2 of the fitted PSF")),
60 _npixelsKey(schema.addField<
int>(name +
"_npixels",
61 "number of pixels that were included in the PSF fit",
"pixel")),
62 _centroidExtractor(schema, name) {
77 fitBBox.clip(exposure.getBBox());
89 badBits |= exposure.getMaskedImage().getMask()->getPlaneBitMask(*
i);
92 ->intersectNot(*exposure.getMaskedImage().getMask(),
badBits)
93 ->clippedTo(exposure.getMaskedImage().getMask()->getBBox()));
104 auto dataNdArray =
fitRegion.getSpans()->flatten(exposure.getMaskedImage().getImage()->getArray(),
139 if (
mapper.getInputSchema().getNames().count(
mapper.getInputSchema().join(name, flag.
name)) == 0)
142 mapper.getInputSchema().find<afw::table::Flag>(name +
"_" + flag.
name).key;
#define LSST_EXCEPT(type,...)
#define LOGL_ERROR(logger, message...)
math::Kernel::Pixel Pixel
std::string getLogName() const
Exception to be thrown when a measurement algorithm experiences a fatal error.
vector-type utility class to build a collection of FlagDefinitions
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.
Exception to be thrown when a measurement algorithm experiences a known failure mode.
Exception to be thrown when a measurement algorithm encounters a NaN or infinite pixel.
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
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.
static FlagDefinition const NO_GOOD_PIXELS
static FlagDefinition const EDGE
static FlagDefinitionList const & getFlagDefinitions()
PsfFluxAlgorithm(Control const &ctrl, std::string const &name, afw::table::Schema &schema, std::string const &logName="")
A C++ control class to handle PsfFluxAlgorithm's configuration.
std::vector< std::string > badMaskPlanes
"Mask planes that indicate pixels that should be excluded from the fit" ;
Simple class used to define and document flags The name and doc constitute the identity of the FlagDe...
A reusable result struct for instFlux measurements.
meas::base::Flux instFlux
Measured instFlux in DN.