LSST Applications  21.0.0+04719a4bac,21.0.0-1-ga51b5d4+f5e6047307,21.0.0-11-g2b59f77+a9c1acf22d,21.0.0-11-ga42c5b2+86977b0b17,21.0.0-12-gf4ce030+76814010d2,21.0.0-13-g1721dae+760e7a6536,21.0.0-13-g3a573fe+768d78a30a,21.0.0-15-g5a7caf0+f21cbc5713,21.0.0-16-g0fb55c1+b60e2d390c,21.0.0-19-g4cded4ca+71a93a33c0,21.0.0-2-g103fe59+bb20972958,21.0.0-2-g45278ab+04719a4bac,21.0.0-2-g5242d73+3ad5d60fb1,21.0.0-2-g7f82c8f+8babb168e8,21.0.0-2-g8f08a60+06509c8b61,21.0.0-2-g8faa9b5+616205b9df,21.0.0-2-ga326454+8babb168e8,21.0.0-2-gde069b7+5e4aea9c2f,21.0.0-2-gecfae73+1d3a86e577,21.0.0-2-gfc62afb+3ad5d60fb1,21.0.0-25-g1d57be3cd+e73869a214,21.0.0-3-g357aad2+ed88757d29,21.0.0-3-g4a4ce7f+3ad5d60fb1,21.0.0-3-g4be5c26+3ad5d60fb1,21.0.0-3-g65f322c+e0b24896a3,21.0.0-3-g7d9da8d+616205b9df,21.0.0-3-ge02ed75+a9c1acf22d,21.0.0-4-g591bb35+a9c1acf22d,21.0.0-4-g65b4814+b60e2d390c,21.0.0-4-gccdca77+0de219a2bc,21.0.0-4-ge8a399c+6c55c39e83,21.0.0-5-gd00fb1e+05fce91b99,21.0.0-6-gc675373+3ad5d60fb1,21.0.0-64-g1122c245+4fb2b8f86e,21.0.0-7-g04766d7+cd19d05db2,21.0.0-7-gdf92d54+04719a4bac,21.0.0-8-g5674e7b+d1bd76f71f,master-gac4afde19b+a9c1acf22d,w.2021.13
LSST Data Management Base Package
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"
28 #include "lsst/afw/math/Kernel.h"
31 
33 
34 namespace lsst {
35 namespace afw {
36 
37 template std::shared_ptr<math::DeltaFunctionKernel> table::io::PersistableFacade<
38  math::DeltaFunctionKernel>::dynamicCast(std::shared_ptr<table::io::Persistable> const&);
39 
40 namespace 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 
98 namespace {
99 
100 struct 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 
116 private:
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 
125 public:
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 
143 namespace {
144 
145 std::string getDeltaFunctionKernelPersistenceName() { return "DeltaFunctionKernel"; }
146 
147 DeltaFunctionKernel::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:430
table::Key< std::string > name
Definition: Amplifier.cc:116
table::PointKey< int > pixel
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
std::ostream * os
Definition: Schema.cc:746
std::string prefix
Definition: SchemaMapper.cc:79
#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
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:671
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:111
int getHeight() const
Return the Kernel's height.
Definition: Kernel.h:230
lsst::geom::Point2I getCtr() const
Return index of kernel's center.
Definition: Kernel.h:235
virtual std::string toString(std::string const &prefix="") const
Return a string representation of the kernel.
Definition: Kernel.cc:195
int getWidth() const
Return the Kernel's width.
Definition: Kernel.h:225
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.
Definition: CatalogVector.h:29
A multi-catalog archive object used to load table::io::Persistable objects.
Definition: InputArchive.h:31
An object passed to Persistable::write to allow it to persist itself.
A base class for factory classes used to reconstruct objects from records.
Definition: Persistable.h:228
PersistableFactory(std::string const &name)
Constructor for the factory.
Definition: Persistable.cc:74
Reports invalid arguments.
Definition: Runtime.h:66
T endl(T... args)
T front(T... args)
Backwards-compatibility support for depersisting the old Calib (FluxMag0/FluxMag0Err) objects.
FilterProperty & operator=(FilterProperty const &)=default
A base class for image defects.
T size(T... args)