LSSTApplications  17.0+124,17.0+14,17.0+73,18.0.0+37,18.0.0+80,18.0.0-4-g68ffd23+4,18.1.0-1-g0001055+12,18.1.0-1-g03d53ef+5,18.1.0-1-g1349e88+55,18.1.0-1-g2505f39+44,18.1.0-1-g5315e5e+4,18.1.0-1-g5e4b7ea+14,18.1.0-1-g7e8fceb+4,18.1.0-1-g85f8cd4+48,18.1.0-1-g8ff0b9f+4,18.1.0-1-ga2c679d+1,18.1.0-1-gd55f500+35,18.1.0-10-gb58edde+2,18.1.0-11-g0997b02+4,18.1.0-13-gfe4edf0b+12,18.1.0-14-g259bd21+21,18.1.0-19-gdb69f3f+2,18.1.0-2-g5f9922c+24,18.1.0-2-gd3b74e5+11,18.1.0-2-gfbf3545+32,18.1.0-26-g728bddb4+5,18.1.0-27-g6ff7ca9+2,18.1.0-3-g52aa583+25,18.1.0-3-g8ea57af+9,18.1.0-3-gb69f684+42,18.1.0-3-gfcaddf3+6,18.1.0-32-gd8786685a,18.1.0-4-gf3f9b77+6,18.1.0-5-g1dd662b+2,18.1.0-5-g6dbcb01+41,18.1.0-6-gae77429+3,18.1.0-7-g9d75d83+9,18.1.0-7-gae09a6d+30,18.1.0-9-gc381ef5+4,w.2019.45
LSSTDataManagementBasePackage
Model.h
Go to the documentation of this file.
1 // -*- lsst-c++ -*-
2 /*
3  * LSST Data Management System
4  * Copyright 2008-2013 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 #ifndef LSST_MEAS_MODELFIT_Model_h_INCLUDED
25 #define LSST_MEAS_MODELFIT_Model_h_INCLUDED
26 
27 #include <vector>
28 
29 #include "lsst/base.h"
33 
34 namespace lsst { namespace meas { namespace modelfit {
35 
36 class Model;
37 class Prior;
38 
40 
42 
56 class Model {
57 public:
58 
59  enum CenterEnum {
60  FIXED_CENTER = 0x0,
63  };
64 
70 
86  static PTR(Model) make(BasisVector basisVector, NameVector const & prefixes, CenterEnum center);
87 
108 
121  static PTR(Model) makeGaussian(CenterEnum center, double radius=1.0);
122 
124  int getNonlinearDim() const { return _nonlinearNames.size(); }
125 
127  int getAmplitudeDim() const { return _amplitudeNames.size(); }
128 
130  int getFixedDim() const { return _fixedNames.size(); }
131 
133  int getBasisCount() const { return _basisVector.size(); }
134 
136  NameVector const & getNonlinearNames() const { return _nonlinearNames; }
137 
139  NameVector const & getAmplitudeNames() const { return _amplitudeNames; }
140 
142  NameVector const & getFixedNames() const { return _fixedNames; }
143 
145  BasisVector const & getBasisVector() const { return _basisVector; }
146 
149  ndarray::Array<Scalar const,1,1> const & nonlinear,
150  ndarray::Array<Scalar const,1,1> const & amplitudes,
151  ndarray::Array<Scalar const,1,1> const & fixed
152  ) const;
153 
155  virtual PTR(Prior) adaptPrior(PTR(Prior) prior) const = 0;
156 
161  virtual EllipseVector makeEllipseVector() const = 0;
162 
179  virtual void writeEllipses(
180  Scalar const * nonlinearIter, Scalar const * fixedIter,
181  EllipseIterator ellipseIter
182  ) const = 0;
183 
200  virtual void readEllipses(
201  EllipseConstIterator ellipseIter,
202  Scalar * nonlinearIter, Scalar * fixedIter
203  ) const = 0;
204 
215  EllipseVector writeEllipses(
216  ndarray::Array<Scalar const,1,1> const & nonlinear,
217  ndarray::Array<Scalar const,1,1> const & fixed
218  ) const;
219 
230  void readEllipses(
231  EllipseVector const & ellipses,
232  ndarray::Array<Scalar,1,1> const & nonlinear,
233  ndarray::Array<Scalar,1,1> const & fixed
234  ) const;
235 
244  virtual void transformParameters(
246  ndarray::Array<Scalar,1,1> const & nonlinear,
247  ndarray::Array<Scalar,1,1> const & amplitudes,
248  ndarray::Array<Scalar,1,1> const & fixed
249  ) const;
250 
251  virtual ~Model() {}
252 
253  // No copying
254  Model (const Model&) = delete;
255  Model& operator=(const Model&) = delete;
256 
257  // No moving
258  Model (Model&&) = delete;
259  Model& operator=(Model&&) = delete;
260 
261 protected:
262 
263  Model(
264  BasisVector basisVector,
265  NameVector nonlinearNames,
266  NameVector amplitudeNames,
267  NameVector fixedNames
268  );
269 
270 private:
271  NameVector _nonlinearNames;
272  NameVector _amplitudeNames;
273  NameVector _fixedNames;
274  BasisVector _basisVector;
275 };
276 
277 }}} // namespace lsst::meas::modelfit
278 
279 #endif // !LSST_MEAS_MODELFIT_Model_h_INCLUDED
NameVector const & getFixedNames() const
Return the names of the fixed nonlinear parameters.
Definition: Model.h:142
virtual EllipseVector makeEllipseVector() const =0
Return an uninitialized vector of afw::geom::ellipses::Ellipse with the parametrization expected by r...
double Scalar
Typedefs to be used for probability and parameter values.
Definition: common.h:44
std::vector< boost::shared_ptr< Model > > ModelVector
Definition: Model.h:39
virtual boost::shared_ptr< Prior > adaptPrior(boost::shared_ptr< Prior > prior) const =0
Given an arbitrary prior, return one compatible with this Model or throw LogicError.
table::Key< table::Array< double > > amplitudes
A basis formed from a linear combination of shapelet bases that differ only in radius.
int getNonlinearDim() const
Return the number of free nonlinear parameters.
Definition: Model.h:124
Abstract base class and concrete factories that define multi-shapelet galaxy models.
Definition: Model.h:56
int getBasisCount() const
Return the number of MultiShapeletBasis objects (equivalently, the number of ellipses) ...
Definition: Model.h:133
A base class for image defects.
A multi-scale shapelet function.
static boost::shared_ptr< Model > make(BasisVector basisVector, NameVector const &prefixes, CenterEnum center)
Construct a concrete Model instance with multiple ellipses and multishapelet bases.
std::vector< std::string > NameVector
Definition: Model.h:65
std::vector< boost::shared_ptr< shapelet::MultiShapeletBasis > > BasisVector
Definition: Model.h:66
table::Key< table::Array< double > > basis
Definition: PsfexPsf.cc:361
shapelet::MultiShapeletFunction makeShapeletFunction(ndarray::Array< Scalar const, 1, 1 > const &nonlinear, ndarray::Array< Scalar const, 1, 1 > const &amplitudes, ndarray::Array< Scalar const, 1, 1 > const &fixed) const
Create a MultiShapeletFunction object from a set of parameter vectors.
Base class for Bayesian priors.
Definition: Prior.h:36
int getFixedDim() const
Return the number of fixed nonlinear parameters.
Definition: Model.h:130
A local mapping between two UnitSystems.
Definition: UnitSystem.h:80
T size(T... args)
STL class.
virtual void transformParameters(LocalUnitTransform const &transform, ndarray::Array< Scalar, 1, 1 > const &nonlinear, ndarray::Array< Scalar, 1, 1 > const &amplitudes, ndarray::Array< Scalar, 1, 1 > const &fixed) const
Transform (in-place) parameter vectors from one unit system to another.
int getAmplitudeDim() const
Return the number of linear parameters.
Definition: Model.h:127
virtual void readEllipses(EllipseConstIterator ellipseIter, Scalar *nonlinearIter, Scalar *fixedIter) const =0
Convert a vector of ellipses to a set of nonlinear+fixed parameter arrays.
Model & operator=(const Model &)=delete
static boost::shared_ptr< Model > makeGaussian(CenterEnum center, double radius=1.0)
Construct a concrete Model instance that represents a single elliptical Gaussian function.
Basic LSST definitions.
std::vector< afw::geom::ellipses::Ellipse >::iterator EllipseIterator
Definition: Model.h:68
virtual void writeEllipses(Scalar const *nonlinearIter, Scalar const *fixedIter, EllipseIterator ellipseIter) const =0
Convert a set of nonlinear+fixed parameter arrays to a vector of ellipses.
Model(const Model &)=delete
NameVector const & getNonlinearNames() const
Return the names of the free nonlinear parameters.
Definition: Model.h:136
BasisVector const & getBasisVector() const
Return the MultiShapeletBasis objects that comprise the Model.
Definition: Model.h:145
std::vector< afw::geom::ellipses::Ellipse > EllipseVector
Definition: Model.h:67
NameVector const & getAmplitudeNames() const
Return the names of the amplitude parameters.
Definition: Model.h:139
std::vector< afw::geom::ellipses::Ellipse >::const_iterator EllipseConstIterator
Definition: Model.h:69
#define PTR(...)
Definition: base.h:41