LSST Applications 24.1.6,g063fba187b+56b85ce14a,g0f08755f38+df8a265115,g12f32b3c4e+891a09f10d,g1524ad2192+7a5d7b3fbd,g1653933729+a8ce1bb630,g168dd56ebc+a8ce1bb630,g28da252d5a+07cb1400be,g2bbee38e9b+ae03bbfc84,g2bc492864f+ae03bbfc84,g3156d2b45e+6e55a43351,g347aa1857d+ae03bbfc84,g35bb328faa+a8ce1bb630,g3a166c0a6a+ae03bbfc84,g3e281a1b8c+c5dd892a6c,g414038480c+6b9177ef31,g41af890bb2+8f257c4c0b,g781aacb6e4+a8ce1bb630,g7af13505b9+7137b3b17d,g80478fca09+6df6903293,g82479be7b0+091ce1d07f,g858d7b2824+df8a265115,g89c8672015+f4add4ffd5,g9125e01d80+a8ce1bb630,g9726552aa6+414189b318,ga5288a1d22+4a2bca08d7,gacef1a1666+c9a8ff65f4,gb58c049af0+d64f4d3760,gbcfae0f0a0+de1d42d831,gc28159a63d+ae03bbfc84,gcf0d15dbbd+72117bf34e,gda6a2b7d83+72117bf34e,gdaeeff99f8+1711a396fd,ge500cccec5+c8c9c9af63,ge79ae78c31+ae03bbfc84,gf0baf85859+c1f95f4921,gfa517265be+df8a265115,gfa999e8aa5+17cd334064,gfb92a5be7c+df8a265115
LSST Data Management Base Package
Loading...
Searching...
No Matches
PsfCandidate.h
Go to the documentation of this file.
1// -*- LSST-C++ -*-
2#if !defined(LSST_MEAS_ALGORITHMS_PSFCANDIDATE_H)
3#define LSST_MEAS_ALGORITHMS_PSFCANDIDATE_H
4
5/*
6 * LSST Data Management System
7 * Copyright 2008, 2009, 2010 LSST Corporation.
8 *
9 * This product includes software developed by the
10 * LSST Project (http://www.lsst.org/).
11 *
12 * This program is free software: you can redistribute it and/or modify
13 * it under the terms of the GNU General Public License as published by
14 * the Free Software Foundation, either version 3 of the License, or
15 * (at your option) any later version.
16 *
17 * This program is distributed in the hope that it will be useful,
18 * but WITHOUT ANY WARRANTY; without even the implied warranty of
19 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
20 * GNU General Public License for more details.
21 *
22 * You should have received a copy of the LSST License Statement and
23 * the GNU General Public License along with this program. If not,
24 * see <http://www.lsstcorp.org/LegalNotices/>.
25 */
26
34#include <memory>
35#include <vector>
36
37#include "lsst/geom/Point.h"
43
44namespace lsst {
45namespace meas {
46namespace algorithms {
54template <typename PixelT>
56public:
60
62
70 parentExposure
71 )
72 : afw::math::SpatialCellImageCandidate(source->getX(), source->getY()),
73 _parentExposure(parentExposure),
74 _offsetImage(),
75 _source(source),
76 _image(nullptr),
77 _amplitude(0.0),
78 _var(1.0) {}
79
85 parentExposure,
86 double xCenter,
87 double yCenter
88 )
89 : afw::math::SpatialCellImageCandidate(xCenter, yCenter),
90 _parentExposure(parentExposure),
91 _offsetImage(),
92 _source(source),
93 _image(nullptr),
94 _amplitude(0.0),
95 _var(1.0) {}
96
98 virtual ~PsfCandidate(){};
99
105 double getCandidateRating() const { return _source->getPsfInstFlux(); }
106
109
111 double getAmplitude() const { return _amplitude; }
112
114 void setAmplitude(double amplitude) { _amplitude = amplitude; }
115
117 double getVar() const { return _var; }
118
120 void setVar(double var) { _var = var; }
121
123 std::shared_ptr<afw::image::MaskedImage<PixelT> const> getMaskedImage(int width, int height) const;
125 getOffsetImage(std::string const algorithm, unsigned int buffer) const;
126
128 static int getBorderWidth();
129
131 static void setBorderWidth(int border);
132
136 static void setPixelThreshold(float threshold);
137
139 static float getPixelThreshold();
140
142 static void setMaskBlends(bool doMaskBlends);
143
145 static bool getMaskBlends();
146
147private:
149 _parentExposure; // the %image that the Sources are found in
150
152 offsetImage(std::shared_ptr<afw::image::MaskedImage<PixelT>> img, std::string const algorithm, unsigned int buffer);
153
155 extractImage(unsigned int width, unsigned int height) const;
156
157 std::shared_ptr<afw::image::MaskedImage<PixelT>> mutable _offsetImage; // %image offset to put center on a pixel
158 std::shared_ptr<afw::table::SourceRecord> _source; // the Source itself
159
160 mutable std::shared_ptr<afw::image::MaskedImage<PixelT>> _image; // cutout image to return (cached)
161 double _amplitude; // best-fit amplitude of current PSF model
162 double _var; // variance to use when fitting this candidate
163 static int _border; // width of border of ignored pixels around _image
164 geom::Point2D _xyCenter;
165 static int _defaultWidth;
166 static float _pixelThreshold;
167 static bool _doMaskBlends;
168};
169
175template <typename PixelT>
177 const & source,
179 image
180 ) {
181 return std::make_shared<PsfCandidate<PixelT>>(source, image);
182}
183
184} // namespace algorithms
185} // namespace meas
186} // namespace lsst
187
188#endif
afw::table::Key< afw::table::Array< ImagePixelT > > image
A class to contain the data, WCS, and other information needed to describe an image of the sky.
Definition Exposure.h:72
A class to manipulate images, masks, and variance as a single object.
Definition MaskedImage.h:74
Base class for candidate objects in a SpatialCell that are able to return an Image of some sort (e....
Class stored in SpatialCells for spatial Psf fitting.
static float getPixelThreshold()
Get threshold for rejecting pixels unconnected with the central footprint.
std::shared_ptr< const PsfCandidate< PixelT > > ConstPtr
double getCandidateRating() const
Return Cell rating.
static bool getMaskBlends()
Get whether blends are masked.
PsfCandidate(std::shared_ptr< afw::table::SourceRecord > const &source, std::shared_ptr< afw::image::Exposure< PixelT > const > parentExposure)
Construct a PsfCandidate from a specified source and image.
static void setBorderWidth(int border)
Set the number of pixels to ignore around the candidate image's edge.
double getVar() const
Return the variance in use when fitting this object.
static int getBorderWidth()
Return the number of pixels being ignored around the candidate image's edge.
static void setPixelThreshold(float threshold)
Set threshold for rejecting pixels unconnected with the central footprint.
void setVar(double var)
Set the variance to use when fitting this object.
std::shared_ptr< afw::image::MaskedImage< PixelT > const > getMaskedImage() const
Return the image at the position of the Source, without any sub-pixel shifts to put the centre of the...
void setAmplitude(double amplitude)
Set the best-fit amplitude.
afw::image::MaskedImage< PixelT > MaskedImageT
PsfCandidate(std::shared_ptr< afw::table::SourceRecord > const &source, std::shared_ptr< afw::image::Exposure< PixelT > const > parentExposure, double xCenter, double yCenter)
Construct a PsfCandidate from a specified source, image and xyCenter.
double getAmplitude() const
Return the best-fit amplitude.
static void setMaskBlends(bool doMaskBlends)
Set whether blends are masked.
std::shared_ptr< PsfCandidate< PixelT > > Ptr
std::shared_ptr< afw::table::SourceRecord > getSource() const
Return the original Source.
std::shared_ptr< afw::image::MaskedImage< PixelT > > getOffsetImage(std::string const algorithm, unsigned int buffer) const
Return an offset version of the image of the source.
std::shared_ptr< PsfCandidate< PixelT > > makePsfCandidate(std::shared_ptr< afw::table::SourceRecord > const &source, std::shared_ptr< afw::image::Exposure< PixelT > > image)
Return a PsfCandidate of the right sort.