LSSTApplications  19.0.0-14-gb0260a2+492360823e,20.0.0+1c463ce1c5,20.0.0+332b39b81b,20.0.0+34a42eae2c,20.0.0+41b643067b,20.0.0+6e3b0e79d6,20.0.0+8558dd3f48,20.0.0+90990fb0e1,20.0.0+b0de6306d8,20.0.0+b290a576ab,20.0.0+b2ea66fa67,20.0.0+c2e041a55f,20.0.0+dcf29472a8,20.0.0+f45b7d88f4,20.0.0+fd13e6d5b6,20.0.0-1-g10df615+0d65e0fd03,20.0.0-1-g253301a+dcf29472a8,20.0.0-1-g498fb60+ff88705a28,20.0.0-1-g4d801e7+20a86399e6,20.0.0-1-g8a53f90+2817c06967,20.0.0-1-gc96f8cb+41b643067b,20.0.0-1-gd1c87d7+2817c06967,20.0.0-1-gdb27ee5+0eaf4f30a3,20.0.0-11-gda4966f+0cec7fe7b8,20.0.0-17-g98646ba4+f780282b58,20.0.0-2-g25c5319+ec537bc764,20.0.0-2-gec03fae+3bc057fb2a,20.0.0-27-g22bab1b+96b59cb0f5,20.0.0-3-gd2e950e+f45b7d88f4,20.0.0-4-g4a2362f+f45b7d88f4,20.0.0-5-gac0d578b1+6c871ee35c,20.0.0-5-gfcebe35+d7a723b4eb,20.0.0-6-g01203fff+106a5fd355,20.0.0-7-g3c4151b+2efbc3ed6b,20.0.0-8-g9573cbd+0eaf4f30a3,20.0.0-8-gc2abeef+41b643067b,w.2020.32
LSSTDataManagementBasePackage
Public Types | Public Member Functions | Static Public Member Functions | Protected Types | Protected Member Functions | Static Protected Member Functions | Friends | Related Functions | List of all members
lsst::meas::extensions::psfex::PsfexPsf Class Reference

Represent a PSF as a linear combination of PSFEX (== Karhunen-Loeve) basis functions. More...

#include <PsfexPsf.h>

Inheritance diagram for lsst::meas::extensions::psfex::PsfexPsf:
lsst::afw::table::io::PersistableFacade< PsfexPsf > 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

Public Types

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

Public Member Functions

 PsfexPsf (lsst::meas::extensions::psfex::Psf const &psf, lsst::geom::Point2D const &averagePosition=lsst::geom::Point2D())
 Constructor for a PsfexPsf. More...
 
virtual ~PsfexPsf ()
 
virtual boost::shared_ptr< lsst::afw::detection::Psfclone () const
 Polymorphic deep copy; should usually be unnecessary as Psfs are immutable.x. More...
 
virtual boost::shared_ptr< afw::detection::Psfresized (int width, int height) const
 Return a clone with specified kernel dimensions. More...
 
virtual lsst::geom::Point2D getAveragePosition () const
 Return average position of stars; used as default position. More...
 
boost::shared_ptr< lsst::afw::math::LinearCombinationKernel const > getKernel (lsst::geom::Point2D=lsst::geom::Point2D(std::numeric_limits< double >::quiet_NaN())) const
 Return the PSF's basis functions as a spatially-invariant LinearCombinationKernel with unit weights. More...
 
virtual bool isPersistable () const noexcept override
 Is this object persistable? More...
 
void write (lsst::afw::table::io::OutputArchiveHandle &handle) const
 Write the object to one or more catalogs. More...
 
std::shared_ptr< typehandling::Storable > cloneStorable () const override
 Polymorphic deep-copy. More...
 
std::shared_ptr< ImagecomputeImage (lsst::geom::Point2D position=makeNullPoint(), 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. More...
 
std::shared_ptr< ImagecomputeKernelImage (lsst::geom::Point2D position=makeNullPoint(), image::Color color=image::Color(), ImageOwnerEnum owner=COPY) const
 Return an Image of the PSF, in a form suitable for convolution. More...
 
double computePeak (lsst::geom::Point2D position=makeNullPoint(), image::Color color=image::Color()) const
 Return the peak value of the PSF image. More...
 
double computeApertureFlux (double radius, lsst::geom::Point2D position=makeNullPoint(), image::Color color=image::Color()) const
 Compute the "flux" of the Psf model within a circular aperture of the given radius. More...
 
geom::ellipses::Quadrupole computeShape (lsst::geom::Point2D position=makeNullPoint(), image::Color color=image::Color()) const
 Compute the ellipse corresponding to the second moments of the Psf. More...
 
std::shared_ptr< math::Kernel const > getLocalKernel (lsst::geom::Point2D position=makeNullPoint(), image::Color color=image::Color()) const
 Return a FixedKernel corresponding to the Psf image at the given point. More...
 
image::Color getAverageColor () const
 Return the average Color of the stars used to construct the Psf. More...
 
lsst::geom::Box2I computeBBox (lsst::geom::Point2D position=makeNullPoint(), image::Color color=image::Color()) const
 Return the bounding box of the image returned by computeKernelImage() More...
 
std::size_t getCacheCapacity () const
 Return the capacity of the caches. More...
 
void setCacheCapacity (std::size_t capacity)
 Set the capacity of the caches. More...
 
virtual std::string toString () const
 Create a string representation of this object (optional operation). More...
 
virtual std::size_t hash_value () const
 Return a hash of this object (optional operation). More...
 
virtual bool equals (Storable const &other) const noexcept
 Compare this object to another Storable. 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< PsfexPsfreadFits (fits::Fits &fitsfile)
 Read an object from an already open FITS object. More...
 
static std::shared_ptr< PsfexPsfreadFits (std::string const &fileName, int hdu=fits::DEFAULT_HDU)
 Read an object from a regular FITS file. More...
 
static std::shared_ptr< PsfexPsfreadFits (fits::MemFileManager &manager, int hdu=fits::DEFAULT_HDU)
 Read an object from a FITS file in memory. More...
 
static std::shared_ptr< PsfexPsfdynamicCast (std::shared_ptr< Persistable > const &ptr)
 Dynamically cast a shared_ptr. More...
 
static std::shared_ptr< ImagePsfreadFits (fits::Fits &fitsfile)
 Read an object from an already open FITS object. More...
 
static std::shared_ptr< ImagePsfreadFits (std::string const &fileName, int hdu=fits::DEFAULT_HDU)
 Read an object from a regular FITS file. More...
 
static std::shared_ptr< ImagePsfreadFits (fits::MemFileManager &manager, int hdu=fits::DEFAULT_HDU)
 Read an object from a FITS file in memory. More...
 
static std::shared_ptr< PsfreadFits (fits::Fits &fitsfile)
 Read an object from an already open FITS object. More...
 
static std::shared_ptr< PsfreadFits (std::string const &fileName, int hdu=fits::DEFAULT_HDU)
 Read an object from a regular FITS file. More...
 
static std::shared_ptr< PsfreadFits (fits::MemFileManager &manager, int hdu=fits::DEFAULT_HDU)
 Read an object from a FITS file in memory. More...
 
static std::shared_ptr< ImagePsfdynamicCast (std::shared_ptr< Persistable > const &ptr)
 Dynamically cast a shared_ptr. More...
 
static std::shared_ptr< PsfdynamicCast (std::shared_ptr< Persistable > const &ptr)
 Dynamically cast a shared_ptr. More...
 
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. More...
 

Protected Types

typedef io::OutputArchiveHandle OutputArchiveHandle
 

Protected Member Functions

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
 

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

Friends

class detail::PsfexPsfFactory
 

Related Functions

(Note that these are not member functions.)

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

Detailed Description

Represent a PSF as a linear combination of PSFEX (== Karhunen-Loeve) basis functions.

Definition at line 39 of file PsfexPsf.h.

Member Typedef Documentation

◆ Image

Image type returned by computeImage.

Definition at line 83 of file Psf.h.

◆ OutputArchiveHandle

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

Definition at line 108 of file Persistable.h.

◆ Pixel

typedef math::Kernel::Pixel lsst::afw::detection::Psf::Pixel
inherited

Pixel type of Image returned by computeImage.

Definition at line 82 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 86 of file Psf.h.

86  {
87  COPY = 0,
88  INTERNAL = 1
92  };

Constructor & Destructor Documentation

◆ PsfexPsf()

lsst::meas::extensions::psfex::PsfexPsf::PsfexPsf ( lsst::meas::extensions::psfex::Psf const &  psf,
lsst::geom::Point2D const &  averagePosition = lsst::geom::Point2D() 
)
explicit

Constructor for a PsfexPsf.

Parameters
[in]psf[in] Psfex PSF model to be wrapped into an LSST Psf
[in]averagePosition[in] Average position of stars used to construct the Psf.

Definition at line 64 of file PsfexPsf.cc.

67  : ImagePsf(), _averagePosition(averagePosition),
68  _size(psf.impl->dim),
69  _comp(psf.impl->npix),
70  _context(psf.impl->poly->ndim)
71 
72 {
73  _poly = poly_copy(psf.impl->poly);
74 
75  _pixstep = psf.impl->pixstep;
76 
77  std::copy(psf.impl->size, psf.impl->size + psf.impl->dim, _size.begin());
78 
79  std::copy(psf.impl->comp, psf.impl->comp + psf.impl->npix, _comp.begin());
80 
81  for (int i = 0; i != psf.impl->poly->ndim; ++i) {
82  _context[i].first = psf.impl->contextoffset[i];
83  _context[i].second = psf.impl->contextscale[i];
84  }
85 }

◆ ~PsfexPsf()

lsst::meas::extensions::psfex::PsfexPsf::~PsfexPsf ( )
virtual

Definition at line 98 of file PsfexPsf.cc.

99 {
100  poly_end(_poly);
101 }

Member Function Documentation

◆ clone()

boost::shared_ptr< afw::detection::Psf > lsst::meas::extensions::psfex::PsfexPsf::clone ( ) const
virtual

Polymorphic deep copy; should usually be unnecessary as Psfs are immutable.x.

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

Definition at line 104 of file PsfexPsf.cc.

104  {
105  return std::make_shared<PsfexPsf>(*this);
106 }

◆ cloneStorable()

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

Polymorphic deep copy; should usually be unnecessary as Psfs are immutable.x.

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 116 of file Psf.h.

116 { return clone(); }

◆ computeApertureFlux()

double Psf::computeApertureFlux ( double  radius,
lsst::geom::Point2D  position = makeNullPoint(),
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 to evaluate the PSF at; defaults to getAveragePosition().
[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 150 of file Psf.cc.

150  {
151  if (isPointNull(position)) position = getAveragePosition();
152  if (color.isIndeterminate()) color = getAverageColor();
153  return doComputeApertureFlux(radius, position, color);
154 }

◆ computeBBox()

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

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

Definition at line 128 of file Psf.cc.

128  {
129  if (isPointNull(position)) position = getAveragePosition();
130  if (color.isIndeterminate()) color = getAverageColor();
131  return doComputeBBox(position, color);
132 }

◆ computeImage()

std::shared_ptr< Psf::Image > Psf::computeImage ( lsst::geom::Point2D  position = makeNullPoint(),
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 to evaluate the PSF at; defaults to getAveragePosition().
[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 102 of file Psf.cc.

103  {
104  if (isPointNull(position)) position = getAveragePosition();
105  if (color.isIndeterminate()) color = getAverageColor();
106  std::shared_ptr<Psf::Image> result = (*_imageCache)(
107  detail::PsfCacheKey(position, color),
108  [this](detail::PsfCacheKey const &key) { return doComputeImage(key.position, key.color); });
109  if (owner == COPY) {
110  result = std::make_shared<Image>(*result, true);
111  }
112  return result;
113 }

◆ computeKernelImage()

std::shared_ptr< Psf::Image > Psf::computeKernelImage ( lsst::geom::Point2D  position = makeNullPoint(),
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 to evaluate the PSF at; defaults to getAveragePosition().
[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 115 of file Psf.cc.

116  {
117  if (_isFixed || isPointNull(position)) position = getAveragePosition();
118  if (_isFixed || color.isIndeterminate()) color = getAverageColor();
119  std::shared_ptr<Psf::Image> result = (*_kernelImageCache)(
120  detail::PsfCacheKey(position, color),
121  [this](detail::PsfCacheKey const &key) { return doComputeKernelImage(key.position, key.color); });
122  if (owner == COPY) {
123  result = std::make_shared<Image>(*result, true);
124  }
125  return result;
126 }

◆ computePeak()

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

Return the peak value of the PSF image.

Parameters
[in]positionPosition to evaluate the PSF at; defaults to getAveragePosition().
[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 143 of file Psf.cc.

143  {
144  if (isPointNull(position)) position = getAveragePosition();
145  if (color.isIndeterminate()) color = getAverageColor();
147  return (*image)(-image->getX0(), -image->getY0());
148 }

◆ computeShape()

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

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

Parameters
[in]positionPosition to evaluate the PSF at; defaults to getAveragePosition().
[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 156 of file Psf.cc.

156  {
157  if (isPointNull(position)) position = getAveragePosition();
158  if (color.isIndeterminate()) color = getAverageColor();
159  return doComputeShape(position, color);
160 }

◆ 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  afw::image::Image<double> const& image(*computeKernelImage(position, color, INTERNAL));
48 
49  geom::Point2D const center(0.0, 0.0);
50  afw::geom::ellipses::Axes const axes(radius, radius);
51  base::ApertureFluxResult result = base::ApertureFluxAlgorithm::computeSincFlux(
52  image, afw::geom::ellipses::Ellipse(axes, center), base::ApertureFluxControl());
53  return result.instFlux;
54 }

◆ 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 56 of file ImagePsf.cc.

57  {
58  PTR(Image) image = computeKernelImage(position, color, INTERNAL);
60  *image, geom::Point2D(0.0, 0.0) // image has origin at the center
61  )
62  .getShape();
63 }

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

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

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

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

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

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

◆ 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::image::PhotoCalib, lsst::afw::geom::SkyWcs, lsst::afw::geom::polygon::Polygon, lsst::afw::image::Filter, lsst::afw::image::VisitInfo, and lsst::afw::typehandling::StorableHelper< Base >.

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 233 of file Psf.h.

233 { return image::Color(); }

◆ getAveragePosition()

virtual lsst::geom::Point2D lsst::meas::extensions::psfex::PsfexPsf::getAveragePosition ( ) const
inlinevirtual

Return average position of stars; used as default position.

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

Definition at line 60 of file PsfexPsf.h.

60 { return _averagePosition; }

◆ 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 170 of file Psf.cc.

170 { return _kernelImageCache->capacity(); }

◆ getKernel()

boost::shared_ptr< afw::math::LinearCombinationKernel const > lsst::meas::extensions::psfex::PsfexPsf::getKernel ( lsst::geom::Point2D  position = lsst::geom::Point2D(std::numeric_limits<double>::quiet_NaN())) const

Return the PSF's basis functions as a spatially-invariant LinearCombinationKernel with unit weights.

Definition at line 114 of file PsfexPsf.cc.

115 {
116  double pos[MAXCONTEXT];
117  int const ndim = _context.size();
118  if (ndim != 2) { // we're only handling spatial variation for now
120  str(boost::format("Only spatial variation (ndim == 2) is supported; saw %d")
121  % ndim));
122 
123  }
124  // where we want to evaluate the basis function's weights
125  if (!std::isfinite(position[0])) {
126  position = _averagePosition;
127  }
128 
129  for (int i = 0; i < ndim; ++i) {
130  pos[i] = (position[i] - _context[i].first)/_context[i].second;
131  }
132 
133  poly_func(_poly, pos); // evaluate polynomial
134 
135  int const w = _size[0], h = _size[1];
136  std::vector<float> fullresIm(w*h); // accumulate full-resolution image into this buffer
137  /*
138  * Create a fixed Kernel out of each component, and then create a LinearCombinationKernel from them
139  */
140  const int nbasis = _size.size() > 2 ? _size[2] : 1; // number of basis functions
141  afw::math::KernelList kernels; kernels.reserve(nbasis); // the insides of the LinearCombinationKernel
142  std::vector<double> weights; weights.reserve(nbasis);
143 
144  float const vigstep = 1/_pixstep;
145  float const dx = 0.0, dy = 0.0;
146 
147  geom::Box2I bbox = _doComputeBBox(position, geom::Point2D(0, 0));
148  afw::detection::Psf::Image kim(bbox); // a basis function image, to be copied into a FixedKernel
149 
150  int sampleW = bbox.getWidth();
151  int sampleH = bbox.getHeight();
152 
153  std::vector<float> sampledBasis(sampleW*sampleH);
154 
155  for (int i = 0; i != nbasis; ++i) {
156  /*
157  * Resample the basis function onto the output resolution (and potentially subpixel offset)
158  */
159  vignet_resample(const_cast<float *>(&_comp[i*w*h]), w, h,
160  &sampledBasis[0], sampleW, sampleH,
161  -dx*vigstep, -dy*vigstep, vigstep, 1.0);
162  //
163  // And copy it into place
164  //
165  {
166  float *pl = &sampledBasis[0];
167  for (int y = 0; y != sampleH; ++y) {
168  for (int x = 0; x != sampleW; ++x) {
169  kim(x, y) = *pl++;
170  }
171  }
172  }
173 
174  kernels.push_back(std::make_shared<afw::math::FixedKernel>(kim));
175  weights.push_back(_poly->basis[i]);
176  }
177 
178  _kernel = std::make_shared<afw::math::LinearCombinationKernel>(kernels, weights);
179 
180  return _kernel;
181 }

◆ getLocalKernel()

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

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

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

Definition at line 134 of file Psf.cc.

135  {
136  if (isPointNull(position)) position = getAveragePosition();
137  if (color.isIndeterminate()) color = getAverageColor();
138  // FixedKernel ctor will deep copy image, so we can use INTERNAL.
140  return std::make_shared<math::FixedKernel>(*image);
141 }

◆ 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::typehandling::StorableHelper< Base >, lsst::afw::image::Filter, lsst::afw::geom::polygon::Polygon, and lsst::afw::image::VisitInfo.

Definition at line 44 of file Storable.cc.

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

◆ isPersistable()

virtual bool lsst::meas::extensions::psfex::PsfexPsf::isPersistable ( ) const
inlineoverridevirtualnoexcept

Is this object persistable?

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

Definition at line 69 of file PsfexPsf.h.

69 { 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.

183  {
184  return dynamicCast(Persistable::_readFits(fitsfile));
185  }

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

183  {
184  return dynamicCast(Persistable::_readFits(fitsfile));
185  }

◆ readFits() [3/9]

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

Read an object from an already open FITS object.

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

Definition at line 183 of file Persistable.h.

183  {
184  return dynamicCast(Persistable::_readFits(fitsfile));
185  }

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

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

◆ readFits() [5/9]

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

Read an object from a FITS file in memory.

Parameters
[in]managerManager for the memory to read from.
[in]hduHDU to read, where 0 is the primary. The special value of afw::fits::DEFAULT_HDU skips the primary HDU if it is empty.

Definition at line 205 of file Persistable.h.

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

◆ readFits() [6/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.

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

◆ readFits() [7/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.

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

◆ readFits() [8/9]

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

Read an object from a regular FITS file.

Parameters
[in]fileNameName of the file to read.
[in]hduHDU to read, where 0 is the primary. The special value of afw::fits::DEFAULT_HDU skips the primary HDU if it is empty.

Definition at line 194 of file Persistable.h.

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

◆ readFits() [9/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.

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

◆ 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 86 of file Psf.cc.

89  {
90  // "ir" : (integer, residual)
91  std::pair<int, double> const irX = image::positionToIndex(position.getX(), true);
92  std::pair<int, double> const irY = image::positionToIndex(position.getY(), true);
93 
94  if (irX.second != 0.0 || irY.second != 0.0) {
95  im = math::offsetImage(*im, irX.second, irY.second, warpAlgorithm, warpBuffer);
96  }
97 
98  im->setXY0(irX.first + im->getX0(), irY.first + im->getY0());
99  return im;
100 }

◆ resized()

boost::shared_ptr< afw::detection::Psf > lsst::meas::extensions::psfex::PsfexPsf::resized ( int  width,
int  height 
) const
virtual

Return a clone with specified kernel dimensions.

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

Definition at line 109 of file PsfexPsf.cc.

109  {
110  throw LSST_EXCEPT(pex::exceptions::LogicError, "Not Implemented");
111 }

◆ 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 172 of file Psf.cc.

172  {
173  _imageCache->reserve(capacity);
174  _kernelImageCache->reserve(capacity);
175 }

◆ 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\n 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::image::PhotoCalib, lsst::afw::geom::SkyWcs, lsst::afw::geom::polygon::Polygon, 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::extensions::psfex::PsfexPsf::write ( lsst::afw::table::io::OutputArchiveHandle handle) const
virtual

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 479 of file PsfexPsf.cc.

479  {
480  // First write the dimensions of the various arrays to an HDU, so we can construct the second
481  // HDU using them when we read the Psf back
482  {
483  PsfexPsfSchema1 const keys;
484  afw::table::BaseCatalog cat = handle.makeCatalog(keys.schema);
485  PTR(afw::table::BaseRecord) record = cat.addNew();
486 
487  // Sizes in _poly
488  record->set(keys.ndim, _poly->ndim);
489  record->set(keys.ngroup, _poly->ngroup);
490  record->set(keys.ncoeff, _poly->ncoeff);
491  // Other scalars
492  record->set(keys.averagePosition, _averagePosition);
493  record->set(keys._pixstep, _pixstep);
494  record->set(keys._size_size, _size.size());
495  record->set(keys._comp_size, _comp.size());
496  record->set(keys._context_size, _context.size());
497 
498  handle.saveCatalog(cat);
499  }
500  // Now we can write the data
501  {
502  PsfexPsfSchema2 const keys(_poly->ndim, _poly->ngroup, _poly->ncoeff,
503  _size.size(), _comp.size(), _context.size());
504  afw::table::BaseCatalog cat = handle.makeCatalog(keys.schema);
505  // _poly
506  PTR(afw::table::BaseRecord) record = cat.addNew();
507  {
508  int *begin = record->getElement(keys.group);
509  std::copy(_poly->group, _poly->group + _poly->ndim, begin);
510  }
511  {
512  int *begin = record->getElement(keys.degree);
513  std::copy(_poly->degree, _poly->degree + _poly->ngroup, begin);
514  }
515  {
516  double *begin = record->getElement(keys.basis);
517  std::copy(_poly->basis, _poly->basis + _poly->ncoeff, begin);
518  }
519  {
520  double *begin = record->getElement(keys.coeff);
521  std::copy(_poly->coeff, _poly->coeff + _poly->ncoeff, begin);
522  }
523  // vectors
524  {
525  int *begin = record->getElement(keys._size);
526  std::copy(_size.begin(), _size.end(), begin);
527  }
528  {
529  float *begin = record->getElement(keys._comp);
530  std::copy(_comp.begin(), _comp.end(), begin);
531  }
532  {
533  double *begin1 = record->getElement(keys._context_first);
534  double *begin2 = record->getElement(keys._context_second);
535  for (unsigned int i = 0; i != _context.size(); ++i) {
536  begin1[i] = _context[i].first;
537  begin2[i] = _context[i].second;
538  }
539  }
540 
541  handle.saveCatalog(cat);
542  }
543 }

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

◆ 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 Function Documentation

◆ detail::PsfexPsfFactory

friend class detail::PsfexPsfFactory
friend

Definition at line 41 of file PsfexPsf.h.

◆ 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:
y
int y
Definition: SpanSet.cc:49
lsst::meas::algorithms::ImagePsf::ImagePsf
ImagePsf(bool isFixed=false)
Definition: ImagePsf.h:42
lsst::afw::detection::Psf::getAveragePosition
virtual lsst::geom::Point2D getAveragePosition() const
Return the average position of the stars used to construct the Psf.
Definition: Psf.cc:168
lsst::afw::image
Backwards-compatibility support for depersisting the old Calib (FluxMag0/FluxMag0Err) objects.
Definition: imageAlgorithm.dox:1
lsst::afw::table::io::Persistable::writeFits
void writeFits(std::string const &fileName, std::string const &mode="w") const
Write the object to a regular FITS file.
Definition: Persistable.cc:24
std::shared_ptr
STL class.
lsst::afw::detection::Psf::INTERNAL
@ INTERNAL
An internal image will be returned without copying.
Definition: Psf.h:88
lsst::afw::table._match.second
second
Definition: _match.py:78
lsst::afw::image::positionToIndex
int positionToIndex(double pos)
Convert image position to nearest integer index.
Definition: ImageUtils.h:69
std::pair
std::vector::reserve
T reserve(T... args)
std::vector< float >
std::vector::size
T size(T... args)
psf
Key< int > psf
Definition: Exposure.cc:65
lsst::afw::table::io::OutputArchiveHandle::saveCatalog
void saveCatalog(BaseCatalog const &catalog)
Save a catalog in the archive.
Definition: OutputArchive.cc:211
lsst::afw::detection::Psf::clone
virtual std::shared_ptr< Psf > clone() const =0
Polymorphic deep-copy.
pex.config.history.format
def format(config, name=None, writeSourceLine=True, prefix="", verbose=False)
Definition: history.py:174
astshim.keyMap.keyMapContinued.keys
def keys(self)
Definition: keyMapContinued.py:6
lsst::afw::image::Color::isIndeterminate
bool isIndeterminate() const noexcept
Whether the color is the special value that indicates that it is unspecified.
Definition: Color.h:37
averagePosition
table::PointKey< double > averagePosition
Definition: PsfexPsf.cc:333
lsst.pipe.tasks.processCcdWithFakes.radius
radius
Definition: processCcdWithFakes.py:301
std::vector::push_back
T push_back(T... args)
lsst::afw::table::io::OutputArchiveHandle::makeCatalog
BaseCatalog makeCatalog(Schema const &schema)
Return a new, empty catalog with the given schema.
Definition: OutputArchive.cc:207
lsst::meas::base::SdssShapeAlgorithm::computeAdaptiveMoments
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.
std::isfinite
T isfinite(T... args)
image
afw::table::Key< afw::table::Array< ImagePixelT > > image
Definition: HeavyFootprint.cc:216
lsst::afw::detection::Psf::Image
image::Image< Pixel > Image
Image type returned by computeImage.
Definition: Psf.h:83
lsst::afw::fits::Fits::AUTO_CHECK
@ AUTO_CHECK
Definition: fits.h:308
lsst::afw::fits::Fits::AUTO_CLOSE
@ AUTO_CLOSE
Definition: fits.h:307
x
double x
Definition: ChebyshevBoundedField.cc:277
other
ItemVariant const * other
Definition: Schema.cc:56
std::copy
T copy(T... args)
ptr
uint64_t * ptr
Definition: RangeSet.cc:88
lsst::afw::detection::Psf::COPY
@ COPY
The image will be copied before returning; caller will own it.
Definition: Psf.h:87
result
py::object result
Definition: _schema.cc:429
LSST_EXCEPT
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
PTR
#define PTR(...)
Definition: base.h:41
lsst::afw::math::KernelList
std::vector< std::shared_ptr< Kernel > > KernelList
Definition: Kernel.h:509
ndim
table::Key< int > ndim
Definition: PsfexPsf.cc:325
lsst::meas::base::ApertureFluxAlgorithm::computeSincFlux
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.
Definition: ApertureFlux.cc:146
os
std::ostream * os
Definition: Schema.cc:746
lsst::afw::image::ImageBase::getX0
int getX0() const
Return the image's column-origin.
Definition: ImageBase.h:343
lsst::afw::table::BaseCatalog
CatalogT< BaseRecord > BaseCatalog
Definition: fwd.h:71
lsst::afw::detection::Psf::getAverageColor
image::Color getAverageColor() const
Return the average Color of the stars used to construct the Psf.
Definition: Psf.h:233
lsst::pex::exceptions::InvalidParameterError
Reports invalid arguments.
Definition: Runtime.h:66
lsst::afw::table::io::PersistableFacade< Psf >::dynamicCast
static std::shared_ptr< Psf > dynamicCast(std::shared_ptr< Persistable > const &ptr)
Dynamically cast a shared_ptr.
Definition: Persistable.cc:18
std::vector::begin
T begin(T... args)
key
Key< U > key
Definition: Schema.cc:281
lsst::geom::Point< double, 2 >
lsst::afw::math::offsetImage
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.
Definition: offsetImage.cc:41
lsst::meas::base::ShapeResult::getShape
Shape const getShape() const
Return an afw::geom::ellipses object corresponding to xx, yy, xy.
Definition: ShapeUtilities.cc:41
lsst::geom::Box2I
An integer coordinate rectangle.
Definition: Box.h:55
lsst::afw::image::ImageBase::setXY0
void setXY0(lsst::geom::Point2I const origin)
Set the ImageBase's origin.
Definition: ImageBase.h:471
w
double w
Definition: CoaddPsf.cc:69
lsst::afw::image::ImageBase::getY0
int getY0() const
Return the image's row-origin.
Definition: ImageBase.h:351
std::vector::end
T end(T... args)
lsst::afw::detection::Psf::computeKernelImage
std::shared_ptr< Image > computeKernelImage(lsst::geom::Point2D position=makeNullPoint(), image::Color color=image::Color(), ImageOwnerEnum owner=COPY) const
Return an Image of the PSF, in a form suitable for convolution.
Definition: Psf.cc:115
lsst::afw::detection::Psf::doComputeImage
virtual std::shared_ptr< Image > doComputeImage(lsst::geom::Point2D const &position, image::Color const &color) const
This virtual member is protected (rather than private) so that python-implemented derived classes may...
Definition: Psf.cc:162
bbox
AmpInfoBoxKey bbox
Definition: Amplifier.cc:117
lsst::afw::image::Color
Describe the colour of a source.
Definition: Color.h:26