LSST Applications g04e9c324dd+8c5ae1fdc5,g134cb467dc+b203dec576,g18429d2f64+358861cd2c,g199a45376c+0ba108daf9,g1fd858c14a+dd066899e3,g262e1987ae+ebfced1d55,g29ae962dfc+72fd90588e,g2cef7863aa+aef1011c0b,g35bb328faa+8c5ae1fdc5,g3fd5ace14f+b668f15bc5,g4595892280+3897dae354,g47891489e3+abcf9c3559,g4d44eb3520+fb4ddce128,g53246c7159+8c5ae1fdc5,g67b6fd64d1+abcf9c3559,g67fd3c3899+1f72b5a9f7,g74acd417e5+cb6b47f07b,g786e29fd12+668abc6043,g87389fa792+8856018cbb,g89139ef638+abcf9c3559,g8d7436a09f+bcf525d20c,g8ea07a8fe4+9f5ccc88ac,g90f42f885a+6054cc57f1,g97be763408+06f794da49,g9dd6db0277+1f72b5a9f7,ga681d05dcb+7e36ad54cd,gabf8522325+735880ea63,gac2eed3f23+abcf9c3559,gb89ab40317+abcf9c3559,gbf99507273+8c5ae1fdc5,gd8ff7fe66e+1f72b5a9f7,gdab6d2f7ff+cb6b47f07b,gdc713202bf+1f72b5a9f7,gdfd2d52018+8225f2b331,ge365c994fd+375fc21c71,ge410e46f29+abcf9c3559,geaed405ab2+562b3308c0,gf9a733ac38+8c5ae1fdc5,w.2025.35
LSST Data Management Base Package
Loading...
Searching...
No Matches
KernelSolution.h
Go to the documentation of this file.
1// -*- lsst-c++ -*-
11
12#ifndef LSST_IP_DIFFIM_KERNELSOLUTION_H
13#define LSST_IP_DIFFIM_KERNELSOLUTION_H
14
15#include <memory>
16#include "Eigen/Core"
17
18#include "lsst/afw/math.h"
19#include "lsst/afw/image.h"
20#include "lsst/geom.h"
21#include "lsst/daf/base.h"
22
23namespace lsst {
24namespace ip {
25namespace diffim {
26
27 /*
28 * @brief Method used to solve for M and B
29 */
30
32 public:
36
44
49
50 explicit KernelSolution(Eigen::MatrixXd mMat,
51 Eigen::VectorXd bVec,
52 bool fitForBackground);
53 explicit KernelSolution(bool fitForBackground);
54 explicit KernelSolution();
55
56 virtual ~KernelSolution() {};
57 virtual void solve();
58 virtual void solve(Eigen::MatrixXd const& mMat,
59 Eigen::VectorXd const& bVec);
61 virtual double getConditionNumber(ConditionNumberType conditionType);
62 virtual double getConditionNumber(Eigen::MatrixXd const& mMat, ConditionNumberType conditionType);
63
64 inline Eigen::MatrixXd const& getM() {return _mMat;}
65 inline Eigen::VectorXd const& getB() {return _bVec;}
69 inline int getId() const { return _id; }
70
71 protected:
72 int _id;
73 Eigen::MatrixXd _mMat;
74 Eigen::VectorXd _bVec;
75 Eigen::VectorXd _aVec;
78 static int _SolutionId;
79
80 };
81
82 template <typename InputT>
84 public:
86
88 bool fitForBackground);
90
91 /* Overrides KernelSolution */
92 void solve();
93
94 /* Used by RegularizedKernelSolution */
95 virtual void build(lsst::afw::image::Image<InputT> const &templateImage,
96 lsst::afw::image::Image<InputT> const &scienceImage,
100 virtual double getBackground();
101 virtual double getKsum();
103
104 protected:
105 Eigen::MatrixXd _cMat;
106 Eigen::VectorXd _iVec;
107 Eigen::VectorXd _ivVec;
108
110 double _background;
111 double _kSum;
112
113 void _setKernel();
114 void _setKernelUncertainty();
115 };
116
117
118 template <typename InputT>
120 public:
122
124 bool fitForBackground);
126 virtual void buildOrig(lsst::afw::image::Image<InputT> const &templateImage,
127 lsst::afw::image::Image<InputT> const &scienceImage,
129 const &varianceEstimate,
131
132 virtual void buildWithMask(lsst::afw::image::Image<InputT> const &templateImage,
133 lsst::afw::image::Image<InputT> const &scienceImage,
135 const &varianceEstimate,
137
138 virtual void buildSingleMaskOrig(lsst::afw::image::Image<InputT> const &templateImage,
139 lsst::afw::image::Image<InputT> const &scienceImage,
141 const &varianceEstimate,
142 lsst::geom::Box2I maskBox);
143 };
144
145
146
147 template <typename InputT>
149 public:
151
153 bool fitForBackground,
154 Eigen::MatrixXd const& hMat,
156 );
158 void solve();
159 double getLambda() {return _lambda;}
160 double estimateRisk(double maxCond);
161
162 /* Include additive term (_lambda * _hMat) in M matrix? */
163 Eigen::MatrixXd getM(bool includeHmat = true);
164
165 private:
166 Eigen::MatrixXd const _hMat;
167 double _lambda;
169
170 std::vector<double> _createLambdaSteps();
171 };
172
173
175 public:
177
178 /* Creates a polynomial SpatialFunction */
183 );
184
186
187 void addConstraint(float xCenter, float yCenter,
188 Eigen::MatrixXd const& qMat,
189 Eigen::VectorXd const& wVec);
190
191 void solve();
195
196 private:
197 lsst::afw::math::Kernel::SpatialFunctionPtr _spatialKernelFunction;
198 bool _constantFirstTerm;
199
202 double _kSum;
203
205 int _nbases;
206 int _nkt;
207 int _nbt;
208 int _nt;
209
210 void _setKernel();
211 void _setKernelUncertainty();
212 };
213
214}}} // end of namespace lsst::ip::diffim
215
216#endif
A class to represent a 2-dimensional array of pixels.
Definition Image.h:51
Represent a 2-dimensional array of bitmask pixels.
Definition Mask.h:82
std::shared_ptr< lsst::afw::math::Function2< double > > SpatialFunctionPtr
Definition Kernel.h:113
Class for storing generic metadata.
Definition PropertySet.h:67
std::shared_ptr< PropertySet > Ptr
Definition PropertySet.h:70
An integer coordinate rectangle.
Definition Box.h:55
bool _fitForBackground
Background terms included in fit.
virtual double getConditionNumber(ConditionNumberType conditionType)
int _id
Unique ID for object.
Eigen::VectorXd _bVec
Derived least squares B vector.
std::shared_ptr< KernelSolution > Ptr
Eigen::MatrixXd const & getM()
Eigen::VectorXd _aVec
Derived least squares solution matrix.
KernelSolvedBy _solvedBy
Type of algorithm used to make solution.
lsst::afw::math::Kernel::Pixel PixelT
Eigen::VectorXd const & getB()
Eigen::MatrixXd _mMat
Derived least squares M matrix.
KernelSolution(Eigen::MatrixXd mMat, Eigen::VectorXd bVec, bool fitForBackground)
static int _SolutionId
Unique identifier for solution.
lsst::afw::image::Image< lsst::afw::math::Kernel::Pixel > ImageT
virtual void buildWithMask(lsst::afw::image::Image< InputT > const &templateImage, lsst::afw::image::Image< InputT > const &scienceImage, lsst::afw::image::Image< lsst::afw::image::VariancePixel > const &varianceEstimate, lsst::afw::image::Mask< lsst::afw::image::MaskPixel > const &pixelMask)
virtual void buildSingleMaskOrig(lsst::afw::image::Image< InputT > const &templateImage, lsst::afw::image::Image< InputT > const &scienceImage, lsst::afw::image::Image< lsst::afw::image::VariancePixel > const &varianceEstimate, lsst::geom::Box2I maskBox)
MaskedKernelSolution(lsst::afw::math::KernelList const &basisList, bool fitForBackground)
std::shared_ptr< MaskedKernelSolution< InputT > > Ptr
virtual void buildOrig(lsst::afw::image::Image< InputT > const &templateImage, lsst::afw::image::Image< InputT > const &scienceImage, lsst::afw::image::Image< lsst::afw::image::VariancePixel > const &varianceEstimate, lsst::afw::image::Mask< lsst::afw::image::MaskPixel > pixelMask)
std::shared_ptr< RegularizedKernelSolution< InputT > > Ptr
RegularizedKernelSolution(lsst::afw::math::KernelList const &basisList, bool fitForBackground, Eigen::MatrixXd const &hMat, lsst::daf::base::PropertySet const &ps)
std::pair< std::shared_ptr< lsst::afw::math::LinearCombinationKernel >, lsst::afw::math::Kernel::SpatialFunctionPtr > getSolutionPair()
std::shared_ptr< lsst::afw::image::Image< lsst::afw::math::Kernel::Pixel > > makeKernelImage(lsst::geom::Point2D const &pos)
SpatialKernelSolution(lsst::afw::math::KernelList const &basisList, lsst::afw::math::Kernel::SpatialFunctionPtr spatialKernelFunction, lsst::afw::math::Kernel::SpatialFunctionPtr background, lsst::daf::base::PropertySet const &ps)
void addConstraint(float xCenter, float yCenter, Eigen::MatrixXd const &qMat, Eigen::VectorXd const &wVec)
std::shared_ptr< SpatialKernelSolution > Ptr
double _background
Derived differential background estimate.
virtual std::pair< std::shared_ptr< lsst::afw::math::Kernel >, double > getSolutionPair()
void _setKernel()
Set kernel after solution.
virtual std::shared_ptr< lsst::afw::math::Kernel > getKernel()
Eigen::VectorXd _iVec
Vectorized I.
std::shared_ptr< lsst::afw::math::Kernel > _kernel
Derived single-object convolution kernel.
void _setKernelUncertainty()
Not implemented.
virtual std::shared_ptr< lsst::afw::image::Image< lsst::afw::math::Kernel::Pixel > > makeKernelImage()
Eigen::VectorXd _ivVec
Inverse variance.
std::shared_ptr< StaticKernelSolution< InputT > > Ptr
StaticKernelSolution(lsst::afw::math::KernelList const &basisList, bool fitForBackground)
virtual void build(lsst::afw::image::Image< InputT > const &templateImage, lsst::afw::image::Image< InputT > const &scienceImage, lsst::afw::image::Image< lsst::afw::image::VariancePixel > const &varianceEstimate)
T endl(T... args)
std::vector< std::shared_ptr< Kernel > > KernelList
Definition Kernel.h:462
Point< double, 2 > Point2D
Definition Point.h:324