QGIS API Documentation  3.23.0-Master (22c16f2067)
qgsgeometrydanglecheck.cpp
Go to the documentation of this file.
1 /***************************************************************************
2  qgsgeometrydanglecheck.cpp
3  ---------------------
4  begin : June 2017
5  copyright : (C) 2017 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 
17 #include "qgsgeometrydanglecheck.h"
18 #include "qgslinestring.h"
19 #include "qgsvectorlayer.h"
20 #include "qgsgeometrycheckerror.h"
21 
22 void QgsGeometryDangleCheck::collectErrors( const QMap<QString, QgsFeaturePool *> &featurePools, QList<QgsGeometryCheckError *> &errors, QStringList &messages, QgsFeedback *feedback, const LayerFeatureIds &ids ) const
23 {
24  Q_UNUSED( messages )
25  const QMap<QString, QgsFeatureIds> featureIds = ids.isEmpty() ? allLayerFeatureIds( featurePools ) : ids.toMap();
26  const QgsGeometryCheckerUtils::LayerFeatures layerFeatures( featurePools, featureIds, compatibleGeometryTypes(), feedback, mContext );
27  for ( const QgsGeometryCheckerUtils::LayerFeature &layerFeature : layerFeatures )
28  {
29  const QgsAbstractGeometry *geom = layerFeature.geometry().constGet();
30  for ( int iPart = 0, nParts = geom->partCount(); iPart < nParts; ++iPart )
31  {
32  const QgsLineString *line = dynamic_cast<const QgsLineString *>( QgsGeometryCheckerUtils::getGeomPart( geom, iPart ) );
33  if ( !line )
34  {
35  // Should not happen
36  continue;
37  }
38  // Check that start and end node lie on a line
39  const int nVerts = geom->vertexCount( iPart, 0 );
40  const QgsPoint &p1 = geom->vertexAt( QgsVertexId( iPart, 0, 0 ) );
41  const QgsPoint &p2 = geom->vertexAt( QgsVertexId( iPart, 0, nVerts - 1 ) );
42 
43  bool p1touches = QgsGeometryCheckerUtils::pointOnLine( p1, line, mContext->tolerance, true );
44  bool p2touches = QgsGeometryCheckerUtils::pointOnLine( p2, line, mContext->tolerance, true );
45 
46  if ( p1touches && p2touches )
47  {
48  // Both endpoints lie on line itself
49  continue;
50  }
51 
52  // Check whether endpoints line on another line in the layer
53  const QgsGeometryCheckerUtils::LayerFeatures checkFeatures( featurePools, QList<QString>() << layerFeature.layer()->id(), line->boundingBox(), {QgsWkbTypes::LineGeometry}, mContext );
54  for ( const QgsGeometryCheckerUtils::LayerFeature &checkFeature : checkFeatures )
55  {
56  const QgsAbstractGeometry *testGeom = checkFeature.geometry().constGet();
57  for ( int jPart = 0, mParts = testGeom->partCount(); jPart < mParts; ++jPart )
58  {
59  if ( checkFeature.feature().id() == layerFeature.feature().id() && iPart == jPart )
60  {
61  // Skip current feature part, it was already checked above
62  continue;
63  }
64  const QgsLineString *testLine = dynamic_cast<const QgsLineString *>( QgsGeometryCheckerUtils::getGeomPart( testGeom, jPart ) );
65  if ( !testLine )
66  {
67  continue;
68  }
69  p1touches = p1touches || QgsGeometryCheckerUtils::pointOnLine( p1, testLine, mContext->tolerance );
70  p2touches = p2touches || QgsGeometryCheckerUtils::pointOnLine( p2, testLine, mContext->tolerance );
71  if ( p1touches && p2touches )
72  {
73  break;
74  }
75  }
76  if ( p1touches && p2touches )
77  {
78  break;
79  }
80  }
81  if ( p1touches && p2touches )
82  {
83  continue;
84  }
85  if ( !p1touches )
86  {
87  errors.append( new QgsGeometryCheckError( this, layerFeature, p1, QgsVertexId( iPart, 0, 0 ) ) );
88  }
89  if ( !p2touches )
90  {
91  errors.append( new QgsGeometryCheckError( this, layerFeature, p2, QgsVertexId( iPart, 0, nVerts - 1 ) ) );
92  }
93  }
94  }
95 }
96 
97 void QgsGeometryDangleCheck::fixError( const QMap<QString, QgsFeaturePool *> &featurePools, QgsGeometryCheckError *error, int method, const QMap<QString, int> & /*mergeAttributeIndices*/, Changes & /*changes*/ ) const
98 {
99  Q_UNUSED( featurePools )
100 
101  if ( method == NoChange )
102  {
103  error->setFixed( method );
104  }
105  else
106  {
107  error->setFixFailed( tr( "Unknown method" ) );
108  }
109 }
110 
112 {
113  static const QStringList methods = QStringList() << tr( "No action" );
114  return methods;
115 }
116 
118 {
120 }
Abstract base class for all geometries.
virtual int vertexCount(int part=0, int ring=0) const =0
Returns the number of vertices of which this geometry is built.
virtual QgsPoint vertexAt(QgsVertexId id) const =0
Returns the point corresponding to a specified vertex id.
virtual int partCount() const =0
Returns count of parts contained in the geometry.
QgsRectangle boundingBox() const override
Returns the minimal bounding box for the geometry.
Definition: qgscurve.cpp:238
Base class for feedback objects to be used for cancellation of something running in a worker thread.
Definition: qgsfeedback.h:45
const double tolerance
The tolerance to allow for in geometry checks.
This represents an error reported by a geometry check.
void setFixed(int method)
Set the status to fixed and specify the method that has been used to fix the error.
void setFixFailed(const QString &reason)
Set the error status to failed and specify the reason for failure.
QMap< QString, QMap< QgsFeatureId, QList< QgsGeometryCheck::Change > > > Changes
A collection of changes.
const QgsGeometryCheckContext * mContext
CheckType
The type of a check.
@ FeatureNodeCheck
The check controls individual nodes.
QMap< QString, QgsFeatureIds > allLayerFeatureIds(const QMap< QString, QgsFeaturePool * > &featurePools) const
Returns all layers and feature ids.
A layer feature combination to uniquely identify and access a feature in a set of layers.
Contains a set of layers and feature ids in those layers to pass to a geometry check.
static QgsAbstractGeometry * getGeomPart(QgsAbstractGeometry *geom, int partIdx)
static bool pointOnLine(const QgsPoint &p, const QgsLineString *line, double tol, bool excludeExtremities=false)
void fixError(const QMap< QString, QgsFeaturePool * > &featurePools, QgsGeometryCheckError *error, int method, const QMap< QString, int > &mergeAttributeIndices, Changes &changes) const override
Fixes the error error with the specified method.
static QgsGeometryCheck::CheckType factoryCheckType()
QList< QgsWkbTypes::GeometryType > compatibleGeometryTypes() const override
A list of geometry types for which this check can be performed.
Q_DECL_DEPRECATED QStringList resolutionMethods() const override
Returns a list of descriptions for available resolutions for errors.
void collectErrors(const QMap< QString, QgsFeaturePool * > &featurePools, QList< QgsGeometryCheckError * > &errors, QStringList &messages, QgsFeedback *feedback, const LayerFeatureIds &ids=LayerFeatureIds()) const override
The main worker method.
Line string geometry type, with support for z-dimension and m-values.
Definition: qgslinestring.h:44
Point geometry type, with support for z-dimension and m-values.
Definition: qgspoint.h:49
A list of layers and feature ids for each of these layers.
QMap< QString, QgsFeatureIds > toMap() const
Utility class for identifying a unique vertex within a geometry.
Definition: qgsvertexid.h:31