LSST Applications g063fba187b+fee0456c91,g0f08755f38+ea96e5a5a3,g1653933729+a8ce1bb630,g168dd56ebc+a8ce1bb630,g1a2382251a+90257ff92a,g20f6ffc8e0+ea96e5a5a3,g217e2c1bcf+937a289c59,g28da252d5a+daa7da44eb,g2bbee38e9b+253935c60e,g2bc492864f+253935c60e,g3156d2b45e+6e55a43351,g32e5bea42b+31359a2a7a,g347aa1857d+253935c60e,g35bb328faa+a8ce1bb630,g3a166c0a6a+253935c60e,g3b1af351f3+a8ce1bb630,g3e281a1b8c+c5dd892a6c,g414038480c+416496e02f,g41af890bb2+afe91b1188,g599934f4f4+0db33f7991,g7af13505b9+e36de7bce6,g80478fca09+da231ba887,g82479be7b0+a4516e59e3,g858d7b2824+ea96e5a5a3,g89c8672015+f4add4ffd5,g9125e01d80+a8ce1bb630,ga5288a1d22+bc6ab8dfbd,gb58c049af0+d64f4d3760,gc28159a63d+253935c60e,gcab2d0539d+3f2b72788c,gcf0d15dbbd+4ea9c45075,gda6a2b7d83+4ea9c45075,gdaeeff99f8+1711a396fd,ge79ae78c31+253935c60e,gef2f8181fd+3031e3cf99,gf0baf85859+c1f95f4921,gfa517265be+ea96e5a5a3,gfa999e8aa5+17cd334064,w.2024.50
LSST Data Management Base Package
|
A kernel that is a linear combination of fixed basis kernels. More...
#include <Kernel.h>
Classes | |
class | Factory |
Public Types | |
using | Pixel = double |
using | SpatialFunctionPtr = std::shared_ptr<lsst::afw::math::Function2<double>> |
using | SpatialFunction = lsst::afw::math::Function2<double> |
using | NullSpatialFunction = lsst::afw::math::NullFunction2<double> |
using | kernel_fill_factor = generic_kernel_tag |
Public Member Functions | |
LinearCombinationKernel () | |
Construct an empty LinearCombinationKernel of size 0x0. | |
LinearCombinationKernel (KernelList const &kernelList, std::vector< double > const &kernelParameters) | |
Construct a spatially invariant LinearCombinationKernel. | |
LinearCombinationKernel (KernelList const &kernelList, Kernel::SpatialFunction const &spatialFunction) | |
Construct a spatially varying LinearCombinationKernel, where the spatial model is described by one function (that is cloned to give one per basis kernel). | |
LinearCombinationKernel (KernelList const &kernelList, std::vector< Kernel::SpatialFunctionPtr > const &spatialFunctionList) | |
Construct a spatially varying LinearCombinationKernel, where the spatial model is described by a list of functions (one per basis kernel). | |
LinearCombinationKernel (const LinearCombinationKernel &)=delete | |
LinearCombinationKernel (LinearCombinationKernel &&)=delete | |
LinearCombinationKernel & | operator= (const LinearCombinationKernel &)=delete |
LinearCombinationKernel & | operator= (LinearCombinationKernel &&)=delete |
~LinearCombinationKernel () override=default | |
std::shared_ptr< Kernel > | clone () const override |
Return a pointer to a deep copy of this kernel. | |
std::shared_ptr< Kernel > | resized (int width, int height) const override |
Return a pointer to a clone with specified kernel dimensions. | |
std::vector< double > | getKernelParameters () const override |
Return the current kernel parameters. | |
virtual KernelList const & | getKernelList () const |
Get the fixed basis kernels. | |
std::vector< double > | getKernelSumList () const |
Get the sum of the pixels of each fixed basis kernel. | |
int | getNBasisKernels () const |
Get the number of basis kernels. | |
void | checkKernelList (const KernelList &kernelList) const |
Check that all kernels have the same size and center and that none are spatially varying. | |
bool | isDeltaFunctionBasis () const |
Return true if all basis kernels are instances of DeltaFunctionKernel. | |
std::shared_ptr< Kernel > | refactor () const |
Refactor the kernel as a linear combination of N bases where N is the number of parameters for the spatial model. | |
std::string | toString (std::string const &prefix="") const override |
Return a string representation of the kernel. | |
bool | isPersistable () const noexcept override |
Return true if this particular object can be persisted using afw::table::io. | |
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. | |
lsst::geom::Extent2I const | getDimensions () const |
Return the Kernel's dimensions (width, height) | |
void | setDimensions (lsst::geom::Extent2I dims) |
void | setWidth (int width) |
void | setHeight (int height) |
int | getWidth () const |
Return the Kernel's width. | |
int | getHeight () const |
Return the Kernel's height. | |
lsst::geom::Point2I | getCtr () const |
Return index of kernel's center. | |
lsst::geom::Box2I | getBBox () const |
return parent bounding box, with XY0 = -center | |
unsigned int | getNKernelParameters () const |
Return the number of kernel parameters (0 if none) | |
int | getNSpatialParameters () const |
Return the number of spatial parameters (0 if not spatially varying) | |
SpatialFunctionPtr | getSpatialFunction (unsigned int index) const |
Return a clone of the specified spatial function (one component of the spatial model) | |
std::vector< SpatialFunctionPtr > | getSpatialFunctionList () const |
Return a list of clones of the spatial functions. | |
virtual double | getKernelParameter (unsigned int i) const |
Return a particular Kernel Parameter (no bounds checking). | |
lsst::geom::Box2I | growBBox (lsst::geom::Box2I const &bbox) const |
Given a bounding box for pixels one wishes to compute by convolving an image with this kernel, return the bounding box of pixels that must be accessed on the image to be convolved. | |
lsst::geom::Box2I | shrinkBBox (lsst::geom::Box2I const &bbox) const |
Given a bounding box for an image one wishes to convolve with this kernel, return the bounding box for the region of pixels that can be computed. | |
void | setCtr (lsst::geom::Point2I ctr) |
Set index of kernel's center. | |
std::vector< std::vector< double > > | getSpatialParameters () const |
Return the spatial parameters parameters (an empty vector if not spatially varying) | |
bool | isSpatiallyVarying () const |
Return true iff the kernel is spatially varying (has a spatial function) | |
void | setKernelParameters (std::vector< double > const ¶ms) |
Set the kernel parameters of a spatially invariant kernel. | |
void | setKernelParameters (std::pair< double, double > const ¶ms) |
Set the kernel parameters of a 2-component spatially invariant kernel. | |
void | setSpatialParameters (const std::vector< std::vector< double > > params) |
Set the parameters of all spatial functions. | |
void | computeKernelParametersFromSpatialModel (std::vector< double > &kernelParams, double x, double y) const |
Compute the kernel parameters at a specified point. | |
virtual void | computeCache (int const) |
Compute a cache of Kernel values, if desired. | |
virtual int | getCacheSize () const |
Get the current size of the kernel cache (0 if none or if caches not supported) | |
void | writeFits (std::string const &fileName, std::string const &mode="w") const |
Write the object to a regular FITS file. | |
void | writeFits (fits::MemFileManager &manager, std::string const &mode="w") const |
Write the object to a FITS image in memory. | |
void | writeFits (fits::Fits &fitsfile) const |
Write the object to an already-open FITS object. | |
Static Public Member Functions | |
static std::shared_ptr< LinearCombinationKernel > | readFits (fits::Fits &fitsfile) |
Read an object from an already open FITS object. | |
static std::shared_ptr< LinearCombinationKernel > | readFits (std::string const &fileName, int hdu=fits::DEFAULT_HDU) |
Read an object from a regular FITS file. | |
static std::shared_ptr< LinearCombinationKernel > | readFits (fits::MemFileManager &manager, int hdu=fits::DEFAULT_HDU) |
Read an object from a FITS file in memory. | |
static std::shared_ptr< LinearCombinationKernel > | dynamicCast (std::shared_ptr< Persistable > const &ptr) |
Dynamically cast a shared_ptr. | |
static std::shared_ptr< Kernel > | readFits (fits::Fits &fitsfile) |
Read an object from an already open FITS object. | |
static std::shared_ptr< Kernel > | readFits (std::string const &fileName, int hdu=fits::DEFAULT_HDU) |
Read an object from a regular FITS file. | |
static std::shared_ptr< Kernel > | readFits (fits::MemFileManager &manager, int hdu=fits::DEFAULT_HDU) |
Read an object from a FITS file in memory. | |
static std::shared_ptr< Kernel > | dynamicCast (std::shared_ptr< Persistable > const &ptr) |
Dynamically cast a shared_ptr. | |
Protected Types | |
using | OutputArchiveHandle = io::OutputArchiveHandle |
Protected Member Functions | |
double | doComputeImage (lsst::afw::image::Image< Pixel > &image, bool doNormalize) const override |
Low-level version of computeImage. | |
std::string | getPersistenceName () const override |
Return the unique name used to persist this object and look up its factory. | |
void | write (OutputArchiveHandle &handle) const override |
Write the object to one or more catalogs. | |
void | setKernelParameter (unsigned int ind, double value) const override |
Set one kernel parameter. | |
std::string | getPythonModule () const override |
Return the fully-qualified Python module that should be imported to guarantee that its factory is registered. | |
void | setKernelParametersFromSpatialModel (double x, double y) const |
Set the kernel parameters from the spatial model (if any). | |
Protected Attributes | |
std::vector< SpatialFunctionPtr > | _spatialFunctionList |
A kernel that is a linear combination of fixed basis kernels.
Convolution may be performed by first convolving the image with each fixed kernel, then adding the resulting images using the (possibly spatially varying) kernel coefficients.
The basis kernels are cloned (deep copied) so you may safely modify your own copies.
Warnings:
|
inherited |
|
inherited |
|
protectedinherited |
Definition at line 108 of file Persistable.h.
|
inherited |
|
inherited |
|
inherited |
|
explicit |
Construct an empty LinearCombinationKernel of size 0x0.
Definition at line 48 of file LinearCombinationKernel.cc.
|
explicit |
Construct a spatially invariant LinearCombinationKernel.
kernelList | list of (shared pointers to const) basis kernels |
kernelParameters | kernel coefficients |
Definition at line 56 of file LinearCombinationKernel.cc.
|
explicit |
Construct a spatially varying LinearCombinationKernel, where the spatial model is described by one function (that is cloned to give one per basis kernel).
kernelList | list of (shared pointers to const) basis kernels |
spatialFunction | spatial function; one deep copy is made for each basis kernel |
Definition at line 74 of file LinearCombinationKernel.cc.
|
explicit |
Construct a spatially varying LinearCombinationKernel, where the spatial model is described by a list of functions (one per basis kernel).
kernelList | list of (shared pointers to const) kernels |
spatialFunctionList | list of spatial functions, one per basis kernel |
lsst::pex::exceptions::InvalidParameterError | if the length of spatialFunctionList != # kernels |
Definition at line 86 of file LinearCombinationKernel.cc.
|
delete |
|
delete |
|
overridedefault |
void lsst::afw::math::LinearCombinationKernel::checkKernelList | ( | const KernelList & | kernelList | ) | const |
Check that all kernels have the same size and center and that none are spatially varying.
lsst::pex::exceptions::InvalidParameterError | if the check fails |
Definition at line 132 of file LinearCombinationKernel.cc.
|
overridevirtual |
Return a pointer to a deep copy of this kernel.
This kernel exists instead of a copy constructor so one can obtain a copy of an actual kernel instead of a useless copy of the base class.
Every kernel subclass must override this method.
Implements lsst::afw::math::Kernel.
Definition at line 104 of file LinearCombinationKernel.cc.
|
inlinevirtualinherited |
Compute a cache of Kernel values, if desired.
Reimplemented in lsst::afw::math::SeparableKernel.
Definition at line 397 of file Kernel.h.
|
inherited |
Compute an image (pixellized representation of the kernel) in place.
image | image whose pixels are to be set (output); xy0 of the image will be set to -kernel.getCtr() |
doNormalize | normalize the image (so sum is 1)? |
x | x (column position) at which to compute spatial function |
y | y (row position) at which to compute spatial function |
lsst::pex::exceptions::InvalidParameterError | if the image is the wrong size |
lsst::pex::exceptions::OverflowError | if doNormalize is true and the kernel sum is exactly 0 |
Definition at line 76 of file Kernel.cc.
|
inherited |
Compute the kernel parameters at a specified point.
Warning: this is a low-level function that assumes kernelParams is the right length. It will fail in unpredictable ways if that condition is not met.
Definition at line 135 of file Kernel.cc.
|
overrideprotectedvirtual |
Low-level version of computeImage.
Before this is called the image dimensions are checked, the image's xy0 is set and the kernel's parameters are set. This routine sets the pixels, including normalization if requested.
image | image whose pixels are to be set (output) |
doNormalize | normalize the image (so sum is 1)? |
Implements lsst::afw::math::Kernel.
Definition at line 247 of file LinearCombinationKernel.cc.
|
staticinherited |
Dynamically cast a shared_ptr.
Dynamically cast a shared pointer and raise on failure.
You must provide an explicit template instantiation in the .cc file for each class that inherits from PersistableFacade. Designed to work around RTTI issues on macOS with hidden symbols;
lsst::pex::exceptions::LogicError | if the cast fails |
param[in] ptr The pointer to be cast.
lsst::pex::exceptions::TypeError | If the dynamic cast fails. |
Definition at line 218 of file Persistable.cc.
|
staticinherited |
Dynamically cast a shared_ptr.
Dynamically cast a shared pointer and raise on failure.
You must provide an explicit template instantiation in the .cc file for each class that inherits from PersistableFacade. Designed to work around RTTI issues on macOS with hidden symbols;
lsst::pex::exceptions::LogicError | if the cast fails |
param[in] ptr The pointer to be cast.
lsst::pex::exceptions::TypeError | If the dynamic cast fails. |
Definition at line 218 of file Persistable.cc.
|
inlineinherited |
return parent bounding box, with XY0 = -center
Definition at line 239 of file Kernel.h.
|
inlinevirtualinherited |
Get the current size of the kernel cache (0 if none or if caches not supported)
Reimplemented in lsst::afw::math::SeparableKernel.
Definition at line 403 of file Kernel.h.
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
virtual |
Get the fixed basis kernels.
Definition at line 156 of file LinearCombinationKernel.cc.
|
inlinevirtualinherited |
Return a particular Kernel Parameter (no bounds checking).
This version is slow, but specialisations may be faster
Reimplemented in lsst::afw::math::SeparableKernel.
Definition at line 277 of file Kernel.h.
|
overridevirtual |
Return the current kernel parameters.
If the kernel is spatially varying then the parameters are those last computed. See also computeKernelParametersFromSpatialModel. If there are no kernel parameters then returns an empty vector.
Reimplemented from lsst::afw::math::Kernel.
Definition at line 160 of file LinearCombinationKernel.cc.
std::vector< double > lsst::afw::math::LinearCombinationKernel::getKernelSumList | ( | ) | const |
Get the sum of the pixels of each fixed basis kernel.
Definition at line 158 of file LinearCombinationKernel.cc.
|
inline |
|
inlineinherited |
|
inlineinherited |
Return the number of spatial parameters (0 if not spatially varying)
|
overrideprotectedvirtual |
Return the unique name used to persist this object and look up its factory.
Must be less than ArchiveIndexSchema::MAX_NAME_LENGTH characters.
Reimplemented from lsst::afw::table::io::Persistable.
Definition at line 365 of file LinearCombinationKernel.cc.
|
overrideprotectedvirtualinherited |
Return the fully-qualified Python module that should be imported to guarantee that its factory is registered.
Must be less than ArchiveIndexSchema::MAX_MODULE_LENGTH characters.
Will be ignored if empty.
Reimplemented from lsst::afw::table::io::Persistable.
Reimplemented in lsst::afw::math::LanczosWarpingKernel, and lsst::afw::math::NearestWarpingKernel.
Definition at line 224 of file Kernel.cc.
|
inherited |
Return a clone of the specified spatial function (one component of the spatial model)
index | index of desired spatial function; must be in range [0, number spatial parameters - 1] |
lsst::pex::exceptions::InvalidParameterError | if kernel not spatially varying |
lsst::pex::exceptions::InvalidParameterError | if index out of range |
|
inherited |
Return a list of clones of the spatial functions.
Definition at line 157 of file Kernel.cc.
|
inlineinherited |
Return the spatial parameters parameters (an empty vector if not spatially varying)
Definition at line 321 of file Kernel.h.
|
inlineinherited |
|
inherited |
Given a bounding box for pixels one wishes to compute by convolving an image with this kernel, return the bounding box of pixels that must be accessed on the image to be convolved.
Thus the box shifted by -kernel.getCtr() and its size is expanded by kernel.getDimensions()-1.
Definition at line 167 of file Kernel.cc.
|
inline |
Return true if all basis kernels are instances of DeltaFunctionKernel.
Definition at line 780 of file Kernel.h.
|
inlineoverridevirtualnoexcept |
Return true if this particular object can be persisted using afw::table::io.
Reimplemented from lsst::afw::table::io::Persistable.
Definition at line 820 of file Kernel.h.
|
inlineinherited |
|
delete |
|
delete |
|
inlinestaticinherited |
Read an object from an already open FITS object.
[in] | fitsfile | FITS object to read from, already positioned at the desired HDU. |
Definition at line 183 of file Persistable.h.
|
inlinestaticinherited |
Read an object from an already open FITS object.
[in] | fitsfile | FITS object to read from, already positioned at the desired HDU. |
Definition at line 183 of file Persistable.h.
|
inlinestaticinherited |
Read an object from a FITS file in memory.
[in] | manager | Manager for the memory to read from. |
[in] | hdu | HDU to read, where 0 is the primary. The special value of afw::fits::DEFAULT_HDU skips the primary HDU if it is empty. |
Definition at line 205 of file Persistable.h.
|
inlinestaticinherited |
Read an object from a FITS file in memory.
[in] | manager | Manager for the memory to read from. |
[in] | hdu | HDU to read, where 0 is the primary. The special value of afw::fits::DEFAULT_HDU skips the primary HDU if it is empty. |
Definition at line 205 of file Persistable.h.
|
inlinestaticinherited |
Read an object from a regular FITS file.
[in] | fileName | Name of the file to read. |
[in] | hdu | HDU to read, where 0 is the primary. The special value of afw::fits::DEFAULT_HDU skips the primary HDU if it is empty. |
Definition at line 194 of file Persistable.h.
|
inlinestaticinherited |
Read an object from a regular FITS file.
[in] | fileName | Name of the file to read. |
[in] | hdu | HDU to read, where 0 is the primary. The special value of afw::fits::DEFAULT_HDU skips the primary HDU if it is empty. |
Definition at line 194 of file Persistable.h.
std::shared_ptr< Kernel > lsst::afw::math::LinearCombinationKernel::refactor | ( | ) | const |
Refactor the kernel as a linear combination of N bases where N is the number of parameters for the spatial model.
Refactoring is only possible if all of the following are true:
Details: A spatially varying LinearCombinationKernel consisting of M basis kernels and using a spatial model that is a linear combination of N coefficients can be expressed as: K(x,y) = K0 (C00 F0(x,y) + C10 F1(x,y) + C20 F2(x,y) + ... + CN0 FN(x,y))
This is equivalent to the following linear combination of N basis kernels:
= K0' F0(x,y) + K1' F1(x,y) + K2' F2(x,y) + ... + KN' FN(x,y) where Ki' = sum over j of Kj Cij
This is what refactor returns provided the required conditions are met. However, the spatial functions for the refactored kernel are the same as those for the original kernel (for generality and simplicity) with all coefficients equal to 0 except one that is set to 1; hence they are not computed optimally.
Thanks to Kresimir Cosic for inventing or reinventing this useful technique.
Definition at line 162 of file LinearCombinationKernel.cc.
|
overridevirtual |
Return a pointer to a clone with specified kernel dimensions.
width | Number of columns in pixels |
height | Number of rows in pixels |
Must be implemented by derived classes.
Implements lsst::afw::math::Kernel.
Definition at line 115 of file LinearCombinationKernel.cc.
|
inlineinherited |
Set index of kernel's center.
|
inlineinherited |
|
inlineinherited |
|
overrideprotectedvirtual |
Set one kernel parameter.
Classes that have kernel parameters must subclass this function.
This function is marked "const", despite modifying unimportant internals, so that computeImage can be const.
lsst::pex::exceptions::InvalidParameterError | always (unless subclassed) |
Reimplemented from lsst::afw::math::Kernel.
Definition at line 270 of file LinearCombinationKernel.cc.
|
inlineinherited |
Set the kernel parameters of a 2-component spatially invariant kernel.
Definition at line 363 of file Kernel.h.
|
inlineinherited |
Set the kernel parameters of a spatially invariant kernel.
lsst::pex::exceptions::RuntimeError | if the kernel has a spatial function |
lsst::pex::exceptions::InvalidParameterError | if the params vector is the wrong length |
Definition at line 341 of file Kernel.h.
|
protectedinherited |
Set the kernel parameters from the spatial model (if any).
This function has no effect if there is no spatial model.
This function is marked "const", despite modifying unimportant internals, so that computeImage can be const.
Definition at line 217 of file Kernel.cc.
|
inherited |
Set the parameters of all spatial functions.
Params is indexed as [kernel parameter][spatial parameter]
lsst::pex::exceptions::InvalidParameterError | if params is the wrong shape (if this exception is thrown then no parameters are changed) |
Definition at line 110 of file Kernel.cc.
|
inlineinherited |
|
inherited |
Given a bounding box for an image one wishes to convolve with this kernel, return the bounding box for the region of pixels that can be computed.
Thus the box shifted by kernel.getCtr() and its size is reduced by kernel.getDimensions()-1.
lsst::pex::exceptions::InvalidParameterError | if the resulting box would have dimension < 1 in either axis |
|
overridevirtual |
Return a string representation of the kernel.
Reimplemented from lsst::afw::math::Kernel.
Definition at line 227 of file LinearCombinationKernel.cc.
|
overrideprotectedvirtual |
Write the object to one or more catalogs.
The handle object passed to this function provides an interface for adding new catalogs and adding nested objects to the same archive (while checking for duplicates). See OutputArchiveHandle for more information.
Reimplemented from lsst::afw::table::io::Persistable.
Definition at line 369 of file LinearCombinationKernel.cc.
|
inherited |
Write the object to an already-open FITS object.
[in] | fitsfile | Open FITS object to write to. |
Definition at line 18 of file Persistable.cc.
|
inherited |
Write the object to a FITS image in memory.
[in] | manager | Name of the file to write to. |
[in] | mode | If "w", any existing file with the given name will be overwritten. If "a", new HDUs will be appended to an existing file. |
Definition at line 29 of file Persistable.cc.
|
inherited |
Write the object to a regular FITS file.
[in] | fileName | Name of the file to write to. |
[in] | mode | If "w", any existing file with the given name will be overwritten. If "a", new HDUs will be appended to an existing file. |
Definition at line 24 of file Persistable.cc.
|
protectedinherited |