LSST Applications 24.1.6,g063fba187b+56b85ce14a,g0f08755f38+df8a265115,g12f32b3c4e+891a09f10d,g1524ad2192+7a5d7b3fbd,g1653933729+a8ce1bb630,g168dd56ebc+a8ce1bb630,g28da252d5a+07cb1400be,g2bbee38e9b+ae03bbfc84,g2bc492864f+ae03bbfc84,g3156d2b45e+6e55a43351,g347aa1857d+ae03bbfc84,g35bb328faa+a8ce1bb630,g3a166c0a6a+ae03bbfc84,g3e281a1b8c+c5dd892a6c,g414038480c+6b9177ef31,g41af890bb2+8f257c4c0b,g781aacb6e4+a8ce1bb630,g7af13505b9+7137b3b17d,g80478fca09+6df6903293,g82479be7b0+091ce1d07f,g858d7b2824+df8a265115,g89c8672015+f4add4ffd5,g9125e01d80+a8ce1bb630,g9726552aa6+414189b318,ga5288a1d22+4a2bca08d7,gacef1a1666+c9a8ff65f4,gb58c049af0+d64f4d3760,gbcfae0f0a0+de1d42d831,gc28159a63d+ae03bbfc84,gcf0d15dbbd+72117bf34e,gda6a2b7d83+72117bf34e,gdaeeff99f8+1711a396fd,ge500cccec5+c8c9c9af63,ge79ae78c31+ae03bbfc84,gf0baf85859+c1f95f4921,gfa517265be+df8a265115,gfa999e8aa5+17cd334064,gfb92a5be7c+df8a265115
LSST Data Management Base Package
Loading...
Searching...
No Matches
Classes | Public Types | Public Member Functions | Static Public Member Functions | Protected Types | Protected Member Functions | Static Protected Member Functions | Related Symbols | List of all members
lsst::meas::algorithms::CoaddPsf Class Reference

CoaddPsf is the Psf derived to be used for non-PSF-matched Coadd images. More...

#include <CoaddPsf.h>

Inheritance diagram for lsst::meas::algorithms::CoaddPsf:
lsst::afw::table::io::PersistableFacade< CoaddPsf > lsst::meas::algorithms::ImagePsf lsst::afw::table::io::PersistableFacade< ImagePsf > lsst::afw::detection::Psf lsst::afw::table::io::PersistableFacade< Psf > lsst::afw::typehandling::Storable lsst::afw::table::io::Persistable

Classes

class  Factory
 

Public Types

enum  ImageOwnerEnum { COPY = 0 , INTERNAL = 1 }
 Enum passed to computeImage and computeKernelImage to determine image ownership. More...
 
using Pixel = math::Kernel::Pixel
 Pixel type of Image returned by computeImage.
 
using Image = image::Image<Pixel>
 Image type returned by computeImage.
 

Public Member Functions

 CoaddPsf (afw::table::ExposureCatalog const &catalog, afw::geom::SkyWcs const &coaddWcs, std::string const &weightFieldName="weight", std::string const &warpingKernelName="lanczos3", int cacheSize=10000)
 Main constructors for CoaddPsf.
 
 CoaddPsf (afw::table::ExposureCatalog const &catalog, afw::geom::SkyWcs const &coaddWcs, CoaddPsfControl const &ctrl, std::string const &weightFieldName="weight")
 Constructor for CoaddPsf.
 
 CoaddPsf (afw::table::ExposureCatalog const &catalog, afw::geom::SkyWcs const &coaddWcs, geom::Point2D const &averagePosition, std::string const &warpingKernelName="lanczos3", int cacheSize=10000)
 
std::shared_ptr< afw::detection::Psfclone () const override
 Polymorphic deep copy. Usually unnecessary, as Psfs are immutable.
 
std::shared_ptr< afw::detection::Psfresized (int width, int height) const override
 Return a clone with specified kernel dimensions.
 
geom::Point2D getAveragePosition () const override
 Return the average of the positions of the stars that went into this Psf.
 
afw::geom::SkyWcs getCoaddWcs ()
 Return the Wcs of the coadd (defines the coordinate system of the Psf).
 
int getComponentCount () const
 Return the number of component Psfs in this CoaddPsf.
 
std::shared_ptr< afw::detection::Psf const > getPsf (int index)
 Get the Psf of the component image at index.
 
afw::geom::SkyWcs getWcs (int index)
 Get the Wcs of the component image at index.
 
double getWeight (int index)
 Get the weight of the component image at index.
 
afw::table::RecordId getId (int index)
 Get the exposure ID of the component image at index.
 
geom::Box2I getBBox (int index)
 Get the bounding box (in component image Pixel coordinates) of the component image at index.
 
std::shared_ptr< afw::geom::polygon::Polygon const > getValidPolygon (int index)
 Get the validPolygon (in component image Pixel coordinates) of the component image at index.
 
bool isPersistable () const noexcept override
 Return true if the CoaddPsf persistable (always true).
 
std::shared_ptr< typehandling::Storable > cloneStorable () const override
 Polymorphic deep-copy.
 
std::shared_ptr< ImagecomputeImage (lsst::geom::Point2D position, image::Color color=image::Color(), ImageOwnerEnum owner=COPY) const
 Return an Image of the PSF, in a form that can be compared directly with star images.
 
std::shared_ptr< ImagecomputeKernelImage (lsst::geom::Point2D position, image::Color color=image::Color(), ImageOwnerEnum owner=COPY) const
 Return an Image of the PSF, in a form suitable for convolution.
 
double computePeak (lsst::geom::Point2D position, image::Color color=image::Color()) const
 Return the peak value of the PSF image.
 
double computeApertureFlux (double radius, lsst::geom::Point2D position, image::Color color=image::Color()) const
 Compute the "flux" of the Psf model within a circular aperture of the given radius.
 
geom::ellipses::Quadrupole computeShape (lsst::geom::Point2D position, image::Color color=image::Color()) const
 Compute the ellipse corresponding to the second moments of the Psf.
 
std::shared_ptr< math::Kernel const > getLocalKernel (lsst::geom::Point2D position, image::Color color=image::Color()) const
 Return a FixedKernel corresponding to the Psf image at the given point.
 
image::Color getAverageColor () const
 Return the average Color of the stars used to construct the Psf.
 
lsst::geom::Box2I computeBBox (lsst::geom::Point2D position, image::Color color=image::Color()) const
 Return the bounding box of the image returned by computeKernelImage()
 
lsst::geom::Box2I computeImageBBox (lsst::geom::Point2D position, image::Color color=image::Color()) const
 Return the bounding box of the image returned by computeImage()
 
lsst::geom::Box2I computeKernelBBox (lsst::geom::Point2D position, image::Color color=image::Color()) const
 Return the bounding box of the image returned by computeImage()
 
std::size_t getCacheCapacity () const
 Return the capacity of the caches.
 
void setCacheCapacity (std::size_t capacity)
 Set the capacity of the caches.
 
virtual std::string toString () const
 Create a string representation of this object (optional operation).
 
virtual std::size_t hash_value () const
 Return a hash of this object (optional operation).
 
virtual bool equals (Storable const &other) const noexcept
 Compare this object to another Storable.
 
void writeFits (std::string const &fileName, std::string const &mode="w") const
 Write the object to a regular FITS file.
 
void writeFits (fits::MemFileManager &manager, std::string const &mode="w") const
 Write the object to a FITS image in memory.
 
void writeFits (fits::Fits &fitsfile) const
 Write the object to an already-open FITS object.
 

Static Public Member Functions

static std::shared_ptr< CoaddPsfreadFits (fits::Fits &fitsfile)
 Read an object from an already open FITS object.
 
static std::shared_ptr< CoaddPsfreadFits (std::string const &fileName, int hdu=fits::DEFAULT_HDU)
 Read an object from a regular FITS file.
 
static std::shared_ptr< CoaddPsfreadFits (fits::MemFileManager &manager, int hdu=fits::DEFAULT_HDU)
 Read an object from a FITS file in memory.
 
static std::shared_ptr< CoaddPsfdynamicCast (std::shared_ptr< Persistable > const &ptr)
 Dynamically cast a shared_ptr.
 
static std::shared_ptr< ImagePsfreadFits (fits::Fits &fitsfile)
 Read an object from an already open FITS object.
 
static std::shared_ptr< ImagePsfreadFits (std::string const &fileName, int hdu=fits::DEFAULT_HDU)
 Read an object from a regular FITS file.
 
static std::shared_ptr< ImagePsfreadFits (fits::MemFileManager &manager, int hdu=fits::DEFAULT_HDU)
 Read an object from a FITS file in memory.
 
static std::shared_ptr< PsfreadFits (fits::Fits &fitsfile)
 Read an object from an already open FITS object.
 
static std::shared_ptr< PsfreadFits (std::string const &fileName, int hdu=fits::DEFAULT_HDU)
 Read an object from a regular FITS file.
 
static std::shared_ptr< PsfreadFits (fits::MemFileManager &manager, int hdu=fits::DEFAULT_HDU)
 Read an object from a FITS file in memory.
 
static std::shared_ptr< ImagePsfdynamicCast (std::shared_ptr< Persistable > const &ptr)
 Dynamically cast a shared_ptr.
 
static std::shared_ptr< PsfdynamicCast (std::shared_ptr< Persistable > const &ptr)
 Dynamically cast a shared_ptr.
 
static std::shared_ptr< ImagerecenterKernelImage (std::shared_ptr< Image > im, lsst::geom::Point2D const &position, std::string const &warpAlgorithm="lanczos5", unsigned int warpBuffer=5)
 Helper function for Psf::doComputeImage(): converts a kernel image (centered at (0,0) when xy0 is taken into account) to an image centered at position when xy0 is taken into account.
 

Protected Types

using OutputArchiveHandle = io::OutputArchiveHandle
 

Protected Member Functions

std::shared_ptr< afw::detection::Psf::ImagedoComputeKernelImage (geom::Point2D const &ccdXY, afw::image::Color const &color) const override
 These virtual member functions are private, not protected, because we only want derived classes to implement them, not call them; they should call the corresponding compute*Image member functions instead so as to let the Psf base class handle caching properly.
 
geom::Box2I doComputeBBox (geom::Point2D const &position, afw::image::Color const &color) const override
 
std::string getPersistenceName () const override
 Return the unique name used to persist this object and look up its factory.
 
std::string getPythonModule () const override
 Return the fully-qualified Python module that should be imported to guarantee that its factory is registered.
 
void write (OutputArchiveHandle &handle) const override
 Write the object to one or more catalogs.
 
virtual double doComputeApertureFlux (double radius, geom::Point2D const &position, afw::image::Color const &color) const
 
virtual afw::geom::ellipses::Quadrupole doComputeShape (geom::Point2D const &position, afw::image::Color const &color) const
 
virtual std::shared_ptr< ImagedoComputeImage (lsst::geom::Point2D const &position, image::Color const &color) const
 These virtual members are protected (rather than private) so that python-implemented derived classes may opt to use the default implementations.
 
virtual lsst::geom::Box2I doComputeImageBBox (lsst::geom::Point2D const &position, image::Color const &color) const
 

Static Protected Member Functions

template<class T >
static bool singleClassEquals (T const &lhs, Storable const &rhs)
 Test if a Storable is of a particular class and equal to another object.
 

Related Symbols

(Note that these are not member symbols.)

std::ostreamoperator<< (std::ostream &os, Storable const &storable)
 Output operator for Storable.
 

Detailed Description

CoaddPsf is the Psf derived to be used for non-PSF-matched Coadd images.

It incorporates the logic of James Jee's Stackfit algorithm for estimating the Psf of coadd by coadding the images of the Psf models of each input exposure.

Definition at line 58 of file CoaddPsf.h.

Member Typedef Documentation

◆ Image

Image type returned by computeImage.

Definition at line 89 of file Psf.h.

◆ OutputArchiveHandle

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

Definition at line 108 of file Persistable.h.

◆ Pixel

using lsst::afw::detection::Psf::Pixel = math::Kernel::Pixel
inherited

Pixel type of Image returned by computeImage.

Definition at line 88 of file Psf.h.

Member Enumeration Documentation

◆ ImageOwnerEnum

Enum passed to computeImage and computeKernelImage to determine image ownership.

Enumerator
COPY 

The image will be copied before returning; caller will own it.

INTERNAL 

An internal image will be returned without copying.

The caller must not modify it, and it may be invalidated the next time a Psf member function is called with different color and/or position.

Definition at line 92 of file Psf.h.

92 {
93 COPY = 0,
94 INTERNAL = 1
98 };
@ COPY
The image will be copied before returning; caller will own it.
Definition Psf.h:93
@ INTERNAL
An internal image will be returned without copying.
Definition Psf.h:94

Constructor & Destructor Documentation

◆ CoaddPsf() [1/3]

lsst::meas::algorithms::CoaddPsf::CoaddPsf ( afw::table::ExposureCatalog const & catalog,
afw::geom::SkyWcs const & coaddWcs,
std::string const & weightFieldName = "weight",
std::string const & warpingKernelName = "lanczos3",
int cacheSize = 10000 )
explicit

Main constructors for CoaddPsf.

The ExposureCatalog contains info about each visit/ccd in Coadd; this must be provided to the constructor, and cannot be changed.

Parameters
[in]catalogExposureCatalog containing the id, bbox, wcs, psf and weight for each ccd/visit. This is usually the same catalog as the "ccds" catalog in the coadd Exposure's CoaddInputs.
[in]coaddWcsWcs for the coadd.
[in]weightFieldNameField name that contains the weight of the exposure in the coadd; defaults to "weight".
[in]warpingKernelNameName of warping kernel
[in]cacheSizeWarping kernel cache size

Definition at line 142 of file CoaddPsf.cc.

144 : _coaddWcs(coaddWcs),
145 _warpingKernelName(warpingKernelName),
146 _warpingControl(std::make_shared<afw::math::WarpingControl>(warpingKernelName, "", cacheSize)) {
147 afw::table::SchemaMapper mapper(catalog.getSchema());
149
150 // copy the field "goodpix", if available, for computeAveragePosition to use
151 try {
152 afw::table::Key<int> goodPixKey = catalog.getSchema()["goodpix"]; // auto does not work
153 mapper.addMapping(goodPixKey, true);
154 } catch (pex::exceptions::NotFoundError &) {
155 }
156
157 // copy the field specified by weightFieldName to field "weight"
158 afw::table::Field<double> weightField = afw::table::Field<double>("weight", "Coadd weight");
159 afw::table::Key<double> weightKey = catalog.getSchema()[weightFieldName];
160 _weightKey = mapper.addMapping(weightKey, weightField);
161
162 _catalog = afw::table::ExposureCatalog(mapper.getOutputSchema());
163 for (afw::table::ExposureCatalog::const_iterator i = catalog.begin(); i != catalog.end(); ++i) {
164 std::shared_ptr<afw::table::ExposureRecord> record = _catalog.getTable()->makeRecord();
165 record->assign(*i, mapper);
166 _catalog.push_back(record);
167 }
168 _averagePosition = computeAveragePosition(_catalog, _coaddWcs, _weightKey);
169}
SchemaMapper * mapper
std::shared_ptr< Table > getTable() const
Return the table associated with the catalog.
Definition Catalog.h:114
void push_back(Record const &r)
Add a copy of the given record to the end of the catalog.
Definition Catalog.h:480
typename Base::const_iterator const_iterator
Definition Exposure.h:318
static Schema makeMinimalSchema()
Return a minimal schema for Exposure tables and records.
Definition Exposure.h:216
ExposureCatalogT< ExposureRecord > ExposureCatalog
Definition Exposure.h:488
afw::table::Key< int > coaddWcs
afw::table::Key< std::string > warpingKernelName
Definition CoaddPsf.cc:355
table::Key< int > cacheSize

◆ CoaddPsf() [2/3]

lsst::meas::algorithms::CoaddPsf::CoaddPsf ( afw::table::ExposureCatalog const & catalog,
afw::geom::SkyWcs const & coaddWcs,
CoaddPsfControl const & ctrl,
std::string const & weightFieldName = "weight" )
inline

Constructor for CoaddPsf.

The ExposureCatalog contains info about each visit/ccd in Coadd; this must be provided to the constructor, and cannot be changed.

Parameters
[in]catalogExposureCatalog containing the id, bbox, wcs, psf and weight for each ccd/visit. This is usually the same catalog as the "ccds" catalog in the coadd Exposure's CoaddInputs.
[in]coaddWcsWcs for the coadd.
[in]ctrlConfiguration options.
[in]weightFieldNameField name that contains the weight of the exposure in the coadd; defaults to "weight".

Definition at line 93 of file CoaddPsf.h.

95 : CoaddPsf(catalog, coaddWcs, weightFieldName, ctrl.warpingKernelName, ctrl.cacheSize) {}
CoaddPsf(afw::table::ExposureCatalog const &catalog, afw::geom::SkyWcs const &coaddWcs, std::string const &weightFieldName="weight", std::string const &warpingKernelName="lanczos3", int cacheSize=10000)
Main constructors for CoaddPsf.
Definition CoaddPsf.cc:142

◆ CoaddPsf() [3/3]

lsst::meas::algorithms::CoaddPsf::CoaddPsf ( afw::table::ExposureCatalog const & catalog,
afw::geom::SkyWcs const & coaddWcs,
geom::Point2D const & averagePosition,
std::string const & warpingKernelName = "lanczos3",
int cacheSize = 10000 )
explicit
Parameters
catalogUnpersisted catalog
coaddWcsWCS for the coadd
averagePositionDefault position for accessors
warpingKernelNameWarping kernel name
cacheSizeKernel cache size

Definition at line 171 of file CoaddPsf.cc.

173 : _catalog(catalog),
174 _coaddWcs(coaddWcs),
175 _weightKey(_catalog.getSchema()["weight"]),
176 _averagePosition(averagePosition),
177 _warpingKernelName(warpingKernelName),
178 _warpingControl(new afw::math::WarpingControl(warpingKernelName, "", cacheSize)) {}
Schema getSchema() const
Return the schema associated with the catalog's table.
Definition Catalog.h:117
afw::table::PointKey< double > averagePosition
Definition CoaddPsf.cc:354

Member Function Documentation

◆ clone()

std::shared_ptr< afw::detection::Psf > lsst::meas::algorithms::CoaddPsf::clone ( ) const
overridevirtual

Polymorphic deep copy. Usually unnecessary, as Psfs are immutable.

Implements lsst::afw::detection::Psf.

Definition at line 180 of file CoaddPsf.cc.

180{ return std::make_shared<CoaddPsf>(*this); }

◆ cloneStorable()

std::shared_ptr< typehandling::Storable > lsst::afw::detection::Psf::cloneStorable ( ) const
inlineoverridevirtualinherited

Polymorphic deep copy. Usually unnecessary, as Psfs are immutable.

This method is an alias of clone that can be called from a reference to Storable.

Reimplemented from lsst::afw::typehandling::Storable.

Definition at line 122 of file Psf.h.

122{ return clone(); }
virtual std::shared_ptr< Psf > clone() const =0
Polymorphic deep-copy.

◆ computeApertureFlux()

double Psf::computeApertureFlux ( double radius,
lsst::geom::Point2D position,
image::Color color = image::Color() ) const
inherited

Compute the "flux" of the Psf model within a circular aperture of the given radius.

Parameters
[in]radiusRadius of the aperture to measure.
[in]positionPosition at which to evaluate the PSF.
[in]colorColor of the source for which to evaluate the PSF; defaults to getAverageColor().

The flux is relative to a Psf image that has been normalized to unit integral, and the radius is in pixels.

Definition at line 165 of file Psf.cc.

165 {
166 if (isPointNull(position)) position = getAveragePosition();
167 if (color.isIndeterminate()) color = getAverageColor();
168 return doComputeApertureFlux(radius, position, color);
169}
image::Color getAverageColor() const
Return the average Color of the stars used to construct the Psf.
Definition Psf.h:245
virtual double doComputeApertureFlux(double radius, lsst::geom::Point2D const &position, image::Color const &color) const =0
virtual lsst::geom::Point2D getAveragePosition() const
Return the average position of the stars used to construct the Psf.
Definition Psf.cc:189
bool isIndeterminate() const noexcept
Whether the color is the special value that indicates that it is unspecified.
Definition Color.h:36

◆ computeBBox()

lsst::geom::Box2I Psf::computeBBox ( lsst::geom::Point2D position,
image::Color color = image::Color() ) const
inherited

Return the bounding box of the image returned by computeKernelImage()

Parameters
[in]positionPosition at which to evaluate the PSF.
[in]colorColor of the source for which to evaluate the PSF; defaults to getAverageColor().

Definition at line 127 of file Psf.cc.

127 {
128 if (isPointNull(position)) position = getAveragePosition();
129 if (color.isIndeterminate()) color = getAverageColor();
130 auto cached_image = _kernelImageCache->get(detail::PsfCacheKey(position, color));
131 if (cached_image.has_value()) {
132 return cached_image.value()->getBBox();
133 } else {
134 return doComputeBBox(position, color);
135 }
136}
virtual lsst::geom::Box2I doComputeBBox(lsst::geom::Point2D const &position, image::Color const &color) const =0

◆ computeImage()

std::shared_ptr< Psf::Image > Psf::computeImage ( lsst::geom::Point2D position,
image::Color color = image::Color(),
ImageOwnerEnum owner = COPY ) const
inherited

Return an Image of the PSF, in a form that can be compared directly with star images.

The specified position is a floating point number, and the resulting image will have a Psf centered on that point when the returned image's xy0 is taken into account.

The returned image is normalized to sum to unity.

Parameters
[in]positionPosition at which to evaluate the PSF.
[in]colorColor of the source for which to evaluate the PSF; defaults to getAverageColor().
[in]ownerWhether to copy the return value or return an internal image that must be handled with care (see ImageOwnerEnum).

The Psf class caches the most recent return value of computeImage, so repeated calls with the same arguments will be highly optimized.

Note
The real work is done in the virtual private member function Psf::doComputeImage; computeImage only handles caching and default arguments.

Definition at line 101 of file Psf.cc.

102 {
103 if (isPointNull(position)) position = getAveragePosition();
104 if (color.isIndeterminate()) color = getAverageColor();
105 std::shared_ptr<Psf::Image> result = (*_imageCache)(
106 detail::PsfCacheKey(position, color),
107 [this](detail::PsfCacheKey const &key) { return doComputeImage(key.position, key.color); });
108 if (owner == COPY) {
109 result = std::make_shared<Image>(*result, true);
110 }
111 return result;
112}
py::object result
Definition _schema.cc:429
virtual std::shared_ptr< Image > doComputeImage(lsst::geom::Point2D const &position, image::Color const &color) const
These virtual members are protected (rather than private) so that python-implemented derived classes ...
Definition Psf.cc:177

◆ computeImageBBox()

lsst::geom::Box2I Psf::computeImageBBox ( lsst::geom::Point2D position,
image::Color color = image::Color() ) const
inherited

Return the bounding box of the image returned by computeImage()

Parameters
[in]positionPosition at which to evaluate the PSF.
[in]colorColor of the source for which to evaluate the PSF; defaults to getAverageColor().

Definition at line 138 of file Psf.cc.

138 {
139 if (isPointNull(position)) position = getAveragePosition();
140 if (color.isIndeterminate()) color = getAverageColor();
141 auto cached_image = _imageCache->get(detail::PsfCacheKey(position, color));
142 if (cached_image.has_value()) {
143 return cached_image.value()->getBBox();
144 } else {
145 return doComputeImageBBox(position, color);
146 }
147}
virtual lsst::geom::Box2I doComputeImageBBox(lsst::geom::Point2D const &position, image::Color const &color) const
Definition Psf.cc:183

◆ computeKernelBBox()

lsst::geom::Box2I lsst::afw::detection::Psf::computeKernelBBox ( lsst::geom::Point2D position,
image::Color color = image::Color() ) const
inlineinherited

Return the bounding box of the image returned by computeImage()

Parameters
[in]positionPosition at which to evaluate the PSF.
[in]colorColor of the source for which to evaluate the PSF; defaults to getAverageColor(). Alias for computeBBox

Definition at line 282 of file Psf.h.

283 {
284 return computeBBox(position, color);
285 }
lsst::geom::Box2I computeBBox(lsst::geom::Point2D position, image::Color color=image::Color()) const
Return the bounding box of the image returned by computeKernelImage()
Definition Psf.cc:127

◆ computeKernelImage()

std::shared_ptr< Psf::Image > Psf::computeKernelImage ( lsst::geom::Point2D position,
image::Color color = image::Color(),
ImageOwnerEnum owner = COPY ) const
inherited

Return an Image of the PSF, in a form suitable for convolution.

While the position need not be an integer, the center of the PSF image returned by computeKernelImage will in the center of the center pixel of the image, which will be (0,0) when the Image's xy0 is taken into account; this is the same behavior as Kernel::computeImage().

The returned image is normalized to sum to unity.

Parameters
[in]positionPosition at which to evaluate the PSF.
[in]colorColor of the source for which to evaluate the PSF; defaults to getAverageColor().
[in]ownerWhether to copy the return value or return an internal image that must be handled with care (see ImageOwnerEnum).

The Psf class caches the most recent return value of computeKernelImage, so repeated calls with the same arguments will be highly optimized.

Note
The real work is done in the virtual private member function Psf::doComputeKernelImage; computeKernelImage only handles caching and default arguments.

Definition at line 114 of file Psf.cc.

115 {
116 if (_isFixed || isPointNull(position)) position = getAveragePosition();
117 if (_isFixed || color.isIndeterminate()) color = getAverageColor();
118 std::shared_ptr<Psf::Image> result = (*_kernelImageCache)(
119 detail::PsfCacheKey(position, color),
120 [this](detail::PsfCacheKey const &key) { return doComputeKernelImage(key.position, key.color); });
121 if (owner == COPY) {
122 result = std::make_shared<Image>(*result, true);
123 }
124 return result;
125}
virtual std::shared_ptr< Image > doComputeKernelImage(lsst::geom::Point2D const &position, image::Color const &color) const =0
These virtual member functions are private, not protected, because we only want derived classes to im...

◆ computePeak()

double Psf::computePeak ( lsst::geom::Point2D position,
image::Color color = image::Color() ) const
inherited

Return the peak value of the PSF image.

Parameters
[in]positionPosition at which to evaluate the PSF.
[in]colorColor of the source for which to evaluate the PSF; defaults to getAverageColor().

This calls computeKernelImage internally, but because this will usually be cached, it shouldn't be expensive (but be careful not to accidentally call it with no arguments when you actually want to call it with the same arguments just used to call computeImage or computeKernelImage).

Definition at line 158 of file Psf.cc.

158 {
159 if (isPointNull(position)) position = getAveragePosition();
160 if (color.isIndeterminate()) color = getAverageColor();
162 return (*image)(-image->getX0(), -image->getY0());
163}
std::shared_ptr< Image > computeKernelImage(lsst::geom::Point2D position, image::Color color=image::Color(), ImageOwnerEnum owner=COPY) const
Return an Image of the PSF, in a form suitable for convolution.
Definition Psf.cc:114

◆ computeShape()

geom::ellipses::Quadrupole Psf::computeShape ( lsst::geom::Point2D position,
image::Color color = image::Color() ) const
inherited

Compute the ellipse corresponding to the second moments of the Psf.

Parameters
[in]positionPosition at which to evaluate the PSF.
[in]colorColor of the source for which to evaluate the PSF; defaults to getAverageColor().

The algorithm used to compute the moments is up to the derived class, and hence this method should not be used when a particular algorithm or weight function is required.

Definition at line 171 of file Psf.cc.

171 {
172 if (isPointNull(position)) position = getAveragePosition();
173 if (color.isIndeterminate()) color = getAverageColor();
174 return doComputeShape(position, color);
175}
virtual geom::ellipses::Quadrupole doComputeShape(lsst::geom::Point2D const &position, image::Color const &color) const =0

◆ doComputeApertureFlux()

double lsst::meas::algorithms::ImagePsf::doComputeApertureFlux ( double radius,
geom::Point2D const & position,
afw::image::Color const & color ) const
protectedvirtualinherited

Implements lsst::afw::detection::Psf.

Definition at line 45 of file ImagePsf.cc.

46 {
47 geom::Point2D const center(0.0, 0.0);
48 afw::geom::ellipses::Axes const axes(radius, radius);
49 afw::geom::ellipses::Ellipse ellipse(axes, center);
50 geom::Box2I box(ellipse.computeBBox());
51 // perhaps needs a grow here
52 afw::image::Image<double> image(*computeKernelImage(position, color, INTERNAL));
53 if (!image.getBBox().contains(box)) {
54 box.include(image.getBBox());
55 afw::image::Image<double> paddedImage(box);
56 paddedImage = 0.0;
57 paddedImage.assign(image, image.getBBox());
58 image = paddedImage;
59 }
61 image, ellipse, base::ApertureFluxControl());
62 return result.instFlux;
63}
afw::table::Key< afw::table::Array< ImagePixelT > > image
An integer coordinate rectangle.
Definition Box.h:55
static Result computeSincFlux(afw::image::Image< T > const &image, afw::geom::ellipses::Ellipse const &ellipse, Control const &ctrl=Control())
Compute the instFlux (and optionally, uncertanties) within an aperture using Sinc photometry.

◆ doComputeBBox()

geom::Box2I lsst::meas::algorithms::CoaddPsf::doComputeBBox ( geom::Point2D const & position,
afw::image::Color const & color ) const
overrideprotectedvirtual

Implements lsst::afw::detection::Psf.

Definition at line 229 of file CoaddPsf.cc.

229 {
230 afw::table::ExposureCatalog subcat = _catalog.subsetContaining(ccdXY, _coaddWcs, true);
231 if (subcat.empty()) {
232 throw LSST_EXCEPT(
234 (boost::format("Cannot compute BBox at point %s; no input images at that point.") % ccdXY)
235 .str());
236 }
237
238 geom::Box2I ret;
239 for (auto const &exposureRecord : subcat) {
240 // compute transform from exposure pixels to coadd pixels
241 auto exposureToCoadd = afw::geom::makeWcsPairTransform(*exposureRecord.getWcs(), _coaddWcs);
242 WarpedPsf warpedPsf = WarpedPsf(exposureRecord.getPsf(), exposureToCoadd, _warpingControl);
243 geom::Box2I componentBBox = warpedPsf.computeBBox(ccdXY, color);
244 ret.include(componentBBox);
245 }
246
247 return ret;
248}
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition Exception.h:48
An exception thrown when we have an invalid PSF.
Definition Psf.h:52
ExposureCatalogT subsetContaining(lsst::geom::SpherePoint const &coord, bool includeValidPolygon=false) const
Return a shallow subset of the catalog with only those records that contain the given point.
Definition Exposure.cc:473
void include(Point2I const &point)
Expand this to ensure that this->contains(point).
Definition Box.cc:152
std::shared_ptr< TransformPoint2ToPoint2 > makeWcsPairTransform(SkyWcs const &src, SkyWcs const &dst)
A Transform obtained by putting two SkyWcs objects "back to back".
Definition SkyWcs.cc:146

◆ doComputeImage()

std::shared_ptr< Psf::Image > Psf::doComputeImage ( lsst::geom::Point2D const & position,
image::Color const & color ) const
protectedvirtualinherited

These virtual members are protected (rather than private) so that python-implemented derived classes may opt to use the default implementations.

C++ derived classes may override these methods, but should not call them. Derived classes should call the corresponding compute* member functions instead so as to let the Psf base class handle caching properly.

Derived classes are responsible for ensuring that returned images sum to one.

Reimplemented in lsst::meas::extensions::psfex::PsfexPsf.

Definition at line 177 of file Psf.cc.

178 {
180 return recenterKernelImage(im, position);
181}
static std::shared_ptr< Image > recenterKernelImage(std::shared_ptr< Image > im, lsst::geom::Point2D const &position, std::string const &warpAlgorithm="lanczos5", unsigned int warpBuffer=5)
Helper function for Psf::doComputeImage(): converts a kernel image (centered at (0,...
Definition Psf.cc:85

◆ doComputeImageBBox()

lsst::geom::Box2I Psf::doComputeImageBBox ( lsst::geom::Point2D const & position,
image::Color const & color ) const
protectedvirtualinherited

Definition at line 183 of file Psf.cc.

184 {
186 return im->getBBox();
187}
std::shared_ptr< Image > computeImage(lsst::geom::Point2D position, image::Color color=image::Color(), ImageOwnerEnum owner=COPY) const
Return an Image of the PSF, in a form that can be compared directly with star images.
Definition Psf.cc:101

◆ doComputeKernelImage()

std::shared_ptr< afw::detection::Psf::Image > lsst::meas::algorithms::CoaddPsf::doComputeKernelImage ( geom::Point2D const & position,
afw::image::Color const & color ) const
overrideprotectedvirtual

These virtual member functions are private, not protected, because we only want derived classes to implement them, not call them; they should call the corresponding compute*Image member functions instead so as to let the Psf base class handle caching properly.

Derived classes are responsible for ensuring that returned images sum to one.

Implements lsst::afw::detection::Psf.

Definition at line 251 of file CoaddPsf.cc.

251 {
252 // Get the subset of exposures which contain our coordinate within their validPolygons.
253 afw::table::ExposureCatalog subcat = _catalog.subsetContaining(ccdXY, _coaddWcs, true);
254 if (subcat.empty()) {
255 throw LSST_EXCEPT(
257 (boost::format("Cannot compute CoaddPsf at point %s; no input images at that point.") % ccdXY)
258 .str());
259 }
260 double weightSum = 0.0;
261
262 // Read all the Psf images into a vector. The code is set up so that this can be done in chunks,
263 // with the image modified to accomodate
264 // However, we currently read all of the images.
266 std::vector<double> weightVector;
267
268 for (auto const &exposureRecord : subcat) {
269 // compute transform from exposure pixels to coadd pixels
270 auto exposureToCoadd = afw::geom::makeWcsPairTransform(*exposureRecord.getWcs(), _coaddWcs);
272 try {
273 WarpedPsf warpedPsf = WarpedPsf(exposureRecord.getPsf(), exposureToCoadd, _warpingControl);
274 componentImg = warpedPsf.computeKernelImage(ccdXY, color);
275 } catch (pex::exceptions::RangeError &exc) {
276 LSST_EXCEPT_ADD(exc, (boost::format("Computing WarpedPsf kernel image for id=%d") %
277 exposureRecord.getId())
278 .str());
279 throw exc;
280 }
281 imgVector.push_back(componentImg);
282 weightSum += exposureRecord.get(_weightKey);
283 weightVector.push_back(exposureRecord.get(_weightKey));
284 }
285
286 geom::Box2I bbox = getOverallBBox(imgVector);
287
288 // create a zero image of the right size to sum into
289 std::shared_ptr<afw::detection::Psf::Image> image = std::make_shared<afw::detection::Psf::Image>(bbox);
290 *image = 0.0;
291 addToImage(image, imgVector, weightVector);
292 *image /= weightSum;
293 return image;
294}
AmpInfoBoxKey bbox
Definition Amplifier.cc:117
#define LSST_EXCEPT_ADD(e, m)
Add the current location and a message to an existing exception before rethrowing it.
Definition Exception.h:54
geom::Box2I getOverallBBox(std::vector< std::shared_ptr< afw::image::Image< double > > > const &imgVector)
Definition CoaddPsf.cc:189
T push_back(T... args)

◆ doComputeShape()

afw::geom::ellipses::Quadrupole lsst::meas::algorithms::ImagePsf::doComputeShape ( geom::Point2D const & position,
afw::image::Color const & color ) const
protectedvirtualinherited

Implements lsst::afw::detection::Psf.

Definition at line 65 of file ImagePsf.cc.

66 {
69 *image, geom::Point2D(0.0, 0.0) // image has origin at the center
70 )
71 .getShape();
72}
static Result computeAdaptiveMoments(ImageT const &image, geom::Point2D const &position, bool negative=false, Control const &ctrl=Control())
Compute the adaptive Gaussian-weighted moments of an image.
Shape const getShape() const
Return an afw::geom::ellipses object corresponding to xx, yy, xy.

◆ dynamicCast() [1/3]

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 218 of file Persistable.cc.

◆ dynamicCast() [2/3]

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 218 of file Persistable.cc.

◆ dynamicCast() [3/3]

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 218 of file Persistable.cc.

◆ equals()

bool lsst::afw::typehandling::Storable::equals ( Storable const & other) const
virtualnoexceptinherited

Compare this object to another Storable.

Subclasses that implement equality comparison must override this method to give results consistent with operator== for all inputs that are accepted by both.

Returns
This implementation returns whether the two objects are the same.
Warning
This method compares an object to any type of Storable, although cross-class comparisons should usually return false. If cross-class comparisons are valid, implementers should take care that they are symmetric and will give the same result no matter what the compile-time types of the left- and right-hand sides are.
See also
singleClassEquals
Note
When called on Python classes, this method delegates to __eq__ if it exists.

Reimplemented in lsst::afw::typehandling::StorableHelper< Base >, lsst::afw::geom::polygon::Polygon, lsst::afw::geom::SkyWcs, lsst::afw::image::PhotoCalib, and lsst::afw::image::VisitInfo.

Definition at line 48 of file Storable.cc.

48{ return this == &other; }

◆ getAverageColor()

image::Color lsst::afw::detection::Psf::getAverageColor ( ) const
inlineinherited

Return the average Color of the stars used to construct the Psf.

This is also the Color used to return an image if you don't specify a Color.

Definition at line 245 of file Psf.h.

245{ return image::Color(); }
Describe the colour of a source.
Definition Color.h:25

◆ getAveragePosition()

geom::Point2D lsst::meas::algorithms::CoaddPsf::getAveragePosition ( ) const
inlineoverridevirtual

Return the average of the positions of the stars that went into this Psf.

For CoaddPsf, this is calculated as the weighted average of the average positions of all the component Psfs.

Reimplemented from lsst::afw::detection::Psf.

Definition at line 115 of file CoaddPsf.h.

115{ return _averagePosition; }

◆ getBBox()

geom::Box2I lsst::meas::algorithms::CoaddPsf::getBBox ( int index)

Get the bounding box (in component image Pixel coordinates) of the component image at index.

Parameters
[in]indexIndex of component for which the bounding box will be fetched.
Returns
Corresponding bounding box.
Exceptions
RangeErrorIndex of component is out of range.

Definition at line 333 of file CoaddPsf.cc.

333 {
334 if (index < 0 || index >= getComponentCount()) {
335 throw LSST_EXCEPT(pex::exceptions::RangeError, "index of CoaddPsf component out of range");
336 }
337 return _catalog[index].getBBox();
338}
int getComponentCount() const
Return the number of component Psfs in this CoaddPsf.
Definition CoaddPsf.cc:296

◆ getCacheCapacity()

std::size_t Psf::getCacheCapacity ( ) const
inherited

Return the capacity of the caches.

Both the image and kernel image caches have the same capacity.

Definition at line 191 of file Psf.cc.

191{ return _kernelImageCache->capacity(); }

◆ getCoaddWcs()

afw::geom::SkyWcs lsst::meas::algorithms::CoaddPsf::getCoaddWcs ( )
inline

Return the Wcs of the coadd (defines the coordinate system of the Psf).

Definition at line 118 of file CoaddPsf.h.

118{ return _coaddWcs; }

◆ getComponentCount()

int lsst::meas::algorithms::CoaddPsf::getComponentCount ( ) const

Return the number of component Psfs in this CoaddPsf.

Definition at line 296 of file CoaddPsf.cc.

296{ return _catalog.size(); }
size_type size() const
Return the number of elements in the catalog.
Definition Catalog.h:412

◆ getId()

afw::table::RecordId lsst::meas::algorithms::CoaddPsf::getId ( int index)

Get the exposure ID of the component image at index.

Parameters
[in]indexIndex of component for which the exposure ID will be fetched.
Returns
Corresponding exposure ID.
Exceptions
RangeErrorIndex of component is out of range.

Definition at line 326 of file CoaddPsf.cc.

326 {
327 if (index < 0 || index >= getComponentCount()) {
328 throw LSST_EXCEPT(pex::exceptions::RangeError, "index of CoaddPsf component out of range");
329 }
330 return _catalog[index].getId();
331}

◆ getLocalKernel()

std::shared_ptr< math::Kernel const > Psf::getLocalKernel ( lsst::geom::Point2D position,
image::Color color = image::Color() ) const
inherited

Return a FixedKernel corresponding to the Psf image at the given point.

Parameters
[in]positionPosition at which to evaluate the PSF.
[in]colorColor of the source for which to evaluate the PSF; defaults to getAverageColor().

This is implemented by calling computeKernelImage, and is simply provided for convenience.

Definition at line 149 of file Psf.cc.

150 {
151 if (isPointNull(position)) position = getAveragePosition();
152 if (color.isIndeterminate()) color = getAverageColor();
153 // FixedKernel ctor will deep copy image, so we can use INTERNAL.
155 return std::make_shared<math::FixedKernel>(*image);
156}

◆ getPersistenceName()

std::string lsst::meas::algorithms::CoaddPsf::getPersistenceName ( ) const
overrideprotectedvirtual

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

Must be less than ArchiveIndexSchema::MAX_NAME_LENGTH characters.

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

Definition at line 430 of file CoaddPsf.cc.

430{ return getCoaddPsfPersistenceName(); }

◆ getPsf()

std::shared_ptr< afw::detection::Psf const > lsst::meas::algorithms::CoaddPsf::getPsf ( int index)

Get the Psf of the component image at index.

Parameters
[in]indexIndex of component for which the Psf will be fetched.
Returns
Corresponding Psf.
Exceptions
RangeErrorIndex of component is out of range.

Definition at line 298 of file CoaddPsf.cc.

298 {
299 if (index < 0 || index >= getComponentCount()) {
300 throw LSST_EXCEPT(pex::exceptions::RangeError, "index of CoaddPsf component out of range");
301 }
302 return _catalog[index].getPsf();
303}

◆ getPythonModule()

std::string lsst::meas::algorithms::CoaddPsf::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 432 of file CoaddPsf.cc.

432{ return "lsst.meas.algorithms"; }

◆ getValidPolygon()

std::shared_ptr< afw::geom::polygon::Polygon const > lsst::meas::algorithms::CoaddPsf::getValidPolygon ( int index)

Get the validPolygon (in component image Pixel coordinates) of the component image at index.

Parameters
[in]indexIndex of component for which the bounding box will be fetched.
Returns
Corresponding validPolygon.
Exceptions
RangeErrorIndex of component is out of range.

Definition at line 312 of file CoaddPsf.cc.

312 {
313 if (index < 0 || index >= getComponentCount()) {
314 throw LSST_EXCEPT(pex::exceptions::RangeError, "index of CoaddPsf component out of range");
315 }
316 return _catalog[index].getValidPolygon();
317}

◆ getWcs()

afw::geom::SkyWcs lsst::meas::algorithms::CoaddPsf::getWcs ( int index)

Get the Wcs of the component image at index.

Parameters
[in]indexIndex of component for which the Wcs will be fetched.
Returns
Corresponding Wcs.
Exceptions
RangeErrorIndex of component is out of range.

Definition at line 305 of file CoaddPsf.cc.

305 {
306 if (index < 0 || index >= getComponentCount()) {
307 throw LSST_EXCEPT(pex::exceptions::RangeError, "index of CoaddPsf component out of range");
308 }
309 return *_catalog[index].getWcs();
310}

◆ getWeight()

double lsst::meas::algorithms::CoaddPsf::getWeight ( int index)

Get the weight of the component image at index.

Parameters
[in]indexIndex of component for which the weight will be fetched.
Returns
Corresponding weight.
Exceptions
RangeErrorIndex of component is out of range.

Definition at line 319 of file CoaddPsf.cc.

319 {
320 if (index < 0 || index >= getComponentCount()) {
321 throw LSST_EXCEPT(pex::exceptions::RangeError, "index of CoaddPsf component out of range");
322 }
323 return _catalog[index].get(_weightKey);
324}
std::shared_ptr< RecordT > const get(size_type i) const
Return a pointer to the record at index i.
Definition Catalog.h:463

◆ hash_value()

std::size_t lsst::afw::typehandling::Storable::hash_value ( ) const
virtualinherited

Return a hash of this object (optional operation).

Exceptions
UnsupportedOperationExceptionThrown if this object is not hashable.
Note
C++ subclass authors are responsible for any associated specializations of std::hash.
When called on Python classes, this method delegates to __hash__ if it exists.

Reimplemented in lsst::afw::geom::polygon::Polygon, lsst::afw::image::FilterLabel, lsst::afw::image::VisitInfo, and lsst::afw::typehandling::StorableHelper< Base >.

Definition at line 44 of file Storable.cc.

44 {
45 throw LSST_EXCEPT(UnsupportedOperationException, "Hashes are not supported.");
46}

◆ isPersistable()

bool lsst::meas::algorithms::CoaddPsf::isPersistable ( ) const
inlineoverridevirtualnoexcept

Return true if the CoaddPsf persistable (always true).

While it's actually possible to construct a CoaddPsf that isn't persistable (because its nested Psfs and Wcss are not persistable) in artificial situations, in realistic situations it's pretty much impossible, because persistence is a necessary part of how CoaddPsfs are built. And it's simpler and much faster if we just always return true, rather than loop over the elements and check each one.

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

Definition at line 186 of file CoaddPsf.h.

186{ return true; }

◆ readFits() [1/9]

static std::shared_ptr< Psf > lsst::afw::table::io::PersistableFacade< Psf >::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.

◆ readFits() [2/9]

static std::shared_ptr< ImagePsf > lsst::afw::table::io::PersistableFacade< ImagePsf >::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.

◆ readFits() [3/9]

static std::shared_ptr< CoaddPsf > lsst::afw::table::io::PersistableFacade< CoaddPsf >::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.

◆ readFits() [4/9]

static std::shared_ptr< Psf > lsst::afw::table::io::PersistableFacade< Psf >::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.

◆ readFits() [5/9]

static std::shared_ptr< ImagePsf > lsst::afw::table::io::PersistableFacade< ImagePsf >::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.

◆ readFits() [6/9]

static std::shared_ptr< CoaddPsf > lsst::afw::table::io::PersistableFacade< CoaddPsf >::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.

◆ readFits() [7/9]

static std::shared_ptr< Psf > lsst::afw::table::io::PersistableFacade< Psf >::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.

◆ readFits() [8/9]

static std::shared_ptr< ImagePsf > lsst::afw::table::io::PersistableFacade< ImagePsf >::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.

◆ readFits() [9/9]

static std::shared_ptr< CoaddPsf > lsst::afw::table::io::PersistableFacade< CoaddPsf >::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.

◆ recenterKernelImage()

std::shared_ptr< image::Image< double > > Psf::recenterKernelImage ( std::shared_ptr< Image > im,
lsst::geom::Point2D const & position,
std::string const & warpAlgorithm = "lanczos5",
unsigned int warpBuffer = 5 )
staticinherited

Helper function for Psf::doComputeImage(): converts a kernel image (centered at (0,0) when xy0 is taken into account) to an image centered at position when xy0 is taken into account.

warpAlgorithm is passed to afw::math::makeWarpingKernel() and can be "nearest", "bilinear", or "lanczosN"

warpBuffer zero-pads the image before recentering. Recommended value is 1 for bilinear, N for lanczosN (note that it would be cleaner to infer this value from the warping algorithm but this would require mild API changes; same issue occurs in e.g. afw::math::offsetImage()).

The point with integer coordinates (0,0) in the source image (with xy0 taken into account) corresponds to the point position in the destination image. If position is not integer-valued then we will need to fractionally shift the image using interpolation.

Note: if fractional recentering is performed, then a new image will be allocated and returned. If not, then the original image will be returned (after setting XY0).

Definition at line 85 of file Psf.cc.

88 {
89 // "ir" : (integer, residual)
90 std::pair<int, double> const irX = image::positionToIndex(position.getX(), true);
91 std::pair<int, double> const irY = image::positionToIndex(position.getY(), true);
92
93 if (irX.second != 0.0 || irY.second != 0.0) {
94 im = math::offsetImage(*im, irX.second, irY.second, warpAlgorithm, warpBuffer);
95 }
96
97 im->setXY0(irX.first + im->getX0(), irY.first + im->getY0());
98 return im;
99}
int positionToIndex(double pos)
Convert image position to nearest integer index.
Definition ImageUtils.h:69
std::shared_ptr< ImageT > offsetImage(ImageT const &image, float dx, float dy, std::string const &algorithmName="lanczos5", unsigned int buffer=0)
Return an image offset by (dx, dy) using the specified algorithm.

◆ resized()

std::shared_ptr< afw::detection::Psf > lsst::meas::algorithms::CoaddPsf::resized ( int width,
int height ) const
overridevirtual

Return a clone with specified kernel dimensions.

Implements lsst::afw::detection::Psf.

Definition at line 182 of file CoaddPsf.cc.

182 {
183 // Not implemented for WarpedPsf
184 throw LSST_EXCEPT(pex::exceptions::LogicError, "Not Implemented");
185}

◆ setCacheCapacity()

void Psf::setCacheCapacity ( std::size_t capacity)
inherited

Set the capacity of the caches.

Both the image and kernel image caches will be set to this capacity.

Definition at line 193 of file Psf.cc.

193 {
194 _imageCache->reserve(capacity);
195 _kernelImageCache->reserve(capacity);
196}

◆ singleClassEquals()

template<class T >
static bool lsst::afw::typehandling::Storable::singleClassEquals ( T const & lhs,
Storable const & rhs )
inlinestaticprotectedinherited

Test if a Storable is of a particular class and equal to another object.

This method template simplifies implementations of equals that delegate to operator== without supporting cross-class comparisons.

Template Parameters
TThe class expected of the two objects to be compared.
Parameters
lhs,rhsThe objects to compare. Note that rhs need not be a T, while lhs must be.
Returns
true if rhs is a T and lhs == rhs; false otherwise.
Exception Safety
Provides the same level of exception safety as operator==. Most implementations of operator== do not throw.
Note
This method template calls operator== with both arguments of compile-time type T const&. Its use is not recommended if there would be any ambiguity as to which operator== gets picked by overload resolution.

This method template is typically called from equals as:

bool MyType::equals(Storable const& other) const noexcept {
    return singleClassEquals(*this, other);
}

Definition at line 151 of file Storable.h.

151 {
152 auto typedRhs = dynamic_cast<T const*>(&rhs);
153 if (typedRhs != nullptr) {
154 return lhs == *typedRhs;
155 } else {
156 return false;
157 }
158 }

◆ toString()

std::string lsst::afw::typehandling::Storable::toString ( ) const
virtualinherited

Create a string representation of this object (optional operation).

Exceptions
UnsupportedOperationExceptionThrown if this object does not have a string representation.
Note
When called on Python classes, this method delegates to __repr__.

Reimplemented in lsst::afw::geom::polygon::Polygon, lsst::afw::geom::SkyWcs, lsst::afw::image::FilterLabel, lsst::afw::image::PhotoCalib, lsst::afw::image::VisitInfo, and lsst::afw::typehandling::StorableHelper< Base >.

Definition at line 40 of file Storable.cc.

40 {
41 throw LSST_EXCEPT(UnsupportedOperationException, "No string representation available.");
42}

◆ write()

void lsst::meas::algorithms::CoaddPsf::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 434 of file CoaddPsf.cc.

434 {
435 CoaddPsfPersistenceHelper const &keys1 = CoaddPsfPersistenceHelper::get();
436 afw::table::BaseCatalog cat1 = handle.makeCatalog(keys1.schema);
437 std::shared_ptr<afw::table::BaseRecord> record1 = cat1.addNew();
438 auto coaddWcsPtr = std::make_shared<afw::geom::SkyWcs>(_coaddWcs);
439 record1->set(keys1.coaddWcs, handle.put(coaddWcsPtr));
440 record1->set(keys1.cacheSize, _warpingControl->getCacheSize());
441 record1->set(keys1.averagePosition, _averagePosition);
442 record1->set(keys1.warpingKernelName, _warpingKernelName);
443 handle.saveCatalog(cat1);
444 _catalog.writeToArchive(handle, false);
445}
void writeToArchive(io::OutputArchiveHandle &handle, bool ignoreUnpersistable=true) const
Convenience output function for Persistables that contain an ExposureCatalog.
Definition Exposure.cc:446
CatalogT< BaseRecord > BaseCatalog
Definition fwd.h:72

◆ 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]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}

◆ 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.

◆ 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]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}

Friends And Related Symbol Documentation

◆ operator<<()

std::ostream & operator<< ( std::ostream & os,
Storable const & storable )
related

Output operator for Storable.

Parameters
osthe desired output stream
storablethe object to print
Returns
a reference to os
Exceptions
UnsupportedOperationExceptionThrown if storable does not have an implementation of Storable::toString.

Definition at line 174 of file Storable.h.

174 {
175 return os << storable.toString();
176}

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