LSST Applications  21.0.0+04719a4bac,21.0.0-1-ga51b5d4+f5e6047307,21.0.0-11-g2b59f77+a9c1acf22d,21.0.0-11-ga42c5b2+86977b0b17,21.0.0-12-gf4ce030+76814010d2,21.0.0-13-g1721dae+760e7a6536,21.0.0-13-g3a573fe+768d78a30a,21.0.0-15-g5a7caf0+f21cbc5713,21.0.0-16-g0fb55c1+b60e2d390c,21.0.0-19-g4cded4ca+71a93a33c0,21.0.0-2-g103fe59+bb20972958,21.0.0-2-g45278ab+04719a4bac,21.0.0-2-g5242d73+3ad5d60fb1,21.0.0-2-g7f82c8f+8babb168e8,21.0.0-2-g8f08a60+06509c8b61,21.0.0-2-g8faa9b5+616205b9df,21.0.0-2-ga326454+8babb168e8,21.0.0-2-gde069b7+5e4aea9c2f,21.0.0-2-gecfae73+1d3a86e577,21.0.0-2-gfc62afb+3ad5d60fb1,21.0.0-25-g1d57be3cd+e73869a214,21.0.0-3-g357aad2+ed88757d29,21.0.0-3-g4a4ce7f+3ad5d60fb1,21.0.0-3-g4be5c26+3ad5d60fb1,21.0.0-3-g65f322c+e0b24896a3,21.0.0-3-g7d9da8d+616205b9df,21.0.0-3-ge02ed75+a9c1acf22d,21.0.0-4-g591bb35+a9c1acf22d,21.0.0-4-g65b4814+b60e2d390c,21.0.0-4-gccdca77+0de219a2bc,21.0.0-4-ge8a399c+6c55c39e83,21.0.0-5-gd00fb1e+05fce91b99,21.0.0-6-gc675373+3ad5d60fb1,21.0.0-64-g1122c245+4fb2b8f86e,21.0.0-7-g04766d7+cd19d05db2,21.0.0-7-gdf92d54+04719a4bac,21.0.0-8-g5674e7b+d1bd76f71f,master-gac4afde19b+a9c1acf22d,w.2021.13
LSST Data Management Base Package
Footprint.h
Go to the documentation of this file.
1 /*
2  * LSST Data Management System
3  * Copyright 2008-2016 AURA/LSST.
4  *
5  * This product includes software developed by the
6  * LSST Project (http://www.lsst.org/).
7  *
8  * This program is free software: you can redistribute it and/or modify
9  * it under the terms of the GNU General Public License as published by
10  * the Free Software Foundation, either version 3 of the License, or
11  * (at your option) any later version.
12  *
13  * This program is distributed in the hope that it will be useful,
14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16  * GNU General Public License for more details.
17  *
18  * You should have received a copy of the LSST License Statement and
19  * the GNU General Public License along with this program. If not,
20  * see <https://www.lsstcorp.org/LegalNotices/>.
21  */
22 #if !defined(LSST_DETECTION_FOOTPRINT_H)
23 #define LSST_DETECTION_FOOTPRINT_H
24 
25 #include <algorithm>
26 #include <list>
27 #include <set>
28 #include <cmath>
29 #include <cstdint>
30 #include <memory>
31 #include "ndarray.h"
32 #include "lsst/base.h"
35 #include "lsst/afw/geom/Span.h"
36 #include "lsst/afw/geom/SpanSet.h"
37 #include "lsst/geom/Box.h"
38 #include "lsst/afw/geom/ellipses.h"
40 #include "lsst/afw/geom/SkyWcs.h"
42 #include "lsst/afw/table/fwd.h"
43 #include "lsst/afw/table/Schema.h"
44 #include "lsst/afw/table/Key.h"
47 
48 namespace lsst {
49 namespace afw {
50 namespace detection {
51 
62 class Footprint : public afw::table::io::PersistableFacade<lsst::afw::detection::Footprint>,
64 public:
65  /*
66  * As a note there is no constructor which accepts a unique_ptr to a SpanSet.
67  * If someone did have a unique_ptr to a SpanSet it is possible to wrap the
68  * unique_ptr in a std::move when calling the constructor of a Footprint and
69  * the compiler will implicitly turn the unique_ptr into a shared_ptr and use
70  * that to construct the Footprint.
71  */
72 
80  explicit Footprint(std::shared_ptr<geom::SpanSet> inputSpans,
81  lsst::geom::Box2I const &region = lsst::geom::Box2I());
82 
91  explicit Footprint(std::shared_ptr<geom::SpanSet> inputSpans, afw::table::Schema const &peakSchema,
92  lsst::geom::Box2I const &region = lsst::geom::Box2I());
93 
96  explicit Footprint()
97  : _spans(std::make_shared<geom::SpanSet>()),
98  _peaks(PeakTable::makeMinimalSchema()),
99  _region(lsst::geom::Box2I()) {}
100 
101  Footprint(Footprint const &other) = default;
102  Footprint(Footprint &&) = default;
103 
104  Footprint &operator=(Footprint const &other) = default;
105  Footprint &operator=(Footprint &&) = default;
106 
107  ~Footprint() override = default;
108 
111  virtual bool isHeavy() const { return false; }
112 
115  std::shared_ptr<geom::SpanSet> getSpans() const { return _spans; }
116 
121  void setSpans(std::shared_ptr<geom::SpanSet> otherSpanSet);
122 
129  PeakCatalog &getPeaks() { return _peaks; }
130  const PeakCatalog &getPeaks() const { return _peaks; }
131 
138  std::shared_ptr<PeakRecord> addPeak(float fx, float fy, float value);
139 
149 
156  void setPeakSchema(afw::table::Schema const &peakSchema);
157 
166  void setPeakCatalog(PeakCatalog const &otherPeaks);
167 
173  std::size_t getArea() const { return _spans->getArea(); }
174 
180  lsst::geom::Point2D getCentroid() const { return _spans->computeCentroid(); }
181 
188  geom::ellipses::Quadrupole getShape() const { return _spans->computeShape(); }
189 
196  void shift(int dx, int dy);
197 
203  void shift(lsst::geom::ExtentI const &d) { shift(d.getX(), d.getY()); }
204 
208  lsst::geom::Box2I getBBox() const { return _spans->getBBox(); }
209 
213  lsst::geom::Box2I getRegion() const { return _region; }
214 
220  void setRegion(lsst::geom::Box2I const &region) { _region = region; }
221 
228  void clipTo(lsst::geom::Box2I const &bbox);
229 
235  bool contains(lsst::geom::Point2I const &pix) const;
236 
248  lsst::geom::Box2I const &region, bool doClip = true) const;
249 
258  lsst::geom::Box2I const &region, bool doClip = true) const;
259 
268  lsst::geom::Box2I const &region, bool doClip = true) const;
269 
278  lsst::geom::Box2I const &region, bool doClip = true) const;
279 
283  bool isPersistable() const noexcept override { return true; }
284 
294  void dilate(int r, geom::Stencil s = geom::Stencil::CIRCLE);
295 
304  void dilate(geom::SpanSet const &other);
305 
315  void erode(int r, geom::Stencil s = geom::Stencil::CIRCLE);
316 
325  void erode(geom::SpanSet const &other);
326 
330  void removeOrphanPeaks();
331 
335  bool isContiguous() const { return getSpans()->isContiguous(); };
336 
345 
351  bool operator==(Footprint const &other) const;
352 
353 protected:
357  std::string getPersistenceName() const override;
358 
362  inline std::string getPythonModule() const override { return "lsst.afw.detection"; }
363 
367  void write(OutputArchiveHandle &handle) const override;
368 
369  friend class FootprintFactory;
370 
379  static void readPeaks(afw::table::BaseCatalog const &, Footprint &);
380 
381 private:
382  friend class FootprintMerge;
383 
385  PeakCatalog _peaks;
386  lsst::geom::Box2I _region;
387 };
388 
394 std::shared_ptr<Footprint> mergeFootprints(Footprint const &footprint1, Footprint const &footprint2);
395 
405 } // namespace detection
406 } // namespace afw
407 } // namespace lsst
408 
409 #endif // LSST_DETECTION_FOOTPRINT_H
AmpInfoBoxKey bbox
Definition: Amplifier.cc:117
Key< Flag > const & target
ItemVariant const * other
Definition: Schema.cc:56
Key< U > key
Definition: Schema.cc:281
Basic LSST definitions.
Class to describe the properties of a detected object from an image.
Definition: Footprint.h:63
Footprint()
Constructor of a empty Footprint object.
Definition: Footprint.h:96
const PeakCatalog & getPeaks() const
Definition: Footprint.h:130
lsst::geom::Box2I getRegion() const
Return the corners of the MaskedImage the footprints live in.
Definition: Footprint.h:213
bool isPersistable() const noexcept override
Report if this object is persistable.
Definition: Footprint.h:283
virtual bool isHeavy() const
Indicates if this object is a HeavyFootprint.
Definition: Footprint.h:111
bool isContiguous() const
Reports if the Footprint is simply connected or has multiple components.
Definition: Footprint.h:335
std::shared_ptr< PeakRecord > addPeak(float fx, float fy, float value)
Convenience function to add a peak.
Definition: Footprint.cc:47
static std::unique_ptr< Footprint > readSpanSet(afw::table::BaseCatalog const &, afw::table::io::InputArchive const &)
Static method used to unpersist the SpanSet member of the Footprint class.
Definition: Footprint.cc:272
void shift(lsst::geom::ExtentI const &d)
Shift a Footprint by a given extent.
Definition: Footprint.h:203
void write(OutputArchiveHandle &handle) const override
Write an instance of a Footprint to an output Archive.
Definition: Footprint.cc:257
Footprint(Footprint const &other)=default
std::shared_ptr< Footprint > transform(std::shared_ptr< geom::SkyWcs > source, std::shared_ptr< geom::SkyWcs > target, lsst::geom::Box2I const &region, bool doClip=true) const
Transform the footprint from one WCS to another.
Definition: Footprint.cc:81
bool contains(lsst::geom::Point2I const &pix) const
Tests if a pixel postion falls inside the Footprint.
Definition: Footprint.cc:79
void setRegion(lsst::geom::Box2I const &region)
Set the corners of the MaskedImage wherein the footprints dwell.
Definition: Footprint.h:220
std::vector< std::shared_ptr< Footprint > > split() const
Split a multi-component Footprint into a vector of contiguous Footprints.
Definition: Footprint.cc:145
std::shared_ptr< geom::SpanSet > getSpans() const
Return a shared pointer to the SpanSet.
Definition: Footprint.h:115
void dilate(int r, geom::Stencil s=geom::Stencil::CIRCLE)
Dilate the Footprint with a defined kernel.
Definition: Footprint.cc:122
lsst::geom::Box2I getBBox() const
Return the Footprint's bounding box.
Definition: Footprint.h:208
void sortPeaks(afw::table::Key< float > const &key=afw::table::Key< float >())
Sort Peaks from most positive value to most negative.
Definition: Footprint.cc:57
void removeOrphanPeaks()
Remove peaks from the PeakCatalog that fall outside the area of the Footprint.
Definition: Footprint.cc:136
void setPeakSchema(afw::table::Schema const &peakSchema)
Set the Schema used by the PeakCatalog (will throw if PeakCatalog is not empty).
Definition: Footprint.cc:417
PeakCatalog & getPeaks()
Return the Peaks contained in this Footprint.
Definition: Footprint.h:129
~Footprint() override=default
std::size_t getArea() const
Return the number of pixels in this Footprint.
Definition: Footprint.h:173
void setSpans(std::shared_ptr< geom::SpanSet > otherSpanSet)
Sets the shared pointer to the SpanSet in the Footprint.
Definition: Footprint.cc:45
std::string getPersistenceName() const override
Return the name correspoinging ot the persistence type.
Definition: Footprint.cc:255
void clipTo(lsst::geom::Box2I const &bbox)
Clip the Footprint such that all values lie inside the supplied Bounding Box.
Definition: Footprint.cc:74
void erode(int r, geom::Stencil s=geom::Stencil::CIRCLE)
Erode the Footprint with a defined kernel.
Definition: Footprint.cc:126
std::string getPythonModule() const override
Return the python module the object will live in.
Definition: Footprint.h:362
lsst::geom::Point2D getCentroid() const
Return the Footprint's centroid.
Definition: Footprint.h:180
void shift(int dx, int dy)
Shift a Footprint by (dx, dy)
Definition: Footprint.cc:64
bool operator==(Footprint const &other) const
equality operator
Definition: Footprint.cc:162
static void readPeaks(afw::table::BaseCatalog const &, Footprint &)
Static method used to unpersist the PeakCatalog member of the Footprint class.
Definition: Footprint.cc:296
Footprint(Footprint &&)=default
geom::ellipses::Quadrupole getShape() const
Return the Footprint's shape (interpreted as an ellipse)
Definition: Footprint.h:188
Footprint & operator=(Footprint &&)=default
Footprint & operator=(Footprint const &other)=default
void setPeakCatalog(PeakCatalog const &otherPeaks)
Set the peakCatalog to a copy of the supplied catalog.
Definition: Footprint.cc:421
Table class for Peaks in Footprints.
Definition: Peak.h:102
A compact representation of a collection of pixels.
Definition: SpanSet.h:78
An ellipse core with quadrupole moments as parameters.
Definition: Quadrupole.h:47
Defines the fields and offsets for a table.
Definition: Schema.h:50
A multi-catalog archive object used to load table::io::Persistable objects.
Definition: InputArchive.h:31
A CRTP facade class for subclasses of Persistable.
Definition: Persistable.h:176
A base class for objects that can be persisted via afw::table::io Archive classes.
Definition: Persistable.h:74
io::OutputArchiveHandle OutputArchiveHandle
Definition: Persistable.h:108
An affine coordinate transformation consisting of a linear transformation and an offset.
An integer coordinate rectangle.
Definition: Box.h:55
A 2D linear coordinate transformation.
const char * source()
Source function that allows astChannel to source from a Stream.
Definition: Stream.h:224
std::vector< lsst::geom::Box2I > footprintToBBoxList(Footprint const &footprint)
Return a list of BBoxs, whose union contains exactly the pixels in the footprint, neither more nor le...
Definition: Footprint.cc:352
std::shared_ptr< Footprint > mergeFootprints(Footprint const &footprint1, Footprint const &footprint2)
Merges two Footprints – appends their peaks, and unions their spans, returning a new Footprint.
Definition: Footprint.cc:327
Stencil
An enumeration class which describes the shapes.
Definition: SpanSet.h:66
Transform< Point2Endpoint, Point2Endpoint > TransformPoint2ToPoint2
Definition: Transform.h:300
A base class for image defects.
STL namespace.