LSST Applications g04c3c9f7ca+2075667efa,g1e125bf412+5f448d5fcf,g2079a07aa2+3e9fd84d81,g2305ad1205+b635cf1488,g2bbee38e9b+6c6beb4891,g337abbeb29+6c6beb4891,g33d1c0ed96+6c6beb4891,g3a166c0a6a+6c6beb4891,g50ff169b8f+96c6868917,g52b1c1532d+585e252eca,g591dd9f2cf+42f171e1e6,g5c3423f6d4+d536b04327,g607f77f49a+d536b04327,g6f43f06aed+ca1339dc19,g858d7b2824+d536b04327,g8ee334c5b4+d7f9608c2f,g9963eaa53e+b3dc1655d3,g998f4353bf+d536b04327,g99cad8db69+8ef2408349,g9ddcbc5298+9a081db1e4,ga1e77700b3+2cbb763275,gadfd92a7e4+aec2f3b930,gae0086650b+585e252eca,gb0e22166c9+0e73c8378f,gb3b7280ab2+cb5fdb229e,gbb8dafda3b+a327199e22,gc120e1dc64+88074880ea,gc28159a63d+6c6beb4891,gcdd4ae20e8+bd241b2308,gcde1bda545+903e937d91,gcf0d15dbbd+bd241b2308,gdaeeff99f8+f9a426f77a,gddc38dedce+585e252eca,ge79ae78c31+6c6beb4891,gfbcc870c63+b310236976,w.2024.23
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)