LSST Applications  21.0.0-172-gfb10e10a+18fedfabac,22.0.0+297cba6710,22.0.0+80564b0ff1,22.0.0+8d77f4f51a,22.0.0+a28f4c53b1,22.0.0+dcf3732eb2,22.0.1-1-g7d6de66+2a20fdde0d,22.0.1-1-g8e32f31+297cba6710,22.0.1-1-geca5380+7fa3b7d9b6,22.0.1-12-g44dc1dc+2a20fdde0d,22.0.1-15-g6a90155+515f58c32b,22.0.1-16-g9282f48+790f5f2caa,22.0.1-2-g92698f7+dcf3732eb2,22.0.1-2-ga9b0f51+7fa3b7d9b6,22.0.1-2-gd1925c9+bf4f0e694f,22.0.1-24-g1ad7a390+a9625a72a8,22.0.1-25-g5bf6245+3ad8ecd50b,22.0.1-25-gb120d7b+8b5510f75f,22.0.1-27-g97737f7+2a20fdde0d,22.0.1-32-gf62ce7b1+aa4237961e,22.0.1-4-g0b3f228+2a20fdde0d,22.0.1-4-g243d05b+871c1b8305,22.0.1-4-g3a563be+32dcf1063f,22.0.1-4-g44f2e3d+9e4ab0f4fa,22.0.1-42-gca6935d93+ba5e5ca3eb,22.0.1-5-g15c806e+85460ae5f3,22.0.1-5-g58711c4+611d128589,22.0.1-5-g75bb458+99c117b92f,22.0.1-6-g1c63a23+7fa3b7d9b6,22.0.1-6-g50866e6+84ff5a128b,22.0.1-6-g8d3140d+720564cf76,22.0.1-6-gd805d02+cc5644f571,22.0.1-8-ge5750ce+85460ae5f3,master-g6e05de7fdc+babf819c66,master-g99da0e417a+8d77f4f51a,w.2021.48
LSST Data Management Base Package
Public Member Functions | Static Public Member Functions | Protected Types | Protected Member Functions | Static Protected Member Functions | Related Functions | List of all members
lsst::afw::geom::SkyWcs Class Referencefinal

A 2-dimensional celestial WCS that transform pixels to ICRS RA/Dec, using the LSST standard for pixels. More...

#include <SkyWcs.h>

Inheritance diagram for lsst::afw::geom::SkyWcs:
lsst::afw::table::io::PersistableFacade< SkyWcs > lsst::afw::typehandling::Storable lsst::afw::table::io::Persistable

Public Member Functions

 SkyWcs (SkyWcs const &)=default
 
 SkyWcs (SkyWcs &&)=default
 
SkyWcsoperator= (SkyWcs const &)=delete
 
SkyWcsoperator= (SkyWcs &&)=delete
 
 ~SkyWcs () override=default
 
bool operator== (SkyWcs const &other) const
 Equality is based on the string representations being equal. More...
 
bool operator!= (SkyWcs const &other) const
 
 SkyWcs (daf::base::PropertySet &metadata, bool strip=false)
 Construct a SkyWcs from FITS keywords. More...
 
 SkyWcs (ast::FrameDict const &frameDict)
 Construct a SkyWcs from an ast::FrameDict. More...
 
std::shared_ptr< SkyWcscopyAtShiftedPixelOrigin (lsst::geom::Extent2D const &shift) const
 Return a copy of this SkyWcs with the pixel origin shifted by the specified amount. More...
 
std::shared_ptr< daf::base::PropertyListgetFitsMetadata (bool precise=false) const
 Return the WCS as FITS WCS metadata. More...
 
lsst::geom::Angle getPixelScale (lsst::geom::Point2D const &pixel) const
 Get the pixel scale at the specified pixel position. More...
 
lsst::geom::Angle getPixelScale () const
 Get the pixel scale at the pixel origin. More...
 
lsst::geom::Point2D getPixelOrigin () const
 Get the pixel origin, in pixels, using the LSST convention. More...
 
lsst::geom::SpherePoint getSkyOrigin () const
 Get the sky origin, the celestial fiducial point. More...
 
Eigen::Matrix2d getCdMatrix (lsst::geom::Point2D const &pixel) const
 Get the 2x2 CD matrix at the specified pixel position. More...
 
Eigen::Matrix2d getCdMatrix () const
 Get the 2x2 CD matrix at the pixel origin. More...
 
std::shared_ptr< SkyWcsgetTanWcs (lsst::geom::Point2D const &pixel) const
 Get a local TAN WCS approximation to this WCS at the specified pixel position. More...
 
std::shared_ptr< const ast::FrameDictgetFrameDict () const
 Get the contained FrameDict. More...
 
std::shared_ptr< const TransformPoint2ToSpherePointgetTransform () const
 Get a TransformPoint2ToSpherePoint that transforms pixels to sky in the forward direction and sky to pixels in the inverse direction. More...
 
bool isFlipped () const
 Does the WCS follow the convention of North=Up, East=Left? More...
 
bool isPersistable () const noexcept override
 Return true if this particular object can be persisted using afw::table::io. More...
 
lsst::geom::AffineTransform linearizePixelToSky (lsst::geom::SpherePoint const &coord, lsst::geom::AngleUnit const &skyUnit) const
 Return the local linear approximation to pixelToSky at a point given in sky coordinates. More...
 
lsst::geom::AffineTransform linearizePixelToSky (lsst::geom::Point2D const &pixel, lsst::geom::AngleUnit const &skyUnit) const
 Return the local linear approximation to pixelToSky at a point given in pixel coordinates. More...
 
lsst::geom::AffineTransform linearizeSkyToPixel (lsst::geom::SpherePoint const &coord, lsst::geom::AngleUnit const &skyUnit) const
 Return the local linear approximation to skyToPixel at a point given in sky coordinates. More...
 
lsst::geom::AffineTransform linearizeSkyToPixel (lsst::geom::Point2D const &pixel, lsst::geom::AngleUnit const &skyUnit) const
 Return the local linear approximation to skyToPixel at a point given in pixel coordinates. More...
 
lsst::geom::SpherePoint pixelToSky (lsst::geom::Point2D const &pixel) const
 Compute sky position(s) from pixel position(s) More...
 
lsst::geom::SpherePoint pixelToSky (double x, double y) const
 
std::vector< lsst::geom::SpherePointpixelToSky (std::vector< lsst::geom::Point2D > const &pixels) const
 
lsst::geom::Point2D skyToPixel (lsst::geom::SpherePoint const &sky) const
 Compute pixel position(s) from sky position(s) More...
 
std::vector< lsst::geom::Point2DskyToPixel (std::vector< lsst::geom::SpherePoint > const &sky) const
 
bool hasFitsApproximation () const
 Does this SkyWcs have an approximate SkyWcs that can be represented as standard FITS WCS? More...
 
bool isFits () const
 Return true getFitsMetadata(true) will succeed, false if not. More...
 
void writeStream (std::ostream &os) const
 Serialize this SkyWcs to an output stream. More...
 
std::string writeString () const
 Serialize this SkyWcs to a string, using the same format as writeStream. More...
 
std::shared_ptr< typehandling::StorablecloneStorable () const override
 Create a new SkyWcs that is a copy of this one. More...
 
std::string toString () const override
 Create a string representation of this object. More...
 
bool equals (typehandling::Storable const &other) const noexcept override
 Compare this object to another Storable. More...
 
virtual std::size_t hash_value () const
 Return a hash of this object (optional operation). 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::string getShortClassName ()
 
static std::shared_ptr< SkyWcsreadStream (std::istream &is)
 Deserialize a SkyWcs from an input stream. More...
 
static std::shared_ptr< SkyWcsreadString (std::string &str)
 Deserialize a SkyWcs from a string, using the same format as readStream. More...
 
static std::shared_ptr< SkyWcsreadFits (fits::Fits &fitsfile)
 Read an object from an already open FITS object. More...
 
static std::shared_ptr< SkyWcsreadFits (std::string const &fileName, int hdu=fits::DEFAULT_HDU)
 Read an object from a regular FITS file. More...
 
static std::shared_ptr< SkyWcsreadFits (fits::MemFileManager &manager, int hdu=fits::DEFAULT_HDU)
 Read an object from a FITS file in memory. More...
 
static std::shared_ptr< SkyWcsdynamicCast (std::shared_ptr< Persistable > const &ptr)
 Dynamically cast a shared_ptr. More...
 

Protected Types

using OutputArchiveHandle = io::OutputArchiveHandle
 

Protected Member Functions

std::string getPersistenceName () const override
 Return the unique name used to persist this object and look up its factory. More...
 
std::string getPythonModule () const override
 Return the fully-qualified Python module that should be imported to guarantee that its factory is registered. More...
 
void write (OutputArchiveHandle &handle) const override
 Write the object to one or more catalogs. More...
 

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

Related Functions

(Note that these are not member functions.)

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

Detailed Description

A 2-dimensional celestial WCS that transform pixels to ICRS RA/Dec, using the LSST standard for pixels.

SkyWcs is an immutable object that can not only represent any standard FITS WCS, but can also contain arbitrary Transforms, e.g. to model optical distortion or pixel imperfections.

In order to make a SkyWcs that models optical distortion, say, it is usually simplest to start with a standard FITS WCS (such as a TAN WCS) as an approximation, then insert a transform that models optical distortion by calling makeModifiedWcs. However, it is also possible to build a SkyWcs entirely from transforms, if you prefer, by building an ast::FrameDict and constructing the SkyWcs from that.

Frames of reference

SkyWcs internally keeps track of the following frames of reference:

Pixel position standards

The LSST standard for pixel position is: 0,0 is the center of the lower left image pixel. The FITS standard for pixel position is: 1,1 is the center of the lower left image pixel.

LSST and FITS also use a different origin for subimages:

Internal details: the contained ast::FrameDict

SkyWcs contains an ast::FrameDict which transforms from pixels to sky (in radians) in the forward direction.

This FrameDict contains the named frames described in frames of reference, e.g. "SKY", "IWC", cameraGeom::PIXELS and possibly cameraGeom::ACTUAL_PIXELS. "SKY" is the current frame. If cameraGeom::ACTUAL_PIXELS is present then it is the base frame, otherwise cameraGeom::PIXELS is the base frame.

The "SKY" frame is of type ast::SkyFrame and has the following attributes:

The other frames are of type ast::Frame and have 2 axes.

Definition at line 117 of file SkyWcs.h.

Member Typedef Documentation

◆ OutputArchiveHandle

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

Definition at line 108 of file Persistable.h.

Constructor & Destructor Documentation

◆ SkyWcs() [1/4]

lsst::afw::geom::SkyWcs::SkyWcs ( SkyWcs const &  )
default

◆ SkyWcs() [2/4]

lsst::afw::geom::SkyWcs::SkyWcs ( SkyWcs &&  )
default

◆ ~SkyWcs()

lsst::afw::geom::SkyWcs::~SkyWcs ( )
overridedefault

◆ SkyWcs() [3/4]

lsst::afw::geom::SkyWcs::SkyWcs ( daf::base::PropertySet metadata,
bool  strip = false 
)
explicit

Construct a SkyWcs from FITS keywords.

Parameters
[in]metadataFITS header metadata
[in]stripIf true: strip items from metadata used to create the WCS, such as RADESYS, EQUINOX, CTYPE12, CRPIX12, CRVAL12, etc. Always keep keywords that might be wanted for other purpposes, including NAXIS12 and date-related keywords such as "DATE-OBS" and "TIMESYS" (but not "EQUINOX").
Exceptions
lsst::pex::exceptions::TypeErrorif the metadata does not describe a celestial WCS.

Definition at line 151 of file SkyWcs.cc.

152  : SkyWcs(detail::readLsstSkyWcs(metadata, strip)) {}
SkyWcs(SkyWcs const &)=default
bool strip
Definition: fits.cc:911
std::shared_ptr< ast::FrameDict > readLsstSkyWcs(daf::base::PropertySet &metadata, bool strip=true)
Read an LSST celestial WCS FrameDict from a FITS header.

◆ SkyWcs() [4/4]

lsst::afw::geom::SkyWcs::SkyWcs ( ast::FrameDict const &  frameDict)
explicit

Construct a SkyWcs from an ast::FrameDict.

This is the most general constructor; it can be used to define any celestial WCS. Note that in many cases the result will not be exactly representable as a FITS WCS.

Parameters
[in]frameDictAn ast::FrameDict that describes the transformation from pixels to sky. It must meet the requirements outlined in the contained ast::FrameDict.
Exceptions
lsst::pex::exceptions::TypeErrorif frameDict is missing any of the required frames of reference.

Definition at line 154 of file SkyWcs.cc.

154 : SkyWcs(_checkFrameDict(frameDict)) {}

Member Function Documentation

◆ cloneStorable()

std::shared_ptr< typehandling::Storable > lsst::afw::geom::SkyWcs::cloneStorable ( ) const
overridevirtual

Create a new SkyWcs that is a copy of this one.

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

Definition at line 348 of file SkyWcs.cc.

348  {
349  return std::make_unique<SkyWcs>(*this);
350 }

◆ copyAtShiftedPixelOrigin()

std::shared_ptr< SkyWcs > lsst::afw::geom::SkyWcs::copyAtShiftedPixelOrigin ( lsst::geom::Extent2D const &  shift) const

Return a copy of this SkyWcs with the pixel origin shifted by the specified amount.

new pixel origin = the old pixel origin + shift
Parameters
[in]shiftThe amount by which to shift the pixel origin (pixels)

Definition at line 208 of file SkyWcs.cc.

208  {
209  auto newToOldPixel = TransformPoint2ToPoint2(ast::ShiftMap({-shift[0], -shift[1]}));
210  return makeModifiedWcs(newToOldPixel, *this, true);
211 }
ShiftMap is a linear Mapping which shifts each axis by a specified constant value.
Definition: ShiftMap.h:40
std::shared_ptr< SkyWcs > makeModifiedWcs(TransformPoint2ToPoint2 const &pixelTransform, SkyWcs const &wcs, bool modifyActualPixels)
Create a new SkyWcs whose pixels are transformed by pixelTransform, as described below.
Definition: SkyWcs.cc:486
Transform< Point2Endpoint, Point2Endpoint > TransformPoint2ToPoint2
Definition: Transform.h:300

◆ dynamicCast()

Dynamically cast a shared_ptr.

Dynamically cast a shared pointer and raise on failure.

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

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

param[in] ptr The pointer to be cast.

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

Definition at line 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 }
py::object result
Definition: _schema.cc:429
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
uint64_t * ptr
Definition: RangeSet.cc:88

◆ equals()

bool lsst::afw::geom::SkyWcs::equals ( typehandling::Storable const &  other) const
overridevirtualnoexcept

Compare this object to another Storable.

Returns
*this == other if other is a SkyWcs; otherwise false.

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

Definition at line 373 of file SkyWcs.cc.

373  {
374  return singleClassEquals(*this, other);
375 }
static bool singleClassEquals(T const &lhs, Storable const &rhs)
Test if a Storable is of a particular class and equal to another object.
Definition: Storable.h:151

◆ getCdMatrix() [1/2]

Eigen::Matrix2d lsst::afw::geom::SkyWcs::getCdMatrix ( ) const

Get the 2x2 CD matrix at the pixel origin.

The elements are in degrees

Definition at line 199 of file SkyWcs.cc.

199 { return getCdMatrix(getPixelOrigin()); }
Eigen::Matrix2d getCdMatrix() const
Get the 2x2 CD matrix at the pixel origin.
Definition: SkyWcs.cc:199
lsst::geom::Point2D getPixelOrigin() const
Get the pixel origin, in pixels, using the LSST convention.
Definition: SkyWcs.h:215

◆ getCdMatrix() [2/2]

Eigen::Matrix2d lsst::afw::geom::SkyWcs::getCdMatrix ( lsst::geom::Point2D const &  pixel) const

Get the 2x2 CD matrix at the specified pixel position.

The elements are in degrees

Definition at line 193 of file SkyWcs.cc.

193  {
194  auto const pixelToIwc = getFrameDict()->getMapping(ast::FrameSet::BASE, "IWC");
195  auto const pixelToIwcTransform = TransformPoint2ToPoint2(*pixelToIwc);
196  return pixelToIwcTransform.getJacobian(pixel);
197 }
table::PointKey< int > pixel
std::shared_ptr< Mapping > getMapping(int from, std::string const &to) const
Variant of FrameSet::getMapping with the second frame specified by domain.
Definition: FrameDict.h:162
static constexpr int BASE
index of base frame
Definition: FrameSet.h:104
std::shared_ptr< const ast::FrameDict > getFrameDict() const
Get the contained FrameDict.
Definition: SkyWcs.cc:262

◆ getFitsMetadata()

std::shared_ptr< daf::base::PropertyList > lsst::afw::geom::SkyWcs::getFitsMetadata ( bool  precise = false) const

Return the WCS as FITS WCS metadata.

Parameters
[in]preciseFail if the WCS cannot be accurately represented as FITS metadata? If False then return an approximation. For now that approximation is pure TAN but as of DM-13170 it will be a fit TAN-SIP. The required precision is set by constant TIGHT_FITS_TOL in SkyWcs.cc

FITS representations of WCS are described in "Representations of World Coordinates in FITS" by Greisen and Calabretta and several related papers.

The required precision is hard-coded as constant TIGHT_FITS_TOL in SkyWcs.cc

Exceptions
lsst::pex::exceptions::RuntimeErrorif precise is true and AST cannot represent this WCS as a FITS WCS to sufficient precision.

Definition at line 213 of file SkyWcs.cc.

213  {
214  // Make a FrameSet that maps from GRID to SKY; GRID = the base frame (PIXELS or ACTUAL_PIXELS) + 1
215  auto const gridToPixel = ast::ShiftMap({-1.0, -1.0});
216  auto thisDict = getFrameDict();
217  auto const pixelToIwc = thisDict->getMapping(ast::FrameSet::BASE, "IWC");
218  auto const iwcToSky = thisDict->getMapping("IWC", "SKY");
219  auto const gridToSky = gridToPixel.then(*pixelToIwc).then(*iwcToSky);
220  ast::FrameSet frameSet(ast::Frame(2, "Domain=GRID"), gridToSky, *thisDict->getFrame("SKY", false));
221 
222  // Write frameSet to a FitsChan and extract the metadata
224  os << "Encoding=FITS-WCS, CDMatrix=1, FitsAxisOrder=<copy>, FitsTol=" << TIGHT_FITS_TOL;
225  ast::StringStream strStream;
226  ast::FitsChan fitsChan(strStream, os.str());
227  int const nObjectsWritten = fitsChan.write(frameSet);
229  if (nObjectsWritten == 0) {
230  if (precise) {
232  "Could not represent this SkyWcs using FITS-WCS metadata");
233  } else {
234  // An exact representation could not be written, so try to write a local TAN approximation;
235  // set precise true to avoid an infinite loop, should something go wrong
236  auto tanWcs = getTanWcs(getPixelOrigin());
237  header = tanWcs->getFitsMetadata(true);
238  }
239  } else {
240  header = detail::getPropertyListFromFitsChan(fitsChan);
241  }
242 
243  // Remove DATE-OBS, MJD-OBS: AST writes these if the EQUINOX is set, but we set them via other mechanisms.
244  header->remove("DATE-OBS");
245  header->remove("MJD-OBS");
246 
247  // If CD matrix is present, explicitly set any missing entries to zero, as a convenience to the user
248  bool const hasCd11 = header->exists("CD1_1");
249  bool const hasCd12 = header->exists("CD1_2");
250  bool const hasCd21 = header->exists("CD2_1");
251  bool const hasCd22 = header->exists("CD2_2");
252  if (hasCd11 || hasCd12 || hasCd21 || hasCd22) {
253  if (!hasCd11) header->set("CD1_1", 0.0, "Transformation matrix element");
254  if (!hasCd12) header->set("CD1_2", 0.0, "Transformation matrix element");
255  if (!hasCd21) header->set("CD2_1", 0.0, "Transformation matrix element");
256  if (!hasCd22) header->set("CD2_2", 0.0, "Transformation matrix element");
257  }
258 
259  return header;
260 }
std::ostream * os
Definition: Schema.cc:557
table::Key< table::Array< std::uint8_t > > frameSet
A specialized form of Channel which reads and writes FITS header cards.
Definition: FitsChan.h:202
Frame is used to represent a coordinate system.
Definition: Frame.h:157
A FrameSet consists of a set of one or more Frames (which describe coordinate systems),...
Definition: FrameSet.h:99
SeriesMap then(Mapping const &next) const
Return a series compound mapping this(first(input)) containing shallow copies of the original.
Definition: Mapping.cc:37
String-based source and sink for channels.
Definition: Stream.h:180
std::shared_ptr< SkyWcs > getTanWcs(lsst::geom::Point2D const &pixel) const
Get a local TAN WCS approximation to this WCS at the specified pixel position.
Definition: SkyWcs.cc:201
Reports errors that are due to events beyond the control of the program.
Definition: Runtime.h:104
std::shared_ptr< daf::base::PropertyList > getPropertyListFromFitsChan(ast::FitsChan &fitsChan)
Copy values from an AST FitsChan into a PropertyList.

◆ getFrameDict()

std::shared_ptr< const ast::FrameDict > lsst::afw::geom::SkyWcs::getFrameDict ( ) const

Get the contained FrameDict.

The base frame will be cameraGeom::PIXELS or cameraGeom::ACTUAL_PIXELS and the current frame will be SKY, so the forward transform goes from pixels (using the LSST zero convention) to sky ICRS RA, Dec (in radians). For more details see the contained ast::FrameDict

Definition at line 262 of file SkyWcs.cc.

262 { return _frameDict; }

◆ getPersistenceName()

std::string lsst::afw::geom::SkyWcs::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 377 of file SkyWcs.cc.

377 { return getSkyWcsPersistenceName(); }

◆ getPixelOrigin()

lsst::geom::Point2D lsst::afw::geom::SkyWcs::getPixelOrigin ( ) const
inline

Get the pixel origin, in pixels, using the LSST convention.

This is CRPIX1 - 1, CRPIX2 -1 in FITS terminology

Definition at line 215 of file SkyWcs.h.

215 { return _pixelOrigin; };

◆ getPixelScale() [1/2]

lsst::geom::Angle lsst::afw::geom::SkyWcs::getPixelScale ( ) const
inline

Get the pixel scale at the pixel origin.

The scale is the square root of the area of the specified pixel on the sky.

The value is cached, so this is a cheap call.

Definition at line 208 of file SkyWcs.h.

208 { return _pixelScaleAtOrigin; };

◆ getPixelScale() [2/2]

lsst::geom::Angle lsst::afw::geom::SkyWcs::getPixelScale ( lsst::geom::Point2D const &  pixel) const

Get the pixel scale at the specified pixel position.

The scale is the square root of the area of the specified pixel on the sky.

Warning
Unlike getPixelScale() the value is not cached, even if pixel = pixel origin.

Definition at line 158 of file SkyWcs.cc.

158  {
159  // Compute pixVec containing the pixel position and two nearby points
160  // (use a vector so all three points can be converted to sky in a single call)
161  double const side = 1.0;
163  pixel,
164  pixel + lsst::geom::Extent2D(side, 0),
165  pixel + lsst::geom::Extent2D(0, side),
166  };
167 
168  auto skyVec = pixelToSky(pixVec);
169 
170  // Work in 3-space to avoid RA wrapping and pole issues
171  auto skyLL = skyVec[0].getVector();
172  auto skyDx = skyVec[1].getVector() - skyLL;
173  auto skyDy = skyVec[2].getVector() - skyLL;
174 
175  // Compute pixel scale in radians = sqrt(pixel area in radians^2)
176  // pixel area in radians^2 = area of parallelogram with sides skyDx, skyDy = |skyDx cross skyDy|
177  // Use squared norm to avoid two square roots
178  double skyAreaSq = skyDx.cross(skyDy).getSquaredNorm();
179  return (std::pow(skyAreaSq, 0.25) / side) * lsst::geom::radians;
180 }
lsst::geom::SpherePoint pixelToSky(lsst::geom::Point2D const &pixel) const
Compute sky position(s) from pixel position(s)
Definition: SkyWcs.h:334
Extent< double, 2 > Extent2D
Definition: Extent.h:400
constexpr AngleUnit radians
constant with units of radians
Definition: Angle.h:108
T pow(T... args)

◆ getPythonModule()

std::string lsst::afw::geom::SkyWcs::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 379 of file SkyWcs.cc.

379 { return "lsst.afw.geom"; }

◆ getShortClassName()

std::string lsst::afw::geom::SkyWcs::getShortClassName ( )
static

Definition at line 294 of file SkyWcs.cc.

294 { return "SkyWcs"; };

◆ getSkyOrigin()

lsst::geom::SpherePoint lsst::afw::geom::SkyWcs::getSkyOrigin ( ) const

Get the sky origin, the celestial fiducial point.

This is CRVAL1, CRVAL2 in FITS terminology

Definition at line 182 of file SkyWcs.cc.

182  {
183  // CRVAL is stored as the SkyRef property of the sky frame (the current frame of the SkyWcs)
184  auto skyFrame = std::dynamic_pointer_cast<ast::SkyFrame>(
185  getFrameDict()->getFrame(ast::FrameDict::CURRENT, false)); // false: do not copy
186  if (!skyFrame) {
187  throw LSST_EXCEPT(pex::exceptions::LogicError, "Current frame is not a SkyFrame");
188  }
189  auto const crvalRad = skyFrame->getSkyRef();
190  return lsst::geom::SpherePoint(crvalRad[0] * lsst::geom::radians, crvalRad[1] * lsst::geom::radians);
191 }
static constexpr int CURRENT
index of current frame
Definition: FrameSet.h:105
lsst::geom::SpherePoint SpherePoint
Definition: misc.h:34

◆ getTanWcs()

std::shared_ptr< SkyWcs > lsst::afw::geom::SkyWcs::getTanWcs ( lsst::geom::Point2D const &  pixel) const

Get a local TAN WCS approximation to this WCS at the specified pixel position.

Definition at line 201 of file SkyWcs.cc.

201  {
202  auto const crval = pixelToSky(pixel);
203  auto const cdMatrix = getCdMatrix(pixel);
204  auto metadata = makeSimpleWcsMetadata(pixel, crval, cdMatrix);
205  return std::make_shared<SkyWcs>(*metadata);
206 }
table::PointKey< double > crval
Definition: OldWcs.cc:128
std::shared_ptr< daf::base::PropertyList > makeSimpleWcsMetadata(lsst::geom::Point2D const &crpix, lsst::geom::SpherePoint const &crval, Eigen::Matrix2d const &cdMatrix, std::string const &projection="TAN")
Make FITS metadata for a simple FITS WCS (one with no distortion).
Definition: wcsUtils.cc:169

◆ getTransform()

std::shared_ptr<const TransformPoint2ToSpherePoint> lsst::afw::geom::SkyWcs::getTransform ( ) const
inline

Get a TransformPoint2ToSpherePoint that transforms pixels to sky in the forward direction and sky to pixels in the inverse direction.

Definition at line 257 of file SkyWcs.h.

257 { return _transform; };

◆ hasFitsApproximation()

bool lsst::afw::geom::SkyWcs::hasFitsApproximation ( ) const
inline

Does this SkyWcs have an approximate SkyWcs that can be represented as standard FITS WCS?

This feature is not yet implemented, so hasFitsApproximation is always false.

Definition at line 364 of file SkyWcs.h.

364 { return false; }

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

Definition at line 44 of file Storable.cc.

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

◆ isFits()

bool lsst::afw::geom::SkyWcs::isFits ( ) const

Return true getFitsMetadata(true) will succeed, false if not.

In other words, true indicates that the WCS can be accurately represented using FITS WCS metadata.

Definition at line 264 of file SkyWcs.cc.

264  {
265  try {
266  getFitsMetadata(true);
267  } catch (const lsst::pex::exceptions::RuntimeError&) {
268  return false;
269  } catch (const std::runtime_error&) {
270  return false;
271  }
272  return true;
273 }
std::shared_ptr< daf::base::PropertyList > getFitsMetadata(bool precise=false) const
Return the WCS as FITS WCS metadata.
Definition: SkyWcs.cc:213

◆ isFlipped()

bool lsst::afw::geom::SkyWcs::isFlipped ( ) const

Does the WCS follow the convention of North=Up, East=Left?

Returns
False/true if E is along -X/+X when the N/E axes are rotated so that N is along image +Y.
Exceptions
lsst::pex::exceptions::RuntimeErrorif the parity cannot be determined because the CD matrix is singular.

Definition at line 296 of file SkyWcs.cc.

296  {
297  double det = getCdMatrix().determinant();
298  if (det == 0) {
299  throw(LSST_EXCEPT(pex::exceptions::RuntimeError, "CD matrix is singular"));
300  }
301  return (det > 0);
302 }

◆ isPersistable()

bool lsst::afw::geom::SkyWcs::isPersistable ( ) const
inlineoverridevirtualnoexcept

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

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

Definition at line 268 of file SkyWcs.h.

268 { return true; }

◆ linearizePixelToSky() [1/2]

lsst::geom::AffineTransform lsst::afw::geom::SkyWcs::linearizePixelToSky ( lsst::geom::Point2D const &  pixel,
lsst::geom::AngleUnit const &  skyUnit 
) const

Return the local linear approximation to pixelToSky at a point given in pixel coordinates.

The local linear approximation is defined such the following is true (ignoring floating-point errors):

wcs.linearizePixelToSky(pixel, skyUnit)(pixel) == wcs.pixelToSky(pixel).getPosition(skyUnit)

(lsst::geom::AffineTransform::operator() applies the transform in the forward direction)

Parameters
[in]pixelPosition in pixel coordinates where transform is desired.
[in]skyUnitUnits to use for sky coordinates; units of matrix elements will be skyUnits/pixel.

Definition at line 279 of file SkyWcs.cc.

280  {
281  return _linearizePixelToSky(pix, pixelToSky(pix), skyUnit);
282 }

◆ linearizePixelToSky() [2/2]

lsst::geom::AffineTransform lsst::afw::geom::SkyWcs::linearizePixelToSky ( lsst::geom::SpherePoint const &  coord,
lsst::geom::AngleUnit const &  skyUnit 
) const

Return the local linear approximation to pixelToSky at a point given in sky coordinates.

The local linear approximation is defined such the following is true (ignoring floating-point errors):

wcs.linearizePixelToSky(sky, skyUnit)(wcs.skyToPixel(sky)) == sky.getPosition(skyUnit);

(lsst::geom::AffineTransform::operator() applies the transform in the forward direction)

Parameters
[in]coordPosition in sky coordinates where transform is desired.
[in]skyUnitUnits to use for sky coordinates; units of matrix elements will be skyUnits/pixel.

Definition at line 275 of file SkyWcs.cc.

276  {
277  return _linearizePixelToSky(skyToPixel(coord), coord, skyUnit);
278 }
lsst::geom::Point2D skyToPixel(lsst::geom::SpherePoint const &sky) const
Compute pixel position(s) from sky position(s)
Definition: SkyWcs.h:349

◆ linearizeSkyToPixel() [1/2]

lsst::geom::AffineTransform lsst::afw::geom::SkyWcs::linearizeSkyToPixel ( lsst::geom::Point2D const &  pixel,
lsst::geom::AngleUnit const &  skyUnit 
) const

Return the local linear approximation to skyToPixel at a point given in pixel coordinates.

The local linear approximation is defined such the following is true (ignoring floating-point errors):

wcs.linearizeSkyToPixel(pixel, skyUnit)(wcs.pixelToSky(pixel).getPosition(skyUnit)) == pixel

(lsst::geom::AffineTransform::operator() applies the transform in the forward direction)

Parameters
[in]pixelPosition in pixel coordinates where transform is desired.
[in]skyUnitUnits to use for sky coordinates; units of matrix elements will be pixels/skyUnit.

Definition at line 289 of file SkyWcs.cc.

290  {
291  return _linearizeSkyToPixel(pix, pixelToSky(pix), skyUnit);
292 }

◆ linearizeSkyToPixel() [2/2]

lsst::geom::AffineTransform lsst::afw::geom::SkyWcs::linearizeSkyToPixel ( lsst::geom::SpherePoint const &  coord,
lsst::geom::AngleUnit const &  skyUnit 
) const

Return the local linear approximation to skyToPixel at a point given in sky coordinates.

The local linear approximation is defined such the following is true (ignoring floating-point errors):

wcs.linearizeSkyToPixel(sky, skyUnit)(sky.getPosition(skyUnit)) == wcs.skyToPixel(sky)

(lsst::geom::AffineTransform::operator() applies the transform in the forward direction)

Parameters
[in]coordPosition in sky coordinates where transform is desired.
[in]skyUnitUnits to use for sky coordinates; units of matrix elements will be pixels/skyUnit.

Definition at line 284 of file SkyWcs.cc.

285  {
286  return _linearizeSkyToPixel(skyToPixel(coord), coord, skyUnit);
287 }

◆ operator!=()

bool lsst::afw::geom::SkyWcs::operator!= ( SkyWcs const &  other) const
inline

Definition at line 136 of file SkyWcs.h.

136 { return !(*this == other); }

◆ operator=() [1/2]

SkyWcs& lsst::afw::geom::SkyWcs::operator= ( SkyWcs &&  )
delete

◆ operator=() [2/2]

SkyWcs& lsst::afw::geom::SkyWcs::operator= ( SkyWcs const &  )
delete

◆ operator==()

bool lsst::afw::geom::SkyWcs::operator== ( SkyWcs const &  other) const

Equality is based on the string representations being equal.

Two SkyWcs constructed the same way will be equal, and a SkyWcs that has been saved and restored will be equal to the original. However, it is possible to construct two SkyWcs that behave identically as far as transforming points go, but will compare as unequal due to subtle internal differences, such as a contained ast::Mapping that has a different ID in one SkyWcs than another.

Thus equality is primarily useful for testing persistence.

Definition at line 156 of file SkyWcs.cc.

156 { return writeString() == other.writeString(); }
std::string writeString() const
Serialize this SkyWcs to a string, using the same format as writeStream.
Definition: SkyWcs.cc:342

◆ pixelToSky() [1/3]

lsst::geom::SpherePoint lsst::afw::geom::SkyWcs::pixelToSky ( double  x,
double  y 
) const
inline

Definition at line 337 of file SkyWcs.h.

337  {
338  return pixelToSky(lsst::geom::Point2D(x, y));
339  }
double x
int y
Definition: SpanSet.cc:48

◆ pixelToSky() [2/3]

lsst::geom::SpherePoint lsst::afw::geom::SkyWcs::pixelToSky ( lsst::geom::Point2D const &  pixel) const
inline

Compute sky position(s) from pixel position(s)

Definition at line 334 of file SkyWcs.h.

334  {
335  return _transform->applyForward(pixel);
336  }

◆ pixelToSky() [3/3]

std::vector<lsst::geom::SpherePoint> lsst::afw::geom::SkyWcs::pixelToSky ( std::vector< lsst::geom::Point2D > const &  pixels) const
inline

Definition at line 340 of file SkyWcs.h.

340  {
341  return _transform->applyForward(pixels);
342  }

◆ readFits() [1/3]

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

Read an object from an already open FITS object.

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

Definition at line 183 of file Persistable.h.

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

◆ readFits() [2/3]

static std::shared_ptr<SkyWcs > lsst::afw::table::io::PersistableFacade< SkyWcs >::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() [3/3]

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

◆ readStream()

std::shared_ptr< SkyWcs > lsst::afw::geom::SkyWcs::readStream ( std::istream is)
static

Deserialize a SkyWcs from an input stream.

Parameters
[in]isinput stream from which to deserialize this SkyWcs

Definition at line 304 of file SkyWcs.cc.

304  {
305  int version;
306  is >> version;
307  if (version != 1) {
308  throw LSST_EXCEPT(pex::exceptions::TypeError, "Unsupported version " + std::to_string(version));
309  }
310  std::string shortClassName;
311  is >> shortClassName;
312  if (shortClassName != SkyWcs::getShortClassName()) {
314  os << "Class name in stream " << shortClassName << " != " << SkyWcs::getShortClassName();
315  throw LSST_EXCEPT(pex::exceptions::TypeError, os.str());
316  }
318  is >> hasFitsApproximation;
319  auto astStream = ast::Stream(&is, nullptr);
320  auto astObjectPtr = ast::Channel(astStream).read();
321  auto frameSet = std::dynamic_pointer_cast<ast::FrameSet>(astObjectPtr);
322  if (!frameSet) {
324  os << "The AST serialization was read as a " << astObjectPtr->getClassName()
325  << " instead of a FrameSet";
326  throw LSST_EXCEPT(pex::exceptions::TypeError, os.str());
327  }
328  ast::FrameDict frameDict(*frameSet);
329  return std::make_shared<SkyWcs>(frameDict);
330 }
Channel provides input/output of AST objects.
Definition: Channel.h:60
std::shared_ptr< Object > read()
Read an object from a channel.
Definition: Channel.cc:52
A FrameSet whose frames can be referenced by domain name.
Definition: FrameDict.h:67
A stream for ast::Channel.
Definition: Stream.h:41
bool hasFitsApproximation() const
Does this SkyWcs have an approximate SkyWcs that can be represented as standard FITS WCS?
Definition: SkyWcs.h:364
static std::string getShortClassName()
Definition: SkyWcs.cc:294
T to_string(T... args)

◆ readString()

std::shared_ptr< SkyWcs > lsst::afw::geom::SkyWcs::readString ( std::string str)
static

Deserialize a SkyWcs from a string, using the same format as readStream.

Definition at line 332 of file SkyWcs.cc.

332  {
333  std::istringstream is(str);
334  return SkyWcs::readStream(is);
335 }
static std::shared_ptr< SkyWcs > readStream(std::istream &is)
Deserialize a SkyWcs from an input stream.
Definition: SkyWcs.cc:304

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

◆ skyToPixel() [1/2]

lsst::geom::Point2D lsst::afw::geom::SkyWcs::skyToPixel ( lsst::geom::SpherePoint const &  sky) const
inline

Compute pixel position(s) from sky position(s)

Definition at line 349 of file SkyWcs.h.

349  {
350  return _transform->applyInverse(sky);
351  }

◆ skyToPixel() [2/2]

std::vector<lsst::geom::Point2D> lsst::afw::geom::SkyWcs::skyToPixel ( std::vector< lsst::geom::SpherePoint > const &  sky) const
inline

Definition at line 352 of file SkyWcs.h.

352  {
353  return _transform->applyInverse(sky);
354  }

◆ toString()

std::string lsst::afw::geom::SkyWcs::toString ( ) const
overridevirtual

Create a string representation of this object.

An example of the output might look like this (note the newlines):

FITS standard SkyWcs:
Sky Origin: (0.000000, +45.000000)
Pixel Origin: (100, 100)
Pixel Scale: 1 arcsec/pixel

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

Definition at line 352 of file SkyWcs.cc.

352  {
354  if (isFits()) {
355  os << "FITS standard SkyWcs:";
356  } else {
357  os << "Non-standard SkyWcs (Frames: ";
358  // Print the frames in index order (frames are numbered from 1).
359  std::string delimiter = "";
360  for (size_t i = 1; i <= getFrameDict()->getAllDomains().size(); ++i) {
361  os << delimiter << getFrameDict()->getFrame(i)->getDomain();
362  delimiter = ", ";
363  }
364  os << "): ";
365  }
366  std::string delimiter = "\n";
367  os << delimiter << "Sky Origin: " << getSkyOrigin();
368  os << delimiter << "Pixel Origin: " << getPixelOrigin();
369  os << delimiter << "Pixel Scale: " << getPixelScale().asArcseconds() << " arcsec/pixel";
370  return os.str();
371 }
std::set< std::string > getAllDomains() const
Get the domain names for all contained Frames (excluding frames with empty or defaulted domain names)...
Definition: FrameDict.cc:45
std::shared_ptr< Frame > getFrame(std::string const &domain, bool copy=true) const
Variant of getFrame(int, bool) where the frame is specified by domain name.
Definition: FrameDict.h:151
lsst::geom::Angle getPixelScale() const
Get the pixel scale at the pixel origin.
Definition: SkyWcs.h:208
lsst::geom::SpherePoint getSkyOrigin() const
Get the sky origin, the celestial fiducial point.
Definition: SkyWcs.cc:182
bool isFits() const
Return true getFitsMetadata(true) will succeed, false if not.
Definition: SkyWcs.cc:264
constexpr double asArcseconds() const noexcept
Return an Angle's value in arcseconds.
Definition: Angle.h:178
T size(T... args)

◆ write()

void lsst::afw::geom::SkyWcs::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 381 of file SkyWcs.cc.

381  {
382  SkyWcsPersistenceHelper const& keys = SkyWcsPersistenceHelper::get();
383  table::BaseCatalog cat = handle.makeCatalog(keys.schema);
384  std::shared_ptr<table::BaseRecord> record = cat.addNew();
385  record->set(keys.wcs, formatters::stringToBytes(writeString()));
386  handle.saveCatalog(cat);
387 }
ndarray::Array< std::uint8_t, 1, 1 > stringToBytes(std::string const &str)
Encode a std::string as a vector of uint8.
Definition: Utils.cc:161
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.
Definition: Persistable.cc:24

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

◆ writeStream()

void lsst::afw::geom::SkyWcs::writeStream ( std::ostream os) const

Serialize this SkyWcs to an output stream.

Version 1 format is as follows:

If and when fits approximation is supported, the approximate WCS will be written as a second FrameDict immediately following the first.

Parameters
[out]osoutput stream to which to serialize this SkyWcs

Definition at line 337 of file SkyWcs.cc.

337  {
338  os << SERIALIZATION_VERSION << " " << SkyWcs::getShortClassName() << " " << hasFitsApproximation() << " ";
339  getFrameDict()->show(os, false); // false = do not write comments
340 }
void show(std::ostream &os, bool showComments=true) const
Print a textual description the object to an ostream.
Definition: Object.cc:158

◆ writeString()

std::string lsst::afw::geom::SkyWcs::writeString ( ) const

Serialize this SkyWcs to a string, using the same format as writeStream.

Definition at line 342 of file SkyWcs.cc.

342  {
344  writeStream(os);
345  return os.str();
346 }
void writeStream(std::ostream &os) const
Serialize this SkyWcs to an output stream.
Definition: SkyWcs.cc:337

Friends And Related Function 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: