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 ) };
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() )
239 index = std::min(
static_cast<int>(
mGeometries.count() ), index );
262 geometry->normalize();
266 return a->compareTo( b ) > 0;
273 QVector< QgsAbstractGeometry * >::const_iterator it =
mGeometries.constBegin();
276 int dim = ( *it )->dimension();
287 return QStringLiteral(
"GeometryCollection" );
294 g->transform( ct, d, transformZ );
303 g->transform( t, zTranslate, zScale, mTranslate, mScale );
310 QVector< QgsAbstractGeometry * >::const_iterator it =
mGeometries.constBegin();
322 QPainterPath partPath = geom->asQPainterPath();
323 if ( !partPath.isEmpty() )
324 p.addPath( partPath );
343 wkbPtr >> nGeometries;
345 QVector<QgsAbstractGeometry *> geometryListBackup =
mGeometries;
348 for (
int i = 0; i < nGeometries; ++i )
361 qDeleteAll( geometryListBackup );
380 int binarySize =
sizeof( char ) +
sizeof( quint32 ) +
sizeof( quint32 );
385 binarySize += geom->wkbSize( flags );
394 int countNonNull = 0;
407 wkb << static_cast<quint32>(
wkbType() );
408 wkb << static_cast<quint32>( countNonNull );
413 wkb << geom->asWkb( flags );
424 wkt += QLatin1String(
" EMPTY" );
427 wkt += QLatin1String(
" (" );
430 QString childWkt = geom->asWkt(
precision );
433 childWkt = childWkt.mid( childWkt.indexOf(
'(' ) );
435 wkt += childWkt +
',';
437 if ( wkt.endsWith(
',' ) )
448 QDomElement elemMultiGeometry = doc.createElementNS( ns, QStringLiteral(
"MultiGeometry" ) );
451 QDomElement elemGeometryMember = doc.createElementNS( ns, QStringLiteral(
"geometryMember" ) );
452 elemGeometryMember.appendChild( geom->asGml2( doc,
precision, ns, axisOrder ) );
453 elemMultiGeometry.appendChild( elemGeometryMember );
455 return elemMultiGeometry;
460 QDomElement elemMultiGeometry = doc.createElementNS( ns, QStringLiteral(
"MultiGeometry" ) );
463 QDomElement elemGeometryMember = doc.createElementNS( ns, QStringLiteral(
"geometryMember" ) );
464 elemGeometryMember.appendChild( geom->asGml3( doc,
precision, ns, axisOrder ) );
465 elemMultiGeometry.appendChild( elemGeometryMember );
467 return elemMultiGeometry;
472 json coordinates( json::array( ) );
475 coordinates.push_back( geom->asJsonObject(
precision ) );
479 {
"type",
"GeometryCollection" },
480 {
"geometries", coordinates }
487 kml.append( QLatin1String(
"<MultiGeometry>" ) );
488 const QVector< QgsAbstractGeometry * > &geometries =
mGeometries;
491 kml.append( geometry->asKml(
precision ) );
493 kml.append( QLatin1String(
"</MultiGeometry>" ) );
499 if ( mBoundingBox.
isNull() )
528 mHasCachedValidity =
false;
529 mValidityFailureReason.clear();
536 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
537 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
541 QgsCoordinateSequence::const_iterator cIt = geomCoords.constBegin();
542 for ( ; cIt != geomCoords.constEnd(); ++cIt )
544 sequence.push_back( *cIt );
555 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
556 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
558 count += ( *geomIt )->nCoordinates();
607 bool success =
mGeometries.at( position.
part )->insertVertex( position, vertex );
622 bool success =
mGeometries.at( position.
part )->moveVertex( position, newPos );
661 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
662 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
664 length += ( *geomIt )->length();
672 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
673 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
675 area += ( *geomIt )->area();
683 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
684 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
699 qDeleteAll( subtypes );
704 QString secondWithoutParentheses =
parts.second;
705 secondWithoutParentheses = secondWithoutParentheses.remove(
'(' ).remove(
')' ).simplified().remove(
' ' );
706 if ( (
parts.second.compare( QLatin1String(
"EMPTY" ), Qt::CaseInsensitive ) == 0 ) ||
707 secondWithoutParentheses.isEmpty() )
709 qDeleteAll( subtypes );
713 QString defChildWkbType = QStringLiteral(
"%1%2%3 " ).arg( defaultChildWkbType,
is3D() ? QStringLiteral(
"Z" ) : QString(),
isMeasure() ? QStringLiteral(
"M" ) : QString() );
716 for (
const QString &childWkt : blocks )
720 bool success =
false;
736 qDeleteAll( subtypes );
740 qDeleteAll( subtypes );
748 hasZ = hasZ || geom->is3D();
749 hasM = hasM || geom->isMeasure();
763 QVector< QgsAbstractGeometry * >::const_iterator it =
mGeometries.constBegin();
766 if ( ( *it )->hasCurvedSegments() )
778 if ( !geomCollection )
784 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
785 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
787 geomCollection->
addGeometry( ( *geomIt )->segmentize( tolerance, toleranceType ) );
789 return geom.release();
851 if (
id.part < 0 || id.part >=
mGeometries.size() )
867 if ( flags == 0 && mHasCachedValidity )
870 error = mValidityFailureReason;
871 return error.isEmpty();
878 mValidityFailureReason = !res ? error : QString();
879 mHasCachedValidity =
true;
893 geom->addZValue( zValue );
908 geom->addMValue( mValue );
948 geom->filterVertices( filter );
968 std::unique_ptr< QgsGeometryCollection > res;
976 if (
const QgsMultiPoint *mp = qgsgeometry_cast< const QgsMultiPoint *>(
this ) )
980 res = std::make_unique< QgsMultiPoint >();
988 if (
const QgsMultiLineString *ml = qgsgeometry_cast< const QgsMultiLineString *>(
this ) )
992 res = std::make_unique< QgsMultiLineString >();
1000 if (
const QgsMultiPolygon *mp = qgsgeometry_cast< const QgsMultiPolygon *>(
this ) )
1004 res = std::make_unique< QgsMultiPolygon>();
1024 if ( thisPartType == filterSinglePartType )
1026 res->addGeometry( part->clone() );
1030 return res.release();
1049 newCollection->addGeometry( geom->toCurveType() );
1051 return newCollection.release();
1071 res = geom->
transform( transformer, feedback );
1103 const QgsGeometryCollection *otherCollection = qgsgeometry_cast<const QgsGeometryCollection *>( other );
1104 if ( !otherCollection )
1113 const int comparison = aGeom->
compareTo( bGeom );
1114 if ( comparison != 0 )
@ AllowSelfTouchingHoles
Indicates that self-touching holes are permitted. OGC validity states that self-touching holes are NO...
WkbType
The WKB type describes the number of dimensions a geometry has.
@ GeometryCollection
GeometryCollection.
TransformDirection
Flags for raster layer temporal capabilities.
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.
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 void transform(const QgsCoordinateTransform &ct, Qgis::TransformDirection d=Qgis::TransformDirection::Forward, bool transformZ=false)=0
Transforms the geometry using a coordinate transform.
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.
QgsGeometryCollection * snappedToGrid(double hSpacing, double vSpacing, double dSpacing=0, double mSpacing=0) const override
Makes a new geometry with all the points or vertices snapped to the closest point of the grid.
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.
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.
QgsGeometryCollection * clone() const override
Clones the geometry by performing a deep copy.
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
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.
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.