QGIS API Documentation  3.26.3-Buenos Aires (65e4edfdad)
qgsmultipoint.cpp
Go to the documentation of this file.
1 /***************************************************************************
2  qgsmultipoint.cpp
3  -------------------------------------------------------------------
4 Date : 29 Oct 2014
5 Copyright : (C) 2014 by Marco Hugentobler
6 email : marco.hugentobler at sourcepole dot com
7  ***************************************************************************
8  * *
9  * This program is free software; you can redistribute it and/or modify *
10  * it under the terms of the GNU General Public License as published by *
11  * the Free Software Foundation; either version 2 of the License, or *
12  * (at your option) any later version. *
13  * *
14  ***************************************************************************/
15 
16 #include "qgsmultipoint.h"
17 #include "qgsapplication.h"
18 #include "qgsgeometryutils.h"
19 #include "qgspoint.h"
20 #include "qgswkbptr.h"
21 
22 #include <QJsonArray>
23 #include <QJsonObject>
24 #include <QRegularExpression>
25 #include <nlohmann/json.hpp>
26 
28 {
30 }
31 
33 {
34  return qgsgeometry_cast< QgsPoint * >( geometryN( index ) );
35 }
36 
37 const QgsPoint *QgsMultiPoint::pointN( int index ) const
38 {
39  return qgsgeometry_cast< const QgsPoint * >( geometryN( index ) );
40 }
41 
43 {
44  return QStringLiteral( "MultiPoint" );
45 }
46 
48 {
49  auto result = std::make_unique< QgsMultiPoint >();
50  result->mWkbType = mWkbType;
51  return result.release();
52 }
53 
55 {
56  return new QgsMultiPoint( *this );
57 }
58 
60 {
61  return clone();
62 }
63 
64 bool QgsMultiPoint::fromWkt( const QString &wkt )
65 {
66  QString collectionWkt( wkt );
67  //test for non-standard MultiPoint(x1 y1, x2 y2) format
68  const thread_local QRegularExpression regex( QStringLiteral( "^\\s*MultiPoint\\s*[ZM]*\\s*\\(\\s*[-\\d]" ), QRegularExpression::CaseInsensitiveOption );
69  const QRegularExpressionMatch match = regex.match( collectionWkt );
70  if ( match.hasMatch() )
71  {
72  //alternate style without extra brackets, upgrade to standard
73  collectionWkt.replace( '(', QLatin1String( "((" ) ).replace( ')', QLatin1String( "))" ) ).replace( ',', QLatin1String( "),(" ) );
74  }
75 
76  return fromCollectionWkt( collectionWkt, QVector<QgsAbstractGeometry *>() << new QgsPoint, QStringLiteral( "Point" ) );
77 }
78 
80 {
83 }
84 
85 QDomElement QgsMultiPoint::asGml2( QDomDocument &doc, int precision, const QString &ns, const AxisOrder axisOrder ) const
86 {
87  QDomElement elemMultiPoint = doc.createElementNS( ns, QStringLiteral( "MultiPoint" ) );
88 
89  if ( isEmpty() )
90  return elemMultiPoint;
91 
92  for ( const QgsAbstractGeometry *geom : mGeometries )
93  {
94  if ( qgsgeometry_cast<const QgsPoint *>( geom ) )
95  {
96  QDomElement elemPointMember = doc.createElementNS( ns, QStringLiteral( "pointMember" ) );
97  elemPointMember.appendChild( geom->asGml2( doc, precision, ns, axisOrder ) );
98  elemMultiPoint.appendChild( elemPointMember );
99  }
100  }
101 
102  return elemMultiPoint;
103 }
104 
105 QDomElement QgsMultiPoint::asGml3( QDomDocument &doc, int precision, const QString &ns, const QgsAbstractGeometry::AxisOrder axisOrder ) const
106 {
107  QDomElement elemMultiPoint = doc.createElementNS( ns, QStringLiteral( "MultiPoint" ) );
108 
109  if ( isEmpty() )
110  return elemMultiPoint;
111 
112  for ( const QgsAbstractGeometry *geom : mGeometries )
113  {
114  if ( qgsgeometry_cast<const QgsPoint *>( geom ) )
115  {
116  QDomElement elemPointMember = doc.createElementNS( ns, QStringLiteral( "pointMember" ) );
117  elemPointMember.appendChild( geom->asGml3( doc, precision, ns, axisOrder ) );
118  elemMultiPoint.appendChild( elemPointMember );
119  }
120  }
121 
122  return elemMultiPoint;
123 }
124 
126 {
127  json j
128  {
129  { "type", "MultiPoint" },
130  { "coordinates", json::array() },
131  };
132  for ( const QgsAbstractGeometry *geom : std::as_const( mGeometries ) )
133  {
134  const QgsPoint *point = static_cast<const QgsPoint *>( geom );
135  if ( point->is3D() )
136  j[ "coordinates" ].push_back( { qgsRound( point->x(), precision ), qgsRound( point->y(), precision ), qgsRound( point->z(), precision ) } );
137  else
138  j[ "coordinates" ].push_back( { qgsRound( point->x(), precision ), qgsRound( point->y(), precision ) } );
139  }
140  return j;
141 }
142 
143 
145 {
146  return mGeometries.size();
147 }
148 
150 {
151  if ( !qgsgeometry_cast<QgsPoint *>( g ) )
152  {
153  delete g;
154  return false;
155  }
156  if ( mGeometries.empty() )
157  {
159  }
160  if ( is3D() && !g->is3D() )
161  g->addZValue();
162  else if ( !is3D() && g->is3D() )
163  g->dropZValue();
164  if ( isMeasure() && !g->isMeasure() )
165  g->addMValue();
166  else if ( !isMeasure() && g->isMeasure() )
167  g->dropMValue();
168 
170 }
171 
173 {
174  if ( !g || QgsWkbTypes::flatType( g->wkbType() ) != QgsWkbTypes::Point )
175  {
176  delete g;
177  return false;
178  }
179 
180  return QgsGeometryCollection::insertGeometry( g, index );
181 }
182 
184 {
185  return nullptr;
186 }
187 
189 {
190  if ( id.part < 0 || id.part >= mGeometries.count() || id.vertex != 0 || id.ring != 0 )
191  return -1;
192 
193  return id.part; // can shortcut the calculation, since each part will have 1 vertex
194 }
195 
197 {
198  return 0.0;
199 }
200 
201 bool QgsMultiPoint::isValid( QString &, Qgis::GeometryValidityFlags ) const
202 {
203  return true;
204 }
205 
206 void QgsMultiPoint::filterVertices( const std::function<bool ( const QgsPoint & )> &filter )
207 {
208  mGeometries.erase( std::remove_if( mGeometries.begin(), mGeometries.end(), // clazy:exclude=detaching-member
209  [&filter]( const QgsAbstractGeometry * part )
210  {
211  if ( const QgsPoint *point = qgsgeometry_cast< const QgsPoint * >( part ) )
212  {
213  if ( !filter( *point ) )
214  {
215  delete point;
216  return true;
217  }
218  else
219  {
220  return false;
221  }
222  }
223  else
224  {
225  delete part;
226  return true;
227  }
228  } ), mGeometries.end() ); // clazy:exclude=detaching-member
229 }
230 
232 {
233  return true;
234 }
QgsAbstractGeometry::dropMValue
virtual bool dropMValue()=0
Drops any measure values which exist in the geometry.
QgsGeometryCollection::insertGeometry
virtual bool insertGeometry(QgsAbstractGeometry *g, int index)
Inserts a geometry before a specified index and takes ownership.
Definition: qgsgeometrycollection.cpp:268
QgsGeometryCollection::isEmpty
bool isEmpty() const override SIP_HOLDGIL
Returns true if the geometry is empty.
Definition: qgsgeometrycollection.cpp:243
QgsMultiPoint::insertGeometry
bool insertGeometry(QgsAbstractGeometry *g, int index) override
Inserts a geometry before a specified index and takes ownership.
Definition: qgsmultipoint.cpp:172
QgsAbstractGeometry::dropZValue
virtual bool dropZValue()=0
Drops any z-dimensions which exist in the geometry.
QgsWkbTypes::Point
@ Point
Definition: qgswkbtypes.h:72
QgsMultiPoint::vertexNumberFromVertexId
int vertexNumberFromVertexId(QgsVertexId id) const override
Returns the vertex number corresponding to a vertex id.
Definition: qgsmultipoint.cpp:188
QgsAbstractGeometry::addZValue
virtual bool addZValue(double zValue=0)=0
Adds a z-dimension to the geometry, initialized to a preset value.
QgsPoint
Point geometry type, with support for z-dimension and m-values.
Definition: qgspoint.h:48
qgswkbptr.h
QgsMultiPoint::nCoordinates
int nCoordinates() const override SIP_HOLDGIL
Returns the number of nodes contained in the geometry.
Definition: qgsmultipoint.cpp:144
QgsAbstractGeometry::addMValue
virtual bool addMValue(double mValue=0)=0
Adds a measure to the geometry, initialized to a preset value.
QgsWkbTypes::flatType
static Type flatType(Type type) SIP_HOLDGIL
Returns the flat type for a WKB type.
Definition: qgswkbtypes.h:732
QgsMultiPoint::addGeometry
bool addGeometry(QgsAbstractGeometry *g) override
Adds a geometry and takes ownership. Returns true in case of success.
Definition: qgsmultipoint.cpp:149
qgsmultipoint.h
QgsPoint::z
double z
Definition: qgspoint.h:71
QgsMultiPoint::clone
QgsMultiPoint * clone() const override
Clones the geometry by performing a deep copy.
Definition: qgsmultipoint.cpp:54
QgsGeometryCollection::mGeometries
QVector< QgsAbstractGeometry * > mGeometries
Definition: qgsgeometrycollection.h:338
qgspoint.h
QgsAbstractGeometry::mWkbType
QgsWkbTypes::Type mWkbType
Definition: qgsabstractgeometry.h:1119
QgsMultiPoint::createEmptyWithSameType
QgsMultiPoint * createEmptyWithSameType() const override
Creates a new geometry with the same class and same WKB type as the original and transfers ownership.
Definition: qgsmultipoint.cpp:47
qgsapplication.h
QgsAbstractGeometry::isMeasure
bool isMeasure() const SIP_HOLDGIL
Returns true if the geometry contains m values.
Definition: qgsabstractgeometry.h:228
QgsPoint::y
double y
Definition: qgspoint.h:70
precision
int precision
Definition: qgswfsgetfeature.cpp:103
QgsMultiPoint::clear
void clear() override
Clears the geometry, ie reset it to a null geometry.
Definition: qgsmultipoint.cpp:79
QgsAbstractGeometry::wkbType
QgsWkbTypes::Type wkbType() const SIP_HOLDGIL
Returns the WKB type of the geometry.
Definition: qgsabstractgeometry.h:206
QgsAbstractGeometry::AxisOrder
AxisOrder
Axis order for GML generation.
Definition: qgsabstractgeometry.h:138
QgsGeometryCollection::addGeometry
virtual bool addGeometry(QgsAbstractGeometry *g)
Adds a geometry and takes ownership. Returns true in case of success.
Definition: qgsgeometrycollection.cpp:256
qgsRound
double qgsRound(double number, int places)
Returns a double number, rounded (as close as possible) to the specified number of places.
Definition: qgis.h:2319
QgsGeometryCollection::fromCollectionWkt
bool fromCollectionWkt(const QString &wkt, const QVector< QgsAbstractGeometry * > &subtypes, const QString &defaultChildWkbType=QString())
Reads a collection from a WKT string.
Definition: qgsgeometrycollection.cpp:744
QgsMultiPoint
Multi point geometry collection.
Definition: qgsmultipoint.h:29
QgsAbstractGeometry
Abstract base class for all geometries.
Definition: qgsabstractgeometry.h:79
qgsgeometryutils.h
QgsAbstractGeometry::is3D
bool is3D() const SIP_HOLDGIL
Returns true if the geometry is 3D and contains a z-value.
Definition: qgsabstractgeometry.h:219
QgsMultiPoint::pointN
QgsPoint * pointN(int index)
Returns the point with the specified index.
Definition: qgsmultipoint.cpp:32
QgsMultiPoint::fromWkt
bool fromWkt(const QString &wkt) override
Sets the geometry from a WKT string.
Definition: qgsmultipoint.cpp:64
QgsMultiPoint::toCurveType
QgsMultiPoint * toCurveType() const override
Returns the geometry converted to the more generic curve type.
Definition: qgsmultipoint.cpp:59
QgsGeometryCollection::geometryN
const QgsAbstractGeometry * geometryN(int n) const
Returns a const reference to a geometry from within the collection.
Definition: qgsgeometrycollection.h:86
QgsMultiPoint::asGml3
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.
Definition: qgsmultipoint.cpp:105
QgsWkbTypes::MultiPoint
@ MultiPoint
Definition: qgswkbtypes.h:76
QgsVertexId
Utility class for identifying a unique vertex within a geometry.
Definition: qgsvertexid.h:30
QgsMultiPoint::segmentLength
double segmentLength(QgsVertexId startVertex) const override
Returns the length of the segment of the geometry which begins at startVertex.
Definition: qgsmultipoint.cpp:196
QgsMultiPoint::filterVertices
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...
Definition: qgsmultipoint.cpp:206
QgsGeometryCollection::clear
void clear() override
Clears the geometry, ie reset it to a null geometry.
Definition: qgsgeometrycollection.cpp:124
QgsMultiPoint::asGml2
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.
Definition: qgsmultipoint.cpp:85
QgsAbstractGeometry::setZMTypeFromSubGeometry
void setZMTypeFromSubGeometry(const QgsAbstractGeometry *subggeom, QgsWkbTypes::Type baseGeomType)
Updates the geometry type based on whether sub geometries contain z or m values.
Definition: qgsabstractgeometry.cpp:77
QgsMultiPoint::isValid
bool isValid(QString &error, Qgis::GeometryValidityFlags flags=Qgis::GeometryValidityFlags()) const override SIP_HOLDGIL
Checks validity of the geometry, and returns true if the geometry is valid.
Definition: qgsmultipoint.cpp:201
QgsMultiPoint::boundary
QgsAbstractGeometry * boundary() const override
Returns the closure of the combinatorial boundary of the geometry (ie the topological boundary of the...
Definition: qgsmultipoint.cpp:183
QgsPoint::x
double x
Definition: qgspoint.h:69
QgsMultiPoint::geometryType
QString geometryType() const override
Returns a unique string representing the geometry type.
Definition: qgsmultipoint.cpp:42
QgsMultiPoint::wktOmitChildType
bool wktOmitChildType() const override
Returns whether child type names are omitted from Wkt representations of the collection.
Definition: qgsmultipoint.cpp:231
QgsMultiPoint::QgsMultiPoint
QgsMultiPoint() SIP_HOLDGIL
Constructor for an empty multipoint geometry.
Definition: qgsmultipoint.cpp:27
QgsMultiPoint::asJsonObject
json asJsonObject(int precision=17) const override
Returns a json object representation of the geometry.
Definition: qgsmultipoint.cpp:125