lsst::afw::math::LanczosWarpingKernel Class Reference

Lanczos warping: accurate but slow and can introduce ringing artifacts. More...

#include <warpExposure.h>

Inheritance diagram for lsst::afw::math::LanczosWarpingKernel:

Public Types
using	KernelFunction = lsst::afw::math::Function1< Pixel >
using	KernelFunctionPtr = std::shared_ptr< KernelFunction >
using	Pixel = double
using	SpatialFunctionPtr = std::shared_ptr< lsst::afw::math::Function2< double > >
using	SpatialFunction = lsst::afw::math::Function2< double >
using	NullSpatialFunction = lsst::afw::math::NullFunction2< double >
using	kernel_fill_factor = generic_kernel_tag

Public Member Functions
	LanczosWarpingKernel (int order)
	LanczosWarpingKernel (const LanczosWarpingKernel &)=delete
	LanczosWarpingKernel (LanczosWarpingKernel &&)=delete
LanczosWarpingKernel &	operator= (const LanczosWarpingKernel &)=delete
LanczosWarpingKernel &	operator= (LanczosWarpingKernel &&)=delete
	~LanczosWarpingKernel () override=default
std::shared_ptr< Kernel >	clone () const override
	Return a pointer to a deep copy of this kernel. More...
int	getOrder () const
	get the order of the kernel More...
bool	isPersistable () const noexcept override
	Return true if this particular object can be persisted using afw::table::io. More...
std::shared_ptr< Kernel >	resized (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...
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< SpatialFunctionPtr >	getSpatialFunctionList () 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...
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...

Static Public Member Functions
static std::shared_ptr< SeparableKernel >	readFits (fits::Fits &fitsfile)
	Read an object from an already open FITS object. More...
static std::shared_ptr< SeparableKernel >	readFits (std::string const &fileName, int hdu=fits::DEFAULT_HDU)
	Read an object from a regular FITS file. More...
static std::shared_ptr< SeparableKernel >	readFits (fits::MemFileManager &manager, int hdu=fits::DEFAULT_HDU)
	Read an object from a FITS file in memory. More...
static std::shared_ptr< Kernel >	readFits (fits::Fits &fitsfile)
	Read an object from an already open FITS object. More...
static std::shared_ptr< Kernel >	readFits (std::string const &fileName, int hdu=fits::DEFAULT_HDU)
	Read an object from a regular FITS file. More...
static std::shared_ptr< Kernel >	readFits (fits::MemFileManager &manager, int hdu=fits::DEFAULT_HDU)
	Read an object from a FITS file in memory. More...
static std::shared_ptr< SeparableKernel >	dynamicCast (std::shared_ptr< Persistable > const &ptr)
	Dynamically cast a shared_ptr. More...
static std::shared_ptr< Kernel >	dynamicCast (std::shared_ptr< Persistable > const &ptr)
	Dynamically cast a shared_ptr. More...

Protected Types
using	OutputArchiveHandle = io::OutputArchiveHandle

Protected Member Functions
void	setKernelParameter (unsigned int ind, double value) const override
	Set one kernel parameter. More...
std::string	getPersistenceName () const override
	Return the unique name used to persist this object and look up its factory. 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	write (OutputArchiveHandle &handle) const override
	Write the object to one or more catalogs. More...
double	doComputeImage (lsst::afw::image::Image< Pixel > &image, bool doNormalize) const override
	Low-level version of computeImage. More...
void	setKernelParametersFromSpatialModel (double x, double y) const
	Set the kernel parameters from the spatial model (if any). More...

Protected Attributes
std::vector< SpatialFunctionPtr >	_spatialFunctionList

Detailed Description

Lanczos warping: accurate but slow and can introduce ringing artifacts.

This kernel is the product of two 1-dimensional Lanczos functions. The number of minima and maxima in the 1-dimensional Lanczos function is 2*order + 1. The kernel has one pixel per function minimum or maximum; but as applied to warping, the first or last pixel is always zero and can be omitted. Thus the kernel size is 2*order x 2*order.

For more information about warping kernels see makeWarpingKernel

Todo:

: make a new class WarpingKernel and make this a subclass.

Definition at line 66 of file warpExposure.h.

Member Typedef Documentation

kernel_fill_factor

using lsst::afw::math::Kernel::kernel_fill_factor = generic_kernel_tag

inherited

Definition at line 118 of file Kernel.h.

KernelFunction

using lsst::afw::math::SeparableKernel::KernelFunction = lsst::afw::math::Function1<Pixel>

inherited

Definition at line 862 of file Kernel.h.

KernelFunctionPtr

using lsst::afw::math::SeparableKernel::KernelFunctionPtr = std::shared_ptr<KernelFunction>

inherited

Definition at line 863 of file Kernel.h.

NullSpatialFunction

using lsst::afw::math::Kernel::NullSpatialFunction = lsst::afw::math::NullFunction2<double>

inherited

Definition at line 115 of file Kernel.h.

OutputArchiveHandle

using lsst::afw::table::io::Persistable::OutputArchiveHandle = io::OutputArchiveHandle

protectedinherited

Definition at line 108 of file Persistable.h.

Pixel

using lsst::afw::math::Kernel::Pixel = double

inherited

Definition at line 112 of file Kernel.h.

SpatialFunction

using lsst::afw::math::Kernel::SpatialFunction = lsst::afw::math::Function2<double>

inherited

Definition at line 114 of file Kernel.h.

SpatialFunctionPtr

using lsst::afw::math::Kernel::SpatialFunctionPtr = std::shared_ptr<lsst::afw::math::Function2<double> >

inherited

Definition at line 113 of file Kernel.h.

Constructor & Destructor Documentation

LanczosWarpingKernel() [1/3]

lsst::afw::math::LanczosWarpingKernel::LanczosWarpingKernel ( int  order )

inlineexplicit

Parameters

order

order of Lanczos function

Definition at line 68 of file warpExposure.h.

            : SeparableKernel(2 * order, 2 * order, LanczosFunction1<Kernel::Pixel>(order),
                              LanczosFunction1<Kernel::Pixel>(order)) {}

LanczosWarpingKernel() [2/3]

lsst::afw::math::LanczosWarpingKernel::LanczosWarpingKernel ( const LanczosWarpingKernel &  )

delete

LanczosWarpingKernel() [3/3]

lsst::afw::math::LanczosWarpingKernel::LanczosWarpingKernel ( LanczosWarpingKernel &&  )

delete

~LanczosWarpingKernel()

lsst::afw::math::LanczosWarpingKernel::~LanczosWarpingKernel ( )

overridedefault

Member Function Documentation

clone()

std::shared_ptr< Kernel > lsst::afw::math::LanczosWarpingKernel::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

Reimplemented from lsst::afw::math::SeparableKernel.

Definition at line 97 of file warpExposure.cc.

                                                        {
    return std::shared_ptr<Kernel>(new LanczosWarpingKernel(this->getOrder()));
}

computeCache()

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

overridevirtualinherited

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 309 of file SeparableKernel.cc.

                                                      {
    SeparableKernel::KernelFunctionPtr func;
 
    func = getKernelColFunction();
    _computeCache(cacheSize, _kernelY, func, &_kernelColCache);
 
    func = getKernelRowFunction();
    _computeCache(cacheSize, _kernelX, func, &_kernelRowCache);
}

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

image	image whose pixels are to be set (output); xy0 of the image will be set to -kernel.getCtr()
doNormalize	normalize the image (so sum is 1)?
x	x (column position) at which to compute spatial function
y	y (row position) at which to compute spatial function

Returns

The kernel sum

Exceptions

lsst::pex::exceptions::InvalidParameterError	if the image is the wrong size
lsst::pex::exceptions::OverflowError	if doNormalize is true and the kernel sum is exactly 0

Note

computeNewImage has been retired; it doesn't need to be a member

Definition at line 76 of file Kernel.cc.

                                                                                                {
    if (image.getDimensions() != this->getDimensions()) {
        std::ostringstream os;
        os << "image dimensions = ( " << image.getWidth() << ", " << image.getHeight() << ") != ("
           << this->getWidth() << ", " << this->getHeight() << ") = kernel dimensions";
        throw LSST_EXCEPT(pexExcept::InvalidParameterError, os.str());
    }
    image.setXY0(-_ctrX, -_ctrY);
    if (this->isSpatiallyVarying()) {
        this->setKernelParametersFromSpatialModel(x, y);
    }
    return doComputeImage(image, doNormalize);
}

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

                                                                     {
    std::vector<double>::iterator paramIter = kernelParams.begin();
    std::vector<SpatialFunctionPtr>::const_iterator funcIter = _spatialFunctionList.begin();
    for (; funcIter != _spatialFunctionList.end(); ++funcIter, ++paramIter) {
        *paramIter = (*(*funcIter))(x, y);
    }
}

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

inherited

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

colList	column vector
rowList	row vector
doNormalize	normalize the image (so sum of each is 1)?
x	x (column position) at which to compute spatial function
y	y (row position) at which to compute spatial function

Returns

the kernel sum (1.0 if doNormalize true)

Exceptions

lsst::pex::exceptions::InvalidParameterError	if colList or rowList is the wrong size
lsst::pex::exceptions::OverflowError	if doNormalize is true and the kernel sum is exactly 0

Definition at line 122 of file SeparableKernel.cc.

                                                                                   {
    if (static_cast<int>(colList.size()) != this->getWidth() ||
        static_cast<int>(rowList.size()) != this->getHeight()) {
        std::ostringstream os;
        os << "colList.size(), rowList.size() = (" << colList.size() << ", " << rowList.size() << ") != ("
           << this->getWidth() << ", " << this->getHeight() << ") = "
           << "kernel dimensions";
        throw LSST_EXCEPT(pexExcept::InvalidParameterError, os.str());
    }
    if (this->isSpatiallyVarying()) {
        this->setKernelParametersFromSpatialModel(x, y);
    }
 
    return basicComputeVectors(colList, rowList, doNormalize);
}

doComputeImage()

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

overrideprotectedvirtualinherited

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

image	image whose pixels are to be set (output)
doNormalize	normalize the image (so sum is 1)?

Returns

The kernel sum

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

Definition at line 171 of file SeparableKernel.cc.

                                                                                       {
    double imSum = basicComputeVectors(_localColList, _localRowList, doNormalize);
 
    for (int y = 0; y != image.getHeight(); ++y) {
        image::Image<Pixel>::x_iterator imPtr = image.row_begin(y);
        for (std::vector<Pixel>::iterator colIter = _localColList.begin(); colIter != _localColList.end();
             ++colIter, ++imPtr) {
            *imPtr = (*colIter) * _localRowList[y];
        }
    }
 
    return imSum;
}

dynamicCast() [1/2]

template std::shared_ptr< meas::extensions::psfex::PsfexPsf > lsst::afw::table::io::PersistableFacade< Kernel >::dynamicCast ( std::shared_ptr< Persistable > const &  ptr )

staticinherited

Dynamically cast a shared_ptr.

Dynamically cast a shared pointer and raise on failure.

You must provide an explicit template instantiation in the .cc file for each class that inherits from PersistableFacade. Designed to work around RTTI issues on macOS with hidden symbols;

Exceptions

lsst::pex::exceptions::LogicError

if the cast fails

param[in] ptr The pointer to be cast.

Returns

The cast pointer.

Exceptions

lsst::pex::exceptions::TypeError

If the dynamic cast fails.

Definition at line 218 of file Persistable.cc.

                                                                                        {
    auto result = std::dynamic_pointer_cast<T>(ptr);
    if (!result) {
        throw LSST_EXCEPT(pex::exceptions::TypeError, "Dynamic pointer cast failed");
    }
    return result;
}

dynamicCast() [2/2]

template std::shared_ptr< meas::extensions::psfex::PsfexPsf > lsst::afw::table::io::PersistableFacade< SeparableKernel >::dynamicCast ( std::shared_ptr< Persistable > const &  ptr )

staticinherited

Dynamically cast a shared_ptr.

Dynamically cast a shared pointer and raise on failure.

You must provide an explicit template instantiation in the .cc file for each class that inherits from PersistableFacade. Designed to work around RTTI issues on macOS with hidden symbols;

Exceptions

lsst::pex::exceptions::LogicError

if the cast fails

param[in] ptr The pointer to be cast.

Returns

The cast pointer.

Exceptions

lsst::pex::exceptions::TypeError

If the dynamic cast fails.

Definition at line 218 of file Persistable.cc.

                                                                                        {
    auto result = std::dynamic_pointer_cast<T>(ptr);
    if (!result) {
        throw LSST_EXCEPT(pex::exceptions::TypeError, "Dynamic pointer cast failed");
    }
    return result;
}

getBBox()

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

inlineinherited

return parent bounding box, with XY0 = -center

Definition at line 239 of file Kernel.h.

                                         {
        return lsst::geom::Box2I(lsst::geom::Point2I(-_ctrX, -_ctrY), lsst::geom::Extent2I(_width, _height));
    }

getCacheSize()

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

overridevirtualinherited

Get the current cache size (0 if none)

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

Definition at line 319 of file SeparableKernel.cc.

{ return _kernelColCache.size(); };

getCtr()

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

inlineinherited

Return index of kernel's center.

Definition at line 234 of file Kernel.h.

{ return lsst::geom::Point2I(_ctrX, _ctrY); }

getDimensions()

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

inlineinherited

Return the Kernel's dimensions (width, height)

Definition at line 212 of file Kernel.h.

{ return lsst::geom::Extent2I(_width, _height); }

getHeight()

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

inlineinherited

Return the Kernel's height.

Definition at line 229 of file Kernel.h.

{ return _height; }

getKernelColFunction()

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

inherited

Get a deep copy of the col kernel function.

Definition at line 139 of file SeparableKernel.cc.

                                                                             {
    return _kernelColFunctionPtr->clone();
}

getKernelParameter()

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

inlineoverridevirtualinherited

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

                                                             {
        unsigned int const ncol = _kernelColFunctionPtr->getNParameters();
        if (i < ncol) {
            return _kernelColFunctionPtr->getParameter(i);
        } else {
            i -= ncol;
            return _kernelRowFunctionPtr->getParameter(i);
        }
    }

getKernelParameters()

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

overridevirtualinherited

Return the current kernel parameters.

If the kernel is spatially varying then the parameters are those last computed. See also computeKernelParametersFromSpatialModel. If there are no kernel parameters then returns an empty vector.

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

Definition at line 160 of file SeparableKernel.cc.

                                                             {
    std::vector<double> allParams = _kernelColFunctionPtr->getParameters();
    std::vector<double> yParams = _kernelRowFunctionPtr->getParameters();
    std::copy(yParams.begin(), yParams.end(), std::back_inserter(allParams));
    return allParams;
}

getKernelRowFunction()

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

inherited

Get a deep copy of the row kernel function.

Definition at line 143 of file SeparableKernel.cc.

                                                                             {
    return _kernelRowFunctionPtr->clone();
}

getNKernelParameters()

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

inlineinherited

Return the number of kernel parameters (0 if none)

Definition at line 246 of file Kernel.h.

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

                                             {
        return this->isSpatiallyVarying() ? _spatialFunctionList[0]->getNParameters() : 0;
    }

getOrder()

int lsst::afw::math::LanczosWarpingKernel::getOrder

(

)

const

get the order of the kernel

Definition at line 101 of file warpExposure.cc.

{ return this->getWidth() / 2; }

getPersistenceName()

std::string lsst::afw::math::LanczosWarpingKernel::getPersistenceName ( ) const

inlineoverrideprotectedvirtual

Return the unique name used to persist this object and look up its factory.

Must be less than ArchiveIndexSchema::MAX_NAME_LENGTH characters.

Reimplemented from lsst::afw::table::io::Persistable.

Definition at line 92 of file warpExposure.h.

{ return "LanczosWarpingKernel"; }

getPythonModule()

std::string lsst::afw::math::LanczosWarpingKernel::getPythonModule ( ) const

inlineoverrideprotectedvirtual

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

Definition at line 94 of file warpExposure.h.

{ 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

index

index 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::InvalidParameterError	if kernel not spatially varying
lsst::pex::exceptions::InvalidParameterError	if index out of range

Definition at line 144 of file Kernel.cc.

                                                                            {
    if (index >= _spatialFunctionList.size()) {
        if (!this->isSpatiallyVarying()) {
            throw LSST_EXCEPT(pexExcept::InvalidParameterError, "kernel is not spatially varying");
        } else {
            std::ostringstream errStream;
            errStream << "index = " << index << "; must be < , " << _spatialFunctionList.size();
            throw LSST_EXCEPT(pexExcept::InvalidParameterError, errStream.str());
        }
    }
    return _spatialFunctionList[index]->clone();
}

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

                                                                         {
    std::vector<SpatialFunctionPtr> spFuncCopyList;
    for (auto const &spFuncIter : _spatialFunctionList) {
        spFuncCopyList.push_back((*spFuncIter).clone());
    }
    return spFuncCopyList;
}

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

                                                                     {
        std::vector<std::vector<double>> spatialParams;
        std::vector<SpatialFunctionPtr>::const_iterator spFuncIter = _spatialFunctionList.begin();
        for (; spFuncIter != _spatialFunctionList.end(); ++spFuncIter) {
            spatialParams.push_back((*spFuncIter)->getParameters());
        }
        return spatialParams;
    }

getWidth()

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

inlineinherited

Return the Kernel's width.

Definition at line 224 of file Kernel.h.

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

                                                              {
    return lsst::geom::Box2I(
            lsst::geom::Point2I(bbox.getMin() - lsst::geom::Extent2I(getCtr())),
            lsst::geom::Extent2I(bbox.getDimensions() + getDimensions() - lsst::geom::Extent2I(1, 1)));
}

isPersistable()

bool lsst::afw::math::LanczosWarpingKernel::isPersistable ( ) const

inlineoverridevirtualnoexcept

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

Reimplemented from lsst::afw::table::io::Persistable.

Definition at line 87 of file warpExposure.h.

{ return true; }

isSpatiallyVarying()

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

inlineinherited

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

Definition at line 333 of file Kernel.h.

{ return _spatialFunctionList.size() != 0; }

operator=() [1/2]

LanczosWarpingKernel& lsst::afw::math::LanczosWarpingKernel::operator= ( const LanczosWarpingKernel &  )

delete

operator=() [2/2]

LanczosWarpingKernel& lsst::afw::math::LanczosWarpingKernel::operator= ( LanczosWarpingKernel &&  )

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]

fitsfile

FITS object to read from, already positioned at the desired HDU.

Definition at line 183 of file Persistable.h.

                                                         {
        return dynamicCast(Persistable::_readFits(fitsfile));
    }

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]

fitsfile

FITS object to read from, already positioned at the desired HDU.

Definition at line 183 of file Persistable.h.

                                                         {
        return dynamicCast(Persistable::_readFits(fitsfile));
    }

readFits() [3/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]	manager	Manager for the memory to read from.
[in]	hdu	HDU to read, where 0 is the primary. The special value of afw::fits::DEFAULT_HDU skips the primary HDU if it is empty.

Definition at line 205 of file Persistable.h.

                                                                                             {
        return dynamicCast(Persistable::_readFits(manager, hdu));
    }

readFits() [4/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]	manager	Manager for the memory to read from.
[in]	hdu	HDU to read, where 0 is the primary. The special value of afw::fits::DEFAULT_HDU skips the primary HDU if it is empty.

Definition at line 205 of file Persistable.h.

                                                                                             {
        return dynamicCast(Persistable::_readFits(manager, hdu));
    }

readFits() [5/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]	fileName	Name of the file to read.
[in]	hdu	HDU to read, where 0 is the primary. The special value of afw::fits::DEFAULT_HDU skips the primary HDU if it is empty.

Definition at line 194 of file Persistable.h.

                                                                                           {
        return dynamicCast(Persistable::_readFits(fileName, hdu));
    }

readFits() [6/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]	fileName	Name of the file to read.
[in]	hdu	HDU to read, where 0 is the primary. The special value of afw::fits::DEFAULT_HDU skips the primary HDU if it is empty.

Definition at line 194 of file Persistable.h.

                                                                                           {
        return dynamicCast(Persistable::_readFits(fileName, hdu));
    }

resized()

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

overridevirtualinherited

Return a pointer to a clone with specified kernel dimensions.

Parameters

width	Number of columns in pixels
height	Number 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 110 of file SeparableKernel.cc.

                                                                          {
    std::shared_ptr<Kernel> retPtr;
    if (isSpatiallyVarying()) {
        retPtr = std::make_shared<SeparableKernel>(width, height, *_kernelColFunctionPtr,
                                                   *_kernelRowFunctionPtr, _spatialFunctionList);
    } else {
        retPtr = std::make_shared<SeparableKernel>(width, height, *_kernelColFunctionPtr,
                                                   *_kernelRowFunctionPtr);
    }
    return retPtr;
}

setCtr()

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

inlineinherited

Set index of kernel's center.

Definition at line 312 of file Kernel.h.

                                            {
        _ctrX = ctr.getX();
        _ctrY = ctr.getY();
        _setKernelXY();
    }

setDimensions()

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

inlineinherited

Definition at line 214 of file Kernel.h.

                                              {
        _width = dims.getX();
        _height = dims.getY();
    }

setHeight()

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

inlineinherited

Definition at line 219 of file Kernel.h.

{ _height = height; }

setKernelParameter()

void lsst::afw::math::LanczosWarpingKernel::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::InvalidParameterError

always (unless subclassed)

Reimplemented from lsst::afw::math::SeparableKernel.

Definition at line 103 of file warpExposure.cc.

                                                                                  {
    checkWarpingKernelParameter(this, ind, value);
    SeparableKernel::setKernelParameter(ind, value);
}

setKernelParameters() [1/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 363 of file Kernel.h.

                                                                           {
        this->setKernelParameter(0, params.first);
        this->setKernelParameter(1, params.second);
    }

setKernelParameters() [2/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::RuntimeError	if the kernel has a spatial function
lsst::pex::exceptions::InvalidParameterError	if the params vector is the wrong length

Definition at line 341 of file Kernel.h.

                                                                     {
        if (this->isSpatiallyVarying()) {
            throw LSST_EXCEPT(lsst::pex::exceptions::RuntimeError, "Kernel is spatially varying");
        }
        const unsigned int nParams = this->getNKernelParameters();
        if (nParams != params.size()) {
            throw LSST_EXCEPT(lsst::pex::exceptions::InvalidParameterError,
                              (boost::format("Number of parameters is wrong, saw %d expected %d") % nParams %
                               params.size())
                                      .str());
        }
        for (unsigned int ii = 0; ii < nParams; ++ii) {
            this->setKernelParameter(ii, params[ii]);
        }
    }

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

                                                                         {
    std::vector<SpatialFunctionPtr>::const_iterator funcIter = _spatialFunctionList.begin();
    for (int ii = 0; funcIter != _spatialFunctionList.end(); ++funcIter, ++ii) {
        this->setKernelParameter(ii, (*(*funcIter))(x, y));
    }
}

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::InvalidParameterError

if params is the wrong shape (if this exception is thrown then no parameters are changed)

Definition at line 110 of file Kernel.cc.

                                                                            {
    // Check params size before changing anything
    unsigned int nKernelParams = this->getNKernelParameters();
    if (params.size() != nKernelParams) {
        throw LSST_EXCEPT(
                pexExcept::InvalidParameterError,
                (boost::format("params has %d entries instead of %d") % params.size() % nKernelParams).str());
    }
    unsigned int nSpatialParams = this->getNSpatialParameters();
    for (unsigned int ii = 0; ii < nKernelParams; ++ii) {
        if (params[ii].size() != nSpatialParams) {
            throw LSST_EXCEPT(pexExcept::InvalidParameterError,
                              (boost::format("params[%d] has %d entries instead of %d") % ii %
                               params[ii].size() % nSpatialParams)
                                      .str());
        }
    }
    // Set parameters
    if (nSpatialParams > 0) {
        for (unsigned int ii = 0; ii < nKernelParams; ++ii) {
            this->_spatialFunctionList[ii]->setParameters(params[ii]);
        }
    }
}

setWidth()

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

inlineinherited

Definition at line 218 of file Kernel.h.

{ _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::InvalidParameterError

if the resulting box would have dimension < 1 in either axis

Definition at line 173 of file Kernel.cc.

                                                                {
    if ((bbox.getWidth() < getWidth()) || ((bbox.getHeight() < getHeight()))) {
        std::ostringstream os;
        os << "bbox dimensions = " << bbox.getDimensions() << " < (" << getWidth() << ", " << getHeight()
           << ") in one or both dimensions";
        throw LSST_EXCEPT(pexExcept::InvalidParameterError, os.str());
    }
    return lsst::geom::Box2I(
            lsst::geom::Point2I(bbox.getMin() + lsst::geom::Extent2I(getCtr())),
            lsst::geom::Extent2I(bbox.getWidth() + 1 - getWidth(), bbox.getHeight() + 1 - getHeight()));
}

toString()

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

overridevirtualinherited

Return a string representation of the kernel.

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

Definition at line 147 of file SeparableKernel.cc.

                                                                 {
    std::ostringstream os;
    os << prefix << "SeparableKernel:" << std::endl;
    os << prefix
       << "..x (width) function: " << (_kernelColFunctionPtr ? _kernelColFunctionPtr->toString() : "None")
       << std::endl;
    os << prefix
       << "..y (rows) function: " << (_kernelRowFunctionPtr ? _kernelRowFunctionPtr->toString() : "None")
       << std::endl;
    os << Kernel::toString(prefix + "\t");
    return os.str();
}

write()

void lsst::afw::math::LanczosWarpingKernel::write ( OutputArchiveHandle &  handle ) const

overrideprotectedvirtual

Write the object to one or more catalogs.

The handle object passed to this function provides an interface for adding new catalogs and adding nested objects to the same archive (while checking for duplicates). See OutputArchiveHandle for more information.

Reimplemented from lsst::afw::table::io::Persistable.

Definition at line 212 of file warpExposure.cc.

                                                                  {
    auto const &keys = LanczosKernelPersistenceHelper::get();
    table::BaseCatalog catalog = handle.makeCatalog(keys.schema);
    std::shared_ptr<table::BaseRecord> record = catalog.addNew();
    record->set(keys.order, getOrder());
    handle.saveCatalog(catalog);
}

writeFits() [1/3]

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

inherited

Write the object to an already-open FITS object.

Parameters

[in]

fitsfile

Open FITS object to write to.

Definition at line 18 of file Persistable.cc.

                                                    {
    OutputArchive archive;
    archive.put(this);
    archive.writeFits(fitsfile);
}

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]	manager	Name of the file to write to.
[in]	mode	If "w", any existing file with the given name will be overwritten. If "a", new HDUs will be appended to an existing file.

Definition at line 29 of file Persistable.cc.

                                                                                    {
    fits::Fits fitsfile(manager, mode, fits::Fits::AUTO_CLOSE | fits::Fits::AUTO_CHECK);
    writeFits(fitsfile);
}

writeFits() [3/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]	fileName	Name of the file to write to.
[in]	mode	If "w", any existing file with the given name will be overwritten. If "a", new HDUs will be appended to an existing file.

Definition at line 24 of file Persistable.cc.

                                                                                  {
    fits::Fits fitsfile(fileName, mode, fits::Fits::AUTO_CLOSE | fits::Fits::AUTO_CHECK);
    writeFits(fitsfile);
}

Member Data Documentation

_spatialFunctionList

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

protectedinherited

Definition at line 449 of file Kernel.h.

---

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

• /j/snowflake/release/lsstsw/stack/lsst-scipipe-0.7.0/Linux64/afw/22.0.1-42-gca6935d93+ba5e5ca3eb/include/lsst/afw/math/warpExposure.h
• /j/snowflake/release/lsstsw/stack/lsst-scipipe-0.7.0/Linux64/afw/22.0.1-42-gca6935d93+ba5e5ca3eb/src/math/warpExposure.cc