LSST Applications g042eb84c57+730a74494b,g04e9c324dd+8c5ae1fdc5,g134cb467dc+1f1e3e7524,g199a45376c+0ba108daf9,g1fd858c14a+fa7d31856b,g210f2d0738+f66ac109ec,g262e1987ae+83a3acc0e5,g29ae962dfc+d856a2cb1f,g2cef7863aa+aef1011c0b,g35bb328faa+8c5ae1fdc5,g3fd5ace14f+a1e0c9f713,g47891489e3+0d594cb711,g4d44eb3520+c57ec8f3ed,g4d7b6aa1c5+f66ac109ec,g53246c7159+8c5ae1fdc5,g56a1a4eaf3+fd7ad03fde,g64539dfbff+f66ac109ec,g67b6fd64d1+0d594cb711,g67fd3c3899+f66ac109ec,g6985122a63+0d594cb711,g74acd417e5+3098891321,g786e29fd12+668abc6043,g81db2e9a8d+98e2ab9f28,g87389fa792+8856018cbb,g89139ef638+0d594cb711,g8d7436a09f+80fda9ce03,g8ea07a8fe4+760ca7c3fc,g90f42f885a+033b1d468d,g97be763408+a8a29bda4b,g99822b682c+e3ec3c61f9,g9d5c6a246b+0d5dac0c3d,ga41d0fce20+9243b26dd2,gbf99507273+8c5ae1fdc5,gd7ef33dd92+0d594cb711,gdab6d2f7ff+3098891321,ge410e46f29+0d594cb711,geaed405ab2+c4bbc419c6,gf9a733ac38+8c5ae1fdc5,w.2025.38
LSST Data Management Base Package
Loading...
Searching...
No Matches
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
13template std::shared_ptr<meas::algorithms::KernelPsf>
14PersistableFacade<meas::algorithms::KernelPsf>::dynamicCast(std::shared_ptr<Persistable> const&);
15
16} // namespace io
17} // namespace table
18} // namespace afw
19namespace meas {
20namespace algorithms {
21
22std::shared_ptr<afw::detection::Psf::Image>
24 std::shared_ptr<Psf::Image> im = std::make_shared<Psf::Image>(_kernel->getDimensions());
25 _kernel->computeImage(*im, true, position.getX(), position.getY());
26 return im;
27}
28
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
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:23
Kernels are used for convolution with MaskedImages and (eventually) Images.
Definition Kernel.h:110
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.
io::OutputArchiveHandle OutputArchiveHandle
An integer coordinate rectangle.
Definition Box.h:55
ImagePsf(bool isFixed=false)
Definition ImagePsf.h:42
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
geom::Box2I doComputeBBox(geom::Point2D const &position, afw::image::Color const &color) const override
Definition KernelPsf.cc:29
std::shared_ptr< Image > doComputeKernelImage(geom::Point2D const &position, afw::image::Color const &color) const override
These virtual member functions are private, not protected, because we only want derived classes to im...
Definition KernelPsf.cc:23
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)
CatalogT< BaseRecord > BaseCatalog
Definition fwd.h:72
Point< double, 2 > Point2D
Definition Point.h:324
A read-only singleton struct containing the schema and key used in persistence for KernelPsf.
KernelPsfPersistenceHelper(const KernelPsfPersistenceHelper &)=delete
static KernelPsfPersistenceHelper const & get()
Definition KernelPsf.cc:55