lsst.meas.algorithms  22.0.0-14-g0da61bd6+74b56c712e
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")),
63  averagePosition(afw::table::PointKey<double>::addFields(
64  schema, "averagePosition", "average position of stars used to make the PSF", "pixel")) {}
65 
66 bool KernelPsf::isPersistable() const noexcept { return _kernel->isPersistable(); }
67 
68 std::string KernelPsf::getPersistenceName() const { return "KernelPsf"; }
69 
70 std::string KernelPsf::getPythonModule() const { return "lsst.meas.algorithms"; }
71 
74  afw::table::BaseCatalog catalog = handle.makeCatalog(keys.schema);
75  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
afw::table::Schema schema
afw::table::PointKey< double > averagePosition
Definition: CoaddPsf.cc:340
Utilities for persisting KernelPsf and subclasses thereof.
#define PTR(...)
image::Image< Pixel > Image
lsst::geom::Extent2I const getDimensions() const
double computeImage(lsst::afw::image::Image< Pixel > &image, bool doNormalize, double x=0.0, double y=0.0) const
lsst::geom::Box2I getBBox() const
std::shared_ptr< BaseRecord > addNew()
static std::shared_ptr< T > dynamicCast(std::shared_ptr< Persistable > const &ptr)
io::OutputArchiveHandle OutputArchiveHandle
An intermediate base class for Psfs that use an image representation.
Definition: ImagePsf.h:40
A PersistableFactory for KernelPsf and its subclasses.
boost::shared_ptr< afw::detection::Psf > clone() const override
Polymorphic deep copy.
Definition: KernelPsf.cc:41
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
boost::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
Definition: KernelPsf.cc:68
std::string getPythonModule() const override
Definition: KernelPsf.cc:70
void write(OutputArchiveHandle &handle) const override
Definition: KernelPsf.cc:72
geom::Point2D getAveragePosition() const override
Return average position of stars; used as default position.
Definition: KernelPsf.cc:47
def keys(self)
clone
Point< double, 2 > Point2D
A read-only singleton struct containing the schema and key used in persistence for KernelPsf.
static KernelPsfPersistenceHelper const & get()
Definition: KernelPsf.cc:55