lsst.meas.algorithms  17.0.1-11-gf0f4e679+2
KernelPsf.cc
Go to the documentation of this file.
1 // -*- LSST-C++ -*-
2 
3 #include "lsst/geom/Box.h"
7 
8 namespace lsst {
9 namespace afw {
10 namespace table {
11 namespace io {
12 
15 
16 } // namespace io
17 } // namespace table
18 } // namespace afw
19 namespace meas {
20 namespace algorithms {
21 
23 KernelPsf::doComputeKernelImage(geom::Point2D const& position, afw::image::Color const& color) const {
24  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 
29 geom::Box2I KernelPsf::doComputeBBox(geom::Point2D const& position, afw::image::Color const& color) const {
30  return _kernel->getBBox();
31 }
32 
33 KernelPsf::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 
41 PTR(afw::detection::Psf) KernelPsf::clone() const { return std::make_shared<KernelPsf>(*this); }
42 
43 PTR(afw::detection::Psf) KernelPsf::resized(int width, int height) const {
44  return std::make_shared<KernelPsf>(*_kernel->resized(width, height), _averagePosition);
45 }
46 
47 geom::Point2D KernelPsf::getAveragePosition() const { return _averagePosition; }
48 
49 namespace {
50 
51 KernelPsfFactory<> registration("KernelPsf");
52 
53 } // namespace
54 
56  static KernelPsfPersistenceHelper instance;
57  return instance;
58 }
59 
60 KernelPsfPersistenceHelper::KernelPsfPersistenceHelper()
61  : schema(),
62  kernel(schema.addField<int>("kernel", "archive ID of nested kernel object")),
64  schema, "averagePosition", "average position of stars used to make the PSF", "pixel")) {
65  schema.getCitizen().markPersistent();
66 }
67 
68 bool KernelPsf::isPersistable() const noexcept { return _kernel->isPersistable(); }
69 
70 std::string KernelPsf::getPersistenceName() const { return "KernelPsf"; }
71 
72 std::string KernelPsf::getPythonModule() const { return "lsst.meas.algorithms"; }
73 
76  afw::table::BaseCatalog catalog = handle.makeCatalog(keys.schema);
77  PTR(afw::table::BaseRecord) record = catalog.addNew();
78  record->set(keys.kernel, handle.put(_kernel));
79  record->set(keys.averagePosition, _averagePosition);
80  handle.saveCatalog(catalog);
81 }
82 
83 } // namespace algorithms
84 } // namespace meas
85 } // namespace lsst
A read-only singleton struct containing the schema and key used in persistence for KernelPsf...
def keys(self)
boost::shared_ptr< afw::detection::Psf > resized(int width, int height) const override
Return a clone with specified kernel dimensions.
Definition: KernelPsf.cc:43
static std::shared_ptr< T > dynamicCast(std::shared_ptr< Persistable > const &ptr)
#define PTR(...)
Point< double, 2 > Point2D
STL class.
void write(OutputArchiveHandle &handle) const override
Definition: KernelPsf.cc:74
afw::table::Schema schema
afw::table::PointKey< double > averagePosition
Definition: CoaddPsf.cc:340
static KernelPsfPersistenceHelper const & get()
Definition: KernelPsf.cc:55
std::string getPythonModule() const override
Definition: KernelPsf.cc:72
boost::shared_ptr< afw::detection::Psf > clone() const override
Polymorphic deep copy.
Definition: KernelPsf.cc:41
io::OutputArchiveHandle OutputArchiveHandle
Utilities for persisting KernelPsf and subclasses thereof.
An intermediate base class for Psfs that use an image representation.
Definition: ImagePsf.h:40
A PersistableFactory for KernelPsf and its subclasses.
bool isPersistable() const noexcept override
Whether this object is persistable; just delegates to the kernel.
Definition: KernelPsf.cc:68
std::string getPersistenceName() const override
Definition: KernelPsf.cc:70
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
clone
image::Image< Pixel > Image
geom::Point2D getAveragePosition() const override
Return average position of stars; used as default position.
Definition: KernelPsf.cc:47
std::shared_ptr< BaseRecord > addNew()