LSST Applications g07dc498a13+5a531fccd6,g1409bbee79+5a531fccd6,g1a7e361dbc+5a531fccd6,g1fd858c14a+bae9e05889,g28da252d5a+b6acab2954,g33399d78f5+749e2df9f6,g35bb328faa+e55fef2c71,g3bd4b5ce2c+9ee8c95c05,g3c79e8cd92+5ffa104366,g43bc871e57+9a9eb91fab,g53246c7159+e55fef2c71,g60b5630c4e+f9e43d3906,g60ed82cc77+7534bcf0a9,g6e5c4a0e23+f441d97430,g78460c75b0+8427c4cc8f,g786e29fd12+307f82e6af,g8534526c7b+af2545e932,g89139ef638+5a531fccd6,g8b49a6ea8e+f9e43d3906,g9125e01d80+e55fef2c71,g989de1cb63+5a531fccd6,g9a9baf55bd+f1bd1a7c26,g9f33ca652e+c963d5c8aa,gaaedd4e678+5a531fccd6,gabe3b4be73+9c0c3c7524,gb1101e3267+ded3a614ca,gb58c049af0+28045f66fd,gc2fcbed0ba+f9e43d3906,gca43fec769+e55fef2c71,gcf25f946ba+749e2df9f6,gd6cbbdb0b4+784e334a77,gde0f65d7ad+500aeadb5f,ge278dab8ac+25667260f6,geab183fbe5+f9e43d3906,gecb8035dfe+0fa5abcb94,gefa07fa684+89734069dd,gf58bf46354+e55fef2c71,gfe7187db8c+5ee9e055d4,w.2025.02
LSST Data Management Base Package
Loading...
Searching...
No Matches
lsst::meas::extensions::photometryKron Namespace Reference

Namespaces

namespace  version
 

Classes

class  BadKronException
 
class  KronAperture
 
class  KronFluxAlgorithm
 A measurement algorithm that estimates flux using Kron photometry. More...
 
class  KronFluxControl
 C++ control object for Kron flux. More...
 

Functions

 PYBIND11_MODULE (photometryKron, mod)
 
template<typename ImageT >
std::pair< double, double > photometer (ImageT const &image, afw::geom::ellipses::Ellipse const &aperture, double const maxSincRadius)
 
double calculatePsfKronRadius (std::shared_ptr< afw::detection::Psf const > const &psf, geom::Point2D const &center, double smoothingSigma=0.0)
 
template std::shared_ptr< KronApertureKronAperture::determineRadius< afw::image::MaskedImage< float > > (afw::image::MaskedImage< float > const &, afw::geom::ellipses::Axes, geom::Point2D const &, KronFluxControl const &)
 
template std::pair< double, double > KronAperture::measureFlux< afw::image::MaskedImage< float > > (afw::image::MaskedImage< float > const &, double const, double const) const
 

Variables

 KronFluxPlugin
 
 KronFluxForcedPlugin
 
 name
 
 Control
 
 executionOrder
 
 shouldApCorr
 
 needsMetadata
 

Function Documentation

◆ calculatePsfKronRadius()

double lsst::meas::extensions::photometryKron::calculatePsfKronRadius ( std::shared_ptr< afw::detection::Psf const > const & psf,
geom::Point2D const & center,
double smoothingSigma = 0.0 )

Definition at line 338 of file KronPhotometry.cc.

343{
344 assert(psf);
345 double const radius = psf->computeShape(center).getDeterminantRadius();
346 // For a Gaussian N(0, sigma^2), the Kron radius is sqrt(pi/2)*sigma
347 return ::sqrt(geom::PI/2)*::hypot(radius, std::max(0.0, smoothingSigma));
348}
T max(T... args)
double constexpr PI
The ratio of a circle's circumference to diameter.
Definition Angle.h:40
Key< int > psf
Definition Exposure.cc:65

◆ KronAperture::determineRadius< afw::image::MaskedImage< float > >()

◆ KronAperture::measureFlux< afw::image::MaskedImage< float > >()

template std::pair< double, double > lsst::meas::extensions::photometryKron::KronAperture::measureFlux< afw::image::MaskedImage< float > > ( afw::image::MaskedImage< float > const & ,
double const ,
double const  ) const

◆ photometer()

template<typename ImageT >
std::pair< double, double > lsst::meas::extensions::photometryKron::photometer ( ImageT const & image,
afw::geom::ellipses::Ellipse const & aperture,
double const maxSincRadius )

Definition at line 311 of file KronPhotometry.cc.

316{
317 afw::geom::ellipses::Axes const& axes = aperture.getCore();
318 if (axes.getB() > maxSincRadius) {
319 FootprintFlux<ImageT> fluxFunctor;
320 auto spans = afw::geom::SpanSet::fromShape(aperture);
321 spans->applyFunctor(
322 fluxFunctor, *(image.getImage()), *(image.getVariance()));
323 return std::make_pair(fluxFunctor.getSum(), ::sqrt(fluxFunctor.getSumVar()));
324 }
325 try {
326 base::ApertureFluxResult fluxResult = base::ApertureFluxAlgorithm::computeSincFlux<float>(image, aperture);
327 return std::make_pair(fluxResult.instFlux, fluxResult.instFluxErr);
328 } catch(pex::exceptions::LengthError &e) {
329 LSST_EXCEPT_ADD(e, (boost::format("Measuring Kron flux for object at (%.3f, %.3f);"
330 " aperture radius %g,%g theta %g")
331 % aperture.getCenter().getX() % aperture.getCenter().getY()
332 % axes.getA() % axes.getB() % geom::radToDeg(axes.getTheta())).str());
333 throw e;
334 }
335}
#define LSST_EXCEPT_ADD(e, m)
Add the current location and a message to an existing exception before rethrowing it.
Definition Exception.h:54
An ellipse core for the semimajor/semiminor axis and position angle parametrization (a,...
Definition Axes.h:47
double const getTheta() const
Definition Axes.h:57
double const getA() const
Definition Axes.h:51
double const getB() const
Definition Axes.h:54
Reports attempts to exceed implementation-defined length limits for some classes.
Definition Runtime.h:76
T make_pair(T... args)
constexpr double radToDeg(double x) noexcept
Definition Angle.h:53
A Result struct for running an aperture flux algorithm with a single radius.
meas::base::Flux instFlux
Measured instFlux in DN.
meas::base::FluxErrElement instFluxErr
Standard deviation of instFlux in DN.

◆ PYBIND11_MODULE()

lsst::meas::extensions::photometryKron::PYBIND11_MODULE ( photometryKron ,
mod  )

Definition at line 132 of file photometryKron.cc.

132 {
133 py::module::import("lsst.afw.geom");
134 py::module::import("lsst.afw.image");
135 py::module::import("lsst.afw.table");
136 py::module::import("lsst.daf.base");
137
138 declareKronFluxControl(mod);
139 declareKronFluxAlgorithm(mod);
140 declareKronAperture(mod);
141}

Variable Documentation

◆ Control

lsst.meas.extensions.photometryKron.Control

Definition at line 31 of file __init__.py.

◆ executionOrder

lsst.meas.extensions.photometryKron.executionOrder

Definition at line 32 of file __init__.py.

◆ KronFluxForcedPlugin

lsst.meas.extensions.photometryKron.KronFluxForcedPlugin

Definition at line 28 of file __init__.py.

◆ KronFluxPlugin

lsst.meas.extensions.photometryKron.KronFluxPlugin

Definition at line 28 of file __init__.py.

◆ name

lsst.meas.extensions.photometryKron.name

Definition at line 30 of file __init__.py.

◆ needsMetadata

lsst.meas.extensions.photometryKron.needsMetadata

Definition at line 34 of file __init__.py.

◆ shouldApCorr

lsst.meas.extensions.photometryKron.shouldApCorr

Definition at line 33 of file __init__.py.