LSST Applications g063fba187b+66a50001ff,g0f08755f38+1a22dc2551,g1653933729+a8ce1bb630,g168dd56ebc+a8ce1bb630,g217e2c1bcf+12e87a5bd8,g246886dfd9+466c7b6c06,g28da252d5a+858b171e04,g2bbee38e9b+404b60ec9b,g2bc492864f+404b60ec9b,g3156d2b45e+6e55a43351,g347aa1857d+404b60ec9b,g35bb328faa+a8ce1bb630,g3a166c0a6a+404b60ec9b,g3e281a1b8c+c5dd892a6c,g414038480c+6b9177ef31,g41af890bb2+70bea58702,g599934f4f4+b8c5400ca5,g781aacb6e4+a8ce1bb630,g7af13505b9+b5b9cefdb8,g80478fca09+c2997882f3,g82479be7b0+8974e6af0f,g858d7b2824+1a22dc2551,g89c8672015+f4add4ffd5,g8f1c07a47a+de51c9b0a5,g9125e01d80+a8ce1bb630,ga5288a1d22+b66f8cf76b,gb58c049af0+d64f4d3760,gc28159a63d+404b60ec9b,gcab2d0539d+66cf1de5d4,gcf0d15dbbd+12cb7e2563,gda6a2b7d83+12cb7e2563,gdaeeff99f8+1711a396fd,ge79ae78c31+404b60ec9b,gef2f8181fd+414189b318,gf0baf85859+c1f95f4921,gf0c06eb49c+1a22dc2551,gfa517265be+1a22dc2551,gfa999e8aa5+17cd334064,v28.0.0.rc2
LSST Data Management Base Package
Loading...
Searching...
No Matches
AnalyticKernel.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
26#include "lsst/pex/exceptions.h"
30
32
33namespace lsst {
34namespace afw {
35
36template std::shared_ptr<math::AnalyticKernel> table::io::PersistableFacade<
37 math::AnalyticKernel>::dynamicCast(std::shared_ptr<table::io::Persistable> const &);
38
39namespace math {
40
41AnalyticKernel::AnalyticKernel() : Kernel(), _kernelFunctionPtr() {}
42
44 Kernel::SpatialFunction const &spatialFunction)
45 : Kernel(width, height, kernelFunction.getNParameters(), spatialFunction),
46 _kernelFunctionPtr(kernelFunction.clone()) {}
47
49 std::vector<Kernel::SpatialFunctionPtr> const &spatialFunctionList)
50 : Kernel(width, height, spatialFunctionList), _kernelFunctionPtr(kernelFunction.clone()) {
51 if (kernelFunction.getNParameters() != spatialFunctionList.size()) {
53 os << "kernelFunction.getNParameters() = " << kernelFunction.getNParameters()
54 << " != " << spatialFunctionList.size() << " = "
55 << "spatialFunctionList.size()";
57 }
58}
59
62 if (this->isSpatiallyVarying()) {
63 retPtr.reset(new AnalyticKernel(this->getWidth(), this->getHeight(), *(this->_kernelFunctionPtr),
65 } else {
66 retPtr.reset(new AnalyticKernel(this->getWidth(), this->getHeight(), *(this->_kernelFunctionPtr)));
67 }
68 retPtr->setCtr(this->getCtr());
69 return retPtr;
70}
71
74 if (isSpatiallyVarying()) {
75 retPtr = std::make_shared<AnalyticKernel>(width, height, *_kernelFunctionPtr, _spatialFunctionList);
76 } else {
77 retPtr = std::make_shared<AnalyticKernel>(width, height, *_kernelFunctionPtr);
78 }
79 return retPtr;
80}
81
82double AnalyticKernel::computeImage(image::Image<Pixel> &image, bool doNormalize, double x, double y) const {
83 lsst::geom::Extent2I llBorder = (image.getDimensions() - getDimensions()) / 2;
85 if (this->isSpatiallyVarying()) {
87 }
88 return doComputeImage(image, doNormalize);
89}
90
94
97 os << prefix << "AnalyticKernel:" << std::endl;
98 os << prefix << "..function: " << (_kernelFunctionPtr ? _kernelFunctionPtr->toString() : "None")
99 << std::endl;
100 os << Kernel::toString(prefix + "\t");
101 return os.str();
102}
103
107
108//
109// Protected Member Functions
110//
112 double imSum = 0;
113 for (int y = 0; y != image.getHeight(); ++y) {
114 double const fy = image.indexToPosition(y, image::Y);
116 for (int x = 0; x != image.getWidth(); ++x, ++ptr) {
117 double const fx = image.indexToPosition(x, image::X);
118 Pixel const pixelVal = (*_kernelFunctionPtr)(fx, fy);
119 *ptr = pixelVal;
120 imSum += pixelVal;
121 }
122 }
123
124 if (doNormalize && (imSum != 1)) {
125 if (imSum == 0) {
126 throw LSST_EXCEPT(pexExcept::OverflowError, "Cannot normalize; kernel sum is 0");
127 }
128 image /= imSum;
129 imSum = 1;
130 }
131
132 return imSum;
133}
134
135void AnalyticKernel::setKernelParameter(unsigned int ind, double value) const {
136 _kernelFunctionPtr->setParameter(ind, value);
137}
138
139// ------ Persistence ---------------------------------------------------------------------------------------
140
141namespace {
142
143struct AnalyticKernelPersistenceHelper : public Kernel::PersistenceHelper {
144 table::Key<int> kernelFunction;
145
146 explicit AnalyticKernelPersistenceHelper(int nSpatialFunctions)
147 : Kernel::PersistenceHelper(nSpatialFunctions),
148 kernelFunction(schema.addField<int>(
149 "kernelfunction", "archive ID for analytic function used to produce kernel images")) {}
150
151 explicit AnalyticKernelPersistenceHelper(table::Schema const &schema_)
152 : Kernel::PersistenceHelper(schema_), kernelFunction(schema["kernelfunction"]) {}
153};
154
155} // namespace
156
158public:
160 CatalogVector const &catalogs) const override {
161 LSST_ARCHIVE_ASSERT(catalogs.size() == 1u);
162 LSST_ARCHIVE_ASSERT(catalogs.front().size() == 1u);
163 AnalyticKernelPersistenceHelper const keys(catalogs.front().getSchema());
164 afw::table::BaseRecord const &record = catalogs.front().front();
166 archive.get<AnalyticKernel::KernelFunction>(record.get(keys.kernelFunction));
168 if (keys.spatialFunctions.isValid()) {
169 result = std::make_shared<AnalyticKernel>(record.get(keys.dimensions.getX()),
170 record.get(keys.dimensions.getY()), *kernelFunction,
171 keys.readSpatialFunctions(archive, record));
172 } else {
173 result = std::make_shared<AnalyticKernel>(record.get(keys.dimensions.getX()),
174 record.get(keys.dimensions.getY()), *kernelFunction);
175 }
176 result->setCtr(record.get(keys.center));
177 return result;
178 }
179
180 explicit Factory(std::string const &name) : afw::table::io::PersistableFactory(name) {}
181};
182
183namespace {
184
185std::string getAnalyticKernelPersistenceName() { return "AnalyticKernel"; }
186
187AnalyticKernel::Factory registration(getAnalyticKernelPersistenceName());
188
189} // namespace
190
191std::string AnalyticKernel::getPersistenceName() const { return getAnalyticKernelPersistenceName(); }
192
194 AnalyticKernelPersistenceHelper const keys(_spatialFunctionList.size());
195 std::shared_ptr<afw::table::BaseRecord> record = keys.write(handle, *this);
196 record->set(keys.kernelFunction, handle.put(_kernelFunctionPtr.get()));
197}
198} // namespace math
199} // namespace afw
200} // namespace lsst
py::object result
Definition _schema.cc:429
table::Key< int > kernelFunction
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition Exception.h:48
std::uint64_t * ptr
Definition RangeSet.cc:95
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
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.
std::vector< double > getKernelParameters() const override
Return the current kernel parameters.
std::string toString(std::string const &prefix="") const override
Return a string representation of the kernel.
std::shared_ptr< Kernel > resized(int width, int height) const override
Return a pointer to a clone with specified kernel dimensions.
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.
std::shared_ptr< Kernel > clone() const override
Return a pointer to a deep copy of this kernel.
void setKernelParameter(unsigned int ind, double value) const override
Set one kernel parameter.
virtual KernelFunctionPtr getKernelFunction() const
Get a deep copy of the kernel function.
double doComputeImage(lsst::afw::image::Image< Pixel > &image, bool doNormalize) const override
Low-level version of computeImage.
void write(OutputArchiveHandle &handle) const override
Write the object to one or more catalogs.
KernelFunctionPtr _kernelFunctionPtr
Definition Kernel.h:633
std::string getPersistenceName() const override
Return the unique name used to persist this object and look up its factory.
AnalyticKernel()
Construct an empty spatially invariant AnalyticKernel of size 0x0.
virtual std::shared_ptr< Function2< ReturnT > > clone() const =0
Return a pointer to a deep copy of this function.
std::string toString(std::string const &prefix="") const override
Return a string representation of the function.
Definition Function.h:299
void setParameter(unsigned int ind, double newValue)
Set one function parameter without range checking.
Definition Function.h:143
std::vector< double > const & getParameters() const noexcept
Return all function parameters.
Definition Function.h:129
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
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition Kernel.h:449
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
bool isSpatiallyVarying() const
Return true iff the kernel is spatially varying (has a spatial function)
Definition Kernel.h:333
void setKernelParametersFromSpatialModel(double x, double y) const
Set the kernel parameters from the spatial model (if any).
Definition Kernel.cc:217
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.
std::shared_ptr< Persistable > get(int id) const
Load the Persistable with the given ID and return it.
An object passed to Persistable::write to allow it to persist itself.
int put(Persistable const *obj, bool permissive=false)
Save an object to the archive and return a unique ID that can be used to retrieve it from an InputArc...
A base class for factory classes used to reconstruct objects from records.
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)
T get(T... args)
T reset(T... args)
T size(T... args)