LSSTApplications  18.0.0+106,18.0.0+50,19.0.0,19.0.0+1,19.0.0+10,19.0.0+11,19.0.0+13,19.0.0+17,19.0.0+2,19.0.0-1-g20d9b18+6,19.0.0-1-g425ff20,19.0.0-1-g5549ca4,19.0.0-1-g580fafe+6,19.0.0-1-g6fe20d0+1,19.0.0-1-g7011481+9,19.0.0-1-g8c57eb9+6,19.0.0-1-gb5175dc+11,19.0.0-1-gdc0e4a7+9,19.0.0-1-ge272bc4+6,19.0.0-1-ge3aa853,19.0.0-10-g448f008b,19.0.0-12-g6990b2c,19.0.0-2-g0d9f9cd+11,19.0.0-2-g3d9e4fb2+11,19.0.0-2-g5037de4,19.0.0-2-gb96a1c4+3,19.0.0-2-gd955cfd+15,19.0.0-3-g2d13df8,19.0.0-3-g6f3c7dc,19.0.0-4-g725f80e+11,19.0.0-4-ga671dab3b+1,19.0.0-4-gad373c5+3,19.0.0-5-ga2acb9c+2,19.0.0-5-gfe96e6c+2,w.2020.01
LSSTDataManagementBasePackage
Public Types | Public Member Functions | Static Public Member Functions | Protected Types | Protected Member Functions | Protected Attributes | List of all members
lsst::afw::math::SeparableKernel Class Reference

A kernel described by a pair of functions: func(x, y) = colFunc(x) * rowFunc(y) More...

#include <Kernel.h>

Inheritance diagram for lsst::afw::math::SeparableKernel:
lsst::afw::table::io::PersistableFacade< SeparableKernel > lsst::afw::math::Kernel lsst::afw::table::io::PersistableFacade< Kernel > lsst::afw::table::io::Persistable lsst::afw::math::BilinearWarpingKernel lsst::afw::math::LanczosWarpingKernel lsst::afw::math::NearestWarpingKernel

Public Types

typedef lsst::afw::math::Function1< PixelKernelFunction
 
typedef std::shared_ptr< KernelFunctionKernelFunctionPtr
 
typedef double Pixel
 
typedef std::shared_ptr< lsst::afw::math::Function2< double > > SpatialFunctionPtr
 
typedef lsst::afw::math::Function2< double > SpatialFunction
 
typedef lsst::afw::math::NullFunction2< double > NullSpatialFunction
 
typedef generic_kernel_tag kernel_fill_factor
 

Public Member Functions

 SeparableKernel ()
 Construct an empty spatially invariant SeparableKernel of size 0x0. More...
 
 SeparableKernel (int width, int height, KernelFunction const &kernelColFunction, KernelFunction const &kernelRowFunction, Kernel::SpatialFunction const &spatialFunction=NullSpatialFunction())
 Construct a spatially invariant SeparableKernel, or a spatially varying SeparableKernel that uses the same functional form to model each function parameter. More...
 
 SeparableKernel (int width, int height, KernelFunction const &kernelColFunction, KernelFunction const &kernelRowFunction, std::vector< Kernel::SpatialFunctionPtr > const &spatialFunctionList)
 Construct a spatially varying SeparableKernel. More...
 
 SeparableKernel (const SeparableKernel &)=delete
 
 SeparableKernel (SeparableKernel &&)=delete
 
SeparableKerneloperator= (const SeparableKernel &)=delete
 
SeparableKerneloperator= (SeparableKernel &&)=delete
 
 ~SeparableKernel () override=default
 
std::shared_ptr< Kernelclone () const override
 Return a pointer to a deep copy of this kernel. More...
 
std::shared_ptr< Kernelresized (int width, int height) const override
 Return a pointer to a clone with specified kernel dimensions. More...
 
double computeVectors (std::vector< Pixel > &colList, std::vector< Pixel > &rowList, bool doNormalize, double x=0.0, double y=0.0) const
 Compute the column and row arrays in place, where kernel(col, row) = colList(col) * rowList(row) More...
 
double getKernelParameter (unsigned int i) const override
 Return a particular Kernel Parameter (no bounds checking). More...
 
std::vector< double > getKernelParameters () const override
 Return the current kernel parameters. More...
 
KernelFunctionPtr getKernelColFunction () const
 Get a deep copy of the col kernel function. More...
 
KernelFunctionPtr getKernelRowFunction () const
 Get a deep copy of the row kernel function. More...
 
std::string toString (std::string const &prefix="") const override
 Return a string representation of the kernel. More...
 
void computeCache (int const cacheSize) override
 Compute a cache of Kernel values, if desired. More...
 
int getCacheSize () const override
 Get the current cache size (0 if none) More...
 
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. More...
 
lsst::geom::Extent2I const getDimensions () const
 Return the Kernel's dimensions (width, height) More...
 
void setDimensions (lsst::geom::Extent2I dims)
 
void setWidth (int width)
 
void setHeight (int height)
 
int getWidth () const
 Return the Kernel's width. More...
 
int getHeight () const
 Return the Kernel's height. More...
 
lsst::geom::Point2I getCtr () const
 Return index of kernel's center. More...
 
int getCtrX () const
 Return x index of kernel's center. More...
 
int getCtrY () const
 Return y index of kernel's center. More...
 
lsst::geom::Box2I getBBox () const
 return parent bounding box, with XY0 = -center More...
 
unsigned int getNKernelParameters () const
 Return the number of kernel parameters (0 if none) More...
 
int getNSpatialParameters () const
 Return the number of spatial parameters (0 if not spatially varying) More...
 
SpatialFunctionPtr getSpatialFunction (unsigned int index) const
 Return a clone of the specified spatial function (one component of the spatial model) More...
 
std::vector< SpatialFunctionPtrgetSpatialFunctionList () const
 Return a list of clones of the spatial functions. More...
 
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. More...
 
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. More...
 
void setCtr (lsst::geom::Point2I ctr)
 Set index of kernel's center. More...
 
void setCtrX (int ctrX)
 Set x index of kernel's center. More...
 
void setCtrY (int ctrY)
 Set y index of kernel's center. More...
 
std::vector< std::vector< double > > getSpatialParameters () const
 Return the spatial parameters parameters (an empty vector if not spatially varying) More...
 
bool isSpatiallyVarying () const
 Return true iff the kernel is spatially varying (has a spatial function) More...
 
void setKernelParameters (std::vector< double > const &params)
 Set the kernel parameters of a spatially invariant kernel. More...
 
void setKernelParameters (std::pair< double, double > const &params)
 Set the kernel parameters of a 2-component spatially invariant kernel. More...
 
void setSpatialParameters (const std::vector< std::vector< double >> params)
 Set the parameters of all spatial functions. More...
 
void computeKernelParametersFromSpatialModel (std::vector< double > &kernelParams, double x, double y) const
 Compute the kernel parameters at a specified point. More...
 
void writeFits (std::string const &fileName, std::string const &mode="w") const
 Write the object to a regular FITS file. More...
 
void writeFits (fits::MemFileManager &manager, std::string const &mode="w") const
 Write the object to a FITS image in memory. More...
 
void writeFits (fits::Fits &fitsfile) const
 Write the object to an already-open FITS object. More...
 
virtual bool isPersistable () const noexcept
 Return true if this particular object can be persisted using afw::table::io. More...
 

Static Public Member Functions

static std::shared_ptr< SeparableKernelreadFits (fits::Fits &fitsfile)
 Read an object from an already open FITS object. More...
 
static std::shared_ptr< SeparableKernelreadFits (std::string const &fileName, int hdu=fits::DEFAULT_HDU)
 Read an object from a regular FITS file. More...
 
static std::shared_ptr< SeparableKernelreadFits (fits::MemFileManager &manager, int hdu=fits::DEFAULT_HDU)
 Read an object from a FITS file in memory. More...
 
static std::shared_ptr< SeparableKerneldynamicCast (std::shared_ptr< Persistable > const &ptr)
 Dynamically cast a shared_ptr. More...
 
static std::shared_ptr< KernelreadFits (fits::Fits &fitsfile)
 Read an object from an already open FITS object. More...
 
static std::shared_ptr< KernelreadFits (std::string const &fileName, int hdu=fits::DEFAULT_HDU)
 Read an object from a regular FITS file. More...
 
static std::shared_ptr< KernelreadFits (fits::MemFileManager &manager, int hdu=fits::DEFAULT_HDU)
 Read an object from a FITS file in memory. More...
 
static std::shared_ptr< KerneldynamicCast (std::shared_ptr< Persistable > const &ptr)
 Dynamically cast a shared_ptr. More...
 

Protected Types

typedef io::OutputArchiveHandle OutputArchiveHandle
 

Protected Member Functions

double doComputeImage (lsst::afw::image::Image< Pixel > &image, bool doNormalize) const override
 Low-level version of computeImage. More...
 
void setKernelParameter (unsigned int ind, double value) const override
 Set one kernel parameter. More...
 
std::string getPythonModule () const override
 Return the fully-qualified Python module that should be imported to guarantee that its factory is registered. More...
 
void setKernelParametersFromSpatialModel (double x, double y) const
 Set the kernel parameters from the spatial model (if any). More...
 
virtual std::string getPersistenceName () const
 Return the unique name used to persist this object and look up its factory. More...
 
virtual void write (OutputArchiveHandle &handle) const
 Write the object to one or more catalogs. More...
 

Protected Attributes

std::vector< SpatialFunctionPtr_spatialFunctionList
 

Detailed Description

A kernel described by a pair of functions: func(x, y) = colFunc(x) * rowFunc(y)

The function's x, y arguments are as follows:

Note: each pixel is set to the value of the kernel function at the center of the pixel (rather than averaging the function over the area of the pixel).

Definition at line 907 of file Kernel.h.

Member Typedef Documentation

◆ kernel_fill_factor

Definition at line 119 of file Kernel.h.

◆ KernelFunction

Definition at line 909 of file Kernel.h.

◆ KernelFunctionPtr

Definition at line 910 of file Kernel.h.

◆ NullSpatialFunction

Definition at line 116 of file Kernel.h.

◆ OutputArchiveHandle

typedef io::OutputArchiveHandle lsst::afw::table::io::Persistable::OutputArchiveHandle
protectedinherited

Definition at line 108 of file Persistable.h.

◆ Pixel

typedef double lsst::afw::math::Kernel::Pixel
inherited

Definition at line 113 of file Kernel.h.

◆ SpatialFunction

Definition at line 115 of file Kernel.h.

◆ SpatialFunctionPtr

Definition at line 114 of file Kernel.h.

Constructor & Destructor Documentation

◆ SeparableKernel() [1/5]

lsst::afw::math::SeparableKernel::SeparableKernel ( )
explicit

Construct an empty spatially invariant SeparableKernel of size 0x0.

Definition at line 43 of file SeparableKernel.cc.

44  : Kernel(),
45  _kernelColFunctionPtr(),
46  _kernelRowFunctionPtr(),
47  _localColList(0),
48  _localRowList(0),
49  _kernelX(0),
50  _kernelY(0),
51  _kernelRowCache(0),
52  _kernelColCache(0) {
53  _setKernelXY();
54 }
Kernel()
Construct a null Kernel of size 0,0.
Definition: Kernel.cc:58

◆ SeparableKernel() [2/5]

lsst::afw::math::SeparableKernel::SeparableKernel ( int  width,
int  height,
KernelFunction const &  kernelColFunction,
KernelFunction const &  kernelRowFunction,
Kernel::SpatialFunction const &  spatialFunction = NullSpatialFunction() 
)
explicit

Construct a spatially invariant SeparableKernel, or a spatially varying SeparableKernel that uses the same functional form to model each function parameter.

Parameters
widthwidth of kernel
heightheight of kernel
kernelColFunctionkernel column function
kernelRowFunctionkernel row function
spatialFunctionspatial function; one deep copy is made for each kernel column and row function parameter; if omitted or set to Kernel::NullSpatialFunction then the kernel is spatially invariant

Definition at line 56 of file SeparableKernel.cc.

59  : Kernel(width, height, kernelColFunction.getNParameters() + kernelRowFunction.getNParameters(),
60  spatialFunction),
61  _kernelColFunctionPtr(kernelColFunction.clone()),
62  _kernelRowFunctionPtr(kernelRowFunction.clone()),
63  _localColList(width),
64  _localRowList(height),
65  _kernelX(width),
66  _kernelY(height),
67  _kernelRowCache(0),
68  _kernelColCache(0) {
69  _setKernelXY();
70 }
Kernel()
Construct a null Kernel of size 0,0.
Definition: Kernel.cc:58

◆ SeparableKernel() [3/5]

lsst::afw::math::SeparableKernel::SeparableKernel ( int  width,
int  height,
KernelFunction const &  kernelColFunction,
KernelFunction const &  kernelRowFunction,
std::vector< Kernel::SpatialFunctionPtr > const &  spatialFunctionList 
)
explicit

Construct a spatially varying SeparableKernel.

Parameters
widthwidth of kernel
heightheight of kernel
kernelColFunctionkernel column function
kernelRowFunctionkernel row function
spatialFunctionListlist of spatial funcs, one per kernel column and row function parameter; a deep copy is made of each function
Exceptions
lsst::pex::exceptions::InvalidParameterErrorif the length of spatialFunctionList != # kernel function parameters.

Definition at line 72 of file SeparableKernel.cc.

75  : Kernel(width, height, spatialFunctionList),
76  _kernelColFunctionPtr(kernelColFunction.clone()),
77  _kernelRowFunctionPtr(kernelRowFunction.clone()),
78  _localColList(width),
79  _localRowList(height),
80  _kernelX(width),
81  _kernelY(height),
82  _kernelRowCache(0),
83  _kernelColCache(0) {
84  if (kernelColFunction.getNParameters() + kernelRowFunction.getNParameters() !=
85  spatialFunctionList.size()) {
87  os << "kernelColFunction.getNParameters() + kernelRowFunction.getNParameters() = "
88  << kernelColFunction.getNParameters() << " + " << kernelRowFunction.getNParameters()
89  << " != " << spatialFunctionList.size() << " = "
90  << "spatialFunctionList.size()";
92  }
93 
94  _setKernelXY();
95 }
T str(T... args)
T size(T... args)
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
Kernel()
Construct a null Kernel of size 0,0.
Definition: Kernel.cc:58
Reports invalid arguments.
Definition: Runtime.h:66
std::ostream * os
Definition: Schema.cc:746

◆ SeparableKernel() [4/5]

lsst::afw::math::SeparableKernel::SeparableKernel ( const SeparableKernel )
delete

◆ SeparableKernel() [5/5]

lsst::afw::math::SeparableKernel::SeparableKernel ( SeparableKernel &&  )
delete

◆ ~SeparableKernel()

lsst::afw::math::SeparableKernel::~SeparableKernel ( )
overridedefault

Member Function Documentation

◆ clone()

std::shared_ptr< Kernel > lsst::afw::math::SeparableKernel::clone ( ) const
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.

Returns
a pointer to a deep copy of the kernel

Implements lsst::afw::math::Kernel.

Reimplemented in lsst::afw::math::NearestWarpingKernel, lsst::afw::math::BilinearWarpingKernel, and lsst::afw::math::LanczosWarpingKernel.

Definition at line 97 of file SeparableKernel.cc.

97  {
99  if (this->isSpatiallyVarying()) {
100  retPtr.reset(new SeparableKernel(this->getWidth(), this->getHeight(), *(this->_kernelColFunctionPtr),
101  *(this->_kernelRowFunctionPtr), this->_spatialFunctionList));
102  } else {
103  retPtr.reset(new SeparableKernel(this->getWidth(), this->getHeight(), *(this->_kernelColFunctionPtr),
104  *(this->_kernelRowFunctionPtr)));
105  }
106  retPtr->setCtr(this->getCtr());
107  retPtr->computeCache(this->getCacheSize());
108  return retPtr;
109 }
int getHeight() const
Return the Kernel&#39;s height.
Definition: Kernel.h:230
bool isSpatiallyVarying() const
Return true iff the kernel is spatially varying (has a spatial function)
Definition: Kernel.h:380
SeparableKernel()
Construct an empty spatially invariant SeparableKernel of size 0x0.
lsst::geom::Point2I getCtr() const
Return index of kernel&#39;s center.
Definition: Kernel.h:235
T reset(T... args)
int getCacheSize() const override
Get the current cache size (0 if none)
int getWidth() const
Return the Kernel&#39;s width.
Definition: Kernel.h:225
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:496

◆ computeCache()

void lsst::afw::math::SeparableKernel::computeCache ( int const  const)
overridevirtual

Compute a cache of Kernel values, if desired.

Warning
: few kernel classes actually support this, in which case this is a no-op and getCacheSize always returns 0.

Reimplemented from lsst::afw::math::Kernel.

Definition at line 310 of file SeparableKernel.cc.

310  {
312 
313  func = getKernelColFunction();
314  _computeCache(cacheSize, _kernelY, func, &_kernelColCache);
315 
316  func = getKernelRowFunction();
317  _computeCache(cacheSize, _kernelX, func, &_kernelRowCache);
318 }
std::shared_ptr< KernelFunction > KernelFunctionPtr
Definition: Kernel.h:910
afw::table::Key< int > cacheSize
Definition: CoaddPsf.cc:339
KernelFunctionPtr getKernelRowFunction() const
Get a deep copy of the row kernel function.
KernelFunctionPtr getKernelColFunction() const
Get a deep copy of the col kernel function.

◆ computeImage()

double lsst::afw::math::Kernel::computeImage ( lsst::afw::image::Image< Pixel > &  image,
bool  doNormalize,
double  x = 0.0,
double  y = 0.0 
) const
inherited

Compute an image (pixellized representation of the kernel) in place.

Parameters
imageimage whose pixels are to be set (output); xy0 of the image will be set to -kernel.getCtr()
doNormalizenormalize the image (so sum is 1)?
xx (column position) at which to compute spatial function
yy (row position) at which to compute spatial function
Returns
The kernel sum
Exceptions
lsst::pex::exceptions::InvalidParameterErrorif the image is the wrong size
lsst::pex::exceptions::OverflowErrorif doNormalize is true and the kernel sum is exactly 0
Note
computeNewImage has been retired; it doesn't need to be a member

Definition at line 85 of file Kernel.cc.

85  {
86  if (image.getDimensions() != this->getDimensions()) {
88  os << "image dimensions = ( " << image.getWidth() << ", " << image.getHeight() << ") != ("
89  << this->getWidth() << ", " << this->getHeight() << ") = kernel dimensions";
91  }
92  image.setXY0(-_ctrX, -_ctrY);
93  if (this->isSpatiallyVarying()) {
95  }
96  return doComputeImage(image, doNormalize);
97 }
int getHeight() const
Return the Kernel&#39;s height.
Definition: Kernel.h:230
int getHeight() const
Return the number of rows in the image.
Definition: ImageBase.h:335
int y
Definition: SpanSet.cc:49
void setXY0(lsst::geom::Point2I const origin)
Set the ImageBase&#39;s origin.
Definition: ImageBase.h:471
bool isSpatiallyVarying() const
Return true iff the kernel is spatially varying (has a spatial function)
Definition: Kernel.h:380
T str(T... args)
double x
void setKernelParametersFromSpatialModel(double x, double y) const
Set the kernel parameters from the spatial model (if any).
Definition: Kernel.cc:228
int getWidth() const
Return the Kernel&#39;s width.
Definition: Kernel.h:225
lsst::geom::Extent2I const getDimensions() const
Return the Kernel&#39;s dimensions (width, height)
Definition: Kernel.h:213
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
int getWidth() const
Return the number of columns in the image.
Definition: ImageBase.h:333
virtual double doComputeImage(lsst::afw::image::Image< Pixel > &image, bool doNormalize) const =0
Low-level version of computeImage.
Reports invalid arguments.
Definition: Runtime.h:66
lsst::geom::Extent2I getDimensions() const
Return the image&#39;s size; useful for passing to constructors.
Definition: ImageBase.h:393
std::ostream * os
Definition: Schema.cc:746

◆ computeKernelParametersFromSpatialModel()

void lsst::afw::math::Kernel::computeKernelParametersFromSpatialModel ( std::vector< double > &  kernelParams,
double  x,
double  y 
) const
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 144 of file Kernel.cc.

145  {
146  std::vector<double>::iterator paramIter = kernelParams.begin();
148  for (; funcIter != _spatialFunctionList.end(); ++funcIter, ++paramIter) {
149  *paramIter = (*(*funcIter))(x, y);
150  }
151 }
int y
Definition: SpanSet.cc:49
double x
STL class.
T begin(T... args)
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:496

◆ computeVectors()

double lsst::afw::math::SeparableKernel::computeVectors ( std::vector< Pixel > &  colList,
std::vector< Pixel > &  rowList,
bool  doNormalize,
double  x = 0.0,
double  y = 0.0 
) const

Compute the column and row arrays in place, where kernel(col, row) = colList(col) * rowList(row)

x, y are ignored if there is no spatial function.

Parameters
colListcolumn vector
rowListrow vector
doNormalizenormalize the image (so sum of each is 1)?
xx (column position) at which to compute spatial function
yy (row position) at which to compute spatial function
Returns
the kernel sum (1.0 if doNormalize true)
Exceptions
lsst::pex::exceptions::InvalidParameterErrorif colList or rowList is the wrong size
lsst::pex::exceptions::OverflowErrorif doNormalize is true and the kernel sum is exactly 0

Definition at line 123 of file SeparableKernel.cc.

124  {
125  if (static_cast<int>(colList.size()) != this->getWidth() ||
126  static_cast<int>(rowList.size()) != this->getHeight()) {
128  os << "colList.size(), rowList.size() = (" << colList.size() << ", " << rowList.size() << ") != ("
129  << this->getWidth() << ", " << this->getHeight() << ") = "
130  << "kernel dimensions";
132  }
133  if (this->isSpatiallyVarying()) {
135  }
136 
137  return basicComputeVectors(colList, rowList, doNormalize);
138 }
int getHeight() const
Return the Kernel&#39;s height.
Definition: Kernel.h:230
int y
Definition: SpanSet.cc:49
bool isSpatiallyVarying() const
Return true iff the kernel is spatially varying (has a spatial function)
Definition: Kernel.h:380
T str(T... args)
double x
void setKernelParametersFromSpatialModel(double x, double y) const
Set the kernel parameters from the spatial model (if any).
Definition: Kernel.cc:228
int getWidth() const
Return the Kernel&#39;s width.
Definition: Kernel.h:225
T size(T... args)
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
Reports invalid arguments.
Definition: Runtime.h:66
std::ostream * os
Definition: Schema.cc:746

◆ doComputeImage()

double lsst::afw::math::SeparableKernel::doComputeImage ( lsst::afw::image::Image< Pixel > &  image,
bool  doNormalize 
) const
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.

Parameters
imageimage whose pixels are to be set (output)
doNormalizenormalize the image (so sum is 1)?
Returns
The kernel sum

Implements lsst::afw::math::Kernel.

Definition at line 172 of file SeparableKernel.cc.

172  {
173  double imSum = basicComputeVectors(_localColList, _localRowList, doNormalize);
174 
175  for (int y = 0; y != image.getHeight(); ++y) {
177  for (std::vector<Pixel>::iterator colIter = _localColList.begin(); colIter != _localColList.end();
178  ++colIter, ++imPtr) {
179  *imPtr = (*colIter) * _localRowList[y];
180  }
181  }
182 
183  return imSum;
184 }
int getHeight() const
Return the number of rows in the image.
Definition: ImageBase.h:335
int y
Definition: SpanSet.cc:49
T end(T... args)
x_iterator row_begin(int y) const
Return an x_iterator to the start of the y&#39;th row.
Definition: ImageBase.h:438
STL class.
T begin(T... args)
A class to represent a 2-dimensional array of pixels.
Definition: Image.h:58

◆ dynamicCast() [1/2]

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;

Exceptions
lsst::pex::exceptions::LogicErrorif the cast fails

param[in] ptr The pointer to be cast.

Returns
The cast pointer.
Exceptions
lsst::pex::exceptions::TypeErrorIf the dynamic cast fails.

Definition at line 18 of file Persistable.cc.

18  {
20  if (!result) {
21  throw LSST_EXCEPT(pex::exceptions::TypeError, "Dynamic pointer cast failed");
22  }
23  return result;
24 }
uint64_t * ptr
Definition: RangeSet.cc:88
T dynamic_pointer_cast(T... args)
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
py::object result
Definition: _schema.cc:429

◆ dynamicCast() [2/2]

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;

Exceptions
lsst::pex::exceptions::LogicErrorif the cast fails

param[in] ptr The pointer to be cast.

Returns
The cast pointer.
Exceptions
lsst::pex::exceptions::TypeErrorIf the dynamic cast fails.

Definition at line 18 of file Persistable.cc.

18  {
20  if (!result) {
21  throw LSST_EXCEPT(pex::exceptions::TypeError, "Dynamic pointer cast failed");
22  }
23  return result;
24 }
uint64_t * ptr
Definition: RangeSet.cc:88
T dynamic_pointer_cast(T... args)
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
py::object result
Definition: _schema.cc:429

◆ getBBox()

lsst::geom::Box2I lsst::afw::math::Kernel::getBBox ( ) const
inlineinherited

return parent bounding box, with XY0 = -center

Definition at line 262 of file Kernel.h.

◆ getCacheSize()

int lsst::afw::math::SeparableKernel::getCacheSize ( ) const
overridevirtual

Get the current cache size (0 if none)

Reimplemented from lsst::afw::math::Kernel.

Definition at line 320 of file SeparableKernel.cc.

320 { return _kernelColCache.size(); };
T size(T... args)

◆ getCtr()

lsst::geom::Point2I lsst::afw::math::Kernel::getCtr ( ) const
inlineinherited

Return index of kernel's center.

Definition at line 235 of file Kernel.h.

235 { return lsst::geom::Point2I(_ctrX, _ctrY); }
Point< int, 2 > Point2I
Definition: Point.h:321

◆ getCtrX()

int lsst::afw::math::Kernel::getCtrX ( ) const
inlineinherited

Return x index of kernel's center.

Deprecated:
Use getCtr instead

Definition at line 244 of file Kernel.h.

244  {
245  return _ctrX;
246  }

◆ getCtrY()

int lsst::afw::math::Kernel::getCtrY ( ) const
inlineinherited

Return y index of kernel's center.

Deprecated:
Use getCtr instead

Definition at line 255 of file Kernel.h.

255  {
256  return _ctrY;
257  }

◆ getDimensions()

lsst::geom::Extent2I const lsst::afw::math::Kernel::getDimensions ( ) const
inlineinherited

Return the Kernel's dimensions (width, height)

Definition at line 213 of file Kernel.h.

213 { return lsst::geom::Extent2I(_width, _height); }
Extent< int, 2 > Extent2I
Definition: Extent.h:397

◆ getHeight()

int lsst::afw::math::Kernel::getHeight ( ) const
inlineinherited

Return the Kernel's height.

Definition at line 230 of file Kernel.h.

230 { return _height; }

◆ getKernelColFunction()

SeparableKernel::KernelFunctionPtr lsst::afw::math::SeparableKernel::getKernelColFunction ( ) const

Get a deep copy of the col kernel function.

Definition at line 140 of file SeparableKernel.cc.

140  {
141  return _kernelColFunctionPtr->clone();
142 }

◆ getKernelParameter()

double lsst::afw::math::SeparableKernel::getKernelParameter ( unsigned int  i) const
inlineoverridevirtual

Return a particular Kernel Parameter (no bounds checking).

This version is slow, but specialisations may be faster

Reimplemented from lsst::afw::math::Kernel.

Definition at line 980 of file Kernel.h.

980  {
981  unsigned int const ncol = _kernelColFunctionPtr->getNParameters();
982  if (i < ncol) {
983  return _kernelColFunctionPtr->getParameter(i);
984  } else {
985  i -= ncol;
986  return _kernelRowFunctionPtr->getParameter(i);
987  }
988  }

◆ getKernelParameters()

std::vector< double > lsst::afw::math::SeparableKernel::getKernelParameters ( ) const
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 161 of file SeparableKernel.cc.

161  {
162  std::vector<double> allParams = _kernelColFunctionPtr->getParameters();
163  std::vector<double> yParams = _kernelRowFunctionPtr->getParameters();
164  std::copy(yParams.begin(), yParams.end(), std::back_inserter(allParams));
165  return allParams;
166 }
T copy(T... args)
T end(T... args)
T begin(T... args)
T back_inserter(T... args)

◆ getKernelRowFunction()

SeparableKernel::KernelFunctionPtr lsst::afw::math::SeparableKernel::getKernelRowFunction ( ) const

Get a deep copy of the row kernel function.

Definition at line 144 of file SeparableKernel.cc.

144  {
145  return _kernelRowFunctionPtr->clone();
146 }

◆ getNKernelParameters()

unsigned int lsst::afw::math::Kernel::getNKernelParameters ( ) const
inlineinherited

Return the number of kernel parameters (0 if none)

Definition at line 269 of file Kernel.h.

269 { return _nKernelParams; }

◆ getNSpatialParameters()

int lsst::afw::math::Kernel::getNSpatialParameters ( ) const
inlineinherited

Return the number of spatial parameters (0 if not spatially varying)

Definition at line 274 of file Kernel.h.

274  {
275  return this->isSpatiallyVarying() ? _spatialFunctionList[0]->getNParameters() : 0;
276  }
bool isSpatiallyVarying() const
Return true iff the kernel is spatially varying (has a spatial function)
Definition: Kernel.h:380
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:496

◆ getPersistenceName()

std::string lsst::afw::table::io::Persistable::getPersistenceName ( ) const
protectedvirtualinherited

◆ getPythonModule()

std::string lsst::afw::math::Kernel::getPythonModule ( ) const
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.

Definition at line 235 of file Kernel.cc.

235 { return "lsst.afw.math"; }

◆ getSpatialFunction()

Kernel::SpatialFunctionPtr lsst::afw::math::Kernel::getSpatialFunction ( unsigned int  index) const
inherited

Return a clone of the specified spatial function (one component of the spatial model)

Parameters
indexindex of desired spatial function; must be in range [0, number spatial parameters - 1]
Returns
a shared pointer to a spatial function. The function is a deep copy, so setting its parameters has no effect on the kernel.
Exceptions
lsst::pex::exceptions::InvalidParameterErrorif kernel not spatially varying
lsst::pex::exceptions::InvalidParameterErrorif index out of range

Definition at line 153 of file Kernel.cc.

153  {
154  if (index >= _spatialFunctionList.size()) {
155  if (!this->isSpatiallyVarying()) {
156  throw LSST_EXCEPT(pexExcept::InvalidParameterError, "kernel is not spatially varying");
157  } else {
158  std::ostringstream errStream;
159  errStream << "index = " << index << "; must be < , " << _spatialFunctionList.size();
161  }
162  }
163  return _spatialFunctionList[index]->clone();
164 }
bool isSpatiallyVarying() const
Return true iff the kernel is spatially varying (has a spatial function)
Definition: Kernel.h:380
T str(T... args)
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:496
Reports invalid arguments.
Definition: Runtime.h:66

◆ getSpatialFunctionList()

std::vector< Kernel::SpatialFunctionPtr > lsst::afw::math::Kernel::getSpatialFunctionList ( ) const
inherited

Return a list of clones of the spatial functions.

Returns
a list of shared pointers to spatial functions. The functions are deep copies, so setting their parameters has no effect on the kernel.

Definition at line 166 of file Kernel.cc.

166  {
167  std::vector<SpatialFunctionPtr> spFuncCopyList;
169  spFuncIter != _spatialFunctionList.end(); ++spFuncIter) {
170  spFuncCopyList.push_back((**spFuncIter).clone());
171  }
172  return spFuncCopyList;
173 }
T push_back(T... args)
STL class.
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:496

◆ getSpatialParameters()

std::vector<std::vector<double> > lsst::afw::math::Kernel::getSpatialParameters ( ) const
inlineinherited

Return the spatial parameters parameters (an empty vector if not spatially varying)

Definition at line 368 of file Kernel.h.

368  {
369  std::vector<std::vector<double>> spatialParams;
371  for (; spFuncIter != _spatialFunctionList.end(); ++spFuncIter) {
372  spatialParams.push_back((*spFuncIter)->getParameters());
373  }
374  return spatialParams;
375  }
T push_back(T... args)
STL class.
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:496

◆ getWidth()

int lsst::afw::math::Kernel::getWidth ( ) const
inlineinherited

Return the Kernel's width.

Definition at line 225 of file Kernel.h.

225 { return _width; }

◆ growBBox()

lsst::geom::Box2I lsst::afw::math::Kernel::growBBox ( lsst::geom::Box2I const &  bbox) const
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.

Returns
the bbox expanded by the kernel.

Definition at line 177 of file Kernel.cc.

177  {
178  return lsst::geom::Box2I(
180  lsst::geom::Extent2I(bbox.getDimensions() + getDimensions() - lsst::geom::Extent2I(1, 1)));
181 }
AmpInfoBoxKey bbox
Definition: Amplifier.cc:117
lsst::geom::Point2I getCtr() const
Return index of kernel&#39;s center.
Definition: Kernel.h:235
lsst::geom::Extent2I const getDimensions() const
Return the Kernel&#39;s dimensions (width, height)
Definition: Kernel.h:213
Extent< int, 2 > Extent2I
Definition: Extent.h:397

◆ isPersistable()

virtual bool lsst::afw::table::io::Persistable::isPersistable ( ) const
inlinevirtualnoexceptinherited

Return true if this particular object can be persisted using afw::table::io.

Reimplemented in lsst::afw::math::LinearCombinationKernel, lsst::afw::math::Chebyshev1Function2< ReturnT >, lsst::afw::math::DeltaFunctionKernel, lsst::afw::geom::SpanSet, lsst::afw::math::AnalyticKernel, lsst::afw::math::FixedKernel, lsst::afw::math::PolynomialFunction2< ReturnT >, lsst::afw::image::PhotoCalib, lsst::meas::modelfit::Mixture, lsst::afw::math::DoubleGaussianFunction2< ReturnT >, lsst::afw::cameraGeom::Detector, lsst::afw::detection::Footprint, lsst::afw::geom::SkyWcs, lsst::afw::geom::Transform< FromEndpoint, ToEndpoint >, lsst::afw::geom::Transform< afw::geom::Point2Endpoint, afw::geom::GenericEndpoint >, lsst::afw::geom::polygon::Polygon, lsst::afw::math::GaussianFunction2< ReturnT >, lsst::afw::image::Filter, lsst::afw::cameraGeom::TransformMap, lsst::afw::math::ChebyshevBoundedField, lsst::meas::algorithms::CoaddPsf, lsst::afw::cameraGeom::DetectorCollection, lsst::afw::image::VisitInfo, lsst::afw::cameraGeom::Camera, lsst::meas::algorithms::CoaddBoundedField, lsst::afw::math::PixelAreaBoundedField, lsst::afw::detection::GaussianPsf, lsst::afw::image::CoaddInputs, lsst::afw::image::ApCorrMap, lsst::afw::math::TransformBoundedField, lsst::afw::math::ProductBoundedField, lsst::meas::extensions::psfex::PsfexPsf, lsst::meas::algorithms::KernelPsf, lsst::meas::algorithms::DoubleGaussianPsf, and lsst::meas::algorithms::SingleGaussianPsf.

Definition at line 102 of file Persistable.h.

102 { return false; }

◆ isSpatiallyVarying()

bool lsst::afw::math::Kernel::isSpatiallyVarying ( ) const
inlineinherited

Return true iff the kernel is spatially varying (has a spatial function)

Definition at line 380 of file Kernel.h.

380 { return _spatialFunctionList.size() != 0; }
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:496

◆ operator=() [1/2]

SeparableKernel& lsst::afw::math::SeparableKernel::operator= ( const SeparableKernel )
delete

◆ operator=() [2/2]

SeparableKernel& lsst::afw::math::SeparableKernel::operator= ( SeparableKernel &&  )
delete

◆ readFits() [1/6]

static std::shared_ptr<Kernel > lsst::afw::table::io::PersistableFacade< Kernel >::readFits ( fits::Fits fitsfile)
inlinestaticinherited

Read an object from an already open FITS object.

Parameters
[in]fitsfileFITS object to read from, already positioned at the desired HDU.

Definition at line 183 of file Persistable.h.

183  {
184  return dynamicCast(Persistable::_readFits(fitsfile));
185  }
static std::shared_ptr< Kernel > dynamicCast(std::shared_ptr< Persistable > const &ptr)
Dynamically cast a shared_ptr.
Definition: Persistable.cc:18

◆ readFits() [2/6]

static std::shared_ptr<SeparableKernel > lsst::afw::table::io::PersistableFacade< SeparableKernel >::readFits ( fits::Fits fitsfile)
inlinestaticinherited

Read an object from an already open FITS object.

Parameters
[in]fitsfileFITS object to read from, already positioned at the desired HDU.

Definition at line 183 of file Persistable.h.

183  {
184  return dynamicCast(Persistable::_readFits(fitsfile));
185  }
static std::shared_ptr< SeparableKernel > dynamicCast(std::shared_ptr< Persistable > const &ptr)
Dynamically cast a shared_ptr.
Definition: Persistable.cc:18

◆ readFits() [3/6]

static std::shared_ptr<SeparableKernel > lsst::afw::table::io::PersistableFacade< SeparableKernel >::readFits ( std::string const &  fileName,
int  hdu = fits::DEFAULT_HDU 
)
inlinestaticinherited

Read an object from a regular FITS file.

Parameters
[in]fileNameName of the file to read.
[in]hduHDU 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.

194  {
195  return dynamicCast(Persistable::_readFits(fileName, hdu));
196  }
static std::shared_ptr< SeparableKernel > dynamicCast(std::shared_ptr< Persistable > const &ptr)
Dynamically cast a shared_ptr.
Definition: Persistable.cc:18

◆ readFits() [4/6]

static std::shared_ptr<Kernel > lsst::afw::table::io::PersistableFacade< Kernel >::readFits ( std::string const &  fileName,
int  hdu = fits::DEFAULT_HDU 
)
inlinestaticinherited

Read an object from a regular FITS file.

Parameters
[in]fileNameName of the file to read.
[in]hduHDU 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.

194  {
195  return dynamicCast(Persistable::_readFits(fileName, hdu));
196  }
static std::shared_ptr< Kernel > dynamicCast(std::shared_ptr< Persistable > const &ptr)
Dynamically cast a shared_ptr.
Definition: Persistable.cc:18

◆ readFits() [5/6]

static std::shared_ptr<SeparableKernel > lsst::afw::table::io::PersistableFacade< SeparableKernel >::readFits ( fits::MemFileManager manager,
int  hdu = fits::DEFAULT_HDU 
)
inlinestaticinherited

Read an object from a FITS file in memory.

Parameters
[in]managerManager for the memory to read from.
[in]hduHDU 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.

205  {
206  return dynamicCast(Persistable::_readFits(manager, hdu));
207  }
static std::shared_ptr< SeparableKernel > dynamicCast(std::shared_ptr< Persistable > const &ptr)
Dynamically cast a shared_ptr.
Definition: Persistable.cc:18

◆ readFits() [6/6]

static std::shared_ptr<Kernel > lsst::afw::table::io::PersistableFacade< Kernel >::readFits ( fits::MemFileManager manager,
int  hdu = fits::DEFAULT_HDU 
)
inlinestaticinherited

Read an object from a FITS file in memory.

Parameters
[in]managerManager for the memory to read from.
[in]hduHDU 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.

205  {
206  return dynamicCast(Persistable::_readFits(manager, hdu));
207  }
static std::shared_ptr< Kernel > dynamicCast(std::shared_ptr< Persistable > const &ptr)
Dynamically cast a shared_ptr.
Definition: Persistable.cc:18

◆ resized()

std::shared_ptr< Kernel > lsst::afw::math::SeparableKernel::resized ( int  width,
int  height 
) const
overridevirtual

Return a pointer to a clone with specified kernel dimensions.

Parameters
widthNumber of columns in pixels
heightNumber of rows in pixels
Returns
a pointer to a clone with new dimensions.

Must be implemented by derived classes.

Implements lsst::afw::math::Kernel.

Definition at line 111 of file SeparableKernel.cc.

111  {
113  if (isSpatiallyVarying()) {
114  retPtr = std::make_shared<SeparableKernel>(width, height, *_kernelColFunctionPtr,
115  *_kernelRowFunctionPtr, _spatialFunctionList);
116  } else {
117  retPtr = std::make_shared<SeparableKernel>(width, height, *_kernelColFunctionPtr,
118  *_kernelRowFunctionPtr);
119  }
120  return retPtr;
121 }
bool isSpatiallyVarying() const
Return true iff the kernel is spatially varying (has a spatial function)
Definition: Kernel.h:380
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:496

◆ setCtr()

void lsst::afw::math::Kernel::setCtr ( lsst::geom::Point2I  ctr)
inlineinherited

Set index of kernel's center.

Definition at line 335 of file Kernel.h.

335  {
336  _ctrX = ctr.getX();
337  _ctrY = ctr.getY();
338  _setKernelXY();
339  }

◆ setCtrX()

void lsst::afw::math::Kernel::setCtrX ( int  ctrX)
inlineinherited

Set x index of kernel's center.

Deprecated:
Use setCtr instead

Definition at line 348 of file Kernel.h.

348  {
349  _ctrX = ctrX;
350  _setKernelXY();
351  }

◆ setCtrY()

void lsst::afw::math::Kernel::setCtrY ( int  ctrY)
inlineinherited

Set y index of kernel's center.

Deprecated:
Use setCtr instead

Definition at line 360 of file Kernel.h.

360  {
361  _ctrY = ctrY;
362  _setKernelXY();
363  }

◆ setDimensions()

void lsst::afw::math::Kernel::setDimensions ( lsst::geom::Extent2I  dims)
inlineinherited

Definition at line 215 of file Kernel.h.

215  {
216  _width = dims.getX();
217  _height = dims.getY();
218  }

◆ setHeight()

void lsst::afw::math::Kernel::setHeight ( int  height)
inlineinherited

Definition at line 220 of file Kernel.h.

220 { _height = height; }

◆ setKernelParameter()

void lsst::afw::math::SeparableKernel::setKernelParameter ( unsigned int  ind,
double  value 
) const
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.

Exceptions
lsst::pex::exceptions::InvalidParameterErroralways (unless subclassed)

Reimplemented from lsst::afw::math::Kernel.

Reimplemented in lsst::afw::math::NearestWarpingKernel, lsst::afw::math::BilinearWarpingKernel, and lsst::afw::math::LanczosWarpingKernel.

Definition at line 186 of file SeparableKernel.cc.

186  {
187  unsigned int const nColParams = _kernelColFunctionPtr->getNParameters();
188  if (ind < nColParams) {
189  _kernelColFunctionPtr->setParameter(ind, value);
190  } else {
191  _kernelRowFunctionPtr->setParameter(ind - nColParams, value);
192  }
193 }

◆ setKernelParameters() [1/2]

void lsst::afw::math::Kernel::setKernelParameters ( std::vector< double > const &  params)
inlineinherited

Set the kernel parameters of a spatially invariant kernel.

Exceptions
lsst::pex::exceptions::RuntimeErrorif the kernel has a spatial function
lsst::pex::exceptions::InvalidParameterErrorif the params vector is the wrong length

Definition at line 388 of file Kernel.h.

388  {
389  if (this->isSpatiallyVarying()) {
390  throw LSST_EXCEPT(lsst::pex::exceptions::RuntimeError, "Kernel is spatially varying");
391  }
392  const unsigned int nParams = this->getNKernelParameters();
393  if (nParams != params.size()) {
395  (boost::format("Number of parameters is wrong, saw %d expected %d") % nParams %
396  params.size())
397  .str());
398  }
399  for (unsigned int ii = 0; ii < nParams; ++ii) {
400  this->setKernelParameter(ii, params[ii]);
401  }
402  }
def format(config, name=None, writeSourceLine=True, prefix="", verbose=False)
Definition: history.py:174
virtual void setKernelParameter(unsigned int ind, double value) const
Set one kernel parameter.
Definition: Kernel.cc:224
bool isSpatiallyVarying() const
Return true iff the kernel is spatially varying (has a spatial function)
Definition: Kernel.h:380
unsigned int getNKernelParameters() const
Return the number of kernel parameters (0 if none)
Definition: Kernel.h:269
T size(T... args)
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
Reports invalid arguments.
Definition: Runtime.h:66
Reports errors that are due to events beyond the control of the program.
Definition: Runtime.h:104

◆ setKernelParameters() [2/2]

void lsst::afw::math::Kernel::setKernelParameters ( std::pair< double, double > const &  params)
inlineinherited

Set the kernel parameters of a 2-component spatially invariant kernel.

Warning
This is a low-level method intended for maximum efficiency when using warping kernels. No error checking is performed. Use the std::vector<double> form if you want safety.

Definition at line 410 of file Kernel.h.

410  {
411  this->setKernelParameter(0, params.first);
412  this->setKernelParameter(1, params.second);
413  }
virtual void setKernelParameter(unsigned int ind, double value) const
Set one kernel parameter.
Definition: Kernel.cc:224

◆ setKernelParametersFromSpatialModel()

void lsst::afw::math::Kernel::setKernelParametersFromSpatialModel ( double  x,
double  y 
) const
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 228 of file Kernel.cc.

228  {
230  for (int ii = 0; funcIter != _spatialFunctionList.end(); ++funcIter, ++ii) {
231  this->setKernelParameter(ii, (*(*funcIter))(x, y));
232  }
233 }
virtual void setKernelParameter(unsigned int ind, double value) const
Set one kernel parameter.
Definition: Kernel.cc:224
int y
Definition: SpanSet.cc:49
double x
STL class.
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:496

◆ setSpatialParameters()

void lsst::afw::math::Kernel::setSpatialParameters ( const std::vector< std::vector< double >>  params)
inherited

Set the parameters of all spatial functions.

Params is indexed as [kernel parameter][spatial parameter]

Exceptions
lsst::pex::exceptions::InvalidParameterErrorif params is the wrong shape (if this exception is thrown then no parameters are changed)

Definition at line 119 of file Kernel.cc.

119  {
120  // Check params size before changing anything
121  unsigned int nKernelParams = this->getNKernelParameters();
122  if (params.size() != nKernelParams) {
123  throw LSST_EXCEPT(
125  (boost::format("params has %d entries instead of %d") % params.size() % nKernelParams).str());
126  }
127  unsigned int nSpatialParams = this->getNSpatialParameters();
128  for (unsigned int ii = 0; ii < nKernelParams; ++ii) {
129  if (params[ii].size() != nSpatialParams) {
131  (boost::format("params[%d] has %d entries instead of %d") % ii %
132  params[ii].size() % nSpatialParams)
133  .str());
134  }
135  }
136  // Set parameters
137  if (nSpatialParams > 0) {
138  for (unsigned int ii = 0; ii < nKernelParams; ++ii) {
139  this->_spatialFunctionList[ii]->setParameters(params[ii]);
140  }
141  }
142 }
def format(config, name=None, writeSourceLine=True, prefix="", verbose=False)
Definition: history.py:174
unsigned int getNKernelParameters() const
Return the number of kernel parameters (0 if none)
Definition: Kernel.h:269
T size(T... args)
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:496
int getNSpatialParameters() const
Return the number of spatial parameters (0 if not spatially varying)
Definition: Kernel.h:274
Reports invalid arguments.
Definition: Runtime.h:66

◆ setWidth()

void lsst::afw::math::Kernel::setWidth ( int  width)
inlineinherited

Definition at line 219 of file Kernel.h.

219 { _width = width; }

◆ shrinkBBox()

lsst::geom::Box2I lsst::afw::math::Kernel::shrinkBBox ( lsst::geom::Box2I const &  bbox) const
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.

Returns
the bbox shrunk by the kernel.
Exceptions
lsst::pex::exceptions::InvalidParameterErrorif the resulting box would have dimension < 1 in either axis

Definition at line 183 of file Kernel.cc.

183  {
184  if ((bbox.getWidth() < getWidth()) || ((bbox.getHeight() < getHeight()))) {
186  os << "bbox dimensions = " << bbox.getDimensions() << " < (" << getWidth() << ", " << getHeight()
187  << ") in one or both dimensions";
189  }
190  return lsst::geom::Box2I(
191  lsst::geom::Point2I(bbox.getMinX() + getCtrX(), bbox.getMinY() + getCtrY()),
192  lsst::geom::Extent2I(bbox.getWidth() + 1 - getWidth(), bbox.getHeight() + 1 - getHeight()));
193 }
AmpInfoBoxKey bbox
Definition: Amplifier.cc:117
int getHeight() const
Return the Kernel&#39;s height.
Definition: Kernel.h:230
int getCtrX() const
Return x index of kernel&#39;s center.
Definition: Kernel.h:244
int getCtrY() const
Return y index of kernel&#39;s center.
Definition: Kernel.h:255
T str(T... args)
int getWidth() const
Return the Kernel&#39;s width.
Definition: Kernel.h:225
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
Reports invalid arguments.
Definition: Runtime.h:66
std::ostream * os
Definition: Schema.cc:746

◆ toString()

std::string lsst::afw::math::SeparableKernel::toString ( std::string const &  prefix = "") const
overridevirtual

Return a string representation of the kernel.

Reimplemented from lsst::afw::math::Kernel.

Definition at line 148 of file SeparableKernel.cc.

148  {
150  os << prefix << "SeparableKernel:" << std::endl;
151  os << prefix
152  << "..x (width) function: " << (_kernelColFunctionPtr ? _kernelColFunctionPtr->toString() : "None")
153  << std::endl;
154  os << prefix
155  << "..y (rows) function: " << (_kernelRowFunctionPtr ? _kernelRowFunctionPtr->toString() : "None")
156  << std::endl;
157  os << Kernel::toString(prefix + "\t");
158  return os.str();
159 }
T endl(T... args)
T str(T... args)
std::string prefix
Definition: SchemaMapper.cc:79
std::ostream * os
Definition: Schema.cc:746
virtual std::string toString(std::string const &prefix="") const
Return a string representation of the kernel.
Definition: Kernel.cc:195

◆ write()

void lsst::afw::table::io::Persistable::write ( OutputArchiveHandle handle) const
protectedvirtualinherited

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 in lsst::afw::math::LinearCombinationKernel, lsst::afw::math::Chebyshev1Function2< ReturnT >, lsst::afw::math::DeltaFunctionKernel, lsst::afw::math::AnalyticKernel, lsst::afw::math::FixedKernel, lsst::afw::math::PolynomialFunction2< ReturnT >, lsst::afw::image::PhotoCalib, lsst::afw::geom::SkyWcs, lsst::meas::modelfit::Mixture, lsst::afw::detection::Footprint, lsst::afw::math::DoubleGaussianFunction2< ReturnT >, lsst::afw::geom::Transform< FromEndpoint, ToEndpoint >, lsst::afw::geom::Transform< afw::geom::Point2Endpoint, afw::geom::GenericEndpoint >, lsst::afw::geom::polygon::Polygon, lsst::afw::math::GaussianFunction2< ReturnT >, lsst::afw::image::Filter, lsst::afw::math::ChebyshevBoundedField, lsst::afw::image::VisitInfo, lsst::meas::algorithms::CoaddPsf, lsst::afw::cameraGeom::DetectorCollection, lsst::afw::cameraGeom::Camera, lsst::afw::detection::HeavyFootprint< ImagePixelT, MaskPixelT, VariancePixelT >, lsst::meas::algorithms::CoaddBoundedField, lsst::afw::math::PixelAreaBoundedField, lsst::afw::image::CoaddInputs, lsst::afw::detection::GaussianPsf, lsst::afw::math::TransformBoundedField, lsst::afw::math::ProductBoundedField, lsst::meas::algorithms::KernelPsf, lsst::meas::extensions::psfex::PsfexPsf, lsst::meas::algorithms::DoubleGaussianPsf, and lsst::meas::algorithms::SingleGaussianPsf.

Definition at line 38 of file Persistable.cc.

38  {
39  assert(!isPersistable());
40  throw LSST_EXCEPT(pex::exceptions::LogicError,
41  "afw::table-based persistence is not supported for this object.");
42 }
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
virtual bool isPersistable() const noexcept
Return true if this particular object can be persisted using afw::table::io.
Definition: Persistable.h:102

◆ writeFits() [1/3]

void lsst::afw::table::io::Persistable::writeFits ( std::string const &  fileName,
std::string const &  mode = "w" 
) const
inherited

Write the object to a regular FITS file.

Parameters
[in]fileNameName of the file to write to.
[in]modeIf "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.

24  {
25  fits::Fits fitsfile(fileName, mode, fits::Fits::AUTO_CLOSE | fits::Fits::AUTO_CHECK);
26  writeFits(fitsfile);
27 }
void writeFits(std::string const &fileName, std::string const &mode="w") const
Write the object to a regular FITS file.
Definition: Persistable.cc:24

◆ writeFits() [2/3]

void lsst::afw::table::io::Persistable::writeFits ( fits::MemFileManager manager,
std::string const &  mode = "w" 
) const
inherited

Write the object to a FITS image in memory.

Parameters
[in]managerName of the file to write to.
[in]modeIf "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.

29  {
30  fits::Fits fitsfile(manager, mode, fits::Fits::AUTO_CLOSE | fits::Fits::AUTO_CHECK);
31  writeFits(fitsfile);
32 }
void writeFits(std::string const &fileName, std::string const &mode="w") const
Write the object to a regular FITS file.
Definition: Persistable.cc:24

◆ writeFits() [3/3]

void lsst::afw::table::io::Persistable::writeFits ( fits::Fits fitsfile) const
inherited

Write the object to an already-open FITS object.

Parameters
[in]fitsfileOpen FITS object to write to.

Definition at line 18 of file Persistable.cc.

18  {
19  OutputArchive archive;
20  archive.put(this);
21  archive.writeFits(fitsfile);
22 }

Member Data Documentation

◆ _spatialFunctionList

std::vector<SpatialFunctionPtr> lsst::afw::math::Kernel::_spatialFunctionList
protectedinherited

Definition at line 496 of file Kernel.h.


The documentation for this class was generated from the following files: