LSST Applications g042eb84c57+730a74494b,g04e9c324dd+8c5ae1fdc5,g134cb467dc+1f1e3e7524,g199a45376c+0ba108daf9,g1fd858c14a+fa7d31856b,g210f2d0738+f66ac109ec,g262e1987ae+83a3acc0e5,g29ae962dfc+d856a2cb1f,g2cef7863aa+aef1011c0b,g35bb328faa+8c5ae1fdc5,g3fd5ace14f+a1e0c9f713,g47891489e3+0d594cb711,g4d44eb3520+c57ec8f3ed,g4d7b6aa1c5+f66ac109ec,g53246c7159+8c5ae1fdc5,g56a1a4eaf3+fd7ad03fde,g64539dfbff+f66ac109ec,g67b6fd64d1+0d594cb711,g67fd3c3899+f66ac109ec,g6985122a63+0d594cb711,g74acd417e5+3098891321,g786e29fd12+668abc6043,g81db2e9a8d+98e2ab9f28,g87389fa792+8856018cbb,g89139ef638+0d594cb711,g8d7436a09f+80fda9ce03,g8ea07a8fe4+760ca7c3fc,g90f42f885a+033b1d468d,g97be763408+a8a29bda4b,g99822b682c+e3ec3c61f9,g9d5c6a246b+0d5dac0c3d,ga41d0fce20+9243b26dd2,gbf99507273+8c5ae1fdc5,gd7ef33dd92+0d594cb711,gdab6d2f7ff+3098891321,ge410e46f29+0d594cb711,geaed405ab2+c4bbc419c6,gf9a733ac38+8c5ae1fdc5,w.2025.38
LSST Data Management Base Package
Loading...
Searching...
No Matches
KernelCandidate.h
Go to the documentation of this file.
1// -*- lsst-c++ -*-
11
12#ifndef LSST_IP_DIFFIM_KERNELCANDIDATE_H
13#define LSST_IP_DIFFIM_KERNELCANDIDATE_H
14
15#include <memory>
16#include "Eigen/Core"
17
18#include "lsst/afw/math.h"
19#include "lsst/afw/image.h"
22#include "lsst/daf/base.h"
23
24namespace lsst {
25namespace ip {
26namespace diffim {
27
28
38 template <typename _PixelT>
40 public:
42 typedef _PixelT PixelT; // _after_ afw::math::Kernel::Pixel
43 public:
48
50 ORIG = 0,
51 PCA = 1,
53 };
54
64 KernelCandidate(float const xCenter,
65 float const yCenter,
66 MaskedImagePtr const& templateMaskedImage,
67 MaskedImagePtr const& scienceMaskedImage,
68 daf::base::PropertySet const& ps);
69
79 KernelCandidate(SourcePtr const& source,
80 MaskedImagePtr const& templateMaskedImage,
81 MaskedImagePtr const& scienceMaskedImage,
82 daf::base::PropertySet const& ps);
84 virtual ~KernelCandidate() {};
85
91 double getCandidateRating() const { return _coreFlux; }
95 SourcePtr getSource() const { return _source; }
99 MaskedImagePtr getTemplateMaskedImage() {return _templateMaskedImage;}
100 MaskedImagePtr getScienceMaskedImage() {return _scienceMaskedImage;}
101
107 double getBackground(CandidateSwitch cand) const;
108 double getKsum(CandidateSwitch cand) const;
110 std::shared_ptr<ImageT const> getImage() const; // For SpatialCellImageCandidate
112
117
125 double background
126 );
127
128 bool isInitialized() const {return _isInitialized;}
129
130
150
164
165 /*
166 * @note This method uses an estimate of the variance which is the
167 * straight difference of the 2 images. If requested in the PropertySet
168 * ("iterateSingleKernel"), the kernel will be rebuilt using the
169 * variance of the difference image resulting from this first
170 * approximate step. This is particularly useful when convolving a
171 * single-depth science image; the variance (and thus resulting kernel)
172 * generally converges after 1 iteration. If
173 * "constantVarianceWeighting" is requested in the PropertySet, no iterations
174 * will be performed even if requested.
175 */
176
177 void build(
178 afw::math::KernelList const& basisList
179 );
180 void build(
181 afw::math::KernelList const& basisList,
182 Eigen::MatrixXd const& hMat
183 );
184
185 private:
186 MaskedImagePtr _templateMaskedImage;
187 MaskedImagePtr _scienceMaskedImage;
188 VariancePtr _varianceEstimate;
190 SourcePtr _source;
191 double _coreFlux;
192 bool _isInitialized;
193 bool _useRegularization;
194 bool _fitForBackground;
195
196 /* best single raw kernel */
198
199 /* with Pca basis */
201
202 void _buildKernelSolution(afw::math::KernelList const& basisList,
203 Eigen::MatrixXd const& hMat);
204 };
205
206
218 template <typename PixelT>
220 makeKernelCandidate(float const xCenter,
221 float const yCenter,
222 std::shared_ptr<afw::image::MaskedImage<PixelT> > const& templateMaskedImage,
223 std::shared_ptr<afw::image::MaskedImage<PixelT> > const& scienceMaskedImage,
224 daf::base::PropertySet const& ps){
225
227 templateMaskedImage,
228 scienceMaskedImage,
229 ps));
230 }
231
243 template <typename PixelT>
246 std::shared_ptr<afw::image::MaskedImage<PixelT> > const& templateMaskedImage,
247 std::shared_ptr<afw::image::MaskedImage<PixelT> > const& scienceMaskedImage,
248 daf::base::PropertySet const& ps){
249
251 templateMaskedImage,
252 scienceMaskedImage,
253 ps));
254 }
255
256
257}}} // end of namespace lsst::ip::diffim
258
259#endif
Declaration of classes to store the solution for convolution kernels.
A class to represent a 2-dimensional array of pixels.
Definition Image.h:51
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 for storing generic metadata.
Definition PropertySet.h:67
std::shared_ptr< PropertySet > Ptr
Definition PropertySet.h:70
Class stored in SpatialCells for spatial Kernel fitting.
std::shared_ptr< ImageT const > getImage() const
std::shared_ptr< afw::image::MaskedImage< PixelT > > MaskedImagePtr
std::shared_ptr< afw::image::Image< afw::image::VariancePixel > > VariancePtr
virtual ~KernelCandidate()
Destructor.
afw::image::Image< afw::math::Kernel::Pixel > ImageT
afw::image::MaskedImage< PixelT > getDifferenceImage(CandidateSwitch cand)
Calculate associated difference image using internal solutions.
double getCandidateRating() const
Return Candidate rating.
std::shared_ptr< KernelCandidate > Ptr
MaskedImagePtr getTemplateMaskedImage()
Return pointers to the image pixels used in kernel determination.
double getBackground(CandidateSwitch cand) const
double getKsum(CandidateSwitch cand) const
std::shared_ptr< afw::table::SourceRecord > SourcePtr
std::shared_ptr< StaticKernelSolution< PixelT > > getKernelSolution(CandidateSwitch cand) const
std::shared_ptr< afw::math::Kernel > getKernel(CandidateSwitch cand) const
Return results of kernel solution.
KernelCandidate(float const xCenter, float const yCenter, MaskedImagePtr const &templateMaskedImage, MaskedImagePtr const &scienceMaskedImage, daf::base::PropertySet const &ps)
Constructor.
SourcePtr getSource() const
Return the original source.
std::shared_ptr< ImageT > getKernelImage(CandidateSwitch cand) const
void build(afw::math::KernelList const &basisList)
Core functionality of KernelCandidate, to build and fill a KernelSolution.
std::vector< std::shared_ptr< Kernel > > KernelList
Definition Kernel.h:462
std::shared_ptr< KernelCandidate< PixelT > > makeKernelCandidate(float const xCenter, float const yCenter, std::shared_ptr< afw::image::MaskedImage< PixelT > > const &templateMaskedImage, std::shared_ptr< afw::image::MaskedImage< PixelT > > const &scienceMaskedImage, daf::base::PropertySet const &ps)
Return a KernelCandidate pointer of the right sort.