28QString QgsGeometryCheckMissingVertexAlgorithm::name()
const
30 return QStringLiteral(
"checkgeometrymissingvertex" );
33QString QgsGeometryCheckMissingVertexAlgorithm::displayName()
const
35 return QObject::tr(
"Missing vertices along borders" );
38QString QgsGeometryCheckMissingVertexAlgorithm::shortDescription()
const
40 return QObject::tr(
"Detects missing vertices along polygon borders for topology compliance." );
43QStringList QgsGeometryCheckMissingVertexAlgorithm::tags()
const
45 return QObject::tr(
"check,geometry,missing,vertex" ).split(
',' );
48QString QgsGeometryCheckMissingVertexAlgorithm::group()
const
50 return QObject::tr(
"Check geometry" );
53QString QgsGeometryCheckMissingVertexAlgorithm::groupId()
const
55 return QStringLiteral(
"checkgeometry" );
58QString QgsGeometryCheckMissingVertexAlgorithm::shortHelpString()
const
60 return QObject::tr(
"This algorithm checks for missing vertices along polygon borders.\n"
61 "To be topologically correct, a vertex at the junction of two polygons must be present on both polygons. "
62 "Missing vertices are errors." );
70QgsGeometryCheckMissingVertexAlgorithm *QgsGeometryCheckMissingVertexAlgorithm::createInstance()
const
72 return new QgsGeometryCheckMissingVertexAlgorithm();
75void QgsGeometryCheckMissingVertexAlgorithm::initAlgorithm(
const QVariantMap &configuration )
77 Q_UNUSED( configuration )
82 QStringLiteral(
"INPUT" ), QObject::tr(
"Input layer" ),
87 QStringLiteral(
"UNIQUE_ID" ), QObject::tr(
"Unique feature identifier" ), QString(), QStringLiteral(
"INPUT" )
98 std::unique_ptr<QgsProcessingParameterNumber> tolerance = std::make_unique<QgsProcessingParameterNumber>(
102 tolerance->setHelp( QObject::tr(
"The \"Tolerance\" advanced parameter defines the numerical precision of geometric operations, "
103 "given as an integer n, meaning that any difference smaller than 10⁻ⁿ (in map units) is considered zero." ) );
104 addParameter( tolerance.release() );
109 mTolerance = parameterAsInt( parameters, QStringLiteral(
"TOLERANCE" ), context );
114QgsFields QgsGeometryCheckMissingVertexAlgorithm::outputFields()
117 fields.
append(
QgsField( QStringLiteral(
"gc_layerid" ), QMetaType::QString ) );
118 fields.
append(
QgsField( QStringLiteral(
"gc_layername" ), QMetaType::QString ) );
119 fields.
append(
QgsField( QStringLiteral(
"gc_partidx" ), QMetaType::Int ) );
120 fields.
append(
QgsField( QStringLiteral(
"gc_ringidx" ), QMetaType::Int ) );
121 fields.
append(
QgsField( QStringLiteral(
"gc_vertidx" ), QMetaType::Int ) );
122 fields.
append(
QgsField( QStringLiteral(
"gc_errorx" ), QMetaType::Double ) );
123 fields.
append(
QgsField( QStringLiteral(
"gc_errory" ), QMetaType::Double ) );
124 fields.
append(
QgsField( QStringLiteral(
"gc_error" ), QMetaType::QString ) );
132 const std::unique_ptr<QgsProcessingFeatureSource> input( parameterAsSource( parameters, QStringLiteral(
"INPUT" ), context ) );
136 QString uniqueIdFieldName( parameterAsString( parameters, QStringLiteral(
"UNIQUE_ID" ), context ) );
137 int uniqueIdFieldIdx = input->fields().indexFromName( uniqueIdFieldName );
138 if ( uniqueIdFieldIdx == -1 )
139 throw QgsProcessingException( QObject::tr(
"Missing field %1 in input layer" ).arg( uniqueIdFieldName ) );
141 const QgsField uniqueIdField = input->fields().at( uniqueIdFieldIdx );
144 fields.
append( uniqueIdField );
146 const std::unique_ptr<QgsFeatureSink> sink_output( parameterAsSink( parameters, QStringLiteral(
"OUTPUT" ), context, dest_output, fields, input->wkbType(), input->sourceCrs() ) );
148 const std::unique_ptr<QgsFeatureSink> sink_errors( parameterAsSink( parameters, QStringLiteral(
"ERRORS" ), context, dest_errors, fields,
Qgis::WkbType::Point, input->sourceCrs() ) );
159 QList<QgsGeometryCheckError *> checkErrors;
160 QStringList messages;
164 multiStepFeedback.setCurrentStep( 1 );
166 QMap<QString, QgsFeaturePool *> featurePools;
167 std::unique_ptr<QgsVectorLayer> inputLayer( input->materialize(
QgsFeatureRequest() ) );
169 featurePools.insert( inputLayer->id(), &featurePool );
171 multiStepFeedback.setCurrentStep( 2 );
173 check.collectErrors( featurePools, checkErrors, messages, feedback );
175 multiStepFeedback.setCurrentStep( 3 );
177 const double step { checkErrors.size() > 0 ? 100.0 / checkErrors.size() : 1 };
190 attrs << error->layerId()
191 << inputLayer->name()
192 << error->vidx().part
193 << error->vidx().ring
194 << error->vidx().vertex
195 << error->location().x()
196 << error->location().y()
197 << error->value().toString()
198 << inputLayer->getFeature( error->featureId() ).attribute( uniqueIdField.
name() );
203 throw QgsProcessingException( writeFeatureError( sink_output.get(), parameters, QStringLiteral(
"OUTPUT" ) ) );
207 throw QgsProcessingException( writeFeatureError( sink_errors.get(), parameters, QStringLiteral(
"ERRORS" ) ) );
210 feedback->
setProgress( 100.0 * step *
static_cast<double>( i ) );
221 outputs.insert( QStringLiteral(
"OUTPUT" ), dest_output );
222 outputs.insert( QStringLiteral(
"ERRORS" ), dest_errors );
@ VectorPoint
Vector point layers.
@ VectorPolygon
Vector polygon layers.
QFlags< ProcessingAlgorithmFlag > ProcessingAlgorithmFlags
Flags indicating how and when an algorithm operates and should be exposed to users.
@ NoThreading
Algorithm is not thread safe and cannot be run in a background thread, e.g. for algorithms which mani...
@ Advanced
Parameter is an advanced parameter which should be hidden from users by default.
Wraps a request for features to a vector layer (or directly its vector data provider).
@ FastInsert
Use faster inserts, at the cost of updating the passed features to reflect changes made at the provid...
The feature class encapsulates a single feature including its unique ID, geometry and a list of field...
void setAttributes(const QgsAttributes &attrs)
Sets the feature's attributes.
void setGeometry(const QgsGeometry &geometry)
Set the feature's geometry.
bool isCanceled() const
Tells whether the operation has been canceled already.
void setProgress(double progress)
Sets the current progress for the feedback object.
Encapsulate a field in an attribute table or data source.
Container of fields for a vector layer.
bool append(const QgsField &field, Qgis::FieldOrigin origin=Qgis::FieldOrigin::Provider, int originIndex=-1)
Appends a field.
Base configuration for geometry checks.
This represents an error reported by a geometry check.
A topology check for missing vertices.
static QgsGeometry fromPoint(const QgsPoint &point)
Creates a new geometry from a QgsPoint object.
Point geometry type, with support for z-dimension and m-values.
virtual Qgis::ProcessingAlgorithmFlags flags() const
Returns the flags indicating how and when the algorithm operates and should be exposed to users.
Contains information about the context in which a processing algorithm is executed.
Custom exception class for processing related exceptions.
Base class for providing feedback from a processing algorithm.
virtual void setProgressText(const QString &text)
Sets a progress report text string.
Processing feedback object for multi-step operations.
A feature sink output for processing algorithms.
An input feature source (such as vector layers) parameter for processing algorithms.
A vector layer or feature source field parameter for processing algorithms.
Encapsulates a QGIS project, including sets of map layers and their styles, layouts,...
static QgsProject * instance()
Returns the QgsProject singleton instance.
QgsCoordinateTransformContext transformContext
A feature pool based on a vector data provider.