LSSTApplications
19.0.0-14-gb0260a2+72efe9b372,20.0.0+7927753e06,20.0.0+8829bf0056,20.0.0+995114c5d2,20.0.0+b6f4b2abd1,20.0.0+bddc4f4cbe,20.0.0-1-g253301a+8829bf0056,20.0.0-1-g2b7511a+0d71a2d77f,20.0.0-1-g5b95a8c+7461dd0434,20.0.0-12-g321c96ea+23efe4bbff,20.0.0-16-gfab17e72e+fdf35455f6,20.0.0-2-g0070d88+ba3ffc8f0b,20.0.0-2-g4dae9ad+ee58a624b3,20.0.0-2-g61b8584+5d3db074ba,20.0.0-2-gb780d76+d529cf1a41,20.0.0-2-ged6426c+226a441f5f,20.0.0-2-gf072044+8829bf0056,20.0.0-2-gf1f7952+ee58a624b3,20.0.0-20-geae50cf+e37fec0aee,20.0.0-25-g3dcad98+544a109665,20.0.0-25-g5eafb0f+ee58a624b3,20.0.0-27-g64178ef+f1f297b00a,20.0.0-3-g4cc78c6+e0676b0dc8,20.0.0-3-g8f21e14+4fd2c12c9a,20.0.0-3-gbd60e8c+187b78b4b8,20.0.0-3-gbecbe05+48431fa087,20.0.0-38-ge4adf513+a12e1f8e37,20.0.0-4-g97dc21a+544a109665,20.0.0-4-gb4befbc+087873070b,20.0.0-4-gf910f65+5d3db074ba,20.0.0-5-gdfe0fee+199202a608,20.0.0-5-gfbfe500+d529cf1a41,20.0.0-6-g64f541c+d529cf1a41,20.0.0-6-g9a5b7a1+a1cd37312e,20.0.0-68-ga3f3dda+5fca18c6a4,20.0.0-9-g4aef684+e18322736b,w.2020.45
LSSTDataManagementBasePackage
|
Go to the documentation of this file.
26 #include "boost/algorithm/string.hpp"
28 #include "ndarray/eigen.h"
40 FlagDefinitionList flagDefinitions;
45 flagDefinitions.add(
"flag_apertureTruncated",
"aperture did not fit within measurement image");
47 "flag_sincCoeffsTruncated",
"full sinc coefficient image did not fit within measurement image");
53 static std::array<double, 10> defaultRadii = {{3.0, 4.5, 6.0, 9.0, 12.0, 17.0, 25.0, 35.0, 50.0, 70.0}};
59 return boost::replace_all_copy(
prefix,
".",
"_");
75 : _ctrl(ctrl), _centroidExtractor(
schema,
name) {
79 boost::to_upper(upperName);
80 metadata.
add(upperName +
"_RADII", ctrl.
radii[i]);
92 _keys[i].flags.handleFailure(measRecord,
error);
98 record.
set(_keys[index].instFluxKey,
result);
114 template <
typename T>
124 if (!
bbox.contains(cImage->getBBox())) {
138 cImage = std::make_shared<afw::image::Image<T> >(*cImage, overlap);
145 template <
typename T>
153 (ndarray::asEigenArray(subImage.
getArray()) * ndarray::asEigenArray(cImage->getArray())).sum();
157 template <
typename T>
165 result.instFlux = (ndarray::asEigenArray(subImage.
getImage()->getArray()) *
166 ndarray::asEigenArray(cImage->getArray()))
170 ndarray::asEigenArray(cImage->getArray()).square())
175 template <
typename T>
187 spanIter != spanEnd; ++spanIter) {
189 image.x_at(spanIter->getBeginX() -
image.getX0(), spanIter->getY() -
image.getY0());
195 template <
typename T>
209 spanIter != spanEnd; ++spanIter) {
211 spanIter->getBeginX() -
image.getX0(), spanIter->getY() -
image.getY0());
213 spanIter->getBeginX() -
image.getX0(), spanIter->getY() -
image.getY0());
222 template <
typename T>
231 template <
typename T>
239 #define INSTANTIATE(T) \
240 template ApertureFluxAlgorithm::Result ApertureFluxAlgorithm::computeFlux( \
241 afw::image::Image<T> const &, afw::geom::ellipses::Ellipse const &, Control const &); \
242 template ApertureFluxAlgorithm::Result ApertureFluxAlgorithm::computeFlux( \
243 afw::image::MaskedImage<T> const &, afw::geom::ellipses::Ellipse const &, Control const &); \
244 template ApertureFluxAlgorithm::Result ApertureFluxAlgorithm::computeSincFlux( \
245 afw::image::Image<T> const &, afw::geom::ellipses::Ellipse const &, Control const &); \
246 template ApertureFluxAlgorithm::Result ApertureFluxAlgorithm::computeSincFlux( \
247 afw::image::MaskedImage<T> const &, afw::geom::ellipses::Ellipse const &, Control const &); \
248 template ApertureFluxAlgorithm::Result ApertureFluxAlgorithm::computeNaiveFlux( \
249 afw::image::Image<T> const &, afw::geom::ellipses::Ellipse const &, Control const &); \
250 template ApertureFluxAlgorithm::Result ApertureFluxAlgorithm::computeNaiveFlux( \
251 afw::image::MaskedImage<T> const &, afw::geom::ellipses::Ellipse const &, Control const &)
268 .find<afw::table::Flag>(
291 for (; inSrc != inputCatalog.end() && outSrc != outputCatalog.
end(); ++inSrc, ++outSrc) {
293 FluxResult instFluxResult = instFluxKeys[i].get(*inSrc);
Backwards-compatibility support for depersisting the old Calib (FluxMag0/FluxMag0Err) objects.
static Result computeFlux(afw::image::Image< T > const &image, afw::geom::ellipses::Ellipse const &ellipse, Control const &ctrl=Control())
Compute the instFlux (and optionally, uncertanties) within an aperture using the algorithm determined...
lsst::geom::Box2I const & getBBox() const
Return the bounding box of the pixel region.
meas::base::Flux instFlux
Measured instFlux in DN.
A reusable result struct for instFlux measurements.
static FlagDefinition const APERTURE_TRUNCATED
Record class that contains measurements made on a single exposure.
std::vector< Span >::const_iterator Iterator
Iterator type used by begin() and end().
double const getB() const
table::Key< table::Array< std::uint8_t > > wcs
lsst::geom::Box2D computeBBox() const
Return the bounding box of the ellipse.
void add(std::string const &name, T const &value)
Append a single value to the vector of values for a property name (possibly hierarchical).
Exception to be thrown when a measurement algorithm experiences a known failure mode.
Defines the fields and offsets for a table.
An ellipse core for the semimajor/semiminor axis and position angle parametrization (a,...
def format(config, name=None, writeSourceLine=True, prefix="", verbose=False)
A 2-dimensional celestial WCS that transform pixels to ICRS RA/Dec, using the LSST standard for pixel...
int getWidth() const
Return the number of columns in the image.
A Result struct for running an aperture flux algorithm with a single radius.
A FunctorKey for FluxResult.
std::string shiftKernel
"Warping kernel used to shift Sinc photometry coefficients to different center positions" ;
static FlagDefinition const FAILURE
static boost::shared_ptr< CoeffT const > get(afw::geom::ellipses::Axes const &outerEllipse, float const innerRadiusFactor=0.0)
Get the coefficients for an aperture.
vector-type utility class to build a collection of FlagDefinitions
meas::base::FluxErrElement instFluxErr
Standard deviation of instFlux in DN.
A class to manipulate images, masks, and variance as a single object.
double maxSincRadius
"Maximum radius (in pixels) for which the sinc algorithm should be used instead of the " "faster naiv...
Utility class for handling flag fields that indicate the failure modes of an algorithm.
Iterator begin() const
Iterator range over Spans whose pixels are within the Ellipse.
void copyResultToRecord(Result const &result, afw::table::SourceRecord &record, int index) const
lsst::geom::Point2D const & getCenter() const
Return the center point.
bool contains(Point2I const &point) const noexcept
Return true if the box contains the point.
A mapping between the keys of two Schemas, used to copy data between them.
static MagResultKey addFields(afw::table::Schema &schema, std::string const &name)
Add a pair of _mag, _magErr fields to a Schema, and return a MagResultKey that points to them.
Iterator class for CatalogT.
The photometric calibration of an exposure.
#define CONST_PTR(...)
A shared pointer to a const object.
_view_t::x_iterator x_iterator
An iterator for traversing the pixels in a row.
static std::string makeFieldPrefix(std::string const &name, double radius)
Construct an appropriate prefix for table fields.
A base class for image defects.
VariancePtr getVariance() const
Return a (shared_ptr to) the MaskedImage's variance.
An ellipse defined by an arbitrary BaseCore and a center point.
A pixelized region containing all pixels whose centers are within an Ellipse.
static Result computeSincFlux(afw::image::Image< T > const &image, afw::geom::ellipses::Ellipse const &ellipse, Control const &ctrl=Control())
Compute the instFlux (and optionally, uncertanties) within an aperture using Sinc photometry.
std::size_t size() const
return the current size (number of defined elements) of the collection
static FlagDefinitionList const & getFlagDefinitions()
Simple class used to define and document flags The name and doc constitute the identity of the FlagDe...
ApertureFluxAlgorithm(Control const &ctrl, std::string const &name, afw::table::Schema &schema, daf::base::PropertySet &metadata)
Construct the algorithm and add its fields to the given Schema.
std::vector< double > radii
"Radius (in pixels) of apertures." ;
void clip(Box2I const &other) noexcept
Shrink this to ensure that other.contains(*this).
std::shared_ptr< ImageT > offsetImage(ImageT const &image, float dx, float dy, std::string const &algorithmName="lanczos5", unsigned int buffer=0)
Return an image offset by (dx, dy) using the specified algorithm.
Class for storing generic metadata.
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.
An integer coordinate rectangle.
ImagePtr getImage() const
Return a (shared_ptr to) the MaskedImage's image.
Base class for multiple-aperture photometry algorithms.
BaseCore const & getCore() const
Return the ellipse core.
void set(Key< T > const &key, U const &value)
Set value of a field for the given key.
Configuration object for multiple-aperture flux algorithms.
A class to represent a 2-dimensional array of pixels.
iterator begin()
Iterator access.
static Result computeNaiveFlux(afw::image::Image< T > const &image, afw::geom::ellipses::Ellipse const &ellipse, Control const &ctrl=Control())
Compute the instFlux (and optionally, uncertanties) within an aperture using naive photometry.
int getWidth() const
Return the number of columns in the image.
static FlagDefinition const SINC_COEFFS_TRUNCATED