evaluate.h File Reference

#include <iostream>
#include <memory>
#include <set>
#include <stdexcept>
#include <string>
#include <vector>
#include "gaussian.h"
#include "image.h"
#include "to_string.h"

Go to the source code of this file.

Classes
struct	lsst::gauss2d::detail::ValuesGauss
struct	lsst::gauss2d::detail::TermsMoments
struct	lsst::gauss2d::detail::Terms
class	lsst::gauss2d::detail::TermsPixel
	Storage for terms common to Gaussians for a single pixel. More...
class	lsst::gauss2d::detail::TermsGradient
class	lsst::gauss2d::detail::GradientsExtra< t, Data, Indices >
class	lsst::gauss2d::GaussianEvaluator< T, Data, Indices >
	A class that evaluates 2D Gaussians and renders them in images. More...

Namespaces
namespace	lsst
namespace	lsst::gauss2d
namespace	lsst::gauss2d::detail

Typedefs
typedef std::vector< double >	lsst::gauss2d::detail::vecd
typedef std::unique_ptr< vecd >	lsst::gauss2d::detail::vecdptr
typedef std::array< double, N_PARAMS_GAUSS2D >	lsst::gauss2d::detail::Weights
typedef std::vector< TermsPixel >	lsst::gauss2d::detail::TermsPixelVec
typedef std::vector< TermsGradient >	lsst::gauss2d::detail::TermsGradientVec

Enumerations
enum class	lsst::gauss2d::BackgroundType : unsigned char { lsst::gauss2d::none = 0 , lsst::gauss2d::constant = 1 }
enum class	lsst::gauss2d::OutputType : unsigned char { lsst::gauss2d::none = 0 , lsst::gauss2d::overwrite = 1 , lsst::gauss2d::add = 2 }
enum class	lsst::gauss2d::GradientType : unsigned char { lsst::gauss2d::none = 0 , lsst::gauss2d::loglike = 1 , lsst::gauss2d::jacobian = 2 }

Functions
TermsMoments	lsst::gauss2d::detail::gaussian_pixel_x_xx (const double cen_x, const double x_min, const double bin_x, const unsigned int dim_x, const double xx_weight, const double xy_weight)
Terms	lsst::gauss2d::detail::terms_from_covar (double weight, const Ellipse &ell)
template<class Data>
void	lsst::gauss2d::detail::gaussian_pixel (Data &image, const double weight, const double cen_x, const double cen_y, const double xx_weight, const double yy_weight, const double xy_weight)
void	lsst::gauss2d::detail::gaussian_pixel_get_jacobian (ValuesGauss &out, const double m, const double m_unweight, const double xmc_norm, const double ymc_weighted, const double xmc, const double ymc, const double norms_yy, const double xmc_t_xy_factor, const double sig_x_inv, const double sig_y_inv, const double xy_norm, const double xx_norm, const double yy_weighted, const double rho_div_one_m_rhosq, const double norms_xy_div_rho, const double dsig_x_conv_src=1, const double dsig_y_conv_src=1, const double drho_c_dsig_x_src=0, const double drho_c_dsig_y_src=0, const double drho_c_drho_s=1)
void	lsst::gauss2d::detail::gaussian_pixel_get_jacobian_from_terms (ValuesGauss &out, const size_t dim1, const TermsPixel &terms_pixel, const TermsGradient &terms_grad, double m, double m_unweighted, double xy_norm)
template<typename t>
void	lsst::gauss2d::detail::gaussian_pixel_add_values (t &cen_x, t &cen_y, t &L, t &sig_x, t &sig_y, t &rho, const ValuesGauss &values, const t weight_cen_x=1, const t weight_cen_y=1, const t weight_L=1, const t weight_sig_x=1, const t weight_sig_y=1, const t weight_rho=1)
template<typename t, class Data, class Indices, bool do_extra>
void	lsst::gauss2d::detail::gaussians_pixel_add_like_grad (Image< t, Data > *output, const Image< idx_type, Indices > &grad_param_map, const Image< t, Data > &grad_param_factor, const size_t N_GAUSS, const std::vector< Weights > &gaussweights, double &chi_pix, double &loglike, const double model, double data, double sigma_inv, unsigned int dim1, unsigned int dim2, const TermsPixelVec &terms_pixel, const TermsGradientVec &terms_grad, ValuesGauss &gradients, const std::shared_ptr< const Image< idx_type, Indices > > extra_param_map, const std::shared_ptr< const Image< t, Data > > extra_param_factor)
template<typename T, class Data, class Indices, GradientType gradient_type, bool do_extra>
T	lsst::gauss2d::detail::gaussian_pixel_add_all (size_t g, size_t j, size_t i, double weight, double sigma_inv, const TermsPixelVec &terms_pixel_vec, ImageArray< T, Data > &output_jac_ref, const Image< idx_type, Indices > &grad_param_map, const Image< T, Data > &grad_param_factor, std::vector< Weights > &gradweights, const TermsGradientVec &terms_grad_vec, ValuesGauss &gradients, GradientsExtra< T, Data, Indices > *grad_extra)
template<typename Indices>
const std::shared_ptr< const Indices >	lsst::gauss2d::detail::_param_map_default (size_t n_gaussians, size_t n_params=N_PARAMS_GAUSS2D, size_t increment=1)
template<typename Data>
const std::shared_ptr< const Data >	lsst::gauss2d::detail::_param_factor_default (size_t n_gaussians, size_t n_params=N_PARAMS_GAUSS2D, double value=1.)
template<typename T, class Data, class Indices>
std::shared_ptr< Data >	lsst::gauss2d::make_gaussians_pixel (const std::shared_ptr< const ConvolvedGaussians > gaussians, std::shared_ptr< Data > output=nullptr, const unsigned int n_rows=0, const unsigned int n_cols=0, const std::shared_ptr< const CoordinateSystem > coordsys=nullptr, bool to_add=false)
	Add gaussians to an image, creating one if needed.

Variables
const size_t	lsst::gauss2d::N_EXTRA_MAP = 2
const size_t	lsst::gauss2d::N_EXTRA_FACTOR = 3
const size_t	lsst::gauss2d::N_PARAMS_GAUSS2D = 6