LSST Applications 27.0.0,g0265f82a02+469cd937ee,g02d81e74bb+21ad69e7e1,g1470d8bcf6+cbe83ee85a,g2079a07aa2+e67c6346a6,g212a7c68fe+04a9158687,g2305ad1205+94392ce272,g295015adf3+81dd352a9d,g2bbee38e9b+469cd937ee,g337abbeb29+469cd937ee,g3939d97d7f+72a9f7b576,g487adcacf7+71499e7cba,g50ff169b8f+5929b3527e,g52b1c1532d+a6fc98d2e7,g591dd9f2cf+df404f777f,g5a732f18d5+be83d3ecdb,g64a986408d+21ad69e7e1,g858d7b2824+21ad69e7e1,g8a8a8dda67+a6fc98d2e7,g99cad8db69+f62e5b0af5,g9ddcbc5298+d4bad12328,ga1e77700b3+9c366c4306,ga8c6da7877+71e4819109,gb0e22166c9+25ba2f69a1,gb6a65358fc+469cd937ee,gbb8dafda3b+69d3c0e320,gc07e1c2157+a98bf949bb,gc120e1dc64+615ec43309,gc28159a63d+469cd937ee,gcf0d15dbbd+72a9f7b576,gdaeeff99f8+a38ce5ea23,ge6526c86ff+3a7c1ac5f1,ge79ae78c31+469cd937ee,gee10cc3b42+a6fc98d2e7,gf1cff7945b+21ad69e7e1,gfbcc870c63+9a11dc8c8f
LSST Data Management Base Package
Loading...
Searching...
No Matches
Exposure.cc
Go to the documentation of this file.
1// -*- LSST-C++ -*- // fixed format comment for emacs
2
3/*
4 * LSST Data Management System
5 * Copyright 2008, 2009, 2010 LSST Corporation.
6 *
7 * This product includes software developed by the
8 * LSST Project (http://www.lsst.org/).
9 *
10 * This program is free software: you can redistribute it and/or modify
11 * it under the terms of the GNU General Public License as published by
12 * the Free Software Foundation, either version 3 of the License, or
13 * (at your option) any later version.
14 *
15 * This program is distributed in the hope that it will be useful,
16 * but WITHOUT ANY WARRANTY; without even the implied warranty of
17 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18 * GNU General Public License for more details.
19 *
20 * You should have received a copy of the LSST License Statement and
21 * the GNU General Public License along with this program. If not,
22 * see <http://www.lsstcorp.org/LegalNotices/>.
23 */
24
25#include <memory>
26#include <stdexcept>
27#include <sstream>
28#include <cstdint>
29
30#include "lsst/pex/exceptions.h"
36#include "lsst/afw/fits.h"
38
39namespace lsst {
40namespace afw {
41namespace image {
42
43// CLASS CONSTRUCTORS and DESTRUCTOR
44
45template <typename ImageT, typename MaskT, typename VarianceT>
46Exposure<ImageT, MaskT, VarianceT>::Exposure(unsigned int width, unsigned int height,
48 : _maskedImage(width, height), _info(new ExposureInfo(wcs)) {}
49
50template <typename ImageT, typename MaskT, typename VarianceT>
53 : _maskedImage(dimensions), _info(new ExposureInfo(wcs)) {}
54
55template <typename ImageT, typename MaskT, typename VarianceT>
59
60template <typename ImageT, typename MaskT, typename VarianceT>
63 : _maskedImage(maskedImage), _info(new ExposureInfo(wcs)) {}
64
65template <typename ImageT, typename MaskT, typename VarianceT>
67 : _maskedImage(maskedImage), _info(info ? info : std::make_shared<ExposureInfo>()) {}
68
69template <typename ImageT, typename MaskT, typename VarianceT>
71 : _maskedImage(src.getMaskedImage(), deep), _info(new ExposureInfo(*src.getInfo(), deep)) {}
72// Delegate to copy-constructor for backwards compatibility
73template <typename ImageT, typename MaskT, typename VarianceT>
75
76template <typename ImageT, typename MaskT, typename VarianceT>
78 ImageOrigin const origin, bool const deep)
79 : _maskedImage(src.getMaskedImage(), bbox, origin, deep),
80 _info(new ExposureInfo(*src.getInfo(), deep)) {}
81
82template <typename ImageT, typename MaskT, typename VarianceT>
84 ImageOrigin origin, bool conformMasks, bool allowUnsafe)
85 : _maskedImage(), _info(new ExposureInfo()) {
86 ExposureFitsReader reader(fileName);
87 *this = reader.read<ImageT, MaskT, VarianceT>(bbox, origin, conformMasks, allowUnsafe);
88}
89
90template <typename ImageT, typename MaskT, typename VarianceT>
92 ImageOrigin origin, bool conformMasks, bool allowUnsafe)
93 : _maskedImage(), _info(new ExposureInfo()) {
94 ExposureFitsReader reader(manager);
95 *this = reader.read<ImageT, MaskT, VarianceT>(bbox, origin, conformMasks, allowUnsafe);
96}
97
98template <typename ImageT, typename MaskT, typename VarianceT>
100 ImageOrigin origin, bool conformMasks, bool allowUnsafe) {
101 ExposureFitsReader reader(&fitsFile);
102 *this = reader.read<ImageT, MaskT, VarianceT>(bbox, origin, conformMasks, allowUnsafe);
103}
104
105template <typename ImageT, typename MaskT, typename VarianceT>
107
108// SET METHODS
110template <typename ImageT, typename MaskT, typename VarianceT>
112 _maskedImage = maskedImage;
113}
114
115template <typename ImageT, typename MaskT, typename VarianceT>
117 lsst::geom::Point2I old(_maskedImage.getXY0());
118 if (_info->hasWcs()) {
119 auto shift = lsst::geom::Extent2D(origin - old);
120 auto newWcs = _info->getWcs()->copyAtShiftedPixelOrigin(shift);
121 _info->setWcs(newWcs);
122 }
123 _maskedImage.setXY0(origin);
124}
125
126template <typename ImageT, typename MaskT, typename VarianceT>
127Exposure<ImageT, MaskT, VarianceT> &Exposure<ImageT, MaskT, VarianceT>::operator=(Exposure const &) = default;
128template <typename ImageT, typename MaskT, typename VarianceT>
129Exposure<ImageT, MaskT, VarianceT> &Exposure<ImageT, MaskT, VarianceT>::operator=(Exposure &&) = default;
130
131// Write FITS
133template <typename ImageT, typename MaskT, typename VarianceT>
135 fits::Fits fitsfile(fileName, "w", fits::Fits::AUTO_CLOSE | fits::Fits::AUTO_CHECK);
136 writeFits(fitsfile);
137}
138
139template <typename ImageT, typename MaskT, typename VarianceT>
141 fits::Fits fitsfile(manager, "w", fits::Fits::AUTO_CLOSE | fits::Fits::AUTO_CHECK);
142 writeFits(fitsfile);
143}
144
145template <typename ImageT, typename MaskT, typename VarianceT>
147 writeFits(fitsfile, fits::ImageWriteOptions(*_maskedImage.getImage()),
148 fits::ImageWriteOptions(*_maskedImage.getMask()),
149 fits::ImageWriteOptions(*_maskedImage.getVariance()));
150}
151
152template <typename ImageT, typename MaskT, typename VarianceT>
154 fits::ImageWriteOptions const &imageOptions,
155 fits::ImageWriteOptions const &maskOptions,
156 fits::ImageWriteOptions const &varianceOptions) const {
157 fits::Fits fitsfile(fileName, "w", fits::Fits::AUTO_CLOSE | fits::Fits::AUTO_CHECK);
158 writeFits(fitsfile, imageOptions, maskOptions, varianceOptions);
159}
161template <typename ImageT, typename MaskT, typename VarianceT>
163 fits::ImageWriteOptions const &imageOptions,
164 fits::ImageWriteOptions const &maskOptions,
165 fits::ImageWriteOptions const &varianceOptions) const {
166 fits::Fits fitsfile(manager, "w", fits::Fits::AUTO_CLOSE | fits::Fits::AUTO_CHECK);
167 writeFits(fitsfile, imageOptions, maskOptions, varianceOptions);
170template <typename ImageT, typename MaskT, typename VarianceT>
172 fits::ImageWriteOptions const &imageOptions,
173 fits::ImageWriteOptions const &maskOptions,
174 fits::ImageWriteOptions const &varianceOptions) const {
175 ExposureInfo::FitsWriteData data = _info->_startWriteFits(getXY0());
176 _maskedImage.writeFits(fitsfile, imageOptions, maskOptions, varianceOptions, data.metadata,
177 data.imageMetadata, data.maskMetadata, data.varianceMetadata);
178 _info->_finishWriteFits(fitsfile, data);
179}
180
181namespace {
190template <class ExposureT>
191void _copyCommonPixels(ExposureT &destination, ExposureT const &source) {
192 lsst::geom::Box2I overlapBox = destination.getBBox();
193 overlapBox.clip(source.getBBox());
194
195 // MaskedImage::assign interprets empty bounding box as "whole image"
196 if (!overlapBox.isEmpty()) {
197 typename ExposureT::MaskedImageT overlapPixels(source.getMaskedImage(), overlapBox);
198 destination.getMaskedImage().assign(overlapPixels, overlapBox);
199 }
200}
201} // namespace
202
203template <typename ImageT, typename MaskT, typename VarianceT>
205 lsst::geom::SpherePoint const &center, lsst::geom::Extent2I const &size) const {
206 if (!hasWcs()) {
207 throw LSST_EXCEPT(pex::exceptions::LogicError, "Cannot look up source position without WCS.");
208 }
209 lsst::geom::Point2D pixelCenter = getWcs()->skyToPixel(center);
210
211 return getCutout(pixelCenter, size);
212}
213
214template <typename ImageT, typename MaskT, typename VarianceT>
216 lsst::geom::Point2D const &center, lsst::geom::Extent2I const &size) const {
217 if (!lsst::geom::Box2D(getBBox()).contains(center)) {
218 std::stringstream buffer;
219 buffer << "Point " << center << " lies at pixel " << center << ", which lies outside Exposure "
220 << getBBox();
222 }
223 if (size[0] <= 0 || size[1] <= 0) {
224 std::stringstream buffer;
225 buffer << "Cannot create bounding box with dimensions " << size;
227 }
229
230 // cutout must have independent ExposureInfo
231 auto copyInfo = std::make_shared<ExposureInfo>(*getInfo());
232 MaskedImageT blank(bbox); // Can't initialize Exposure with a temporary
233 blank = math::edgePixel<MaskedImageT>(
235 Exposure cutout(blank, copyInfo);
236
237 _copyCommonPixels(cutout, *this);
238 return cutout;
239}
240
241template <typename ImageT, typename MaskT, typename VarianceT>
246
247// Explicit instantiations
249template class Exposure<std::uint16_t>;
250template class Exposure<int>;
251template class Exposure<float>;
252template class Exposure<double>;
253template class Exposure<std::uint64_t>;
255} // namespace image
256} // namespace afw
257} // namespace lsst
AmpInfoBoxKey bbox
Definition Amplifier.cc:117
char * data
Definition BaseRecord.cc:61
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition Exception.h:48
std::shared_ptr< RecordT > src
Definition Match.cc:48
A simple struct that combines the two arguments that must be passed to most cfitsio routines and cont...
Definition fits.h:308
Lifetime-management for memory that goes into FITS memory files.
Definition fits.h:125
A FITS reader class for Exposures and their components.
Exposure< ImagePixelT, MaskPixelT, VariancePixelT > read(lsst::geom::Box2I const &bbox=lsst::geom::Box2I(), ImageOrigin origin=PARENT, bool conformMasks=false, bool allowUnsafe=false)
Read the full Exposure.
A class to contain the data, WCS, and other information needed to describe an image of the sky.
Definition Exposure.h:72
void writeFits(std::string const &fileName) const
Write an Exposure to a regular multi-extension FITS file.
Definition Exposure.cc:134
Exposure(unsigned int width, unsigned int height, std::shared_ptr< geom::SkyWcs const > wcs=std::shared_ptr< geom::SkyWcs const >())
Construct an Exposure with a blank MaskedImage of specified size (default 0x0) and a SkyWcs (which ma...
Definition Exposure.cc:46
A collection of all the things that make an Exposure different from a MaskedImage.
A class to manipulate images, masks, and variance as a single object.
Definition MaskedImage.h:74
A floating-point coordinate rectangle geometry.
Definition Box.h:413
An integer coordinate rectangle.
Definition Box.h:55
void clip(Box2I const &other) noexcept
Shrink this to ensure that other.contains(*this).
Definition Box.cc:189
bool isEmpty() const noexcept
Return true if the box contains no points.
Definition Box.h:213
static Box2I makeCenteredBox(Point2D const &center, Extent const &size)
Create a box centered as closely as possible on a particular point.
Definition Box.cc:97
Extent2I const getDimensions() const noexcept
Definition Box.h:186
Point2D const getCenter() const noexcept
1-d interval accessors
Definition Box.cc:93
Point in an unspecified spherical coordinate system.
Definition SpherePoint.h:57
Reports invalid arguments.
Definition Runtime.h:66
Reports errors in the logical structure of the program.
Definition Runtime.h:46
Extent< double, 2 > Extent2D
Definition Extent.h:400
STL namespace.
T str(T... args)
Options for writing an image to FITS.
Definition fits.h:223
typename ImageT::image_category image_category
Definition ImageBase.h:67