LSST Applications g063fba187b+fee0456c91,g0f08755f38+ea96e5a5a3,g1653933729+a8ce1bb630,g168dd56ebc+a8ce1bb630,g1a2382251a+90257ff92a,g20f6ffc8e0+ea96e5a5a3,g217e2c1bcf+937a289c59,g28da252d5a+daa7da44eb,g2bbee38e9b+253935c60e,g2bc492864f+253935c60e,g3156d2b45e+6e55a43351,g32e5bea42b+31359a2a7a,g347aa1857d+253935c60e,g35bb328faa+a8ce1bb630,g3a166c0a6a+253935c60e,g3b1af351f3+a8ce1bb630,g3e281a1b8c+c5dd892a6c,g414038480c+416496e02f,g41af890bb2+afe91b1188,g599934f4f4+0db33f7991,g7af13505b9+e36de7bce6,g80478fca09+da231ba887,g82479be7b0+a4516e59e3,g858d7b2824+ea96e5a5a3,g89c8672015+f4add4ffd5,g9125e01d80+a8ce1bb630,ga5288a1d22+bc6ab8dfbd,gb58c049af0+d64f4d3760,gc28159a63d+253935c60e,gcab2d0539d+3f2b72788c,gcf0d15dbbd+4ea9c45075,gda6a2b7d83+4ea9c45075,gdaeeff99f8+1711a396fd,ge79ae78c31+253935c60e,gef2f8181fd+3031e3cf99,gf0baf85859+c1f95f4921,gfa517265be+ea96e5a5a3,gfa999e8aa5+17cd334064,w.2024.50
LSST Data Management Base Package
Loading...
Searching...
No Matches
DeltaFunctionKernel.cc
Go to the documentation of this file.
1// -*- LSST-C++ -*-
2
3/*
4 * LSST Data Management System
5 * Copyright 2008, 2009, 2010 LSST Corporation.
6 *
7 * This product includes software developed by the
8 * LSST Project (http://www.lsst.org/).
9 *
10 * This program is free software: you can redistribute it and/or modify
11 * it under the terms of the GNU General Public License as published by
12 * the Free Software Foundation, either version 3 of the License, or
13 * (at your option) any later version.
14 *
15 * This program is distributed in the hope that it will be useful,
16 * but WITHOUT ANY WARRANTY; without even the implied warranty of
17 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18 * GNU General Public License for more details.
19 *
20 * You should have received a copy of the LSST License Statement and
21 * the GNU General Public License along with this program. If not,
22 * see <http://www.lsstcorp.org/LegalNotices/>.
23 */
24#include <sstream>
25#include <vector>
26
27#include "lsst/pex/exceptions.h"
31
33
34namespace lsst {
35namespace afw {
36
37template std::shared_ptr<math::DeltaFunctionKernel> table::io::PersistableFacade<
38 math::DeltaFunctionKernel>::dynamicCast(std::shared_ptr<table::io::Persistable> const&);
39
40namespace math {
41
43 : Kernel(width, height, 0), _pixel(point) {
44 if (point.getX() < 0 || point.getX() >= width || point.getY() < 0 || point.getY() >= height) {
46 os << "point (" << point.getX() << ", " << point.getY() << ") lies outside " << width << "x" << height
47 << " sized kernel";
49 }
50}
51
54 new DeltaFunctionKernel(this->getWidth(), this->getHeight(), this->_pixel));
55 retPtr->setCtr(this->getCtr());
56 return retPtr;
57}
58
60 int padX = width - getWidth();
61 int padY = height - getHeight();
62 if ((padX % 2) || (padY % 2)) {
64 os << "Cannot resize DeltaFunctionKernel from (" << getWidth() << ", " << getHeight() << ") to ("
65 << width << ", " << height << "), because at least one dimension would change by an odd value.";
67 }
68 int newPixelX = getPixel().getX() + padX / 2;
69 int newPixelY = getPixel().getY() + padY / 2;
71 std::make_shared<DeltaFunctionKernel>(width, height, lsst::geom::Point2I(newPixelX, newPixelY));
72 return retPtr;
73}
74
76 const int pixelX = getPixel().getX(); // active pixel in Kernel
77 const int pixelY = getPixel().getY();
78
80 os << prefix << "DeltaFunctionKernel:" << std::endl;
81 os << prefix << "Pixel (c,r) " << pixelX << "," << pixelY << ")" << std::endl;
82 os << Kernel::toString(prefix + "\t");
83 return os.str();
84}
85
87 const int pixelX = getPixel().getX(); // active pixel in Kernel
88 const int pixelY = getPixel().getY();
89
90 image = 0;
91 *image.xy_at(pixelX, pixelY) = 1;
92
93 return 1;
94}
95
96// ------ Persistence ---------------------------------------------------------------------------------------
97
98namespace {
99
100struct DeltaFunctionKernelPersistenceHelper : public Kernel::PersistenceHelper {
101 table::PointKey<int> pixel;
102
103 static DeltaFunctionKernelPersistenceHelper const& get() {
104 static DeltaFunctionKernelPersistenceHelper const instance;
105 return instance;
106 }
107
108 // No copying
109 DeltaFunctionKernelPersistenceHelper(const DeltaFunctionKernelPersistenceHelper&) = delete;
110 DeltaFunctionKernelPersistenceHelper& operator=(const DeltaFunctionKernelPersistenceHelper&) = delete;
111
112 // No moving
113 DeltaFunctionKernelPersistenceHelper(DeltaFunctionKernelPersistenceHelper&&) = delete;
114 DeltaFunctionKernelPersistenceHelper& operator=(DeltaFunctionKernelPersistenceHelper&&) = delete;
115
116private:
117 explicit DeltaFunctionKernelPersistenceHelper()
118 : Kernel::PersistenceHelper(0),
119 pixel(table::PointKey<int>::addFields(schema, "pixel", "position of nonzero pixel", "pixel")) {}
120};
121
122} // namespace
123
125public:
127 CatalogVector const& catalogs) const override {
128 LSST_ARCHIVE_ASSERT(catalogs.size() == 1u);
129 LSST_ARCHIVE_ASSERT(catalogs.front().size() == 1u);
130 DeltaFunctionKernelPersistenceHelper const& keys = DeltaFunctionKernelPersistenceHelper::get();
131 LSST_ARCHIVE_ASSERT(catalogs.front().getSchema() == keys.schema);
132 afw::table::BaseRecord const& record = catalogs.front().front();
134 new DeltaFunctionKernel(record.get(keys.dimensions.getX()),
135 record.get(keys.dimensions.getY()), record.get(keys.pixel)));
136 result->setCtr(record.get(keys.center));
137 return result;
138 }
139
140 explicit Factory(std::string const& name) : afw::table::io::PersistableFactory(name) {}
141};
142
143namespace {
144
145std::string getDeltaFunctionKernelPersistenceName() { return "DeltaFunctionKernel"; }
146
147DeltaFunctionKernel::Factory registration(getDeltaFunctionKernelPersistenceName());
148
149} // namespace
150
152 return getDeltaFunctionKernelPersistenceName();
153}
154
156 DeltaFunctionKernelPersistenceHelper const& keys = DeltaFunctionKernelPersistenceHelper::get();
157 std::shared_ptr<afw::table::BaseRecord> record = keys.write(handle, *this);
158 record->set(keys.pixel, _pixel);
159}
160} // namespace math
161} // namespace afw
162} // namespace lsst
py::object result
Definition _schema.cc:429
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition Exception.h:48
std::ostream * os
Definition Schema.cc:557
std::string prefix
#define LSST_ARCHIVE_ASSERT(EXPR)
An assertion macro used to validate the structure of an InputArchive.
Definition Persistable.h:48
table::Schema schema
Definition python.h:134
A class to represent a 2-dimensional array of pixels.
Definition Image.h:51
std::shared_ptr< afw::table::io::Persistable > read(InputArchive const &archive, CatalogVector const &catalogs) const override
Construct a new object from the given InputArchive and vector of catalogs.
DeltaFunctionKernel(int width, int height, lsst::geom::Point2I const &point)
Construct a spatially invariant DeltaFunctionKernel.
std::shared_ptr< Kernel > resized(int width, int height) const override
Return a pointer to a clone with specified kernel dimensions.
std::shared_ptr< Kernel > clone() const override
Return a pointer to a deep copy of this kernel.
double doComputeImage(lsst::afw::image::Image< Pixel > &image, bool doNormalize) const override
Low-level version of computeImage.
lsst::geom::Point2I getPixel() const
Definition Kernel.h:670
std::string getPersistenceName() const override
Return the unique name used to persist this object and look up its factory.
std::string toString(std::string const &prefix="") const override
Return a string representation of the kernel.
void write(OutputArchiveHandle &handle) const override
Write the object to one or more catalogs.
Kernels are used for convolution with MaskedImages and (eventually) Images.
Definition Kernel.h:110
int getHeight() const
Return the Kernel's height.
Definition Kernel.h:229
lsst::geom::Point2I getCtr() const
Return index of kernel's center.
Definition Kernel.h:234
virtual std::string toString(std::string const &prefix="") const
Return a string representation of the kernel.
Definition Kernel.cc:185
int getWidth() const
Return the Kernel's width.
Definition Kernel.h:224
Base class for all records.
Definition BaseRecord.h:31
Field< T >::Value get(Key< T > const &key) const
Return the value of a field for the given key.
Definition BaseRecord.h:151
A vector of catalogs used by Persistable.
A multi-catalog archive object used to load table::io::Persistable objects.
An object passed to Persistable::write to allow it to persist itself.
A base class for factory classes used to reconstruct objects from records.
Reports invalid arguments.
Definition Runtime.h:66
T endl(T... args)