LSST Applications g0f08755f38+82efc23009,g12f32b3c4e+e7bdf1200e,g1653933729+a8ce1bb630,g1a0ca8cf93+50eff2b06f,g28da252d5a+52db39f6a5,g2bbee38e9b+37c5a29d61,g2bc492864f+37c5a29d61,g2cdde0e794+c05ff076ad,g3156d2b45e+41e33cbcdc,g347aa1857d+37c5a29d61,g35bb328faa+a8ce1bb630,g3a166c0a6a+37c5a29d61,g3e281a1b8c+fb992f5633,g414038480c+7f03dfc1b0,g41af890bb2+11b950c980,g5fbc88fb19+17cd334064,g6b1c1869cb+12dd639c9a,g781aacb6e4+a8ce1bb630,g80478fca09+72e9651da0,g82479be7b0+04c31367b4,g858d7b2824+82efc23009,g9125e01d80+a8ce1bb630,g9726552aa6+8047e3811d,ga5288a1d22+e532dc0a0b,gae0086650b+a8ce1bb630,gb58c049af0+d64f4d3760,gc28159a63d+37c5a29d61,gcf0d15dbbd+2acd6d4d48,gd7358e8bfb+778a810b6e,gda3e153d99+82efc23009,gda6a2b7d83+2acd6d4d48,gdaeeff99f8+1711a396fd,ge2409df99d+6b12de1076,ge79ae78c31+37c5a29d61,gf0baf85859+d0a5978c5a,gf3967379c6+4954f8c433,gfb92a5be7c+82efc23009,gfec2e1e490+2aaed99252,w.2024.46
LSST Data Management Base Package
Loading...
Searching...
No Matches
FixedKernel.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 <stdexcept>
25#include <memory>
26#include <string>
27
28#include "lsst/pex/exceptions.h"
32
34
35namespace lsst {
36namespace afw {
37
40
41namespace math {
42
43FixedKernel::FixedKernel() : Kernel(), _image(), _sum(0) {}
44
46 : Kernel(image.getWidth(), image.getHeight(), 0), _image(image, true), _sum(0) {
48 double imSum = 0.0;
49 for (int y = 0; y != image.getHeight(); ++y) {
50 for (XIter imPtr = image.row_begin(y), imEnd = image.row_end(y); imPtr != imEnd; ++imPtr) {
51 imSum += *imPtr;
52 }
53 }
54 this->_sum = imSum;
55}
56
58 : Kernel(kernel.getWidth(), kernel.getHeight(), 0), _image(kernel.getDimensions()), _sum(0) {
59 _sum = kernel.computeImage(_image, false, pos[0], pos[1]);
60}
61
63 std::shared_ptr<Kernel> retPtr(new FixedKernel(_image));
64 retPtr->setCtr(this->getCtr());
65 return retPtr;
66}
67
68std::shared_ptr<Kernel> FixedKernel::resized(int width, int height) const {
69 if ((width <= 0) || (height <= 0)) {
71 os << "Cannot create FixedKernel with dimensions (" << width << ", " << height << "). ";
73 }
74 if ((width - getWidth()) % 2 || (height - getHeight()) % 2) {
76 os << "Cannot resize FixedKernel from (" << getWidth() << ", " << getHeight() << ") to (" << width
77 << ", " << height << "), because at least one dimension would change by an odd value.";
79 }
80
81 lsst::geom::Box2I bboxNew(lsst::geom::Point2I((1 - width) / 2, (1 - height) / 2),
82 lsst::geom::Extent2I(width, height));
83 std::shared_ptr<image::Image<Pixel>> imNew = std::make_shared<image::Image<Pixel>>(bboxNew);
84
85 // getBBox() instantiates a new BBox from member data _width, _height, _ctrX, _ctrY
86 // so modifying it is OK
87 lsst::geom::Box2I bboxIntersect = getBBox();
88 bboxIntersect.clip(bboxNew);
89
90 // Kernel (this) and member image (this->_image) do not always have same XY0.
91 // Member image of resized kernel will not have same BBox as orig,
92 // but BBox of member image is ignored by the kernel.
93 int offsetX = _image.getX0() - getBBox().getMinX();
94 int offsetY = _image.getY0() - getBBox().getMinY();
95 lsst::geom::Box2I bboxIntersectShifted = lsst::geom::Box2I(
96 lsst::geom::Point2I(bboxIntersect.getMinX() + offsetX, bboxIntersect.getMinY() + offsetY),
97 bboxIntersect.getDimensions());
98 image::Image<Pixel> imIntersect = image::Image<Pixel>(_image, bboxIntersectShifted);
99
100 imNew->assign(imIntersect, bboxIntersect);
101 std::shared_ptr<Kernel> retPtr = std::make_shared<FixedKernel>(*imNew);
102 return retPtr;
103}
104
105double FixedKernel::doComputeImage(image::Image<Pixel>& image, bool doNormalize) const {
106 double multFactor = 1.0;
107 double imSum = this->_sum;
108 if (doNormalize) {
109 if (imSum == 0) {
110 throw LSST_EXCEPT(pexExcept::OverflowError, "Cannot normalize; kernel sum is 0");
111 }
112 multFactor = 1.0 / static_cast<double>(this->_sum);
113 imSum = 1.0;
114 }
115
117
118 for (int y = 0; y != this->getHeight(); ++y) {
119 for (XIter imPtr = image.row_begin(y), imEnd = image.row_end(y), kPtr = this->_image.row_begin(y);
120 imPtr != imEnd; ++imPtr, ++kPtr) {
121 imPtr[0] = multFactor * kPtr[0];
122 }
123 }
124
125 return imSum;
126}
127
130 os << prefix << "FixedKernel:" << std::endl;
131 os << prefix << "..sum: " << _sum << std::endl;
132 os << Kernel::toString(prefix + "\t");
133 return os.str();
134}
135
136// ------ Persistence ---------------------------------------------------------------------------------------
137
138namespace {
139
140struct FixedKernelPersistenceHelper : public Kernel::PersistenceHelper {
141 table::Key<table::Array<Kernel::Pixel>> image;
142
143 explicit FixedKernelPersistenceHelper(lsst::geom::Extent2I const& dimensions)
144 : Kernel::PersistenceHelper(0),
145 image(schema.addField<table::Array<Kernel::Pixel>>("image", "pixel values (row-major)",
146 dimensions.getX() * dimensions.getY())) {}
147
148 explicit FixedKernelPersistenceHelper(table::Schema const& schema_)
149 : Kernel::PersistenceHelper(schema_), image(schema["image"]) {}
150};
151
152} // namespace
153
155public:
157 CatalogVector const& catalogs) const override {
158 LSST_ARCHIVE_ASSERT(catalogs.size() == 1u);
159 LSST_ARCHIVE_ASSERT(catalogs.front().size() == 1u);
160 FixedKernelPersistenceHelper const keys(catalogs.front().getSchema());
161 afw::table::BaseRecord const& record = catalogs.front().front();
162 image::Image<Pixel> image(lsst::geom::Extent2I(record.get(keys.dimensions)));
163 ndarray::flatten<1>(ndarray::static_dimension_cast<2>(image.getArray())) = record[keys.image];
164 std::shared_ptr<FixedKernel> result = std::make_shared<FixedKernel>(image);
165 result->setCtr(record.get(keys.center));
166 return result;
167 }
168
169 explicit Factory(std::string const& name) : afw::table::io::PersistableFactory(name) {}
170};
171
172namespace {
173
174std::string getFixedKernelPersistenceName() { return "FixedKernel"; }
175
176FixedKernel::Factory registration(getFixedKernelPersistenceName());
177
178} // namespace
179
180std::string FixedKernel::getPersistenceName() const { return getFixedKernelPersistenceName(); }
181
183 FixedKernelPersistenceHelper const keys(getDimensions());
184 std::shared_ptr<afw::table::BaseRecord> record = keys.write(handle, *this);
185 (*record)[keys.image] = ndarray::flatten<1>(ndarray::copy(_image.getArray()));
186}
187} // namespace math
188} // namespace afw
189} // 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
int y
Definition SpanSet.cc:48
#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
int getX0() const
Return the image's column-origin.
Definition ImageBase.h:306
int getY0() const
Return the image's row-origin.
Definition ImageBase.h:314
x_iterator row_begin(int y) const
Return an x_iterator to the start of the y'th row.
Definition ImageBase.h:401
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.
Factory(std::string const &name)
std::string getPersistenceName() const override
Return the unique name used to persist this object and look up its factory.
std::shared_ptr< Kernel > clone() const override
Return a pointer to a deep copy of this kernel.
FixedKernel()
Construct an empty FixedKernel of size 0x0.
double doComputeImage(lsst::afw::image::Image< Pixel > &image, bool doNormalize) const override
Low-level version of computeImage.
std::shared_ptr< Kernel > resized(int width, int height) const override
Return a pointer to a clone with specified kernel dimensions.
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
lsst::geom::Extent2I const getDimensions() const
Return the Kernel's dimensions (width, height)
Definition Kernel.h:212
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
double computeImage(lsst::afw::image::Image< Pixel > &image, bool doNormalize, double x=0.0, double y=0.0) const
Compute an image (pixellized representation of the kernel) in place.
Definition Kernel.cc:76
lsst::geom::Box2I getBBox() const
return parent bounding box, with XY0 = -center
Definition Kernel.h:239
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.
static std::shared_ptr< T > dynamicCast(std::shared_ptr< Persistable > const &ptr)
Dynamically cast a shared_ptr.
A base class for factory classes used to reconstruct objects from records.
An integer coordinate rectangle.
Definition Box.h:55
int getMinY() const noexcept
Definition Box.h:158
int getMinX() const noexcept
Definition Box.h:157
Reports invalid arguments.
Definition Runtime.h:66
Reports when the result of an arithmetic operation is too large for the destination type.
Definition Runtime.h:124
T endl(T... args)