LSST Applications g013ef56533+b8d55c8942,g083dd6704c+a047e97985,g199a45376c+0ba108daf9,g1fd858c14a+7a3b874d60,g210f2d0738+7416ca6900,g262e1987ae+1d557ba9a3,g29ae962dfc+519d34895e,g2cef7863aa+aef1011c0b,g30d7c61c20+36d16ea71a,g35bb328faa+8c5ae1fdc5,g3fd5ace14f+cb326ad149,g47891489e3+f459a6810c,g53246c7159+8c5ae1fdc5,g54cd7ddccb+890c8e1e5d,g5a60e81ecd+6240c63dbc,g64539dfbff+7416ca6900,g67b6fd64d1+f459a6810c,g6ebf1fc0d4+8c5ae1fdc5,g74acd417e5+0bae3c876a,g786e29fd12+668abc6043,g87389fa792+8856018cbb,g89139ef638+f459a6810c,g8d7436a09f+dee7680868,g8ea07a8fe4+81eaaadc04,g90f42f885a+34c0557caf,g97be763408+14b8164b5b,g98a1a72a9c+8389601a76,g98df359435+fff771c62d,gb8cb2b794d+6728931916,gbf99507273+8c5ae1fdc5,gc2a301910b+7416ca6900,gca7fc764a6+f459a6810c,gd7ef33dd92+f459a6810c,gdab6d2f7ff+0bae3c876a,ge410e46f29+f459a6810c,ge41e95a9f2+7416ca6900,geaed405ab2+e3b4b2a692,gf9a733ac38+8c5ae1fdc5,w.2025.43
LSST Data Management Base Package
Loading...
Searching...
No Matches
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< vecdlsst::gauss2d::detail::vecdptr
 
typedef std::array< double, N_PARAMS_GAUSS2Dlsst::gauss2d::detail::Weights
 
typedef std::vector< TermsPixellsst::gauss2d::detail::TermsPixelVec
 
typedef std::vector< TermsGradientlsst::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>
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