LSSTApplications  10.0+286,10.0+36,10.0+46,10.0-2-g4f67435,10.1+152,10.1+37,11.0,11.0+1,11.0-1-g47edd16,11.0-1-g60db491,11.0-1-g7418c06,11.0-2-g04d2804,11.0-2-g68503cd,11.0-2-g818369d,11.0-2-gb8b8ce7
LSSTDataManagementBasePackage
Public Types | Public Member Functions | Static Public Member Functions | Private Attributes | List of all members
lsst::meas::base::GaussianCentroidAlgorithm Class Reference

A class that calculates a centroid by fitting a circular Gaussian to the image. More...

#include <GaussianCentroid.h>

Inheritance diagram for lsst::meas::base::GaussianCentroidAlgorithm:
lsst::meas::base::SimpleAlgorithm lsst::meas::base::SingleFrameAlgorithm lsst::meas::base::ForcedAlgorithm lsst::meas::base::BaseAlgorithm lsst::meas::base::BaseAlgorithm

Public Types

enum  { FAILURE =FlagHandler::FAILURE, NO_PEAK, N_FLAGS }
 Flag bits to be used with the 'flags' data member of the Result object. More...
 
typedef GaussianCentroidControl Control
 

Public Member Functions

 GaussianCentroidAlgorithm (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
 
virtual void fail (afw::table::SourceRecord &measRecord, MeasurementError *error=NULL) const
 
- Public Member Functions inherited from lsst::meas::base::SimpleAlgorithm
virtual void measureForced (afw::table::SourceRecord &measRecord, afw::image::Exposure< float > const &exposure, afw::table::SourceRecord const &refRecord, afw::image::Wcs const &refWcs) const
 
virtual void measureNForced (afw::table::SourceCatalog const &measCat, afw::image::Exposure< float > const &exposure, afw::table::SourceCatalog const &refRecord, afw::image::Wcs const &refWcs) const
 
- Public Member Functions inherited from lsst::meas::base::SingleFrameAlgorithm
virtual void measureN (afw::table::SourceCatalog const &measCat, afw::image::Exposure< float > const &exposure) const
 
- Public Member Functions inherited from lsst::meas::base::BaseAlgorithm
virtual ~BaseAlgorithm ()
 

Static Public Member Functions

template<typename PixelT >
static afw::geom::Point2D fitCentroid (afw::image::Image< PixelT > const &im, double x0, double y0)
 x0, y0 is an initial guess for position, column More...
 

Private Attributes

Control _ctrl
 
CentroidResultKey _centroidKey
 
FlagHandler _flagHandler
 
SafeCentroidExtractor _centroidExtractor
 

Detailed Description

A class that calculates a centroid by fitting a circular Gaussian to the image.

Definition at line 82 of file GaussianCentroid.h.

Member Typedef Documentation

A typedef to the Control object for this algorithm, defined above. The control object contains the configuration parameters for this algorithm.

Definition at line 98 of file GaussianCentroid.h.

Member Enumeration Documentation

anonymous enum

Flag bits to be used with the 'flags' data member of the Result object.

Inspect getFlagDefinitions() for more detailed explanations of each flag.

Enumerator
FAILURE 
NO_PEAK 
N_FLAGS 

Definition at line 90 of file GaussianCentroid.h.

Constructor & Destructor Documentation

lsst::meas::base::GaussianCentroidAlgorithm::GaussianCentroidAlgorithm ( Control const &  ctrl,
std::string const &  name,
afw::table::Schema schema 
)

Definition at line 512 of file GaussianCentroid.cc.

516  : _ctrl(ctrl),
517  _centroidKey(
518  CentroidResultKey::addFields(schema, name, "centroid from Gaussian Centroid algorithm", NO_UNCERTAINTY)
519  ),
521 {
523  getFlagDefinitions().begin(), getFlagDefinitions().end());
524 }
table::Key< std::string > name
Definition: ApCorrMap.cc:71
afw::table::Schema schema
Definition: GaussianPsf.cc:41
static CentroidResultKey addFields(afw::table::Schema &schema, std::string const &name, std::string const &doc, UncertaintyEnum uncertainty)
Add the appropriate fields to a Schema, and return a CentroidResultKey that manages them...
Algorithm provides no uncertainy information at all.
Definition: constants.h:42
static FlagHandler addFields(afw::table::Schema &schema, std::string const &prefix, FlagDefinition const *begin, FlagDefinition const *end)
Definition: FlagHandler.cc:28

Member Function Documentation

void lsst::meas::base::GaussianCentroidAlgorithm::fail ( afw::table::SourceRecord measRecord,
MeasurementError error = NULL 
) const
virtual

Handle an exception thrown by the current algorithm by setting flags in the given record.

fail() is called by the measurement framework when an exception is allowed to propagate out of one the algorithm's measure() methods. It should generally set both a general failure flag for the algorithm as well as a specific flag indicating the error condition, if possible. To aid in this, if the exception was an instance of MeasurementError, it will be passed in, carrying information about what flag to set.

An algorithm can also to chose to set flags within its own measure() methods, and then just return, rather than throw an exception. However, fail() should be implemented even when all known failure modes do not throw exceptions, to ensure that unexpected exceptions thrown in lower-level code are properly handled.

Implements lsst::meas::base::BaseAlgorithm.

Definition at line 573 of file GaussianCentroid.cc.

573  {
574  _flagHandler.handleFailure(measRecord, error);
575 }
def error
Definition: log.py:108
void handleFailure(afw::table::BaseRecord &record, MeasurementError const *error=NULL) const
Definition: FlagHandler.cc:59
template<typename PixelT >
template afw::geom::Point2D lsst::meas::base::GaussianCentroidAlgorithm::fitCentroid ( afw::image::Image< PixelT > const &  im,
double  x0,
double  y0 
)
static

x0, y0 is an initial guess for position, column

Compute centroids with 2-D Gaussian fitter

Definition at line 528 of file GaussianCentroid.cc.

531  {
532  base::FittedModel fit = base::twodg(image, x, y);
533  double const xCen = fit.params[base::FittedModel::X0];
534  double const yCen = fit.params[base::FittedModel::Y0];
535  return afw::geom::Point2D(xCen, yCen);
536 }
int y
Point< double, 2 > Point2D
Definition: Point.h:286
table::Key< table::Array< Kernel::Pixel > > image
Definition: FixedKernel.cc:117
double x
void lsst::meas::base::GaussianCentroidAlgorithm::measure ( afw::table::SourceRecord measRecord,
afw::image::Exposure< float > const &  exposure 
) const
virtual

Called to measure a single child source in an image.

Before this method is called, all neighbors will be replaced with noise, using the outputs of the deblender. Outputs should be saved in the given SourceRecord, which can also be used to obtain centroid (see SafeCentroidExtractor) and shape (see SafeShapeExtractor) information.

Implements lsst::meas::base::SingleFrameAlgorithm.

Definition at line 539 of file GaussianCentroid.cc.

542  {
543  // get our current best guess about the centroid: either a centroider measurement or peak.
544  afw::geom::Point2D center = _centroidExtractor(measRecord, _flagHandler);
545  CentroidResult result;
546  result.x = center.getX();
547  result.y = center.getY();
548  measRecord.set(_centroidKey, result); // better than NaN
549  typedef afw::image::Image<float> ImageT;
550  ImageT const& image = *exposure.getMaskedImage().getImage();
551 
552  int x = static_cast<int>(center.getX() + 0.5);
553  int y = static_cast<int>(center.getY() + 0.5);
554 
555  x -= image.getX0(); // work in image Pixel coordinates
556  y -= image.getY0();
557 
558  FittedModel fit = twodg(image, x, y); // here's the fitter
559  if (fit.params[FittedModel::PEAK] <= 0) {
560  throw LSST_EXCEPT(
561  MeasurementError,
563  NO_PEAK
564  );
565  }
566  result.x = lsst::afw::image::indexToPosition(image.getX0()) + fit.params[FittedModel::X0];
567  result.y = lsst::afw::image::indexToPosition(image.getY0()) + fit.params[FittedModel::Y0];
568  _flagHandler.setValue(measRecord, FAILURE, false); // if we had a suspect flag, we'd set that instead
569  measRecord.set(_centroidKey, result);
570 }
int y
double indexToPosition(double ind)
Convert image index to image position.
Definition: ImageUtils.h:54
Point< double, 2 > Point2D
Definition: Point.h:286
table::Key< table::Array< Kernel::Pixel > > image
Definition: FixedKernel.cc:117
FlagDefinition getDefinition(int i) const
Definition: FlagHandler.h:66
void setValue(afw::table::BaseRecord &record, int i, bool value) const
Definition: FlagHandler.h:72
double x
#define LSST_EXCEPT(type,...)
Definition: Exception.h:46

Member Data Documentation

SafeCentroidExtractor lsst::meas::base::GaussianCentroidAlgorithm::_centroidExtractor
private

Definition at line 127 of file GaussianCentroid.h.

CentroidResultKey lsst::meas::base::GaussianCentroidAlgorithm::_centroidKey
private

Definition at line 125 of file GaussianCentroid.h.

Control lsst::meas::base::GaussianCentroidAlgorithm::_ctrl
private

Definition at line 124 of file GaussianCentroid.h.

FlagHandler lsst::meas::base::GaussianCentroidAlgorithm::_flagHandler
private

Definition at line 126 of file GaussianCentroid.h.


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