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 const QgsGeometryCollection *otherCollection = qgsgeometry_cast< const QgsGeometryCollection * >( &other );
79 if ( !otherCollection )
115 auto result = std::make_unique< QgsGeometryCollection >();
117 return result.release();
134 std::unique_ptr<QgsGeometryCollection> result;
138 std::unique_ptr<QgsAbstractGeometry> gridified { geom->
snappedToGrid( hSpacing, vSpacing, dSpacing, mSpacing ) };
148 return result.release();
156 if ( geom->removeDuplicateNodes( epsilon, useZValues ) ) result =
true;
180 if (
id.part < 0 || id.part >=
mGeometries.count() )
187 if ( part ==
id.part )
189 int partNumber = geometry->vertexNumberFromVertexId(
QgsVertexId( 0,
id.ring,
id.vertex ) );
190 if ( partNumber == -1 )
192 return number + partNumber;
196 number += geometry->nCoordinates();
210 if ( !mBoundingBox.
isNull() )
222 if ( geometry->boundingBoxIntersects( box3d ) )
251 if ( !geometry->isEmpty() )
276 index = std::min(
static_cast<int>(
mGeometries.count() ), index );
299 geometry->normalize();
303 return a->compareTo( b ) > 0;
310 QVector< QgsAbstractGeometry * >::const_iterator it =
mGeometries.constBegin();
313 int dim = ( *it )->dimension();
324 return QStringLiteral(
"GeometryCollection" );
331 g->transform( ct, d, transformZ );
340 g->transform( t, zTranslate, zScale, mTranslate, mScale );
347 QVector< QgsAbstractGeometry * >::const_iterator it =
mGeometries.constBegin();
359 QPainterPath partPath = geom->asQPainterPath();
360 if ( !partPath.isEmpty() )
361 p.addPath( partPath );
380 wkbPtr >> nGeometries;
382 QVector<QgsAbstractGeometry *> geometryListBackup =
mGeometries;
385 for (
int i = 0; i < nGeometries; ++i )
398 qDeleteAll( geometryListBackup );
417 int binarySize =
sizeof( char ) +
sizeof( quint32 ) +
sizeof( quint32 );
422 binarySize += geom->wkbSize( flags );
431 int countNonNull = 0;
444 wkb << static_cast<quint32>(
wkbType() );
445 wkb << static_cast<quint32>( countNonNull );
450 wkb << geom->asWkb( flags );
461 wkt += QLatin1String(
" EMPTY" );
464 wkt += QLatin1String(
" (" );
467 QString childWkt = geom->asWkt(
precision );
470 childWkt = childWkt.mid( childWkt.indexOf(
'(' ) );
472 wkt += childWkt +
',';
474 if ( wkt.endsWith(
',' ) )
485 QDomElement elemMultiGeometry = doc.createElementNS( ns, QStringLiteral(
"MultiGeometry" ) );
488 QDomElement elemGeometryMember = doc.createElementNS( ns, QStringLiteral(
"geometryMember" ) );
489 elemGeometryMember.appendChild( geom->asGml2( doc,
precision, ns, axisOrder ) );
490 elemMultiGeometry.appendChild( elemGeometryMember );
492 return elemMultiGeometry;
497 QDomElement elemMultiGeometry = doc.createElementNS( ns, QStringLiteral(
"MultiGeometry" ) );
500 QDomElement elemGeometryMember = doc.createElementNS( ns, QStringLiteral(
"geometryMember" ) );
501 elemGeometryMember.appendChild( geom->asGml3( doc,
precision, ns, axisOrder ) );
502 elemMultiGeometry.appendChild( elemGeometryMember );
504 return elemMultiGeometry;
509 json coordinates( json::array( ) );
512 coordinates.push_back( geom->asJsonObject(
precision ) );
516 {
"type",
"GeometryCollection" },
517 {
"geometries", coordinates }
524 kml.append( QLatin1String(
"<MultiGeometry>" ) );
525 const QVector< QgsAbstractGeometry * > &geometries =
mGeometries;
528 kml.append( geometry->asKml(
precision ) );
530 kml.append( QLatin1String(
"</MultiGeometry>" ) );
536 if ( mBoundingBox.
isNull() )
565 mHasCachedValidity =
false;
566 mValidityFailureReason.clear();
573 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
574 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
578 QgsCoordinateSequence::const_iterator cIt = geomCoords.constBegin();
579 for ( ; cIt != geomCoords.constEnd(); ++cIt )
581 sequence.push_back( *cIt );
592 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
593 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
595 count += ( *geomIt )->nCoordinates();
644 bool success =
mGeometries.at( position.
part )->insertVertex( position, vertex );
659 bool success =
mGeometries.at( position.
part )->moveVertex( position, newPos );
698 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
699 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
701 length += ( *geomIt )->length();
709 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
710 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
712 area += ( *geomIt )->area();
720 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
721 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
736 qDeleteAll( subtypes );
741 QString secondWithoutParentheses =
parts.second;
742 secondWithoutParentheses = secondWithoutParentheses.remove(
'(' ).remove(
')' ).simplified().remove(
' ' );
743 if ( (
parts.second.compare( QLatin1String(
"EMPTY" ), Qt::CaseInsensitive ) == 0 ) ||
744 secondWithoutParentheses.isEmpty() )
746 qDeleteAll( subtypes );
750 QString defChildWkbType = QStringLiteral(
"%1%2%3 " ).arg( defaultChildWkbType,
is3D() ? QStringLiteral(
"Z" ) : QString(),
isMeasure() ? QStringLiteral(
"M" ) : QString() );
753 for (
const QString &childWkt : blocks )
757 bool success =
false;
773 qDeleteAll( subtypes );
777 qDeleteAll( subtypes );
785 hasZ = hasZ || geom->is3D();
786 hasM = hasM || geom->isMeasure();
800 QVector< QgsAbstractGeometry * >::const_iterator it =
mGeometries.constBegin();
803 if ( ( *it )->hasCurvedSegments() )
815 if ( !geomCollection )
821 QVector< QgsAbstractGeometry * >::const_iterator geomIt =
mGeometries.constBegin();
822 for ( ; geomIt !=
mGeometries.constEnd(); ++geomIt )
824 geomCollection->
addGeometry( ( *geomIt )->segmentize( tolerance, toleranceType ) );
826 return geom.release();
888 if (
id.part < 0 || id.part >=
mGeometries.size() )
904 if ( flags == 0 && mHasCachedValidity )
907 error = mValidityFailureReason;
908 return error.isEmpty();
915 mValidityFailureReason = !res ? error : QString();
916 mHasCachedValidity =
true;
930 geom->addZValue( zValue );
945 geom->addMValue( mValue );
985 geom->filterVertices( filter );
1016 newCollection->addGeometry( geom->toCurveType() );
1018 return newCollection.release();
1038 res = geom->
transform( transformer, feedback );
1070 const QgsGeometryCollection *otherCollection = qgsgeometry_cast<const QgsGeometryCollection *>( other );
1071 if ( !otherCollection )
1080 const int comparison = aGeom->
compareTo( bGeom );
1081 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.
QgsBox3D boundingBox3D() const override
Returns the 3D bounding box for the geometry.
bool operator==(const QgsAbstractGeometry &other) const override
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.
bool operator!=(const QgsAbstractGeometry &other) const override
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::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 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.