QGIS API Documentation  3.8.0-Zanzibar (11aff65)
qgsgeometryduplicatenodescheck.cpp
Go to the documentation of this file.
1 /***************************************************************************
2  qgsgeometryduplicatenodescheck.cpp
3  ---------------------
4  begin : September 2015
5  copyright : (C) 2014 by Sandro Mani / Sourcepole AG
6  email : smani at sourcepole dot ch
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 
18 #include "qgsgeometryutils.h"
19 #include "qgsfeaturepool.h"
20 #include "qgsgeometrycheckerror.h"
21 
22 void QgsGeometryDuplicateNodesCheck::collectErrors( const QMap<QString, QgsFeaturePool *> &featurePools, QList<QgsGeometryCheckError *> &errors, QStringList &messages, QgsFeedback *feedback, const LayerFeatureIds &ids ) const
23 {
24  Q_UNUSED( messages )
25 
26  QMap<QString, QgsFeatureIds> featureIds = ids.isEmpty() ? allLayerFeatureIds( featurePools ) : ids.toMap();
27  QgsGeometryCheckerUtils::LayerFeatures layerFeatures( featurePools, featureIds, compatibleGeometryTypes(), feedback, mContext );
28  for ( const QgsGeometryCheckerUtils::LayerFeature &layerFeature : layerFeatures )
29  {
30  const QgsAbstractGeometry *geom = layerFeature.geometry().constGet();
31  for ( int iPart = 0, nParts = geom->partCount(); iPart < nParts; ++iPart )
32  {
33  for ( int iRing = 0, nRings = geom->ringCount( iPart ); iRing < nRings; ++iRing )
34  {
35  int nVerts = QgsGeometryCheckerUtils::polyLineSize( geom, iPart, iRing );
36  if ( nVerts < 2 )
37  continue;
38  for ( int iVert = nVerts - 1, jVert = 0; jVert < nVerts; iVert = jVert++ )
39  {
40  QgsPoint pi = geom->vertexAt( QgsVertexId( iPart, iRing, iVert ) );
41  QgsPoint pj = geom->vertexAt( QgsVertexId( iPart, iRing, jVert ) );
43  {
44  errors.append( new QgsGeometryCheckError( this, layerFeature, pj, QgsVertexId( iPart, iRing, jVert ) ) );
45  }
46  }
47  }
48  }
49  }
50 }
51 
52 void QgsGeometryDuplicateNodesCheck::fixError( const QMap<QString, QgsFeaturePool *> &featurePools, QgsGeometryCheckError *error, int method, const QMap<QString, int> & /*mergeAttributeIndices*/, Changes &changes ) const
53 {
54  QgsFeaturePool *featurePool = featurePools[ error->layerId() ];
55  QgsFeature feature;
56  if ( !featurePool->getFeature( error->featureId(), feature ) )
57  {
58  error->setObsolete();
59  return;
60  }
61  QgsGeometry featureGeom = feature.geometry();
62  QgsAbstractGeometry *geom = featureGeom.get();
63  QgsVertexId vidx = error->vidx();
64 
65  // Check if point still exists
66  if ( !vidx.isValid( geom ) )
67  {
68  error->setObsolete();
69  return;
70  }
71 
72  // Check if error still applies
73  int nVerts = QgsGeometryCheckerUtils::polyLineSize( geom, vidx.part, vidx.ring );
74  QgsPoint pi = geom->vertexAt( QgsVertexId( vidx.part, vidx.ring, ( vidx.vertex + nVerts - 1 ) % nVerts ) );
75  QgsPoint pj = geom->vertexAt( error->vidx() );
77  {
78  error->setObsolete();
79  return;
80  }
81 
82  // Fix error
83  if ( method == NoChange )
84  {
85  error->setFixed( method );
86  }
87  else if ( method == RemoveDuplicates )
88  {
89  if ( !QgsGeometryCheckerUtils::canDeleteVertex( geom, vidx.part, vidx.ring ) )
90  {
91  error->setFixFailed( tr( "Resulting geometry is degenerate" ) );
92  }
93  else if ( !geom->deleteVertex( error->vidx() ) )
94  {
95  error->setFixFailed( tr( "Failed to delete vertex" ) );
96  }
97  else
98  {
99  feature.setGeometry( featureGeom );
100  featurePool->updateFeature( feature );
101  error->setFixed( method );
102  changes[error->layerId()][error->featureId()].append( Change( ChangeNode, ChangeRemoved, error->vidx() ) );
103  }
104  }
105  else
106  {
107  error->setFixFailed( tr( "Unknown method" ) );
108  }
109 }
110 
112 {
113  static QStringList methods = QStringList() << tr( "Delete duplicate node" ) << tr( "No action" );
114  return methods;
115 }
116 
118 {
120 }
virtual bool deleteVertex(QgsVertexId position)=0
Deletes a vertex within the geometry.
bool isValid() const
Returns true if the vertex id is valid.
QMap< QString, QgsFeatureIds > toMap() const
Contains a set of layers and feature ids in those layers to pass to a geometry check.
CheckType
The type of a check.
static int polyLineSize(const QgsAbstractGeometry *geom, int iPart, int iRing, bool *isClosed=nullptr)
Returns the number of points in a polyline, accounting for duplicate start and end point if the polyl...
A geometry is the spatial representation of a feature.
Definition: qgsgeometry.h:111
void setObsolete()
Set the error status to obsolete.
The feature class encapsulates a single feature including its id, geometry and a list of field/values...
Definition: qgsfeature.h:55
void fixError(const QMap< QString, QgsFeaturePool *> &featurePools, QgsGeometryCheckError *error, int method, const QMap< QString, int > &mergeAttributeIndices, Changes &changes) const override
Fix the error error with the specified method.
virtual void updateFeature(QgsFeature &feature)=0
Updates a feature in this pool.
void setFixFailed(const QString &reason)
Set the error status to failed and specify the reason for failure.
Base class for feedback objects to be used for cancellation of something running in a worker thread...
Definition: qgsfeedback.h:44
Something has been removed.
Utility class for identifying a unique vertex within a geometry.
bool getFeature(QgsFeatureId id, QgsFeature &feature)
Retrieves the feature with the specified id into feature.
A layer feature combination to uniquely identify and access a feature in a set of layers...
virtual int ringCount(int part=0) const =0
Returns the number of rings of which this geometry is built.
Abstract base class for all geometries.
QMap< QString, QgsFeatureIds > allLayerFeatureIds(const QMap< QString, QgsFeaturePool *> &featurePools) const
Returns all layers and feature ids.
const QString & layerId() const
The id of the layer on which this error has been detected.
const double tolerance
The tolerance to allow for in geometry checks.
const QgsGeometryCheckContext * mContext
Point geometry type, with support for z-dimension and m-values.
Definition: qgspoint.h:37
A list of layers and feature ids for each of these layers.
QgsAbstractGeometry * get()
Returns a modifiable (non-const) reference to the underlying abstract geometry primitive.
static bool canDeleteVertex(const QgsAbstractGeometry *geom, int iPart, int iRing)
static double sqrDistance2D(const QgsPoint &pt1, const QgsPoint &pt2)
Returns the squared 2D distance between two points.
QMap< QString, QMap< QgsFeatureId, QList< QgsGeometryCheck::Change > > > Changes
A collection of changes.
A feature pool is based on a vector layer and caches features.
Descripts a change to fix a geometry.
void setFixed(int method)
Set the status to fixed and specify the method that has been used to fix the error.
void setGeometry(const QgsGeometry &geometry)
Set the feature&#39;s geometry.
Definition: qgsfeature.cpp:137
QStringList resolutionMethods() const override
Returns a list of descriptions for available resolutions for errors.
static QgsGeometryCheck::CheckType factoryCheckType()
QgsFeatureId featureId() const
The id of the feature on which this error has been detected.
QList< QgsWkbTypes::GeometryType > compatibleGeometryTypes() const override
A list of geometry types for which this check can be performed.
QgsGeometry geometry
Definition: qgsfeature.h:67
This represents an error reported by a geometry check.
This change happens on node level.
virtual QgsPoint vertexAt(QgsVertexId id) const =0
Returns the point corresponding to a specified vertex id.
const QgsVertexId & vidx() const
The id of the affected vertex.
void collectErrors(const QMap< QString, QgsFeaturePool *> &featurePools, QList< QgsGeometryCheckError *> &errors, QStringList &messages, QgsFeedback *feedback, const LayerFeatureIds &ids=LayerFeatureIds()) const override
The main worker method.
virtual int partCount() const =0
Returns count of parts contained in the geometry.
The check controls individual nodes.