lsst::gauss2d::EllipseValues Class Reference

An EllipseData storing sigma_x, sigma_y, rho values as shared_ptrs. More...

#include <ellipse.h>

Inheritance diagram for lsst::gauss2d::EllipseValues:

Public Member Functions
	EllipseValues (std::shared_ptr< double > sigma_x, std::shared_ptr< double > sigma_y, std::shared_ptr< double > rho=nullptr)
	Construct a new EllipseValues object.
	EllipseValues (double sigma_x=0, double sigma_y=0, double rho=0)
	Construct a new EllipseValues object, creating new pointers for every value.
double	get_sigma_x () const override
	Get sigma_x.
double	get_sigma_y () const override
	Get sigma_y.
double	get_rho () const override
	Get rho.
std::array< double, 3 >	get_xyr () const override
	Get sigma_x, sigma_y, rho.
void	set (double sigma_x, double sigma_y, double rho) override
	Set sigma_x, sigma_y, rho.
void	set_h (double hwhm_x, double hwhm_y, double rho) override
	Set hwhm_x, hwhm_y, rho (half-width at half-max)
void	set_sigma_x (double sigma_x) override
	Set the x-axis dispersion (sigma)
void	set_sigma_y (double sigma_y) override
	Set the y-axis dispersion (sigma)
void	set_rho (double rho) override
	Set the correlation parameter (rho)
std::string	repr (bool name_keywords=false, std::string_view namespace_separator=Object::CC_NAMESPACE_SEPARATOR) const override
	Return a full, callable string representation of this.
std::string	str () const override
	Return a brief, human-readable string representation of this.
virtual void	convolve (const Ellipse &ell)
	Convolve this ellipse with another.
virtual double	get_area () const
	Return the area of this ellipse, equal to pi*sigma_major*sigma_minor.
virtual double	get_cov_xy () const
	Return the covariance, equal to sigma_x*sigma_y*rho.
virtual double	get_hwhm_x () const
	Get the x-axis half-width at half-maximum.
virtual double	get_hwhm_y () const
	Get the y-axis half-width at half-maximum.
virtual double	get_radius_trace () const
	Return the trace radius, equal to sqrt(sigma_x^2 + sigma_y^2)
virtual double	get_sigma_x_sq () const
	Get the square of sigma_x.
virtual double	get_sigma_y_sq () const
	Get the square of sigma_y.
virtual double	get_sigma_xy () const
	Return sigma_x*sigma_y.
virtual std::array< double, 3 >	get_hxyr () const
	Get hwhm_x, hwhm_y, rho.
virtual void	set (const Covariance &covar)
	Set values from a Covariance object.
virtual void	set (const EllipseMajor &ellipse)
	Set values from an EllipseMajor object.
virtual void	set_hwhm_x (double hwhm_x)
	Set the x-axis half-width at half-max (FWHM/2)
virtual void	set_hwhm_y (double hwhm_y)
	Set the y-axis half-width at half-max (FWHM/2)
virtual void	set_hxyr (const std::array< double, 3 > &hxyr)
	Set hwhm_x, hwhm_y, rho from an array.
virtual void	set_xyr (const std::array< double, 3 > &xyr)
	Set sigma_x, sigma_y, rho from an array.
bool	operator== (const EllipseData &other) const
bool	operator!= (const EllipseData &other) const

Static Public Member Functions
static void	check (double size_x, double size_y, double rho, std::string_view error_suffix="")
	Check whether Ellipse parameter values are valid, throwing if not.
static void	check_size (double size, std::string_view error_suffix="")
	Check whether an x- or x-yaxis size value is valid.
static void	check_rho (double rho, std::string_view error_suffix="")
	Check whether a rho value is valid.
static std::string_view	null_str (const std::string_view &namespace_separator)

Static Public Attributes
static constexpr std::string_view	CC_NAMESPACE_SEPARATOR = "::"
	The C++ namespace separator.
static constexpr std::string_view	NULL_STR_GENERAL = "None"
static constexpr std::string_view	PY_NAMESPACE_SEPARATOR = "."

Detailed Description

An EllipseData storing sigma_x, sigma_y, rho values as shared_ptrs.

This implementation stores values in shared_ptrs, allowing sharing of values with other instances.

Definition at line 232 of file ellipse.h.

Constructor & Destructor Documentation

EllipseValues() [1/2]

lsst::gauss2d::EllipseValues::EllipseValues ( std::shared_ptr< double > sigma_x, std::shared_ptr< double > sigma_y, std::shared_ptr< double > rho = nullptr )

inlineexplicit

Construct a new EllipseValues object.

Parameters

sigma_x	The sigma_x shared_ptr
sigma_y	The sigma_y shared_ptr
rho	The rho shared_ptr, defaulting to a new zero-valued double

Definition at line 241 of file ellipse.h.

            : _sigma_x(sigma_x == nullptr ? std::make_shared<double>(0) : std::move(sigma_x)),
              _sigma_y(sigma_y == nullptr ? std::make_shared<double>(0) : std::move(sigma_y)),
              _rho(rho == nullptr ? std::make_shared<double>(0) : std::move(rho)) {};

EllipseValues() [2/2]

lsst::gauss2d::EllipseValues::EllipseValues ( double sigma_x = 0, double sigma_y = 0, double rho = 0 )

inlineexplicit

Construct a new EllipseValues object, creating new pointers for every value.

Definition at line 247 of file ellipse.h.

            : _sigma_x(std::make_shared<double>(sigma_x)),
              _sigma_y(std::make_shared<double>(sigma_y)),
              _rho(std::make_shared<double>(rho)) {};

Member Function Documentation

check()

static void lsst::gauss2d::EllipseData::check ( double size_x, double size_y, double rho, std::string_view error_suffix = "" )

inlinestaticinherited

Check whether Ellipse parameter values are valid, throwing if not.

Definition at line 137 of file ellipse.h.

                                                                                                  {
        if (!(size_x >= 0) || !(size_y >= 0) || !(rho >= -1 && rho <= 1)) {
            throw std::invalid_argument("Invalid size_x, size_y, rho=" + to_string_float(size_x) + ","
                                        + to_string_float(size_y) + "," + to_string_float(rho)
                                        + "; sigma_x,y >= 0 and 1 >= rho >= -1 required."
                                        + std::string(error_suffix));
        }
    }

check_rho()

static void lsst::gauss2d::EllipseData::check_rho ( double rho, std::string_view error_suffix = "" )

inlinestaticinherited

Check whether a rho value is valid.

Definition at line 153 of file ellipse.h.

                                                                        {
        if (!(rho >= -1 && rho <= 1)) {
            throw std::invalid_argument("Invalid rho=" + to_string_float(rho) + "; 1 >= rho >= -1 required."
                                        + std::string(error_suffix));
        }
    }

check_size()

static void lsst::gauss2d::EllipseData::check_size ( double size, std::string_view error_suffix = "" )

inlinestaticinherited

Check whether an x- or x-yaxis size value is valid.

Definition at line 146 of file ellipse.h.

                                                                          {
        if (!(size >= 0)) {
            throw std::invalid_argument("Invalid size=" + to_string_float(size) + "; size >= 0 required."
                                        + std::string(error_suffix));
        }
    }

convolve()

void lsst::gauss2d::EllipseData::convolve ( const Ellipse & ell )

virtualinherited

Convolve this ellipse with another.

Definition at line 149 of file ellipse.cc.

                                             {
    double sigma_x_ell = ell.get_sigma_x();
    double sigma_y_ell = ell.get_sigma_y();
    double sigma_x = this->get_sigma_x();
    sigma_x = sqrt(sigma_x * sigma_x + sigma_x_ell * sigma_x_ell);
    double sigma_y = this->get_sigma_y();
    sigma_y = sqrt(sigma_y * sigma_y + sigma_y_ell * sigma_y_ell);
    double rho = (this->get_cov_xy() + ell.get_cov_xy()) / (sigma_x * sigma_y);
    this->set(sigma_x, sigma_y, rho);
}

get_area()

double lsst::gauss2d::EllipseData::get_area ( ) const

virtualinherited

Return the area of this ellipse, equal to pi*sigma_major*sigma_minor.

Definition at line 142 of file ellipse.cc.

                                   {
    double rho = this->get_rho();
    return M_PI * this->get_sigma_xy() * sqrt(1 - rho * rho);
}

get_cov_xy()

double lsst::gauss2d::EllipseData::get_cov_xy ( ) const

virtualinherited

Return the covariance, equal to sigma_x*sigma_y*rho.

Definition at line 160 of file ellipse.cc.

{ return this->get_sigma_x() * this->get_sigma_y() * this->get_rho(); }

get_hwhm_x()

double lsst::gauss2d::EllipseData::get_hwhm_x ( ) const

virtualinherited

Get the x-axis half-width at half-maximum.

Reimplemented in lsst::gauss2d::fit::GaussianParametricEllipse.

Definition at line 162 of file ellipse.cc.

{ return M_SIGMA_HWHM * this->get_sigma_x(); };

get_hwhm_y()

double lsst::gauss2d::EllipseData::get_hwhm_y ( ) const

virtualinherited

Get the y-axis half-width at half-maximum.

Reimplemented in lsst::gauss2d::fit::GaussianParametricEllipse.

Definition at line 164 of file ellipse.cc.

{ return M_SIGMA_HWHM * this->get_sigma_y(); };

get_hxyr()

std::array< double, 3 > lsst::gauss2d::EllipseData::get_hxyr ( ) const

virtualinherited

Get hwhm_x, hwhm_y, rho.

Reimplemented in lsst::gauss2d::fit::GaussianParametricEllipse.

Definition at line 166 of file ellipse.cc.

                                                {
    return {this->get_hwhm_x(), this->get_hwhm_y(), this->get_rho()};
}

get_radius_trace()

double lsst::gauss2d::EllipseData::get_radius_trace ( ) const

virtualinherited

Return the trace radius, equal to sqrt(sigma_x^2 + sigma_y^2)

Definition at line 173 of file ellipse.cc.

{ return sqrt(this->get_sigma_x_sq() + this->get_sigma_y_sq()); }

get_rho()

double lsst::gauss2d::EllipseValues::get_rho ( ) const

overridevirtual

Get rho.

Implements lsst::gauss2d::EllipseData.

Definition at line 187 of file ellipse.cc.

{ return *_rho; }

get_sigma_x()

double lsst::gauss2d::EllipseValues::get_sigma_x ( ) const

overridevirtual

Get sigma_x.

Implements lsst::gauss2d::EllipseData.

Definition at line 185 of file ellipse.cc.

{ return *_sigma_x; }

get_sigma_x_sq()

double lsst::gauss2d::EllipseData::get_sigma_x_sq ( ) const

virtualinherited

Get the square of sigma_x.

Definition at line 175 of file ellipse.cc.

                                         {
    double sigma = this->get_sigma_x();
    return sigma * sigma;
}

get_sigma_xy()

double lsst::gauss2d::EllipseData::get_sigma_xy ( ) const

virtualinherited

Return sigma_x*sigma_y.

Definition at line 147 of file ellipse.cc.

{ return this->get_sigma_x() * this->get_sigma_y(); }

get_sigma_y()

double lsst::gauss2d::EllipseValues::get_sigma_y ( ) const

overridevirtual

Get sigma_y.

Implements lsst::gauss2d::EllipseData.

Definition at line 186 of file ellipse.cc.

{ return *_sigma_y; }

get_sigma_y_sq()

double lsst::gauss2d::EllipseData::get_sigma_y_sq ( ) const

virtualinherited

Get the square of sigma_y.

Definition at line 180 of file ellipse.cc.

                                         {
    double sigma = this->get_sigma_y();
    return sigma * sigma;
}

get_xyr()

std::array< double, 3 > lsst::gauss2d::EllipseValues::get_xyr ( ) const

overridevirtual

Get sigma_x, sigma_y, rho.

Reimplemented from lsst::gauss2d::EllipseData.

Definition at line 188 of file ellipse.cc.

{ return {*_sigma_x, *_sigma_y, *_rho}; }

null_str()

static std::string_view lsst::gauss2d::Object::null_str ( const std::string_view & namespace_separator )

inlinestaticinherited

Definition at line 49 of file object.h.

                                                                              {
        return namespace_separator == CC_NAMESPACE_SEPARATOR ? "nullptr" : NULL_STR_GENERAL;
    }

operator!=()

bool lsst::gauss2d::EllipseData::operator!= ( const EllipseData & other ) const

inlineinherited

Definition at line 217 of file ellipse.h.

{ return !(*this == other); };

operator==()

bool lsst::gauss2d::EllipseData::operator== ( const EllipseData & other ) const

inlineinherited

Definition at line 216 of file ellipse.h.

{ return get_xyr() == other.get_xyr(); };

repr()

std::string lsst::gauss2d::EllipseValues::repr ( bool name_keywords = false, std::string_view namespace_separator = Object::CC_NAMESPACE_SEPARATOR ) const

overridevirtual

Return a full, callable string representation of this.

Parameters

name_keywords	Whether to prefix arguments with "{name}=", where name is the arg name in the header (as with keyword arguments in Python).
namespace_separator	The string to use to delimit namespaces, i.e. :: in C++ and . in Python.

Returns

A callable string representation of this, which should return an an identical object to this.

Note

The representation with name_keywords=false must be callable in C++. The representation with name_keywords=true should be callable in Python, if there are any bindings.

Implements lsst::gauss2d::EllipseData.

Definition at line 219 of file ellipse.cc.

                                                                                          {
    return type_name_str<EllipseValues>(false, namespace_separator) + "(" + (name_keywords ? "sigma_x=" : "")
           + to_string_float(this->get_sigma_x()) + ", " + (name_keywords ? "sigma_y=" : "")
           + to_string_float(this->get_sigma_y()) + ", " + (name_keywords ? "rho=" : "")
           + to_string_float(this->get_rho()) + ")";
}

set() [1/3]

void lsst::gauss2d::EllipseData::set ( const Covariance & covar )

virtualinherited

Set values from a Covariance object.

Definition at line 263 of file ellipse.cc.

                                             {
    double sigma_x = covar.get_sigma_x_sq();
    double sigma_y = covar.get_sigma_y_sq();
    sigma_x = sqrt(sigma_x);
    sigma_y = sqrt(sigma_y);
    double rho = (sigma_x == 0 || sigma_y == 0) ? 0 : covar.get_cov_xy() / (sigma_x * sigma_y);
    this->set(sigma_x, sigma_y, rho);
}

set() [2/3]

void lsst::gauss2d::EllipseData::set ( const EllipseMajor & ellipse )

virtualinherited

Set values from an EllipseMajor object.

Definition at line 272 of file ellipse.cc.

                                                 {
    const double r_major = ellipse.get_r_major();
    if (r_major == 0) {
        this->set(0, 0, 0);
        return;
    }
    const double axrat = ellipse.get_axrat();
    if (axrat == 1) {
        this->set(r_major, r_major, 0);
        return;
    }
    const auto [sin_th, cos_th] = sincos(ellipse.get_angle_radians());
    const double sin_th_sq = sin_th * sin_th;
    const double cos_th_sq = cos_th * cos_th;
 
    const double r_major_sq = r_major * r_major;
    const double r_minor_sq = r_major_sq * axrat * axrat;
    const double sigma_x = sqrt(cos_th_sq * r_major_sq + sin_th_sq * r_minor_sq);
    const double sigma_y = sqrt(sin_th_sq * r_major_sq + cos_th_sq * r_minor_sq);
    double rho = (sigma_x == 0 || sigma_y == 0)
                         ? 0
                         : sin_th * cos_th * (r_major_sq - r_minor_sq) / (sigma_x * sigma_y);
    this->set(sigma_x, sigma_y, rho);
}

set() [3/3]

void lsst::gauss2d::EllipseValues::set ( double sigma_x, double sigma_y, double rho )

overridevirtual

Set sigma_x, sigma_y, rho.

Reimplemented from lsst::gauss2d::EllipseData.

Definition at line 205 of file ellipse.cc.

                                                                  {
    Ellipse::check(sigma_x, sigma_y, rho);
    *_sigma_x = sigma_x;
    *_sigma_y = sigma_y;
    *_rho = rho;
}

set_h()

void lsst::gauss2d::EllipseValues::set_h ( double hwhm_x, double hwhm_y, double rho )

overridevirtual

Set hwhm_x, hwhm_y, rho (half-width at half-max)

Reimplemented from lsst::gauss2d::EllipseData.

Definition at line 212 of file ellipse.cc.

                                                                  {
    Ellipse::check(hwhm_x, hwhm_y, rho);
    *_sigma_x = M_HWHM_SIGMA * hwhm_x;
    *_sigma_y = M_HWHM_SIGMA * hwhm_y;
    *_rho = rho;
}

set_hwhm_x()

void lsst::gauss2d::EllipseData::set_hwhm_x ( double hwhm_x )

virtualinherited

Set the x-axis half-width at half-max (FWHM/2)

Reimplemented in lsst::gauss2d::fit::GaussianParametricEllipse.

Definition at line 303 of file ellipse.cc.

{ this->set_sigma_x(M_HWHM_SIGMA * hwhm_x); }

set_hwhm_y()

void lsst::gauss2d::EllipseData::set_hwhm_y ( double hwhm_y )

virtualinherited

Set the y-axis half-width at half-max (FWHM/2)

Reimplemented in lsst::gauss2d::fit::GaussianParametricEllipse.

Definition at line 304 of file ellipse.cc.

{ this->set_sigma_y(M_HWHM_SIGMA * hwhm_y); }

set_hxyr()

void lsst::gauss2d::EllipseData::set_hxyr ( const std::array< double, 3 > & hxyr )

virtualinherited

Set hwhm_x, hwhm_y, rho from an array.

Reimplemented in lsst::gauss2d::fit::GaussianParametricEllipse, and lsst::gauss2d::fit::SersicEllipseData.

Definition at line 305 of file ellipse.cc.

{ this->set_h(hxyr[0], hxyr[1], hxyr[2]); }

set_rho()

void lsst::gauss2d::EllipseValues::set_rho ( double rho )

overridevirtual

Set the correlation parameter (rho)

Implements lsst::gauss2d::EllipseData.

Definition at line 200 of file ellipse.cc.

                                      {
    EllipseData::check_rho(rho);
    *_rho = rho;
}

set_sigma_x()

void lsst::gauss2d::EllipseValues::set_sigma_x ( double sigma_x )

overridevirtual

Set the x-axis dispersion (sigma)

Implements lsst::gauss2d::EllipseData.

Definition at line 190 of file ellipse.cc.

                                              {
    EllipseData::check_size(sigma_x, "(size=sigma_x)");
    *_sigma_x = sigma_x;
}

set_sigma_y()

void lsst::gauss2d::EllipseValues::set_sigma_y ( double sigma_y )

overridevirtual

Set the y-axis dispersion (sigma)

Implements lsst::gauss2d::EllipseData.

Definition at line 195 of file ellipse.cc.

                                              {
    EllipseData::check_size(sigma_y, "(size=sigma_y)");
    *_sigma_y = sigma_y;
}

set_xyr()

void lsst::gauss2d::EllipseData::set_xyr ( const std::array< double, 3 > & xyr )

virtualinherited

Set sigma_x, sigma_y, rho from an array.

Reimplemented in lsst::gauss2d::fit::GaussianParametricEllipse, and lsst::gauss2d::fit::SersicEllipseData.

Definition at line 306 of file ellipse.cc.

{ this->set(xyr[0], xyr[1], xyr[2]); }

str()

std::string lsst::gauss2d::EllipseValues::str ( ) const

overridevirtual

Return a brief, human-readable string representation of this.

Implements lsst::gauss2d::EllipseData.

Definition at line 226 of file ellipse.cc.

                                   {
    return type_name_str<EllipseValues>(true) + "(sigma_x=" + to_string_float(this->get_sigma_x())
           + ", sigma_y=" + to_string_float(this->get_sigma_y()) + ", rho=" + to_string_float(this->get_rho())
           + ")";
}

Member Data Documentation

CC_NAMESPACE_SEPARATOR

std::string_view lsst::gauss2d::Object::CC_NAMESPACE_SEPARATOR = "::"

staticconstexprinherited

The C++ namespace separator.

Definition at line 45 of file object.h.

NULL_STR_GENERAL

std::string_view lsst::gauss2d::Object::NULL_STR_GENERAL = "None"

staticconstexprinherited

Definition at line 46 of file object.h.

PY_NAMESPACE_SEPARATOR

std::string_view lsst::gauss2d::Object::PY_NAMESPACE_SEPARATOR = "."

staticconstexprinherited

Definition at line 47 of file object.h.

---

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

• /j/snowflake/release/lsstsw/stack/lsst-scipipe-10.0.0/Linux64/gauss2d/g9125e01d80+fcb1d3bbc8/include/lsst/gauss2d/ellipse.h
• /j/snowflake/release/lsstsw/stack/lsst-scipipe-10.0.0/Linux64/gauss2d/g9125e01d80+fcb1d3bbc8/src/ellipse.cc