LSSTApplications  10.0-2-g4f67435,11.0.rc2+1,11.0.rc2+12,11.0.rc2+3,11.0.rc2+4,11.0.rc2+5,11.0.rc2+6,11.0.rc2+7,11.0.rc2+8
LSSTDataManagementBasePackage
CircularApertureFlux.cc
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 #include "lsst/afw/table/Source.h"
28 
29 namespace lsst { namespace meas { namespace base {
30 
32  Control const & ctrl,
33  std::string const & name,
35  daf::base::PropertySet & metadata
36 ) : ApertureFluxAlgorithm(ctrl, name, schema, metadata)
37 {
38  for (std::size_t i = 0; i < ctrl.radii.size(); ++i) {
39  if (ctrl.radii[i] > ctrl.maxSincRadius) break;
40  SincCoeffs<float>::cache(0.0, ctrl.radii[i]);
41  }
42 }
43 
45  afw::table::SourceRecord & measRecord,
46  afw::image::Exposure<float> const & exposure
47 ) const {
50  = boost::static_pointer_cast<afw::geom::ellipses::Axes>(ellipse.getCorePtr());
51  for (std::size_t i = 0; i < _ctrl.radii.size(); ++i) {
52  // Each call to _centroidExtractor within this loop goes through exactly the same error-checking
53  // logic and returns the same result, but it's not expensive logic, so we just call it repeatedly
54  // instead of coming up with a new interface that would allow us to move it outside the loop.
55  ellipse.setCenter(_centroidExtractor(measRecord, getFlagHandler(i)));
56  axes->setA(_ctrl.radii[i]);
57  axes->setB(_ctrl.radii[i]);
58  ApertureFluxAlgorithm::Result result = computeFlux(exposure.getMaskedImage(), ellipse, _ctrl);
59  copyResultToRecord(result, measRecord, i);
60  }
61 }
62 
63 }}} // namespace lsst::meas::base
Defines the fields and offsets for a table.
Definition: Schema.h:46
MaskedImageT getMaskedImage()
Return the MaskedImage.
Definition: Exposure.h:150
#define PTR(...)
Definition: base.h:41
table::Key< std::string > name
Definition: ApCorrMap.cc:71
void setCenter(Point2D const &center)
Set the center point.
Definition: Ellipse.h:72
void copyResultToRecord(Result const &result, afw::table::SourceRecord &record, int index) const
static Result computeFlux(afw::image::Image< T > const &image, afw::geom::ellipses::Ellipse const &ellipse, Control const &ctrl=Control())
An ellipse defined by an arbitrary BaseCore and a center point.
Definition: Ellipse.h:50
double maxSincRadius
&quot;Maximum radius (in pixels) for which the sinc algorithm should be used instead of the &quot; &quot;faster naiv...
Definition: ApertureFlux.h:56
FlagHandler const & getFlagHandler(int index) const
Definition: ApertureFlux.h:225
tbl::Schema schema
static void cache(float rInner, float rOuter)
Definition: SincCoeffs.cc:515
An ellipse core for the semimajor/semiminor axis and position angle parametrization (a...
Definition: Axes.h:45
Class for storing generic metadata.
Definition: PropertySet.h:82
Record class that contains measurements made on a single exposure.
Definition: Source.h:81
virtual void measure(afw::table::SourceRecord &record, afw::image::Exposure< float > const &exposure) const
CircularApertureFluxAlgorithm(Control const &ctrl, std::string const &name, afw::table::Schema &schema, daf::base::PropertySet &metadata)
std::vector< double > radii
&quot;Radius (in pixels) of apertures.&quot; ;
Definition: ApertureFlux.h:50
BaseCore::ConstPtr getCorePtr() const
Return the ellipse core.
Definition: Ellipse.h:81
SafeCentroidExtractor _centroidExtractor
Definition: ApertureFlux.h:228