LSSTApplications  18.1.0
LSSTDataManagementBasePackage
Classes | Public Types | Public Member Functions | Static Public Member Functions | Protected Types | Protected Member Functions | Protected Attributes | List of all members
lsst::afw::math::Kernel Class Referenceabstract

Kernels are used for convolution with MaskedImages and (eventually) Images. More...

#include <Kernel.h>

Inheritance diagram for lsst::afw::math::Kernel:
lsst::daf::base::Citizen lsst::afw::table::io::PersistableFacade< Kernel > lsst::afw::table::io::Persistable lsst::afw::math::AnalyticKernel lsst::afw::math::DeltaFunctionKernel lsst::afw::math::FixedKernel lsst::afw::math::LinearCombinationKernel lsst::afw::math::SeparableKernel lsst::afw::math::BilinearWarpingKernel lsst::afw::math::LanczosWarpingKernel lsst::afw::math::NearestWarpingKernel

Classes

struct  PersistenceHelper
 

Public Types

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
 
enum  { magicSentinel = 0xdeadbeef }
 
typedef unsigned long memId
 Type of the block's ID. More...
 
typedef memId(* memNewCallback) (const memId cid)
 A function used to register a callback. More...
 
typedef memId(* memCallback) (const Citizen *ptr)
 

Public Member Functions

 Kernel ()
 Construct a null Kernel of size 0,0. More...
 
 Kernel (int width, int height, unsigned int nKernelParams, SpatialFunction const &spatialFunction=NullSpatialFunction())
 Construct a spatially invariant Kernel or a spatially varying Kernel with one spatial function that is duplicated as needed. More...
 
 Kernel (int width, int height, const std::vector< SpatialFunctionPtr > spatialFunctionList)
 Construct a spatially varying Kernel with a list of spatial functions (one per kernel parameter) More...
 
 Kernel (const Kernel &)=delete
 
 Kernel (Kernel &&)=delete
 
Kerneloperator= (const Kernel &)=delete
 
Kerneloperator= (Kernel &&)=delete
 
 ~Kernel () override=default
 
virtual std::shared_ptr< Kernelclone () const =0
 Return a pointer to a deep copy of this kernel. More...
 
virtual std::shared_ptr< Kernelresized (int width, int height) const =0
 Return a pointer to a clone with specified kernel dimensions. 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...
 
virtual double getKernelParameter (unsigned int i) const
 Return a particular Kernel Parameter (no bounds checking). More...
 
virtual std::vector< double > getKernelParameters () const
 Return the current kernel parameters. 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...
 
virtual std::string toString (std::string const &prefix="") const
 Return a string representation of the kernel. More...
 
virtual void computeCache (int const)
 Compute a cache of Kernel values, if desired. More...
 
virtual int getCacheSize () const
 Get the current size of the kernel cache (0 if none or if caches not supported) More...
 
std::string repr () const
 Return a string representation of a Citizen. More...
 
void markPersistent (void)
 Mark a Citizen as persistent and not destroyed until process end. More...
 
memId getId () const
 Return the Citizen's ID. 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 bool hasBeenCorrupted ()
 Check all allocated blocks for corruption. More...
 
static memId getNextMemId ()
 Return the memId of the next object to be allocated. More...
 
static int init ()
 Called once when the memory system is being initialised. 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...
 
Census

Provide a list of current Citizens

static int census (int, memId startingMemId=0)
 How many active Citizens are there? More...
 
static void census (std::ostream &stream, memId startingMemId=0)
 Print a list of all active Citizens to stream, sorted by ID. More...
 
static const std::vector< const Citizen * > * census ()
 Return a (newly allocated) std::vector of active Citizens sorted by ID. More...
 
callbackIDs

Set callback Ids.

The old Id is returned

static memId setNewCallbackId (memId id)
 Call the NewCallback when block is allocated. More...
 
static memId setDeleteCallbackId (memId id)
 Call the current DeleteCallback when block is deleted. More...
 
callbacks

Set the New/Delete callback functions; in each case the previously installed callback is returned.

These callback functions return a value which is Added to the previously registered id.

The default callback functions are called default{New,Delete}Callback; you may want to set a break point in these callbacks from your favourite debugger

static memNewCallback setNewCallback (memNewCallback func)
 Set the NewCallback function. More...
 
static memCallback setDeleteCallback (memCallback func)
 Set the DeleteCallback function. More...
 
static memCallback setCorruptionCallback (memCallback func)
 Set the CorruptionCallback function. More...
 

Protected Types

typedef io::OutputArchiveHandle OutputArchiveHandle
 

Protected Member Functions

std::string getPythonModule () const override
 Return the fully-qualified Python module that should be imported to guarantee that its factory is registered. More...
 
virtual void setKernelParameter (unsigned int ind, double value) const
 Set one kernel parameter. More...
 
void setKernelParametersFromSpatialModel (double x, double y) const
 Set the kernel parameters from the spatial model (if any). More...
 
virtual double doComputeImage (lsst::afw::image::Image< Pixel > &image, bool doNormalize) const =0
 Low-level version of computeImage. 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

Kernels are used for convolution with MaskedImages and (eventually) Images.

Kernel is a virtual base class; it cannot be instantiated. The following notes apply to Kernel and to its subclasses.

The template type should usually be float or double; integer kernels should be used with caution because they do not normalize well.

The center pixel of a Kernel is at index: (width-1)/2, (height-1)/2. Thus it is centered along columns/rows if the kernel has an odd number of columns/rows and shifted 1/2 pixel towards 0 otherwise. A kernel should have an odd number of columns and rows unless it is intended to shift an image.

Spatially Varying Kernels

Kernels may optionally vary spatially (so long as they have any kernel parameters). To make a spatially varying kernel, specify a spatial function at construction (you cannot change your mind once the kernel is constructed). You must also specify a set of spatial parameters, and you may do this at construction and/or later by calling setSpatialParameters. The spatial parameters are a vector (one per kernel function parameter) of spatial function parameters. In other words the spatial parameters are a vector of vectors indexed as [kernel parameter][spatial parameter]. The one spatial function is used to compute the kernel parameters at a given spatial position by computing each kernel parameter using its associated vector of spatial function parameters.

The convolve function computes the spatial function at the pixel position (not index) of the image. See the convolve function for details.

Note that if a kernel is spatially varying then you may not set the kernel parameters directly; that is the job of the spatial function! However, you may change the spatial parameters at any time.

Design Notes

The basic design is to use the same kernel class for both spatially varying and spatially invariant kernels. The user either does or does not supply a function describing the spatial variation at creation time. In addition, analytic kernels are described by a user-supplied function of the same basic type as the spatial variation function.

Several other designs were considered, including: A) Use different classes for spatially varying and spatially invariant versions of each kernel. Thus instead of three basic kernel classes (FixedKernel, AnalyticKernel and LinearCombinationKernel) we would have five (since FixedKernel cannot be spatially varying). Robert Lupton argued that was a needless expansion of the class hiearchy and I agreed. B) Construct analytic kernels by defining a subclass of AnalyticKernel that is specific to the desired functional (e.g. GaussianAnalyticKernel). If spatial models are handled the same way then this creates a serious proliferation of kernel classes (even if we only have two different spatial models, e.g. polynomial and Chebyshev polynomial). I felt it made more sense to define the spatial model by some kind of function class (often called a "functor"), and since we needed such a class, I chose to use it for the analytic kernel as well.

However, having a separate function class does introduce some potential inefficiencies. If a function is part of the class it can potentially be evaluated more quickly than calling a function for each pixel or spatial position.

A possible variant on the current design is to define the spatial model and analytic kernel by specifying the functions as template parameters. This has the potential to regain some efficiency in evaluating the functions. However, it would be difficult or impossible to pre-instantiate the desired template classes, a requirement of the LSST coding standards.

Definition at line 112 of file Kernel.h.

Member Typedef Documentation

◆ kernel_fill_factor

Definition at line 122 of file Kernel.h.

◆ memCallback

typedef memId(* lsst::daf::base::Citizen::memCallback) (const Citizen *ptr)
inherited

Definition at line 61 of file Citizen.h.

◆ memId

typedef unsigned long lsst::daf::base::Citizen::memId
inherited

Type of the block's ID.

Definition at line 58 of file Citizen.h.

◆ memNewCallback

typedef memId(* lsst::daf::base::Citizen::memNewCallback) (const memId cid)
inherited

A function used to register a callback.

Definition at line 60 of file Citizen.h.

◆ NullSpatialFunction

Definition at line 119 of file Kernel.h.

◆ OutputArchiveHandle

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

Definition at line 108 of file Persistable.h.

◆ Pixel

Definition at line 116 of file Kernel.h.

◆ SpatialFunction

Definition at line 118 of file Kernel.h.

◆ SpatialFunctionPtr

Definition at line 117 of file Kernel.h.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum
inherited
Enumerator
magicSentinel 

Definition at line 88 of file Citizen.h.

Constructor & Destructor Documentation

◆ Kernel() [1/5]

lsst::afw::math::Kernel::Kernel ( )
explicit

Construct a null Kernel of size 0,0.

A null constructor is primarily intended for persistence.

Definition at line 58 of file Kernel.cc.

59  : daf::base::Citizen(typeid(this)),
61  _width(0),
62  _height(0),
63  _ctrX(0),
64  _ctrY(0),
65  _nKernelParams(0) {}
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:487

◆ Kernel() [2/5]

lsst::afw::math::Kernel::Kernel ( int  width,
int  height,
unsigned int  nKernelParams,
SpatialFunction const &  spatialFunction = NullSpatialFunction() 
)
explicit

Construct a spatially invariant Kernel or a spatially varying Kernel with one spatial function that is duplicated as needed.

Parameters
widthnumber of columns
heightnumber of height
nKernelParamsnumber of kernel parameters
spatialFunctionspatial function, or NullSpatialFunction() if none specified
Exceptions
lsst::pex::exceptions::InvalidParameterErrorif a spatial function is specified and the kernel has no parameters.
lsst::pex::exceptions::InvalidParameterErrorif a width or height < 1

Definition at line 67 of file Kernel.cc.

68  : daf::base::Citizen(typeid(this)),
70  _width(width),
71  _height(height),
72  _ctrX((width - 1) / 2),
73  _ctrY((height - 1) / 2),
74  _nKernelParams(nKernelParams) {
75  if ((width < 1) || (height < 1)) {
77  os << "kernel height = " << height << " and/or width = " << width << " < 1";
79  }
80  if (dynamic_cast<const NullSpatialFunction *>(&spatialFunction)) {
81  // spatialFunction is not really present
82  } else {
83  if (nKernelParams == 0) {
84  throw LSST_EXCEPT(pexExcept::InvalidParameterError, "Kernel function has no parameters");
85  }
86  for (unsigned int ii = 0; ii < nKernelParams; ++ii) {
87  SpatialFunctionPtr spatialFunctionCopy = spatialFunction.clone();
88  this->_spatialFunctionList.push_back(spatialFunctionCopy);
89  }
90  }
91 }
std::shared_ptr< lsst::afw::math::Function2< double > > SpatialFunctionPtr
Definition: Kernel.h:117
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:487
Reports invalid arguments.
Definition: Runtime.h:66
std::ostream * os
Definition: Schema.cc:746

◆ Kernel() [3/5]

lsst::afw::math::Kernel::Kernel ( int  width,
int  height,
const std::vector< SpatialFunctionPtr spatialFunctionList 
)
explicit

Construct a spatially varying Kernel with a list of spatial functions (one per kernel parameter)

Note: if the list of spatial functions is empty then the kernel is not spatially varying.

Parameters
widthnumber of columns
heightnumber of height
spatialFunctionListlist of spatial function, one per kernel parameter
Exceptions
lsst::pex::exceptions::InvalidParameterErrorif a width or height < 1

Definition at line 107 of file Kernel.cc.

108  : daf::base::Citizen(typeid(this)),
109  _width(width),
110  _height(height),
111  _ctrX(width / 2),
112  _ctrY(height / 2),
113  _nKernelParams(spatialFunctionList.size()) {
114  if ((width < 1) || (height < 1)) {
116  os << "kernel height = " << height << " and/or width = " << width << " < 1";
118  }
119  for (unsigned int ii = 0; ii < spatialFunctionList.size(); ++ii) {
120  SpatialFunctionPtr spatialFunctionCopy = spatialFunctionList[ii]->clone();
121  this->_spatialFunctionList.push_back(spatialFunctionCopy);
122  }
123 }
std::shared_ptr< lsst::afw::math::Function2< double > > SpatialFunctionPtr
Definition: Kernel.h:117
T str(T... args)
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:487
Reports invalid arguments.
Definition: Runtime.h:66
std::ostream * os
Definition: Schema.cc:746

◆ Kernel() [4/5]

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

◆ Kernel() [5/5]

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

◆ ~Kernel()

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

Member Function Documentation

◆ census() [1/3]

int lsst::daf::base::Citizen::census ( int  ,
memId  startingMemId = 0 
)
staticinherited

How many active Citizens are there?

Parameters
startingMemIdDon't print Citizens with lower IDs

Definition at line 238 of file Citizen.cc.

240  {
241  if (startingMemId == 0) { // easy
242  ReadGuard guard(citizenLock);
243  return _activeCitizens.size();
244  }
245 
246  int n = 0;
247  ReadGuard guard(citizenLock);
248  for (table::iterator cur = _activeCitizens.begin(); cur != _activeCitizens.end(); cur++) {
249  if (cur->first->_CitizenId >= startingMemId) {
250  n++;
251  }
252  }
253 
254  return n;
255 }
T end(T... args)
T size(T... args)
T begin(T... args)

◆ census() [2/3]

void lsst::daf::base::Citizen::census ( std::ostream stream,
memId  startingMemId = 0 
)
staticinherited

Print a list of all active Citizens to stream, sorted by ID.

Parameters
streamstream to print to
startingMemIdDon't print Citizens with lower IDs

Definition at line 259 of file Citizen.cc.

261  {
262  ReadGuard guard(citizenLock);
263 
265 
266  for (std::vector<Citizen const*>::const_iterator citizen = leaks->begin(), end = leaks->end();
267  citizen != end; ++citizen) {
268  if ((*citizen)->getId() >= startingMemId) {
269  stream << (*citizen)->repr() << "\n";
270  }
271  }
272 }
STL class.
STL class.
T begin(T... args)
int end
static const std::vector< const Citizen * > * census()
Return a (newly allocated) std::vector of active Citizens sorted by ID.
Definition: Citizen.cc:287

◆ census() [3/3]

std::vector< dafBase::Citizen const * > const * lsst::daf::base::Citizen::census ( )
staticinherited

Return a (newly allocated) std::vector of active Citizens sorted by ID.

You are responsible for deleting it; or you can say std::unique_ptr<std::vector<Citizen const*> const> leaks(Citizen::census()); and not bother (that becomes std::unique_ptr in C++11)

Definition at line 287 of file Citizen.cc.

287  {
289  ReadGuard guard(citizenLock);
290  vec->reserve(_activeCitizens.size());
291 
292  for (table::iterator cur = _activeCitizens.begin(); cur != _activeCitizens.end(); cur++) {
293  vec->push_back(dynamic_cast<Citizen const*>(cur->first));
294  }
295 
296  std::sort(vec->begin(), vec->end(), cmpId);
297 
298  return vec;
299 }
T end(T... args)
T push_back(T... args)
T size(T... args)
STL class.
T begin(T... args)
T sort(T... args)
T reserve(T... args)

◆ clone()

virtual std::shared_ptr<Kernel> lsst::afw::math::Kernel::clone ( ) const
pure virtual

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

Implemented in lsst::afw::math::SeparableKernel, lsst::afw::math::LinearCombinationKernel, lsst::afw::math::DeltaFunctionKernel, lsst::afw::math::AnalyticKernel, lsst::afw::math::FixedKernel, lsst::afw::math::NearestWarpingKernel, lsst::afw::math::BilinearWarpingKernel, and lsst::afw::math::LanczosWarpingKernel.

◆ computeCache()

virtual void lsst::afw::math::Kernel::computeCache ( int  const)
inlinevirtual

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.
Parameters
constdesired cache size

Reimplemented in lsst::afw::math::SeparableKernel.

Definition at line 435 of file Kernel.h.

436  {}

◆ computeImage()

double lsst::afw::math::Kernel::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.

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 93 of file Kernel.cc.

93  {
94  if (image.getDimensions() != this->getDimensions()) {
96  os << "image dimensions = ( " << image.getWidth() << ", " << image.getHeight() << ") != ("
97  << this->getWidth() << ", " << this->getHeight() << ") = kernel dimensions";
99  }
100  image.setXY0(-_ctrX, -_ctrY);
101  if (this->isSpatiallyVarying()) {
103  }
104  return doComputeImage(image, doNormalize);
105 }
int getHeight() const
Return the Kernel&#39;s height.
Definition: Kernel.h:233
int getHeight() const
Return the number of rows in the image.
Definition: ImageBase.h:316
int y
Definition: SpanSet.cc:49
void setXY0(lsst::geom::Point2I const origin)
Set the ImageBase&#39;s origin.
Definition: ImageBase.h:452
bool isSpatiallyVarying() const
Return true iff the kernel is spatially varying (has a spatial function)
Definition: Kernel.h:371
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:237
int getWidth() const
Return the Kernel&#39;s width.
Definition: Kernel.h:228
lsst::geom::Extent2I const getDimensions() const
Return the Kernel&#39;s dimensions (width, height)
Definition: Kernel.h:216
#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:314
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:374
std::ostream * os
Definition: Schema.cc:746

◆ computeKernelParametersFromSpatialModel()

void lsst::afw::math::Kernel::computeKernelParametersFromSpatialModel ( std::vector< double > &  kernelParams,
double  x,
double  y 
) const

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 153 of file Kernel.cc.

154  {
155  std::vector<double>::iterator paramIter = kernelParams.begin();
157  for (; funcIter != _spatialFunctionList.end(); ++funcIter, ++paramIter) {
158  *paramIter = (*(*funcIter))(x, y);
159  }
160 }
int y
Definition: SpanSet.cc:49
double x
STL class.
T begin(T... args)
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:487

◆ doComputeImage()

virtual double lsst::afw::math::Kernel::doComputeImage ( lsst::afw::image::Image< Pixel > &  image,
bool  doNormalize 
) const
protectedpure virtual

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

Implemented in lsst::afw::math::SeparableKernel, lsst::afw::math::LinearCombinationKernel, lsst::afw::math::DeltaFunctionKernel, lsst::afw::math::AnalyticKernel, and lsst::afw::math::FixedKernel.

◆ dynamicCast()

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
py::object result
Definition: schema.cc:418
T dynamic_pointer_cast(T... args)
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48

◆ getBBox()

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

return parent bounding box, with XY0 = -center

Definition at line 257 of file Kernel.h.

257  {
258  return lsst::geom::Box2I(lsst::geom::Point2I(-_ctrX, -_ctrY), lsst::geom::Extent2I(_width, _height));
259  }
An integer coordinate rectangle.
Definition: Box.h:54

◆ getCacheSize()

virtual int lsst::afw::math::Kernel::getCacheSize ( ) const
inlinevirtual

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 441 of file Kernel.h.

441 { return 0; };

◆ getCtr()

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

Return index of kernel's center.

Definition at line 238 of file Kernel.h.

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

◆ getCtrX()

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

Return x index of kernel's center.

Deprecated:
Use getCtr instead

Definition at line 245 of file Kernel.h.

245 { return _ctrX; }

◆ getCtrY()

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

Return y index of kernel's center.

Deprecated:
Use getCtr instead

Definition at line 252 of file Kernel.h.

252 { return _ctrY; }

◆ getDimensions()

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

Return the Kernel's dimensions (width, height)

Definition at line 216 of file Kernel.h.

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

◆ getHeight()

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

Return the Kernel's height.

Definition at line 233 of file Kernel.h.

233 { return _height; }

◆ getId()

dafBase::Citizen::memId lsst::daf::base::Citizen::getId ( ) const
inherited

Return the Citizen's ID.

Definition at line 206 of file Citizen.cc.

206 { return _CitizenId; }

◆ getKernelParameter()

virtual double lsst::afw::math::Kernel::getKernelParameter ( unsigned int  i) const
inlinevirtual

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 295 of file Kernel.h.

295 { return getKernelParameters()[i]; }
virtual std::vector< double > getKernelParameters() const
Return the current kernel parameters.
Definition: Kernel.cc:184

◆ getKernelParameters()

std::vector< double > lsst::afw::math::Kernel::getKernelParameters ( ) const
virtual

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 in lsst::afw::math::SeparableKernel, lsst::afw::math::LinearCombinationKernel, and lsst::afw::math::AnalyticKernel.

Definition at line 184 of file Kernel.cc.

◆ getNextMemId()

dafBase::Citizen::memId lsst::daf::base::Citizen::getNextMemId ( )
staticinherited

Return the memId of the next object to be allocated.

Definition at line 209 of file Citizen.cc.

209 { return _nextMemId(); }

◆ getNKernelParameters()

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

Return the number of kernel parameters (0 if none)

Definition at line 264 of file Kernel.h.

264 { return _nKernelParams; }

◆ getNSpatialParameters()

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

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

Definition at line 269 of file Kernel.h.

269  {
270  return this->isSpatiallyVarying() ? _spatialFunctionList[0]->getNParameters() : 0;
271  }
bool isSpatiallyVarying() const
Return true iff the kernel is spatially varying (has a spatial function)
Definition: Kernel.h:371
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:487

◆ getPersistenceName()

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

◆ getPythonModule()

std::string lsst::afw::math::Kernel::getPythonModule ( ) const
overrideprotectedvirtual

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 244 of file Kernel.cc.

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

◆ getSpatialFunction()

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

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 162 of file Kernel.cc.

162  {
163  if (index >= _spatialFunctionList.size()) {
164  if (!this->isSpatiallyVarying()) {
165  throw LSST_EXCEPT(pexExcept::InvalidParameterError, "kernel is not spatially varying");
166  } else {
167  std::ostringstream errStream;
168  errStream << "index = " << index << "; must be < , " << _spatialFunctionList.size();
170  }
171  }
172  return _spatialFunctionList[index]->clone();
173 }
bool isSpatiallyVarying() const
Return true iff the kernel is spatially varying (has a spatial function)
Definition: Kernel.h:371
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:487
Reports invalid arguments.
Definition: Runtime.h:66

◆ getSpatialFunctionList()

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

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 175 of file Kernel.cc.

175  {
176  std::vector<SpatialFunctionPtr> spFuncCopyList;
178  spFuncIter != _spatialFunctionList.end(); ++spFuncIter) {
179  spFuncCopyList.push_back((**spFuncIter).clone());
180  }
181  return spFuncCopyList;
182 }
T push_back(T... args)
STL class.
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:487

◆ getSpatialParameters()

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

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

Definition at line 359 of file Kernel.h.

359  {
360  std::vector<std::vector<double>> spatialParams;
362  for (; spFuncIter != _spatialFunctionList.end(); ++spFuncIter) {
363  spatialParams.push_back((*spFuncIter)->getParameters());
364  }
365  return spatialParams;
366  }
T push_back(T... args)
STL class.
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:487

◆ getWidth()

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

Return the Kernel's width.

Definition at line 228 of file Kernel.h.

228 { return _width; }

◆ growBBox()

lsst::geom::Box2I lsst::afw::math::Kernel::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.

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 186 of file Kernel.cc.

186  {
187  return lsst::geom::Box2I(
189  lsst::geom::Extent2I(bbox.getDimensions() + getDimensions() - lsst::geom::Extent2I(1, 1)));
190 }
lsst::geom::Point2I getCtr() const
Return index of kernel&#39;s center.
Definition: Kernel.h:238
table::Box2IKey bbox
Definition: Detector.cc:169
lsst::geom::Extent2I const getDimensions() const
Return the Kernel&#39;s dimensions (width, height)
Definition: Kernel.h:216
Extent< int, 2 > Extent2I
Definition: Extent.h:397
An integer coordinate rectangle.
Definition: Box.h:54

◆ hasBeenCorrupted()

bool lsst::daf::base::Citizen::hasBeenCorrupted ( )
staticinherited

Check all allocated blocks for corruption.

Definition at line 316 of file Citizen.cc.

316  {
317  ReadGuard guard(citizenLock);
318  for (table::iterator cur = _activeCitizens.begin(); cur != _activeCitizens.end(); cur++) {
319  if (cur->first->_hasBeenCorrupted()) {
320  return true;
321  }
322  }
323  for (table::iterator cur = _persistentCitizens.begin(); cur != _persistentCitizens.end(); cur++) {
324  if (cur->first->_hasBeenCorrupted()) {
325  return true;
326  }
327  }
328 
329  return false;
330 }
T end(T... args)
T begin(T... args)

◆ init()

int lsst::daf::base::Citizen::init ( )
staticinherited

Called once when the memory system is being initialised.

Definition at line 196 of file Citizen.cc.

196  {
197  volatile int dummy = 1;
198  return dummy;
199 }

◆ isPersistable()

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

◆ isSpatiallyVarying()

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

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

Definition at line 371 of file Kernel.h.

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

◆ markPersistent()

void lsst::daf::base::Citizen::markPersistent ( void  )
inherited

Mark a Citizen as persistent and not destroyed until process end.

Definition at line 225 of file Citizen.cc.

225  {
226  WriteGuard guard(citizenLock);
227  _persistentCitizens[this] = _activeCitizens[this];
228  _activeCitizens.erase(this);
229 }
T erase(T... args)

◆ operator=() [1/2]

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

◆ operator=() [2/2]

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

◆ readFits() [1/3]

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/3]

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() [3/3]

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

◆ repr()

std::string lsst::daf::base::Citizen::repr ( ) const
inherited

Return a string representation of a Citizen.

Definition at line 219 of file Citizen.cc.

219  {
220  return boost::str(boost::format("%d: %08x %s") % _CitizenId % this %
221  lsst::utils::demangleType(_typeName));
222 }
std::string demangleType(std::string const _typeName)
Definition: Demangle.cc:113
def format(config, name=None, writeSourceLine=True, prefix="", verbose=False)
Definition: history.py:168

◆ resized()

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

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.

Implemented in lsst::afw::math::SeparableKernel, lsst::afw::math::LinearCombinationKernel, lsst::afw::math::DeltaFunctionKernel, lsst::afw::math::AnalyticKernel, and lsst::afw::math::FixedKernel.

◆ setCorruptionCallback()

dafBase::Citizen::memCallback lsst::daf::base::Citizen::setCorruptionCallback ( Citizen::memCallback  func)
staticinherited

Set the CorruptionCallback function.

Parameters
funcfunction be called when block is found to be corrupted

Definition at line 391 of file Citizen.cc.

393  {
394  Citizen::memCallback old = _corruptionCallback;
395  _corruptionCallback = func;
396 
397  return old;
398 }
memId(* memCallback)(const Citizen *ptr)
Definition: Citizen.h:61

◆ setCtr()

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

Set index of kernel's center.

Definition at line 330 of file Kernel.h.

330  {
331  _ctrX = ctr.getX();
332  _ctrY = ctr.getY();
333  _setKernelXY();
334  }

◆ setCtrX()

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

Set x index of kernel's center.

Deprecated:
Use setCtr instead

Definition at line 341 of file Kernel.h.

341  {
342  _ctrX = ctrX;
343  _setKernelXY();
344  }

◆ setCtrY()

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

Set y index of kernel's center.

Deprecated:
Use setCtr instead

Definition at line 351 of file Kernel.h.

351  {
352  _ctrY = ctrY;
353  _setKernelXY();
354  }

◆ setDeleteCallback()

dafBase::Citizen::memCallback lsst::daf::base::Citizen::setDeleteCallback ( Citizen::memCallback  func)
staticinherited

Set the DeleteCallback function.

Parameters
funcfunction be called when desired block is deleted

Definition at line 381 of file Citizen.cc.

383  {
384  Citizen::memCallback old = _deleteCallback;
385  _deleteCallback = func;
386 
387  return old;
388 }
memId(* memCallback)(const Citizen *ptr)
Definition: Citizen.h:61

◆ setDeleteCallbackId()

dafBase::Citizen::memId lsst::daf::base::Citizen::setDeleteCallbackId ( Citizen::memId  id)
staticinherited

Call the current DeleteCallback when block is deleted.

Parameters
idDesired ID

Definition at line 347 of file Citizen.cc.

348  {
349  WriteGuard guard(citizenLock);
350  Citizen::memId oldId = _deleteId;
351  _deleteId = id;
352 
353  return oldId;
354 }
unsigned long memId
Type of the block&#39;s ID.
Definition: Citizen.h:58
table::Key< int > id
Definition: Detector.cc:166

◆ setDimensions()

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

Definition at line 218 of file Kernel.h.

218  {
219  _width = dims.getX();
220  _height = dims.getY();
221  }

◆ setHeight()

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

Definition at line 223 of file Kernel.h.

223 { _height = height; }

◆ setKernelParameter()

void lsst::afw::math::Kernel::setKernelParameter ( unsigned int  ind,
double  value 
) const
protectedvirtual

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 in lsst::afw::math::SeparableKernel, lsst::afw::math::LinearCombinationKernel, lsst::afw::math::AnalyticKernel, lsst::afw::math::NearestWarpingKernel, lsst::afw::math::BilinearWarpingKernel, and lsst::afw::math::LanczosWarpingKernel.

Definition at line 233 of file Kernel.cc.

233  {
234  throw LSST_EXCEPT(pexExcept::InvalidParameterError, "Kernel has no kernel parameters");
235 }
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
Reports invalid arguments.
Definition: Runtime.h:66

◆ setKernelParameters() [1/2]

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

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 379 of file Kernel.h.

379  {
380  if (this->isSpatiallyVarying()) {
381  throw LSST_EXCEPT(lsst::pex::exceptions::RuntimeError, "Kernel is spatially varying");
382  }
383  const unsigned int nParams = this->getNKernelParameters();
384  if (nParams != params.size()) {
386  (boost::format("Number of parameters is wrong, saw %d expected %d") % nParams %
387  params.size())
388  .str());
389  }
390  for (unsigned int ii = 0; ii < nParams; ++ii) {
391  this->setKernelParameter(ii, params[ii]);
392  }
393  }
virtual void setKernelParameter(unsigned int ind, double value) const
Set one kernel parameter.
Definition: Kernel.cc:233
bool isSpatiallyVarying() const
Return true iff the kernel is spatially varying (has a spatial function)
Definition: Kernel.h:371
def format(config, name=None, writeSourceLine=True, prefix="", verbose=False)
Definition: history.py:168
unsigned int getNKernelParameters() const
Return the number of kernel parameters (0 if none)
Definition: Kernel.h:264
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)
inline

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 401 of file Kernel.h.

401  {
402  this->setKernelParameter(0, params.first);
403  this->setKernelParameter(1, params.second);
404  }
virtual void setKernelParameter(unsigned int ind, double value) const
Set one kernel parameter.
Definition: Kernel.cc:233

◆ setKernelParametersFromSpatialModel()

void lsst::afw::math::Kernel::setKernelParametersFromSpatialModel ( double  x,
double  y 
) const
protected

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 237 of file Kernel.cc.

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

◆ setNewCallback()

dafBase::Citizen::memNewCallback lsst::daf::base::Citizen::setNewCallback ( Citizen::memNewCallback  func)
staticinherited

Set the NewCallback function.

Parameters
funcThe new function to be called when a designated block is allocated

Definition at line 371 of file Citizen.cc.

373  {
374  Citizen::memNewCallback old = _newCallback;
375  _newCallback = func;
376 
377  return old;
378 }
memId(* memNewCallback)(const memId cid)
A function used to register a callback.
Definition: Citizen.h:60

◆ setNewCallbackId()

dafBase::Citizen::memId lsst::daf::base::Citizen::setNewCallbackId ( Citizen::memId  id)
staticinherited

Call the NewCallback when block is allocated.

Parameters
idDesired ID

Definition at line 337 of file Citizen.cc.

338  {
339  WriteGuard guard(citizenLock);
340  Citizen::memId oldId = _newId;
341  _newId = id;
342 
343  return oldId;
344 }
unsigned long memId
Type of the block&#39;s ID.
Definition: Citizen.h:58
table::Key< int > id
Definition: Detector.cc:166

◆ setSpatialParameters()

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

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 128 of file Kernel.cc.

128  {
129  // Check params size before changing anything
130  unsigned int nKernelParams = this->getNKernelParameters();
131  if (params.size() != nKernelParams) {
132  throw LSST_EXCEPT(
134  (boost::format("params has %d entries instead of %d") % params.size() % nKernelParams).str());
135  }
136  unsigned int nSpatialParams = this->getNSpatialParameters();
137  for (unsigned int ii = 0; ii < nKernelParams; ++ii) {
138  if (params[ii].size() != nSpatialParams) {
140  (boost::format("params[%d] has %d entries instead of %d") % ii %
141  params[ii].size() % nSpatialParams)
142  .str());
143  }
144  }
145  // Set parameters
146  if (nSpatialParams > 0) {
147  for (unsigned int ii = 0; ii < nKernelParams; ++ii) {
148  this->_spatialFunctionList[ii]->setParameters(params[ii]);
149  }
150  }
151 }
def format(config, name=None, writeSourceLine=True, prefix="", verbose=False)
Definition: history.py:168
unsigned int getNKernelParameters() const
Return the number of kernel parameters (0 if none)
Definition: Kernel.h:264
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:487
int getNSpatialParameters() const
Return the number of spatial parameters (0 if not spatially varying)
Definition: Kernel.h:269
Reports invalid arguments.
Definition: Runtime.h:66

◆ setWidth()

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

Definition at line 222 of file Kernel.h.

222 { _width = width; }

◆ shrinkBBox()

lsst::geom::Box2I lsst::afw::math::Kernel::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.

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 192 of file Kernel.cc.

192  {
193  if ((bbox.getWidth() < getWidth()) || ((bbox.getHeight() < getHeight()))) {
195  os << "bbox dimensions = " << bbox.getDimensions() << " < (" << getWidth() << ", " << getHeight()
196  << ") in one or both dimensions";
198  }
199  return lsst::geom::Box2I(
200  lsst::geom::Point2I(bbox.getMinX() + getCtrX(), bbox.getMinY() + getCtrY()),
201  lsst::geom::Extent2I(bbox.getWidth() + 1 - getWidth(), bbox.getHeight() + 1 - getHeight()));
202 }
int getHeight() const
Return the Kernel&#39;s height.
Definition: Kernel.h:233
int getCtrX() const
Return x index of kernel&#39;s center.
Definition: Kernel.h:245
int getCtrY() const
Return y index of kernel&#39;s center.
Definition: Kernel.h:252
T str(T... args)
table::Box2IKey bbox
Definition: Detector.cc:169
int getWidth() const
Return the Kernel&#39;s width.
Definition: Kernel.h:228
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
Reports invalid arguments.
Definition: Runtime.h:66
An integer coordinate rectangle.
Definition: Box.h:54
std::ostream * os
Definition: Schema.cc:746

◆ toString()

std::string lsst::afw::math::Kernel::toString ( std::string const &  prefix = "") const
virtual

Return a string representation of the kernel.

Reimplemented in lsst::afw::math::SeparableKernel, lsst::afw::math::LinearCombinationKernel, lsst::afw::math::DeltaFunctionKernel, lsst::afw::math::AnalyticKernel, and lsst::afw::math::FixedKernel.

Definition at line 204 of file Kernel.cc.

204  {
206  os << prefix << "Kernel:" << std::endl;
207  os << prefix << "..height, width: " << _height << ", " << _width << std::endl;
208  os << prefix << "..ctr (X, Y): " << _ctrX << ", " << _ctrY << std::endl;
209  os << prefix << "..nKernelParams: " << _nKernelParams << std::endl;
210  os << prefix << "..isSpatiallyVarying: " << (this->isSpatiallyVarying() ? "True" : "False") << std::endl;
211  if (this->isSpatiallyVarying()) {
212  os << prefix << "..spatialFunctions:" << std::endl;
214  spFuncPtr != _spatialFunctionList.end(); ++spFuncPtr) {
215  os << prefix << "...." << (*spFuncPtr)->toString() << std::endl;
216  }
217  }
218  return os.str();
219 }
T endl(T... args)
std::string prefix
Definition: SchemaMapper.cc:79
bool isSpatiallyVarying() const
Return true iff the kernel is spatially varying (has a spatial function)
Definition: Kernel.h:371
T str(T... args)
STL class.
std::vector< SpatialFunctionPtr > _spatialFunctionList
Definition: Kernel.h:487
std::ostream * os
Definition: Schema.cc:746

◆ 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::Chebyshev1Function2< ReturnT >, lsst::afw::math::LinearCombinationKernel, lsst::afw::math::DeltaFunctionKernel, lsst::afw::math::AnalyticKernel, lsst::afw::image::PhotoCalib, lsst::afw::math::PolynomialFunction2< ReturnT >, lsst::afw::math::FixedKernel, 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::Camera, lsst::afw::detection::HeavyFootprint< ImagePixelT, MaskPixelT, VariancePixelT >, lsst::afw::cameraGeom::DetectorCollection, lsst::meas::algorithms::CoaddBoundedField, lsst::afw::image::CoaddInputs, lsst::afw::detection::GaussianPsf, lsst::afw::math::TransformBoundedField, 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
protected

Definition at line 487 of file Kernel.h.


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