33#include <nlohmann/json.hpp>
43 mBoundingBox(
c.mBoundingBox ),
44 mHasCachedValidity(
c.mHasCachedValidity ),
45 mValidityFailureReason(
c.mValidityFailureReason )
47 int nGeoms =
c.mGeometries.size();
49 for (
int i = 0; i < nGeoms; ++i )
61 int nGeoms =
c.mGeometries.size();
63 for (
int i = 0; i < nGeoms; ++i )
78 auto result = std::make_unique< QgsGeometryCollection >();
80 return result.release();
97 std::unique_ptr<QgsGeometryCollection> result;
101 std::unique_ptr<QgsAbstractGeometry> gridified { geom->
snappedToGrid( hSpacing, vSpacing, dSpacing, mSpacing, removeRedundantPoints ) };
111 return result.release();
119 if ( geom->removeDuplicateNodes( epsilon, useZValues ) ) result =
true;
143 if (
id.part < 0 || id.part >=
mGeometries.count() )
150 if ( part ==
id.part )
152 int partNumber = geometry->vertexNumberFromVertexId(
QgsVertexId( 0,
id.ring,
id.vertex ) );
153 if ( partNumber == -1 )
155 return number + partNumber;
159 number += geometry->nCoordinates();
173 if ( !mBoundingBox.
isNull() )
185 if ( geometry->boundingBoxIntersects( box3d ) )
214 if ( !geometry->isEmpty() )
246 index = std::min(
static_cast<int>(
mGeometries.count() ), index );
267 QVector< QgsAbstractGeometry * > results =
mGeometries;
277 geometry->normalize();
281 return a->compareTo( b ) > 0;
288 QVector< QgsAbstractGeometry * >::const_iterator it =
mGeometries.constBegin();
291 int dim = ( *it )->dimension();
302 return QStringLiteral(
"GeometryCollection" );
309 g->transform( ct, d, transformZ );
318 g->transform( t, zTranslate, zScale, mTranslate, mScale );
325 QVector< QgsAbstractGeometry * >::const_iterator it =
mGeometries.constBegin();
337 QPainterPath partPath = geom->asQPainterPath();
338 if ( !partPath.isEmpty() )
339 p.addPath( partPath );
358 wkbPtr >> nGeometries;
360 QVector<QgsAbstractGeometry *> geometryListBackup =
mGeometries;
363 for (
int i = 0; i < nGeometries; ++i )
376 qDeleteAll( geometryListBackup );
395 int binarySize =
sizeof( char ) +
sizeof( quint32 ) +
sizeof( quint32 );
400 binarySize += geom->wkbSize( flags );
409 int countNonNull = 0;
422 wkb << static_cast<quint32>(
wkbType() );
423 wkb << static_cast<quint32>( countNonNull );
428 wkb << geom->asWkb( flags );
439 wkt += QLatin1String(
" EMPTY" );
442 wkt += QLatin1String(
" (" );
445 QString childWkt = geom->asWkt(
precision );
448 childWkt = childWkt.mid( childWkt.indexOf(
'(' ) );
450 wkt += childWkt +
',';
452 if ( wkt.endsWith(
',' ) )
463 QDomElement elemMultiGeometry = doc.createElementNS( ns, QStringLiteral(
"MultiGeometry" ) );
466 QDomElement elemGeometryMember = doc.createElementNS( ns, QStringLiteral(
"geometryMember" ) );
467 elemGeometryMember.appendChild( geom->asGml2( doc,
precision, ns, axisOrder ) );
468 elemMultiGeometry.appendChild( elemGeometryMember );
470 return elemMultiGeometry;
475 QDomElement elemMultiGeometry = doc.createElementNS( ns, QStringLiteral(
"MultiGeometry" ) );
478 QDomElement elemGeometryMember = doc.createElementNS( ns, QStringLiteral(
"geometryMember" ) );
479 elemGeometryMember.appendChild( geom->asGml3( doc,
precision, ns, axisOrder ) );
480 elemMultiGeometry.appendChild( elemGeometryMember );
482 return elemMultiGeometry;
487 json coordinates( json::array( ) );
490 coordinates.push_back( geom->asJsonObject(
precision ) );
494 {
"type",
"GeometryCollection" },
495 {
"geometries", coordinates }
502 kml.append( QLatin1String(
"<MultiGeometry>" ) );
503 const QVector< QgsAbstractGeometry * > &geometries =
mGeometries;
506 kml.append( geometry->asKml(
precision ) );
508 kml.append( QLatin1String(
"</MultiGeometry>" ) );
514 if ( mBoundingBox.
isNull() )
543 mHasCachedValidity =
false;
544 mValidityFailureReason.clear();
551 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
552 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
556 QgsCoordinateSequence::const_iterator cIt = geomCoords.constBegin();
557 for ( ; cIt != geomCoords.constEnd(); ++cIt )
559 sequence.push_back( *cIt );
570 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
571 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
573 count += ( *geomIt )->nCoordinates();
622 bool success =
mGeometries.at( position.
part )->insertVertex( position, vertex );
637 bool success =
mGeometries.at( position.
part )->moveVertex( position, newPos );
676 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
677 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
679 length += ( *geomIt )->length();
687 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
688 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
690 area += ( *geomIt )->area();
698 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
699 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
714 qDeleteAll( subtypes );
719 QString secondWithoutParentheses =
parts.second;
720 secondWithoutParentheses = secondWithoutParentheses.remove(
'(' ).remove(
')' ).simplified().remove(
' ' );
721 if ( (
parts.second.compare( QLatin1String(
"EMPTY" ), Qt::CaseInsensitive ) == 0 ) ||
722 secondWithoutParentheses.isEmpty() )
724 qDeleteAll( subtypes );
728 QString defChildWkbType = QStringLiteral(
"%1%2%3 " ).arg( defaultChildWkbType,
is3D() ? QStringLiteral(
"Z" ) : QString(),
isMeasure() ? QStringLiteral(
"M" ) : QString() );
731 for (
const QString &childWkt : blocks )
735 bool success =
false;
751 qDeleteAll( subtypes );
755 qDeleteAll( subtypes );
763 hasZ = hasZ || geom->is3D();
764 hasM = hasM || geom->isMeasure();
778 QVector< QgsAbstractGeometry * >::const_iterator it =
mGeometries.constBegin();
781 if ( ( *it )->hasCurvedSegments() )
793 if ( !geomCollection )
799 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
800 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
802 geomCollection->
addGeometry( ( *geomIt )->segmentize( tolerance, toleranceType ) );
804 return geom.release();
866 if (
id.part < 0 || id.part >=
mGeometries.size() )
882 if ( flags == 0 && mHasCachedValidity )
885 error = mValidityFailureReason;
886 return error.isEmpty();
893 mValidityFailureReason = !res ? error : QString();
894 mHasCachedValidity =
true;
908 geom->addZValue( zValue );
923 geom->addMValue( mValue );
963 geom->filterVertices( filter );
983 std::unique_ptr< QgsGeometryCollection > res;
991 if (
const QgsMultiPoint *mp = qgsgeometry_cast< const QgsMultiPoint *>(
this ) )
995 res = std::make_unique< QgsMultiPoint >();
1003 if (
const QgsMultiLineString *ml = qgsgeometry_cast< const QgsMultiLineString *>(
this ) )
1007 res = std::make_unique< QgsMultiLineString >();
1015 if (
const QgsMultiPolygon *mp = qgsgeometry_cast< const QgsMultiPolygon *>(
this ) )
1019 res = std::make_unique< QgsMultiPolygon>();
1039 if ( thisPartType == filterSinglePartType )
1041 res->addGeometry( part->clone() );
1045 return res.release();
1064 newCollection->addGeometry( geom->toCurveType() );
1066 return newCollection.release();
1079 std::unique_ptr< QgsGeometryCollection > res = std::make_unique< QgsGeometryCollection >();
1083 res->addGeometry(
mGeometries.at( i )->simplifyByDistance( tolerance ) );
1085 return res.release();
1097 res = geom->transform( transformer, feedback );
1129 const QgsGeometryCollection *otherCollection = qgsgeometry_cast<const QgsGeometryCollection *>( other );
1130 if ( !otherCollection )
1139 const int comparison = aGeom->
compareTo( bGeom );
1140 if ( comparison != 0 )
@ AllowSelfTouchingHoles
Indicates that self-touching holes are permitted. OGC validity states that self-touching holes are NO...
QFlags< GeometryValidityFlag > GeometryValidityFlags
Geometry validity flags.
@ RejectOnInvalidSubGeometry
Don't allow geometries with invalid sub-geometries to be created.
WkbType
The WKB type describes the number of dimensions a geometry has.
@ GeometryCollection
GeometryCollection.
TransformDirection
Indicates the direction (forward or inverse) of a transform.
Abstract base class for all geometries.
SegmentationToleranceType
Segmentation tolerance as maximum angle or maximum difference between approximation and circle.
virtual double vertexAngle(QgsVertexId vertex) const =0
Returns approximate angle at a vertex.
virtual const QgsAbstractGeometry * simplifiedTypeRef() const
Returns a reference to the simplest lossless representation of this geometry, e.g.
bool isMeasure() const
Returns true if the geometry contains m values.
QFlags< WkbFlag > WkbFlags
bool is3D() const
Returns true if the geometry is 3D and contains a z-value.
AxisOrder
Axis order for GML generation.
QString wktTypeStr() const
Returns the WKT type string of the geometry.
virtual QgsPoint vertexAt(QgsVertexId id) const =0
Returns the point corresponding to a specified vertex id.
virtual void clearCache() const
Clears any cached parameters associated with the geometry, e.g., bounding boxes.
QgsAbstractGeometry & operator=(const QgsAbstractGeometry &geom)
Qgis::WkbType wkbType() const
Returns the WKB type of the geometry.
virtual bool isEmpty() const
Returns true if the geometry is empty.
virtual bool boundingBoxIntersects(const QgsRectangle &rectangle) const
Returns true if the bounding box of this geometry intersects with a rectangle.
virtual bool deleteVertex(QgsVertexId position)=0
Deletes a vertex within the geometry.
virtual double segmentLength(QgsVertexId startVertex) const =0
Returns the length of the segment of the geometry which begins at startVertex.
virtual int compareTo(const QgsAbstractGeometry *other) const
Comparator for sorting of geometry.
QgsGeometryConstPartIterator parts() const
Returns Java-style iterator for traversal of parts of the geometry.
virtual bool nextVertex(QgsVertexId &id, QgsPoint &vertex) const =0
Returns next vertex id and coordinates.
virtual QgsAbstractGeometry * clone() const =0
Clones the geometry by performing a deep copy.
static endian_t endian()
Returns whether this machine uses big or little endian.
A 3-dimensional box composed of x, y, z coordinates.
bool intersects(const QgsBox3D &other) const
Returns true if box intersects with another box.
void combineWith(const QgsBox3D &box)
Expands the bbox so that it covers both the original rectangle and the given rectangle.
bool isNull() const
Test if the box is null (holding no spatial information).
Circular string geometry type.
Compound curve geometry type.
Qgis::WkbType readHeader() const
readHeader
Curve polygon geometry type.
Base class for feedback objects to be used for cancellation of something running in a worker thread.
bool isCanceled() const
Tells whether the operation has been canceled already.
QgsGeometryCollection * toCurveType() const override
Returns the geometry converted to the more generic curve type.
bool dropMValue() override
Drops any measure values which exist in the geometry.
QString asKml(int precision=17) const override
Returns a KML representation of the geometry.
double perimeter() const override
Returns the planar, 2-dimensional perimeter of the geometry.
QVector< QgsAbstractGeometry * > mGeometries
void draw(QPainter &p) const override
Draws the geometry using the specified QPainter.
json asJsonObject(int precision=17) const override
Returns a json object representation of the geometry.
bool fromWkt(const QString &wkt) override
Sets the geometry from a WKT string.
void transform(const QgsCoordinateTransform &ct, Qgis::TransformDirection d=Qgis::TransformDirection::Forward, bool transformZ=false) override
Transforms the geometry using a coordinate transform.
bool addMValue(double mValue=0) override
Adds a measure to the geometry, initialized to a preset value.
void clear() override
Clears the geometry, ie reset it to a null geometry.
void reserve(int size)
Attempts to allocate memory for at least size geometries.
QPainterPath asQPainterPath() const override
Returns the geometry represented as a QPainterPath.
QgsAbstractGeometry * segmentize(double tolerance=M_PI_2/90, SegmentationToleranceType toleranceType=MaximumAngle) const override
Returns a geometry without curves.
void clearCache() const override
Clears any cached parameters associated with the geometry, e.g., bounding boxes.
int vertexCount(int part=0, int ring=0) const override
Returns the number of vertices of which this geometry is built.
QString geometryType() const override
Returns a unique string representing the geometry type.
QgsGeometryCollection & operator=(const QgsGeometryCollection &c)
bool dropZValue() override
Drops any z-dimensions which exist in the geometry.
void adjacentVertices(QgsVertexId vertex, QgsVertexId &previousVertex, QgsVertexId &nextVertex) const override
Returns the vertices adjacent to a specified vertex within a geometry.
int wkbSize(QgsAbstractGeometry::WkbFlags flags=QgsAbstractGeometry::WkbFlags()) const override
Returns the length of the QByteArray returned by asWkb()
virtual bool insertGeometry(QgsAbstractGeometry *g, int index)
Inserts a geometry before a specified index and takes ownership.
int vertexNumberFromVertexId(QgsVertexId id) const override
Returns the vertex number corresponding to a vertex id.
virtual QgsGeometryCollection * simplifyByDistance(double tolerance) const override
Simplifies the geometry by applying the Douglas Peucker simplification by distance algorithm.
bool isEmpty() const override
Returns true if the geometry is empty.
QgsGeometryCollection()
Constructor for an empty geometry collection.
bool nextVertex(QgsVertexId &id, QgsPoint &vertex) const override
Returns next vertex id and coordinates.
QString asWkt(int precision=17) const override
Returns a WKT representation of the geometry.
bool fromWkb(QgsConstWkbPtr &wkb) override
Sets the geometry from a WKB string.
QgsAbstractGeometry * boundary() const override
Returns the closure of the combinatorial boundary of the geometry (ie the topological boundary of the...
bool isValid(QString &error, Qgis::GeometryValidityFlags flags=Qgis::GeometryValidityFlags()) const override
Checks validity of the geometry, and returns true if the geometry is valid.
double area() const override
Returns the planar, 2-dimensional area of the geometry.
bool removeDuplicateNodes(double epsilon=4 *std::numeric_limits< double >::epsilon(), bool useZValues=false) override
Removes duplicate nodes from the geometry, wherever removing the nodes does not result in a degenerat...
int ringCount(int part=0) const override
Returns the number of rings of which this geometry is built.
void normalize() final
Reorganizes the geometry into a normalized form (or "canonical" form).
QgsAbstractGeometry * childGeometry(int index) const override
Returns pointer to child geometry (for geometries with child geometries - i.e.
void transformVertices(const std::function< QgsPoint(const QgsPoint &) > &transform) override
Transforms the vertices from the geometry in place, applying the transform function to every vertex.
void filterVertices(const std::function< bool(const QgsPoint &) > &filter) override
Filters the vertices from the geometry in place, removing any which do not return true for the filter...
QDomElement asGml2(QDomDocument &doc, int precision=17, const QString &ns="gml", QgsAbstractGeometry::AxisOrder axisOrder=QgsAbstractGeometry::AxisOrder::XY) const override
Returns a GML2 representation of the geometry.
bool moveVertex(QgsVertexId position, const QgsPoint &newPos) override
Moves a vertex within the geometry.
int nCoordinates() const override
Returns the number of nodes contained in the geometry.
QgsGeometryCollection * extractPartsByType(Qgis::WkbType type, bool useFlatType=true) const
Returns a new QgsGeometryCollection subclass which consists of the parts of this collection which mat...
QgsBox3D boundingBox3D() const override
Returns the 3D bounding box for the geometry.
QVector< QgsAbstractGeometry * > takeGeometries()
Removes all geometries from the collection, returning them and their ownership to the caller.
QgsGeometryCollection * clone() const override
Clones the geometry by performing a deep copy.
virtual bool addGeometries(const QVector< QgsAbstractGeometry * > &geometries)
Adds a list of geometries to the collection, transferring ownership to the collection.
int childCount() const override
Returns number of child geometries (for geometries with child geometries) or child points (for geomet...
bool hasCurvedSegments() const override
Returns true if the geometry contains curved segments.
bool insertVertex(QgsVertexId position, const QgsPoint &vertex) override
Inserts a vertex into the geometry.
double length() const override
Returns the planar, 2-dimensional length of the geometry.
int dimension() const override
Returns the inherent dimension of the geometry.
bool deleteVertex(QgsVertexId position) override
Deletes a vertex within the geometry.
QgsPoint vertexAt(QgsVertexId id) const override
Returns the point corresponding to a specified vertex id.
virtual bool removeGeometry(int nr)
Removes a geometry from the collection.
void swapXy() override
Swaps the x and y coordinates from the geometry.
bool boundingBoxIntersects(const QgsBox3D &box3d) const override
Returns true if the bounding box of this geometry intersects with a box3d.
QgsGeometryCollection * createEmptyWithSameType() const override
Creates a new geometry with the same class and same WKB type as the original and transfers ownership.
virtual bool addGeometry(QgsAbstractGeometry *g)
Adds a geometry and takes ownership. Returns true in case of success.
QgsBox3D calculateBoundingBox3D() const override
Calculates the minimal 3D bounding box for the geometry.
bool addZValue(double zValue=0) override
Adds a z-dimension to the geometry, initialized to a preset value.
QByteArray asWkb(QgsAbstractGeometry::WkbFlags flags=QgsAbstractGeometry::WkbFlags()) const override
Returns a WKB representation of the geometry.
int partCount() const override
Returns count of parts contained in the geometry.
~QgsGeometryCollection() override
virtual bool wktOmitChildType() const
Returns whether child type names are omitted from Wkt representations of the collection.
bool fromCollectionWkt(const QString &wkt, const QVector< QgsAbstractGeometry * > &subtypes, const QString &defaultChildWkbType=QString())
Reads a collection from a WKT string.
QgsCoordinateSequence coordinateSequence() const override
Retrieves the sequence of geometries, rings and nodes.
int numGeometries() const
Returns the number of geometries within the collection.
QgsGeometryCollection * snappedToGrid(double hSpacing, double vSpacing, double dSpacing=0, double mSpacing=0, bool removeRedundantPoints=false) const override
Makes a new geometry with all the points or vertices snapped to the closest point of the grid.
double segmentLength(QgsVertexId startVertex) const override
Returns the length of the segment of the geometry which begins at startVertex.
QDomElement asGml3(QDomDocument &doc, int precision=17, const QString &ns="gml", QgsAbstractGeometry::AxisOrder axisOrder=QgsAbstractGeometry::AxisOrder::XY) const override
Returns a GML3 representation of the geometry.
const QgsAbstractGeometry * geometryN(int n) const
Returns a const reference to a geometry from within the collection.
double vertexAngle(QgsVertexId vertex) const override
Returns approximate angle at a vertex.
int compareToSameClass(const QgsAbstractGeometry *other) const final
Compares to an other geometry of the same class, and returns a integer for sorting of the two geometr...
double closestSegment(const QgsPoint &pt, QgsPoint &segmentPt, QgsVertexId &vertexAfter, int *leftOf=nullptr, double epsilon=4 *std::numeric_limits< double >::epsilon()) const override
Searches for the closest segment of the geometry to a given point.
const QgsAbstractGeometry * simplifiedTypeRef() const override
Returns a reference to the simplest lossless representation of this geometry, e.g.
static std::unique_ptr< QgsAbstractGeometry > geomFromWkb(QgsConstWkbPtr &wkb)
Construct geometry from a WKB string.
static std::unique_ptr< QgsAbstractGeometry > geomFromWkbType(Qgis::WkbType t)
Returns empty geometry from wkb type.
static QStringList wktGetChildBlocks(const QString &wkt, const QString &defaultType=QString())
Parses a WKT string and returns of list of blocks contained in the WKT.
static QPair< Qgis::WkbType, QString > wktReadBlock(const QString &wkt)
Parses a WKT block of the format "TYPE( contents )" and returns a pair of geometry type to contents (...
static double closestSegmentFromComponents(T &container, ComponentType ctype, const QgsPoint &pt, QgsPoint &segmentPt, QgsVertexId &vertexAfter, int *leftOf, double epsilon)
Does vector analysis using the geos library and handles import, export, exception handling*.
Line string geometry type, with support for z-dimension and m-values.
Multi curve geometry collection.
Multi line string geometry collection.
Multi point geometry collection.
Multi polygon geometry collection.
Multi surface geometry collection.
Point geometry type, with support for z-dimension and m-values.
static Qgis::WkbType dropM(Qgis::WkbType type)
Drops the m dimension (if present) for a WKB type and returns the new type.
static Qgis::GeometryType geometryType(Qgis::WkbType type)
Returns the geometry type for a WKB type, e.g., both MultiPolygon and CurvePolygon would have a Polyg...
static Qgis::WkbType linearType(Qgis::WkbType type)
Returns the linear type for a WKB type.
static Qgis::WkbType dropZ(Qgis::WkbType type)
Drops the z dimension (if present) for a WKB type and returns the new type.
static Qgis::WkbType addM(Qgis::WkbType type)
Adds the m dimension to a WKB type and returns the new type.
static Qgis::WkbType addZ(Qgis::WkbType type)
Adds the z dimension to a WKB type and returns the new type.
static Qgis::WkbType singleType(Qgis::WkbType type)
Returns the single type for a WKB type.
static bool hasZ(Qgis::WkbType type)
Tests whether a WKB type contains the z-dimension.
static bool hasM(Qgis::WkbType type)
Tests whether a WKB type contains m values.
static Qgis::WkbType flatType(Qgis::WkbType type)
Returns the flat type for a WKB type.
Contains geos related utilities and functions.
As part of the API refactoring and improvements which landed in the Processing API was substantially reworked from the x version This was done in order to allow much of the underlying Processing framework to be ported into c
QVector< QgsRingSequence > QgsCoordinateSequence
Utility class for identifying a unique vertex within a geometry.