LSST Applications 27.0.0,g0265f82a02+469cd937ee,g02d81e74bb+21ad69e7e1,g1470d8bcf6+cbe83ee85a,g2079a07aa2+e67c6346a6,g212a7c68fe+04a9158687,g2305ad1205+94392ce272,g295015adf3+81dd352a9d,g2bbee38e9b+469cd937ee,g337abbeb29+469cd937ee,g3939d97d7f+72a9f7b576,g487adcacf7+71499e7cba,g50ff169b8f+5929b3527e,g52b1c1532d+a6fc98d2e7,g591dd9f2cf+df404f777f,g5a732f18d5+be83d3ecdb,g64a986408d+21ad69e7e1,g858d7b2824+21ad69e7e1,g8a8a8dda67+a6fc98d2e7,g99cad8db69+f62e5b0af5,g9ddcbc5298+d4bad12328,ga1e77700b3+9c366c4306,ga8c6da7877+71e4819109,gb0e22166c9+25ba2f69a1,gb6a65358fc+469cd937ee,gbb8dafda3b+69d3c0e320,gc07e1c2157+a98bf949bb,gc120e1dc64+615ec43309,gc28159a63d+469cd937ee,gcf0d15dbbd+72a9f7b576,gdaeeff99f8+a38ce5ea23,ge6526c86ff+3a7c1ac5f1,ge79ae78c31+469cd937ee,gee10cc3b42+a6fc98d2e7,gf1cff7945b+21ad69e7e1,gfbcc870c63+9a11dc8c8f
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
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)