lsst::sphgeom::ConvexPolygon Class Reference

ConvexPolygon is a closed convex polygon on the unit sphere. More...

#include <ConvexPolygon.h>

Inheritance diagram for lsst::sphgeom::ConvexPolygon:

Public Member Functions
	ConvexPolygon (std::vector< UnitVector3d > const &points)
	This constructor creates a convex polygon that is the convex hull of the given set of points.
	ConvexPolygon (UnitVector3d const &v0, UnitVector3d const &v1, UnitVector3d const &v2)
	This constructor creates a triangle with the given vertices.
	ConvexPolygon (UnitVector3d const &v0, UnitVector3d const &v1, UnitVector3d const &v2, UnitVector3d const &v3)
	This constructor creates a quadrilateral with the given vertices.
bool	operator== (ConvexPolygon const &p) const
	Two convex polygons are equal iff they contain the same points.
bool	operator!= (ConvexPolygon const &p) const
std::vector< UnitVector3d > const &	getVertices () const
UnitVector3d	getCentroid () const
	The centroid of a polygon is its center of mass projected onto S², assuming a uniform mass distribution over the polygon surface.
std::unique_ptr< Region >	clone () const override
	clone returns a deep copy of this region.
Box	getBoundingBox () const override
	getBoundingBox returns a bounding-box for this region.
Box3d	getBoundingBox3d () const override
	getBoundingBox3d returns a 3-dimensional bounding-box for this region.
Circle	getBoundingCircle () const override
	getBoundingCircle returns a bounding-circle for this region.
Relationship	relate (Region const &r) const override
Relationship	relate (Box const &) const override
Relationship	relate (Circle const &) const override
Relationship	relate (ConvexPolygon const &) const override
Relationship	relate (Ellipse const &) const override
std::vector< std::uint8_t >	encode () const override
	encode serializes this region into an opaque byte string.
virtual bool	contains (UnitVector3d const &) const=0
	contains tests whether the given unit vector is inside this region.
bool	contains (double x, double y, double z) const
	contains tests whether the unit vector defined by the given (not necessarily normalized) coordinates is inside this region.
bool	contains (double lon, double lat) const
	contains tests whether the unit vector defined by the given longitude and latitude coordinates (in radians) is inside this region.
bool	contains (UnitVector3d const &v) const override
bool	contains (Region const &r) const
bool	isDisjointFrom (Region const &r) const
bool	intersects (Region const &r) const
bool	isWithin (Region const &r) const

Static Public Member Functions
static ConvexPolygon	convexHull (std::vector< UnitVector3d > const &points)
	convexHull returns the convex hull of the given set of points if it exists and throws an exception otherwise.
static std::vector< std::unique_ptr< Region > >	getRegions (Region const &region)
	getRegions returns a vector of Region.
static std::unique_ptr< ConvexPolygon >	decode (std::vector< std::uint8_t > const &s)
static std::unique_ptr< ConvexPolygon >	decode (std::uint8_t const *buffer, size_t n)

Static Public Attributes
static constexpr std::uint8_t	TYPE_CODE = 'p'

Detailed Description

ConvexPolygon is a closed convex polygon on the unit sphere.

Its edges are great circles (geodesics), and the shorter of the two great circle segments between any two points on the polygon boundary is contained in the polygon.

The vertices of a convex polygon are distinct and have counter-clockwise orientation when viewed from outside the unit sphere. No three consecutive vertices are coplanar and edges do not intersect except at vertices.

Furthermore, if a convex polygon contains a point p of S², then we require that it be disjoint from point -p. This guarantees the existence of a unique shortest great circle segment between any 2 points contained in the polygon, but means e.g. that hemispheres and lunes cannot be represented by convex polygons.

Currently, the only way to construct a convex polygon is to compute the convex hull of a point set.

Definition at line 65 of file ConvexPolygon.h.

Constructor & Destructor Documentation

ConvexPolygon() [1/3]

lsst::sphgeom::ConvexPolygon::ConvexPolygon ( std::vector< UnitVector3d > const & points )

explicit

This constructor creates a convex polygon that is the convex hull of the given set of points.

Definition at line 262 of file ConvexPolygon.cc.

                                                                   :
    _vertices(points)
{
    computeHull(_vertices);
}

ConvexPolygon() [2/3]

lsst::sphgeom::ConvexPolygon::ConvexPolygon ( UnitVector3d const & v0, UnitVector3d const & v1, UnitVector3d const & v2 )

inline

This constructor creates a triangle with the given vertices.

It is assumed that orientation(v0, v1, v2) = 1. Use with caution - for performance reasons, this is not verified!

Definition at line 85 of file ConvexPolygon.h.

                                           :
        _vertices{v0, v1, v2}
    {}

ConvexPolygon() [3/3]

lsst::sphgeom::ConvexPolygon::ConvexPolygon ( UnitVector3d const & v0, UnitVector3d const & v1, UnitVector3d const & v2, UnitVector3d const & v3 )

inline

This constructor creates a quadrilateral with the given vertices.

It is assumed that orientation(v0, v1, v2), orientation(v1, v2, v3), orientation(v2, v3, v0), and orientation (v3, v0, v1) are all 1. Use with caution - for performance reasons, this is not verified!

Definition at line 96 of file ConvexPolygon.h.

                                           :
        _vertices{v0, v1, v2, v3}
    {}

Member Function Documentation

clone()

std::unique_ptr< Region > lsst::sphgeom::ConvexPolygon::clone ( ) const

inlineoverridevirtual

clone returns a deep copy of this region.

Implements lsst::sphgeom::Region.

Definition at line 116 of file ConvexPolygon.h.

                                                 {
        return std::unique_ptr<ConvexPolygon>(new ConvexPolygon(*this));
    }

contains() [1/5]

bool lsst::sphgeom::Region::contains	(	double	lon,
		double	lat ) const

contains tests whether the unit vector defined by the given longitude and latitude coordinates (in radians) is inside this region.

Definition at line 112 of file Region.cc.

                                                  {
    return contains(UnitVector3d(LonLat::fromRadians(lon, lat)));
}

contains() [2/5]

bool lsst::sphgeom::Region::contains	(	double	x,
		double	y,
		double	z ) const

contains tests whether the unit vector defined by the given (not necessarily normalized) coordinates is inside this region.

Definition at line 108 of file Region.cc.

                                                        {
    return contains(UnitVector3d(x, y, z));
}

contains() [3/5]

bool lsst::sphgeom::ConvexPolygon::contains

(

Region const &

r

)

const

contains returns true if the intersection of this convex polygon and x is equal to x.

Definition at line 323 of file ConvexPolygon.cc.

                                                   {
    return (relate(r) & CONTAINS) != 0;
}

contains() [4/5]

virtual bool lsst::sphgeom::Region::contains ( UnitVector3d const & ) const

virtual

contains tests whether the given unit vector is inside this region.

Implements lsst::sphgeom::Region.

contains() [5/5]

bool lsst::sphgeom::ConvexPolygon::contains ( UnitVector3d const & v ) const

overridevirtual

contains returns true if the intersection of this convex polygon and x is equal to x.

Implements lsst::sphgeom::Region.

Definition at line 319 of file ConvexPolygon.cc.

                                                         {
    return detail::contains(_vertices.begin(), _vertices.end(), v);
}

convexHull()

static ConvexPolygon lsst::sphgeom::ConvexPolygon::convexHull ( std::vector< UnitVector3d > const & points )

inlinestatic

convexHull returns the convex hull of the given set of points if it exists and throws an exception otherwise.

Though points are supplied in a vector, they really are conceptually a set - the ConvexPolygon returned is invariant under permutation of the input array.

Definition at line 73 of file ConvexPolygon.h.

                                                                            {
        return ConvexPolygon(points);
    }

decode() [1/2]

std::unique_ptr< ConvexPolygon > lsst::sphgeom::ConvexPolygon::decode ( std::uint8_t const * buffer, size_t n )

static

decode deserializes a ConvexPolygon from a byte string produced by encode.

Definition at line 368 of file ConvexPolygon.cc.

{
    if (buffer == nullptr || *buffer != TYPE_CODE ||
        n < 1 + 24*3 || (n - 1) % 24 != 0) {
        throw std::runtime_error("Byte-string is not an encoded ConvexPolygon");
    }
    std::unique_ptr<ConvexPolygon> poly(new ConvexPolygon);
    ++buffer;
    size_t nv = (n - 1) / 24;
    poly->_vertices.reserve(nv);
    for (size_t i = 0; i < nv; ++i, buffer += 24) {
        poly->_vertices.push_back(UnitVector3d::fromNormalized(
            decodeDouble(buffer),
            decodeDouble(buffer + 8),
            decodeDouble(buffer + 16)
        ));
    }
    return poly;
}

decode() [2/2]

static std::unique_ptr< ConvexPolygon > lsst::sphgeom::ConvexPolygon::decode ( std::vector< std::uint8_t > const & s )

inlinestatic

decode deserializes a ConvexPolygon from a byte string produced by encode.

Definition at line 165 of file ConvexPolygon.h.

                                                                                {
        return decode(s.data(), s.size());
    }

encode()

std::vector< std::uint8_t > lsst::sphgeom::ConvexPolygon::encode ( ) const

overridevirtual

encode serializes this region into an opaque byte string.

Byte strings emitted by encode can be deserialized with decode.

Implements lsst::sphgeom::Region.

Definition at line 355 of file ConvexPolygon.cc.

                                                  {
    std::vector<std::uint8_t> buffer;
    std::uint8_t tc = TYPE_CODE;
    buffer.reserve(1 + 24 * _vertices.size());
    buffer.push_back(tc);
    for (UnitVector3d const & v: _vertices) {
        encodeDouble(v.x(), buffer);
        encodeDouble(v.y(), buffer);
        encodeDouble(v.z(), buffer);
    }
    return buffer;
}

getBoundingBox()

Box lsst::sphgeom::ConvexPolygon::getBoundingBox ( ) const

overridevirtual

getBoundingBox returns a bounding-box for this region.

Implements lsst::sphgeom::Region.

Definition at line 311 of file ConvexPolygon.cc.

                                        {
    return detail::boundingBox(_vertices.begin(), _vertices.end());
}

getBoundingBox3d()

Box3d lsst::sphgeom::ConvexPolygon::getBoundingBox3d ( ) const

overridevirtual

getBoundingBox3d returns a 3-dimensional bounding-box for this region.

Implements lsst::sphgeom::Region.

Definition at line 315 of file ConvexPolygon.cc.

                                            {
    return detail::boundingBox3d(_vertices.begin(), _vertices.end());
}

getBoundingCircle()

Circle lsst::sphgeom::ConvexPolygon::getBoundingCircle ( ) const

overridevirtual

getBoundingCircle returns a bounding-circle for this region.

Implements lsst::sphgeom::Region.

Definition at line 307 of file ConvexPolygon.cc.

                                              {
    return detail::boundingCircle(_vertices.begin(), _vertices.end());
}

getCentroid()

UnitVector3d lsst::sphgeom::ConvexPolygon::getCentroid

(

)

const

The centroid of a polygon is its center of mass projected onto S², assuming a uniform mass distribution over the polygon surface.

Definition at line 303 of file ConvexPolygon.cc.

                                              {
    return detail::centroid(_vertices.begin(), _vertices.end());
}

getRegions()

std::vector< std::unique_ptr< Region > > lsst::sphgeom::Region::getRegions ( Region const & region )

staticinherited

getRegions returns a vector of Region.

Definition at line 76 of file Region.cc.

                                                                        {
    std::vector<std::unique_ptr<Region>> result;
    if (auto union_region = dynamic_cast<UnionRegion const *>(&region)) {
        for(int i = 0; i < 2; ++i) {
            result.emplace_back(union_region->getOperand(i).clone());
        }
    } else if(auto intersection_region = dynamic_cast<IntersectionRegion const *>(&region)) {
        for(int i = 0; i < 2; ++i) {
            intersection_region->getOperand(i);
            result.emplace_back(intersection_region->getOperand(i).clone());
        }
    } else {
        result.emplace_back(region.clone());
    }
    return result;
}

getVertices()

std::vector< UnitVector3d > const & lsst::sphgeom::ConvexPolygon::getVertices ( ) const

inline

Definition at line 107 of file ConvexPolygon.h.

                                                        {
        return _vertices;
    }

intersects()

bool lsst::sphgeom::ConvexPolygon::intersects

(

Region const &

r

)

const

intersects returns true if the intersection of this convex polygon and x is non-empty.

Definition at line 331 of file ConvexPolygon.cc.

                                                     {
    return !isDisjointFrom(r);
}

isDisjointFrom()

bool lsst::sphgeom::ConvexPolygon::isDisjointFrom

(

Region const &

r

)

const

isDisjointFrom returns true if the intersection of this convex polygon and x is empty.

Definition at line 327 of file ConvexPolygon.cc.

                                                         {
    return (relate(r) & DISJOINT) != 0;
}

isWithin()

bool lsst::sphgeom::ConvexPolygon::isWithin

(

Region const &

r

)

const

isWithin returns true if the intersection of this convex polygon and x is this convex polygon.

Definition at line 335 of file ConvexPolygon.cc.

                                                   {
    return (relate(r) & WITHIN) != 0;
}

operator!=()

bool lsst::sphgeom::ConvexPolygon::operator!= ( ConvexPolygon const & p ) const

inline

Definition at line 105 of file ConvexPolygon.h.

{ return !(*this == p); }

operator==()

bool lsst::sphgeom::ConvexPolygon::operator==

(

ConvexPolygon const &

p

)

const

Two convex polygons are equal iff they contain the same points.

Definition at line 268 of file ConvexPolygon.cc.

                                                            {
    if (this == &p) {
        return true;
    }
    if (_vertices.size() != p._vertices.size()) {
        return false;
    }
    VertexIterator i = _vertices.begin();
    VertexIterator f = p._vertices.begin();
    VertexIterator const ep = p._vertices.end();
    // Find vertex f of p equal to the first vertex of this polygon.
    for (; f != ep; ++f) {
        if (*i == *f) {
            break;
        }
    }
    if (f == ep) {
        // No vertex of p is equal to the first vertex of this polygon.
        return false;
    }
    // Now, compare all vertices.
    ++i;
    for (VertexIterator j = f + 1; j != ep; ++i, ++j) {
        if (*i != *j) {
            return false;
        }
    }
    for (VertexIterator j = p._vertices.begin(); j != f; ++i, ++j) {
        if (*i != *j) {
            return false;
        }
    }
    return true;
}

relate() [1/5]

Relationship lsst::sphgeom::ConvexPolygon::relate ( Box const & ) const

overridevirtual

relate computes the spatial relationships between this region A and another region B. The return value S is a bitset with the following properties:

• Bit S & DISJOINT is set only if A and B do not have any points in common.
• Bit S & CONTAINS is set only if A contains all points in B.
• Bit S & WITHIN is set only if B contains all points in A.

Said another way: if the CONTAINS, WITHIN or DISJOINT bit is set, then the corresponding spatial relationship between the two regions holds conclusively. If it is not set, the relationship may or may not hold.

These semantics allow for conservative relationship computations. In particular, a Region may choose to implement relate by replacing itself and/or the argument with a simplified bounding region.

Implements lsst::sphgeom::Region.

Definition at line 339 of file ConvexPolygon.cc.

                                                      {
    return detail::relate(_vertices.begin(), _vertices.end(), b);
}

relate() [2/5]

Relationship lsst::sphgeom::ConvexPolygon::relate ( Circle const & ) const

overridevirtual

relate computes the spatial relationships between this region A and another region B. The return value S is a bitset with the following properties:

• Bit S & DISJOINT is set only if A and B do not have any points in common.
• Bit S & CONTAINS is set only if A contains all points in B.
• Bit S & WITHIN is set only if B contains all points in A.

Said another way: if the CONTAINS, WITHIN or DISJOINT bit is set, then the corresponding spatial relationship between the two regions holds conclusively. If it is not set, the relationship may or may not hold.

These semantics allow for conservative relationship computations. In particular, a Region may choose to implement relate by replacing itself and/or the argument with a simplified bounding region.

Implements lsst::sphgeom::Region.

Definition at line 343 of file ConvexPolygon.cc.

                                                         {
    return detail::relate(_vertices.begin(), _vertices.end(), c);
}

relate() [3/5]

Relationship lsst::sphgeom::ConvexPolygon::relate ( ConvexPolygon const & ) const

overridevirtual

relate computes the spatial relationships between this region A and another region B. The return value S is a bitset with the following properties:

• Bit S & DISJOINT is set only if A and B do not have any points in common.
• Bit S & CONTAINS is set only if A contains all points in B.
• Bit S & WITHIN is set only if B contains all points in A.

Said another way: if the CONTAINS, WITHIN or DISJOINT bit is set, then the corresponding spatial relationship between the two regions holds conclusively. If it is not set, the relationship may or may not hold.

These semantics allow for conservative relationship computations. In particular, a Region may choose to implement relate by replacing itself and/or the argument with a simplified bounding region.

Implements lsst::sphgeom::Region.

Definition at line 347 of file ConvexPolygon.cc.

                                                                {
    return detail::relate(_vertices.begin(), _vertices.end(), p);
}

relate() [4/5]

Relationship lsst::sphgeom::ConvexPolygon::relate ( Ellipse const & ) const

overridevirtual

relate computes the spatial relationships between this region A and another region B. The return value S is a bitset with the following properties:

• Bit S & DISJOINT is set only if A and B do not have any points in common.
• Bit S & CONTAINS is set only if A contains all points in B.
• Bit S & WITHIN is set only if B contains all points in A.

Said another way: if the CONTAINS, WITHIN or DISJOINT bit is set, then the corresponding spatial relationship between the two regions holds conclusively. If it is not set, the relationship may or may not hold.

These semantics allow for conservative relationship computations. In particular, a Region may choose to implement relate by replacing itself and/or the argument with a simplified bounding region.

Implements lsst::sphgeom::Region.

Definition at line 351 of file ConvexPolygon.cc.

                                                          {
    return detail::relate(_vertices.begin(), _vertices.end(), e);
}

relate() [5/5]

Relationship lsst::sphgeom::ConvexPolygon::relate ( Region const & ) const

inlineoverridevirtual

relate computes the spatial relationships between this region A and another region B. The return value S is a bitset with the following properties:

• Bit S & DISJOINT is set only if A and B do not have any points in common.
• Bit S & CONTAINS is set only if A contains all points in B.
• Bit S & WITHIN is set only if B contains all points in A.

Said another way: if the CONTAINS, WITHIN or DISJOINT bit is set, then the corresponding spatial relationship between the two regions holds conclusively. If it is not set, the relationship may or may not hold.

These semantics allow for conservative relationship computations. In particular, a Region may choose to implement relate by replacing itself and/or the argument with a simplified bounding region.

Implements lsst::sphgeom::Region.

Definition at line 151 of file ConvexPolygon.h.

                                                         {
        // Dispatch on the type of r.
        return invert(r.relate(*this));
    }

Member Data Documentation

TYPE_CODE

constexpr std::uint8_t lsst::sphgeom::ConvexPolygon::TYPE_CODE = 'p'

staticconstexpr

Definition at line 67 of file ConvexPolygon.h.

---

The documentation for this class was generated from the following files:

• /j/snowflake/release/lsstsw/stack/lsst-scipipe-9.0.0/Linux64/sphgeom/g781aacb6e4+a8ce1bb630/include/lsst/sphgeom/ConvexPolygon.h
• /j/snowflake/release/lsstsw/stack/lsst-scipipe-9.0.0/Linux64/sphgeom/g781aacb6e4+a8ce1bb630/src/ConvexPolygon.cc