lsst::jointcal::PhotometryTransfoChebyshev Class Reference

nth-order 2d Chebyshev photometry transfo. More...

#include <PhotometryTransfo.h>

Inheritance diagram for lsst::jointcal::PhotometryTransfoChebyshev:

Public Member Functions
	PhotometryTransfoChebyshev (size_t order, afw::geom::Box2D const &bbox, bool identity)
	Create a Chebyshev transfo with terms up to order in (x*y). More...
	PhotometryTransfoChebyshev (ndarray::Array< double, 2, 2 > const &coefficients, afw::geom::Box2D const &bbox)
	Create a Chebyshev transfo with the specified coefficients. More...
double	transformError (double x, double y, double value, double valueErr) const override
	Return the transformed valueErr at Point(x,y). More...
void	dump (std::ostream &stream=std::cout) const override
	dumps the transfo coefficients to stream. More...
int	getNpar () const override
	Return the number of parameters (used to compute chisq) More...
void	offsetParams (Eigen::VectorXd const &delta) override
	Offset the parameters by some (negative) amount during fitting. More...
ndarray::Array< double, 2, 2 >	getCoefficients () const
	Get a copy of the coefficients of the polynomials, as a 2d array (NOTE: layout is [y][x]) More...
Eigen::VectorXd	getParameters () const override
	Get a copy of the parameters of this model, in the same order as offsetParams. More...
ndarray::Size	getOrder () const
afw::geom::Box2D	getBBox () const
double	mean () const
virtual double	transform (double x, double y, double value) const =0
	Return the transform of value at (x,y). More...
double	transform (Point const &in, double value) const
	Return the transformed value at Point(x,y). More...
double	transformError (Point const &in, double value, double valueErr) const
	Return the transformed valueErr at Point(x,y). More...
virtual std::shared_ptr< PhotometryTransfo >	clone () const =0
	return a copy (allocated by new) of the transformation. More...
virtual void	computeParameterDerivatives (double x, double y, double value, Eigen::Ref< Eigen::VectorXd > derivatives) const =0
	Compute the derivatives with respect to the parameters (i.e. More...

Protected Member Functions
double	computeChebyshev (double x, double y) const
	Return the value of this polynomial at x,y. More...
void	computeChebyshevDerivatives (double x, double y, Eigen::Ref< Eigen::VectorXd > derivatives) const
	Set the derivatives of this polynomial at x,y. More...

Detailed Description

nth-order 2d Chebyshev photometry transfo.

The 2-d Chebyshev polynomial used here is defined as:

\[ f(x,y) = \sum_i \sum_j a_{i,j} T_i(x) T_j(y) \]

where \(T_n(x)\) is the n-th order Chebyshev polynomial of \(x\) and \(a_{i,j}\) is the corresponding coefficient of the (i,j) polynomial term.

Note that the polynomial order=n means that the highest terms will be of the form:

\[ a_{0,n}*x^n*y^0, a_{n-1,1}*x^(n-1)*y^1, ..., a_{1,n-1}*x^1*y^(n-1), a_{n,0}*x^0*y^n \]

Definition at line 218 of file PhotometryTransfo.h.

Constructor & Destructor Documentation

PhotometryTransfoChebyshev() [1/2]

lsst::jointcal::PhotometryTransfoChebyshev::PhotometryTransfoChebyshev	(	size_t	order,
		afw::geom::Box2D const &	bbox,
		bool	identity
	)

Create a Chebyshev transfo with terms up to order in (x*y).

Parameters

[in]	order	The maximum order in (x*y).
[in]	bbox	The bounding box it is valid within, to rescale it to [-1,1].
[in]	identity	If true, set a_0,0==1, otherwise all coefficients are 0.

Definition at line 91 of file PhotometryTransfo.cc.

        : _bbox(bbox),
          _toChebyshevRange(makeChebyshevRangeTransform(bbox)),
          _coefficients(_initializeChebyshev(order, identity)),
          _order(order),
          _nParameters((order + 1) * (order + 2) / 2) {}

PhotometryTransfoChebyshev() [2/2]

lsst::jointcal::PhotometryTransfoChebyshev::PhotometryTransfoChebyshev	(	ndarray::Array< double, 2, 2 > const &	coefficients,
		afw::geom::Box2D const &	bbox
	)

Create a Chebyshev transfo with the specified coefficients.

The polynomial order is determined from the number of coefficients, taking only the anti-diagonal upper triangle portion of the passed-in coefficients

Parameters

	coefficients	The polynomial coefficients.
[in]	bbox	The bounding box it is valid within, to rescale it to [-1,1].

Definition at line 99 of file PhotometryTransfo.cc.

        : _bbox(bbox),
          _toChebyshevRange(makeChebyshevRangeTransform(bbox)),
          _coefficients(coefficients),
          _order(coefficients.size() - 1),
          _nParameters((_order + 1) * (_order + 2) / 2) {}

Member Function Documentation

clone()

virtual std::shared_ptr<PhotometryTransfo> lsst::jointcal::PhotometryTransfo::clone ( ) const

pure virtualinherited

return a copy (allocated by new) of the transformation.

Implemented in lsst::jointcal::MagnitudeTransfoChebyshev, lsst::jointcal::FluxTransfoChebyshev, lsst::jointcal::MagnitudeTransfoSpatiallyInvariant, and lsst::jointcal::FluxTransfoSpatiallyInvariant.

computeChebyshev()

double lsst::jointcal::PhotometryTransfoChebyshev::computeChebyshev ( double  x, double  y  ) const

protected

Return the value of this polynomial at x,y.

For use in the sublcass transform() methods.

Definition at line 156 of file PhotometryTransfo.cc.

                                                                            {
    afw::geom::Point2D p = _toChebyshevRange(afw::geom::Point2D(x, y));
    return evaluateFunction1d(RecursionArrayImitator(_coefficients, p.getX()), p.getY(),
                              _coefficients.getSize<0>());
}

computeChebyshevDerivatives()

void lsst::jointcal::PhotometryTransfoChebyshev::computeChebyshevDerivatives ( double  x, double  y, Eigen::Ref< Eigen::VectorXd >  derivatives  ) const

protected

Set the derivatives of this polynomial at x,y.

For use in the sublcass computeParameterDerivatives() methods.

Definition at line 162 of file PhotometryTransfo.cc.

                                                                                                        {
    afw::geom::Point2D p = _toChebyshevRange(afw::geom::Point2D(x, y));
    // Algorithm: compute all the individual components recursively (since we'll need them anyway),
    // then combine them into the final answer vectors.
    Eigen::VectorXd Tnx(_order + 1);
    Eigen::VectorXd Tmy(_order + 1);
    Tnx[0] = 1;
    Tmy[0] = 1;
    if (_order >= 1) {
        Tnx[1] = p.getX();
        Tmy[1] = p.getY();
    }
    for (ndarray::Size i = 2; i <= _order; ++i) {
        Tnx[i] = 2 * p.getX() * Tnx[i - 1] - Tnx[i - 2];
        Tmy[i] = 2 * p.getY() * Tmy[i - 1] - Tmy[i - 2];
    }

    // NOTE: the indexing in this method and offsetParams must be kept consistent!
    Eigen::VectorXd::Index k = 0;
    for (ndarray::Size j = 0; j <= _order; ++j) {
        ndarray::Size const iMax = _order - j;  // to save re-computing `i+j <= order` every inner step.
        for (ndarray::Size i = 0; i <= iMax; ++i, ++k) {
            derivatives[k] = Tmy[j] * Tnx[i];
        }
    }
}

computeParameterDerivatives()

virtual void lsst::jointcal::PhotometryTransfo::computeParameterDerivatives ( double  x, double  y, double  value, Eigen::Ref< Eigen::VectorXd >  derivatives  ) const

pure virtualinherited

Compute the derivatives with respect to the parameters (i.e.

the coefficients).

Parameters

[in]	x	The x coordinate to compute at (in the appropriate units for this transfo).
[in]	y	The y coordinate to compute at (in the appropriate units for this transfo).
[in]	value	The instrument flux or magnitude to compute the derivative at.
[out]	derivatives	The computed derivatives, in the same order as the deltas in offsetParams.

Implemented in lsst::jointcal::MagnitudeTransfoChebyshev, lsst::jointcal::FluxTransfoChebyshev, lsst::jointcal::MagnitudeTransfoSpatiallyInvariant, and lsst::jointcal::FluxTransfoSpatiallyInvariant.

dump()

void lsst::jointcal::PhotometryTransfoChebyshev::dump ( std::ostream &  stream = std::cout ) const

inlineoverridevirtual

dumps the transfo coefficients to stream.

Implements lsst::jointcal::PhotometryTransfo.

Definition at line 245 of file PhotometryTransfo.h.

{ stream << _coefficients; }

getBBox()

afw::geom::Box2D lsst::jointcal::PhotometryTransfoChebyshev::getBBox ( ) const

inline

Definition at line 261 of file PhotometryTransfo.h.

{ return _bbox; }

getCoefficients()

ndarray::Array<double, 2, 2> lsst::jointcal::PhotometryTransfoChebyshev::getCoefficients ( ) const

inline

Get a copy of the coefficients of the polynomials, as a 2d array (NOTE: layout is [y][x])

Definition at line 254 of file PhotometryTransfo.h.

{ return ndarray::copy(_coefficients); }

getNpar()

int lsst::jointcal::PhotometryTransfoChebyshev::getNpar ( ) const

inlineoverridevirtual

Return the number of parameters (used to compute chisq)

Implements lsst::jointcal::PhotometryTransfo.

Definition at line 248 of file PhotometryTransfo.h.

{ return _nParameters; }

getOrder()

ndarray::Size lsst::jointcal::PhotometryTransfoChebyshev::getOrder ( ) const

inline

Definition at line 259 of file PhotometryTransfo.h.

{ return _order; }

getParameters()

Eigen::VectorXd lsst::jointcal::PhotometryTransfoChebyshev::getParameters ( ) const

overridevirtual

Get a copy of the parameters of this model, in the same order as offsetParams.

Implements lsst::jointcal::PhotometryTransfo.

Definition at line 142 of file PhotometryTransfo.cc.

                                                              {
    Eigen::VectorXd parameters(_nParameters);
    // NOTE: the indexing in this method and offsetParams must be kept consistent!
    Eigen::VectorXd::Index k = 0;
    for (ndarray::Size j = 0; j <= _order; ++j) {
        ndarray::Size const iMax = _order - j;  // to save re-computing `i+j <= order` every inner step.
        for (ndarray::Size i = 0; i <= iMax; ++i, ++k) {
            parameters[k] = _coefficients[j][i];
        }
    }

    return parameters;
}

mean()

double lsst::jointcal::PhotometryTransfoChebyshev::mean

(

)

const

Definition at line 140 of file PhotometryTransfo.cc.

{ return integrate() / _bbox.getArea(); }

offsetParams()

void lsst::jointcal::PhotometryTransfoChebyshev::offsetParams ( Eigen::VectorXd const &  delta )

overridevirtual

Offset the parameters by some (negative) amount during fitting.

Equivalent to flatten(parameters) -= delta

Ordering of delta is the same as the ordering of the derivatives returned from computeParameterDerivatives.

Implements lsst::jointcal::PhotometryTransfo.

Definition at line 107 of file PhotometryTransfo.cc.

                                                                        {
    // NOTE: the indexing in this method and computeParameterDerivatives must be kept consistent!
    Eigen::VectorXd::Index k = 0;
    for (ndarray::Size j = 0; j <= _order; ++j) {
        ndarray::Size const iMax = _order - j;  // to save re-computing `i+j <= order` every inner step.
        for (ndarray::Size i = 0; i <= iMax; ++i, ++k) {
            _coefficients[j][i] -= delta[k];
        }
    }
}

transform() [1/2]

virtual double lsst::jointcal::PhotometryTransfo::transform ( double  x, double  y, double  value  ) const

pure virtualinherited

Return the transform of value at (x,y).

Implemented in lsst::jointcal::MagnitudeTransfoChebyshev, lsst::jointcal::FluxTransfoChebyshev, lsst::jointcal::MagnitudeTransfoSpatiallyInvariant, and lsst::jointcal::FluxTransfoSpatiallyInvariant.

transform() [2/2]

double lsst::jointcal::PhotometryTransfo::transform ( Point const &  in, double  value  ) const

inlineinherited

Return the transformed value at Point(x,y).

Definition at line 58 of file PhotometryTransfo.h.

{ return transform(in.x, in.y, value); }

transformError() [1/2]

double lsst::jointcal::PhotometryTransfo::transformError ( Point const &  in, double  value, double  valueErr  ) const

inlineinherited

Return the transformed valueErr at Point(x,y).

Definition at line 64 of file PhotometryTransfo.h.

                                                                                {
        return transformError(in.x, in.y, value, valueErr);
    }

transformError() [2/2]

double lsst::jointcal::PhotometryTransfoChebyshev::transformError ( double  x, double  y, double  value, double  valueErr  ) const

inlineoverridevirtual

Return the transformed valueErr at Point(x,y).

Implements lsst::jointcal::PhotometryTransfo.

Definition at line 242 of file PhotometryTransfo.h.

{ return 0; }

---

The documentation for this class was generated from the following files:

• /j/snowflake/release/lsstsw/stack/Linux64/jointcal/16.0-25-g07af9f2+1/include/lsst/jointcal/PhotometryTransfo.h
• /j/snowflake/release/lsstsw/stack/Linux64/jointcal/16.0-25-g07af9f2+1/src/PhotometryTransfo.cc