LSST Applications  21.0.0-147-g0e635eb1+1acddb5be5,22.0.0+052faf71bd,22.0.0+1ea9a8b2b2,22.0.0+6312710a6c,22.0.0+729191ecac,22.0.0+7589c3a021,22.0.0+9f079a9461,22.0.1-1-g7d6de66+b8044ec9de,22.0.1-1-g87000a6+536b1ee016,22.0.1-1-g8e32f31+6312710a6c,22.0.1-10-gd060f87+016f7cdc03,22.0.1-12-g9c3108e+df145f6f68,22.0.1-16-g314fa6d+c825727ab8,22.0.1-19-g93a5c75+d23f2fb6d8,22.0.1-19-gb93eaa13+aab3ef7709,22.0.1-2-g8ef0a89+b8044ec9de,22.0.1-2-g92698f7+9f079a9461,22.0.1-2-ga9b0f51+052faf71bd,22.0.1-2-gac51dbf+052faf71bd,22.0.1-2-gb66926d+6312710a6c,22.0.1-2-gcb770ba+09e3807989,22.0.1-20-g32debb5+b8044ec9de,22.0.1-23-gc2439a9a+fb0756638e,22.0.1-3-g496fd5d+09117f784f,22.0.1-3-g59f966b+1e6ba2c031,22.0.1-3-g849a1b8+f8b568069f,22.0.1-3-gaaec9c0+c5c846a8b1,22.0.1-32-g5ddfab5d3+60ce4897b0,22.0.1-4-g037fbe1+64e601228d,22.0.1-4-g8623105+b8044ec9de,22.0.1-5-g096abc9+d18c45d440,22.0.1-5-g15c806e+57f5c03693,22.0.1-7-gba73697+57f5c03693,master-g6e05de7fdc+c1283a92b8,master-g72cdda8301+729191ecac,w.2021.39
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
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:51
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.