LSSTApplications  17.0+11,17.0+34,17.0+56,17.0+57,17.0+59,17.0+7,17.0-1-g377950a+33,17.0.1-1-g114240f+2,17.0.1-1-g4d4fbc4+28,17.0.1-1-g55520dc+49,17.0.1-1-g5f4ed7e+52,17.0.1-1-g6dd7d69+17,17.0.1-1-g8de6c91+11,17.0.1-1-gb9095d2+7,17.0.1-1-ge9fec5e+5,17.0.1-1-gf4e0155+55,17.0.1-1-gfc65f5f+50,17.0.1-1-gfc6fb1f+20,17.0.1-10-g87f9f3f+1,17.0.1-11-ge9de802+16,17.0.1-16-ga14f7d5c+4,17.0.1-17-gc79d625+1,17.0.1-17-gdae4c4a+8,17.0.1-2-g26618f5+29,17.0.1-2-g54f2ebc+9,17.0.1-2-gf403422+1,17.0.1-20-g2ca2f74+6,17.0.1-23-gf3eadeb7+1,17.0.1-3-g7e86b59+39,17.0.1-3-gb5ca14a,17.0.1-3-gd08d533+40,17.0.1-30-g596af8797,17.0.1-4-g59d126d+4,17.0.1-4-gc69c472+5,17.0.1-6-g5afd9b9+4,17.0.1-7-g35889ee+1,17.0.1-7-gc7c8782+18,17.0.1-9-gc4bbfb2+3,w.2019.22
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...
Basic LSST definitions.
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.
#define PTR(...)
Definition: base.h:41
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