Loading [MathJax]/extensions/tex2jax.js
LSST Applications 28.0.0,g1653933729+a8ce1bb630,g1a997c3884+a8ce1bb630,g28da252d5a+5bd70b7e6d,g2bbee38e9b+638fca75ac,g2bc492864f+638fca75ac,g3156d2b45e+07302053f8,g347aa1857d+638fca75ac,g35bb328faa+a8ce1bb630,g3a166c0a6a+638fca75ac,g3e281a1b8c+7bbb0b2507,g4005a62e65+17cd334064,g414038480c+5b5cd4fff3,g41af890bb2+4ffae9de63,g4e1a3235cc+0f1912dca3,g6249c6f860+3c3976f90c,g80478fca09+46aba80bd6,g82479be7b0+77990446f6,g858d7b2824+78ba4d1ce1,g89c8672015+f667a5183b,g9125e01d80+a8ce1bb630,ga5288a1d22+2a6264e9ca,gae0086650b+a8ce1bb630,gb58c049af0+d64f4d3760,gc22bb204ba+78ba4d1ce1,gc28159a63d+638fca75ac,gcf0d15dbbd+32ddb6096f,gd6b7c0dfd1+3e339405e9,gda3e153d99+78ba4d1ce1,gda6a2b7d83+32ddb6096f,gdaeeff99f8+1711a396fd,gdd5a9049c5+b18c39e5e3,ge2409df99d+a5e4577cdc,ge33fd446bb+78ba4d1ce1,ge79ae78c31+638fca75ac,gf0baf85859+64e8883e75,gf5289d68f6+e1b046a8d7,gfa443fc69c+91d9ed1ecf,gfda6b12a05+8419469a56
LSST Data Management Base Package
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
KernelPsf.cc
Go to the documentation of this file.
1// -*- LSST-C++ -*-
2
3#include "lsst/geom/Box.h"
7
8namespace lsst {
9namespace afw {
10namespace table {
11namespace io {
12
15
16} // namespace io
17} // namespace table
18} // namespace afw
19namespace meas {
20namespace algorithms {
21
23KernelPsf::doComputeKernelImage(geom::Point2D const& position, afw::image::Color const& color) const {
25 _kernel->computeImage(*im, true, position.getX(), position.getY());
26 return im;
27}
28
29geom::Box2I KernelPsf::doComputeBBox(geom::Point2D const& position, afw::image::Color const& color) const {
30 return _kernel->getBBox();
31}
32
33KernelPsf::KernelPsf(afw::math::Kernel const& kernel, geom::Point2D const& averagePosition)
34 : ImagePsf(!kernel.isSpatiallyVarying()),
35 _kernel(kernel.clone()),
36 _averagePosition(averagePosition) {}
37
39 : ImagePsf(!kernel->isSpatiallyVarying()), _kernel(kernel), _averagePosition(averagePosition) {}
40
42
44 return std::make_shared<KernelPsf>(*_kernel->resized(width, height), _averagePosition);
45}
46
47geom::Point2D KernelPsf::getAveragePosition() const { return _averagePosition; }
48
49namespace {
50
51KernelPsfFactory<> registration("KernelPsf");
52
53} // namespace
54
56 static KernelPsfPersistenceHelper instance;
57 return instance;
58}
59
60KernelPsfPersistenceHelper::KernelPsfPersistenceHelper()
61 : schema(),
62 kernel(schema.addField<int>("kernel", "archive ID of nested kernel object")),
63 averagePosition(afw::table::PointKey<double>::addFields(
64 schema, "averagePosition", "average position of stars used to make the PSF", "pixel")) {}
65
66bool KernelPsf::isPersistable() const noexcept { return _kernel->isPersistable(); }
67
68std::string KernelPsf::getPersistenceName() const { return "KernelPsf"; }
69
70std::string KernelPsf::getPythonModule() const { return "lsst.meas.algorithms"; }
71
74 afw::table::BaseCatalog catalog = handle.makeCatalog(keys.schema);
75 std::shared_ptr<afw::table::BaseRecord> record = catalog.addNew();
76 record->set(keys.kernel, handle.put(_kernel));
77 record->set(keys.averagePosition, _averagePosition);
78 handle.saveCatalog(catalog);
79}
80
81} // namespace algorithms
82} // namespace meas
83} // namespace lsst
Utilities for persisting KernelPsf and subclasses thereof.
Describe the colour of a source.
Definition Color.h:25
Kernels are used for convolution with MaskedImages and (eventually) Images.
Definition Kernel.h:110
lsst::geom::Extent2I const getDimensions() const
Return the Kernel's dimensions (width, height)
Definition Kernel.h:212
double computeImage(lsst::afw::image::Image< Pixel > &image, bool doNormalize, double x=0.0, double y=0.0) const
Compute an image (pixellized representation of the kernel) in place.
Definition Kernel.cc:76
lsst::geom::Box2I getBBox() const
return parent bounding box, with XY0 = -center
Definition Kernel.h:239
An object passed to Persistable::write to allow it to persist itself.
void saveCatalog(BaseCatalog const &catalog)
Save a catalog in the archive.
BaseCatalog makeCatalog(Schema const &schema)
Return a new, empty catalog with the given schema.
int put(Persistable const *obj, bool permissive=false)
Save an object to the archive and return a unique ID that can be used to retrieve it from an InputArc...
static std::shared_ptr< T > dynamicCast(std::shared_ptr< Persistable > const &ptr)
Dynamically cast a shared_ptr.
An integer coordinate rectangle.
Definition Box.h:55
An intermediate base class for Psfs that use an image representation.
Definition ImagePsf.h:40
A PersistableFactory for KernelPsf and its subclasses.
KernelPsf(afw::math::Kernel const &kernel, geom::Point2D const &averagePosition=geom::Point2D())
Construct a KernelPsf with a clone of the given kernel.
Definition KernelPsf.cc:33
bool isPersistable() const noexcept override
Whether this object is persistable; just delegates to the kernel.
Definition KernelPsf.cc:66
std::shared_ptr< afw::detection::Psf > clone() const override
Polymorphic deep copy.
Definition KernelPsf.cc:41
std::shared_ptr< afw::detection::Psf > resized(int width, int height) const override
Return a clone with specified kernel dimensions.
Definition KernelPsf.cc:43
std::string getPersistenceName() const override
Return the unique name used to persist this object and look up its factory.
Definition KernelPsf.cc:68
std::string getPythonModule() const override
Return the fully-qualified Python module that should be imported to guarantee that its factory is reg...
Definition KernelPsf.cc:70
void write(OutputArchiveHandle &handle) const override
Write the object to one or more catalogs.
Definition KernelPsf.cc:72
geom::Point2D getAveragePosition() const override
Return average position of stars; used as default position.
Definition KernelPsf.cc:47
T make_shared(T... args)
afw::table::PointKey< double > averagePosition
Definition CoaddPsf.cc:354
A read-only singleton struct containing the schema and key used in persistence for KernelPsf.
static KernelPsfPersistenceHelper const & get()
Definition KernelPsf.cc:55