LSSTApplications  20.0.0
LSSTDataManagementBasePackage
CentroidUtilities.h
Go to the documentation of this file.
1 // -*- lsst-c++ -*-
2 /*
3  * LSST Data Management System
4  * Copyright 2008-2014 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 #ifndef LSST_MEAS_BASE_CentroidUtilities_h_INCLUDED
25 #define LSST_MEAS_BASE_CentroidUtilities_h_INCLUDED
26 
28 #include "lsst/geom/Point.h"
29 #include "lsst/afw/geom/SkyWcs.h"
33 
34 namespace lsst {
35 namespace meas {
36 namespace base {
37 
47 
50 
53  : x(x_),
54  y(y_),
55  xErr(std::sqrt(matrix(0, 0))),
56  yErr(std::sqrt(matrix(1, 1))),
57  x_y_Cov(matrix(0, 1)) {}
58 
61  : x(x_), y(y_), xErr(xErr_), yErr(yErr_), x_y_Cov(0.0) {}
62 
64  Centroid const getCentroid() const;
65 
67  void setCentroid(Centroid const& centroid);
68 
71 
73  CentroidCov const getCentroidErr() const;
74 
76  void setCentroidErr(CentroidCov const& matrix);
77 
79  void setCentroidErr(ErrElement _xErr, ErrElement _yErr);
80 };
81 
88 class CentroidResultKey : public afw::table::FunctorKey<CentroidResult> {
89 public:
102  std::string const& doc, UncertaintyEnum uncertainty);
103 
105  CentroidResultKey() : _centroid(), _centroidErr() {}
106 
110  : _centroid(centroid), _centroidErr(centroidErr) {}
111 
122 
124  virtual CentroidResult get(afw::table::BaseRecord const& record) const;
125 
127  virtual void set(afw::table::BaseRecord& record, CentroidResult const& value) const;
128 
130  bool operator==(CentroidResultKey const& other) const {
132  return _centroid == other._centroid && _centroidErr == other._centroidErr;
133  }
134  bool operator!=(CentroidResultKey const& other) const { return !(*this == other); }
136 
138  bool isValid() const { return _centroid.isValid() && _centroidErr.isValid(); }
139 
142 
145 
147  afw::table::Key<CentroidElement> getX() const { return _centroid.getX(); }
148 
150  afw::table::Key<CentroidElement> getY() const { return _centroid.getY(); }
151 
152 private:
155 };
156 
169 public:
171 
172  /*
173  * @brief Perform transformation from inputCatalog to outputCatalog.
174  *
175  * @param[in] inputCatalog Source of data to be transformed
176  * @param[in,out] outputCatalog Container for transformed results
177  * @param[in] wcs World coordinate system under which transformation will take place
178  * @param[in] photoCalib Photometric calibration under which transformation will take place
179  * @throws LengthError Catalog sizes do not match
180  */
181  virtual void operator()(afw::table::SourceCatalog const& inputCatalog,
182  afw::table::BaseCatalog& outputCatalog, afw::geom::SkyWcs const& wcs,
183  afw::image::PhotoCalib const& photoCalib) const;
184 
185 private:
186  afw::table::CoordKey _coordKey;
188 };
189 
191 public:
213  CentroidChecker(afw::table::Schema& schema, std::string const& name, bool inside = true,
214  double maxDistFromPeak = -1.0);
215 
222  bool operator()(afw::table::SourceRecord& record) const;
223 
224 private:
225  bool _doFootprintCheck;
226  double _maxDistFromPeak;
234 };
235 } // namespace base
236 } // namespace meas
237 } // namespace lsst
238 
239 #endif // !LSST_MEAS_BASE_CentroidUtilities_h_INCLUDED
lsst::meas::base::CentroidResult::setCentroidErr
void setCentroidErr(CentroidCov const &matrix)
Set the struct uncertainty fields from the given matrix, with rows and columns ordered (x,...
Definition: CentroidUtilities.cc:54
schema
table::Schema schema
Definition: Amplifier.cc:115
Transform.h
lsst::afw::table::PointKey::getX
Key< T > getX() const noexcept
Return the underlying x Key.
Definition: aggregates.h:107
lsst::meas::base::CentroidResultKey::set
virtual void set(afw::table::BaseRecord &record, CentroidResult const &value) const
Set a CentroidResult in the given record.
Definition: CentroidUtilities.cc:114
std::string
STL class.
lsst::meas::base::UncertaintyEnum
UncertaintyEnum
An enum used to specify how much uncertainty information measurement algorithms provide.
Definition: constants.h:43
lsst::meas::base::CentroidResult::getPoint
geom::Point< CentroidElement > getPoint()
Return the 2D point type corresponding to this result.
Definition: CentroidUtilities.h:70
lsst::meas::base::CentroidResultKey::getCentroid
afw::table::PointKey< CentroidElement > getCentroid() const
Return a FunctorKey to just the centroid value.
Definition: CentroidUtilities.h:141
lsst::meas::base::CentroidChecker
Definition: CentroidUtilities.h:190
lsst::afw::table::SourceRecord
Record class that contains measurements made on a single exposure.
Definition: Source.h:80
base
Definition: __init__.py:1
wcs
table::Key< table::Array< std::uint8_t > > wcs
Definition: SkyWcs.cc:71
lsst::meas::base::CentroidResult::setCentroid
void setCentroid(Centroid const &centroid)
Set the struct fields from the given Point object.
Definition: CentroidUtilities.cc:43
lsst::afw::table::FunctorKey
Convenience base class that combines the OutputFunctorKey and InputFunctorKey.
Definition: FunctorKey.h:74
lsst::afw::table::CoordKey
A FunctorKey used to get or set celestial coordinates from a pair of lsst::geom::Angle keys.
Definition: aggregates.h:210
lsst::meas::base::CentroidTransform::CentroidTransform
CentroidTransform(std::string const &name, afw::table::SchemaMapper &mapper)
Definition: CentroidUtilities.cc:121
SkyWcs.h
lsst::meas::base::CentroidResult::CentroidResult
CentroidResult(CentroidElement x_, CentroidElement y_, CentroidCov const &matrix)
Constructor; initializes everything from values.
Definition: CentroidUtilities.h:52
lsst::afw::table::Schema
Defines the fields and offsets for a table.
Definition: Schema.h:50
lsst::afw::geom::SkyWcs
A 2-dimensional celestial WCS that transform pixels to ICRS RA/Dec, using the LSST standard for pixel...
Definition: SkyWcs.h:117
lsst::meas::base::CentroidResultKey
A FunctorKey for CentroidResult.
Definition: CentroidUtilities.h:88
lsst::meas::base::CentroidResultKey::getCentroidErr
afw::table::CovarianceMatrixKey< ErrElement, 2 > getCentroidErr() const
Return a FunctorKey to just the uncertainty matrix.
Definition: CentroidUtilities.h:144
lsst::afw::geom.transform.transformContinued.name
string name
Definition: transformContinued.py:32
aggregates.h
lsst::meas::base::CentroidTransform
Base for centroid measurement transformations.
Definition: CentroidUtilities.h:168
lsst::meas::base::CentroidResult::x
CentroidElement x
x (column) coordinate of the measured position
Definition: CentroidUtilities.h:42
lsst::meas::base::CentroidResult
A reusable struct for centroid measurements.
Definition: CentroidUtilities.h:41
lsst::afw::table::CovarianceMatrixKey< ErrElement, 2 >
lsst::meas::base::CentroidResult::CentroidResult
CentroidResult()
Constructor; initializes everything to NaN.
Definition: CentroidUtilities.cc:34
lsst::meas::base::CentroidResult::x_y_Cov
ErrElement x_y_Cov
x,y term in the uncertainty convariance matrix
Definition: CentroidUtilities.h:46
lsst::meas::base::CentroidElement
double CentroidElement
Definition: constants.h:56
lsst::meas::base::CentroidResultKey::addFields
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.
Definition: CentroidUtilities.cc:66
lsst::meas::base::CentroidResultKey::get
virtual CentroidResult get(afw::table::BaseRecord const &record) const
Get a CentroidResult from the given record.
Definition: CentroidUtilities.cc:105
lsst::meas::base::CentroidResult::xErr
ErrElement xErr
standard deviation of x
Definition: CentroidUtilities.h:44
lsst::afw::table::CovarianceMatrixKey::isValid
bool isValid() const noexcept
Return True if all the constituent error Keys are valid.
Definition: aggregates.cc:294
lsst::afw::table::PointKey< CentroidElement >
lsst::afw::table::PointKey::isValid
bool isValid() const noexcept
Return True if both the x and y Keys are valid.
Definition: aggregates.h:104
lsst::sphgeom::detail::centroid
UnitVector3d centroid(VertexIterator const begin, VertexIterator const end)
Definition: ConvexPolygonImpl.h:48
lsst::meas::base::CentroidChecker::CentroidChecker
CentroidChecker(afw::table::Schema &schema, std::string const &name, bool inside=true, double maxDistFromPeak=-1.0)
Check source record produced by a centroid algorithm called "name".
Definition: CentroidUtilities.cc:174
lsst::meas::base::CentroidChecker::operator()
bool operator()(afw::table::SourceRecord &record) const
Set the centroid to the first footprint if the centroid is either more than _dist pixels from the foo...
Definition: CentroidUtilities.cc:197
lsst::afw::table._source.SourceCatalog
Definition: _source.py:33
other
ItemVariant const * other
Definition: Schema.cc:56
lsst::afw::table::BaseRecord
Base class for all records.
Definition: BaseRecord.h:31
lsst::meas::base::CentroidResultKey::operator!=
bool operator!=(CentroidResultKey const &other) const
Definition: CentroidUtilities.h:134
lsst::afw::table::SchemaMapper
A mapping between the keys of two Schemas, used to copy data between them.
Definition: SchemaMapper.h:21
lsst::afw::table::Key< CentroidElement >
lsst::meas::base::CentroidResultKey::isValid
bool isValid() const
Return True if the centroid key is valid.
Definition: CentroidUtilities.h:138
lsst::meas::base::CentroidCov
Eigen::Matrix< ErrElement, 2, 2, Eigen::DontAlign > CentroidCov
Definition: constants.h:59
lsst::afw::image::PhotoCalib
The photometric calibration of an exposure.
Definition: PhotoCalib.h:114
lsst
A base class for image defects.
Definition: imageAlgorithm.dox:1
lsst::meas::base::CentroidResult::y
CentroidElement y
y (row) coordinate of the measured position
Definition: CentroidUtilities.h:43
photoCalib
Key< int > photoCalib
Definition: Exposure.cc:67
lsst::meas::base::CentroidResult::getCentroidErr
CentroidCov const getCentroidErr() const
Return the 2x2 symmetric covariance matrix, with rows and columns ordered (x, y)
Definition: CentroidUtilities.cc:48
PhotoCalib.h
Implementation of the Photometric Calibration class.
lsst::meas::base::ErrElement
float ErrElement
Definition: constants.h:55
lsst::afw::table::SubSchema
A proxy type for name lookups in a Schema.
Definition: Schema.h:357
std
STL namespace.
lsst::meas::base::CentroidResultKey::getX
afw::table::Key< CentroidElement > getX() const
Return a Key for the x coordinate.
Definition: CentroidUtilities.h:147
lsst::geom::Point
A coordinate class intended to represent absolute positions.
Definition: CoordinateBase.h:39
lsst::meas::base::CentroidResultKey::CentroidResultKey
CentroidResultKey()
Default constructor; instance will not be usuable unless subsequently assigned to.
Definition: CentroidUtilities.h:105
mapper
SchemaMapper * mapper
Definition: SchemaMapper.cc:78
constants.h
lsst::meas::base::CentroidResultKey::CentroidResultKey
CentroidResultKey(afw::table::PointKey< CentroidElement > const &centroid, afw::table::CovarianceMatrixKey< ErrElement, 2 > const &centroidErr)
Construct from a pair of Keys.
Definition: CentroidUtilities.h:108
Point.h
lsst::meas::base::CentroidResult::CentroidResult
CentroidResult(CentroidElement x_, CentroidElement y_, ErrElement xErr_, ErrElement yErr_)
Constructor; initializes everything from values.
Definition: CentroidUtilities.h:60
lsst::meas::base::BaseTransform
Abstract base class for all C++ measurement transformations.
Definition: Transform.h:86
lsst::meas::base::CentroidResult::getCentroid
Centroid const getCentroid() const
Return a Point object containing the measured x and y.
Definition: CentroidUtilities.cc:41
lsst::meas::base::CentroidResultKey::getY
afw::table::Key< CentroidElement > getY() const
Return a Key for the y coordinate.
Definition: CentroidUtilities.h:150
lsst::afw::table::CatalogT< BaseRecord >
lsst::meas::base::CentroidResult::yErr
ErrElement yErr
standard deviation of y
Definition: CentroidUtilities.h:45
lsst::meas::base::CentroidResultKey::operator==
bool operator==(CentroidResultKey const &other) const
Compare the FunctorKey for equality with another, using the underlying Keys.
Definition: CentroidUtilities.h:131
lsst::meas::base::CentroidTransform::operator()
virtual void operator()(afw::table::SourceCatalog const &inputCatalog, afw::table::BaseCatalog &outputCatalog, afw::geom::SkyWcs const &wcs, afw::image::PhotoCalib const &photoCalib) const
Definition: CentroidUtilities.cc:144
lsst::afw::table::PointKey::getY
Key< T > getY() const noexcept
Return the underlying y Key.
Definition: aggregates.h:110