QGIS API Documentation 3.40.0-Bratislava (b56115d8743)
Loading...
Searching...
No Matches
qgsalgorithmsplitwithlines.cpp
Go to the documentation of this file.
1/***************************************************************************
2 qgsalgorithmsplitwithlines.cpp
3 ---------------------
4 begin : April 2017
5 copyright : (C) 2017 by Nyall Dawson
6 email : nyall dot dawson at gmail dot com
7 ***************************************************************************/
8
9/***************************************************************************
10 * *
11 * This program is free software; you can redistribute it and/or modify *
12 * it under the terms of the GNU General Public License as published by *
13 * the Free Software Foundation; either version 2 of the License, or *
14 * (at your option) any later version. *
15 * *
16 ***************************************************************************/
17
19#include "qgsgeometryengine.h"
20#include "qgsvectorlayer.h"
21#include "qgsspatialindex.h"
22
24
25QString QgsSplitWithLinesAlgorithm::name() const
26{
27 return QStringLiteral( "splitwithlines" );
28}
29
30QString QgsSplitWithLinesAlgorithm::displayName() const
31{
32 return QObject::tr( "Split with lines" );
33}
34
35QStringList QgsSplitWithLinesAlgorithm::tags() const
36{
37 return QObject::tr( "split,cut,lines" ).split( ',' );
38}
39
40QString QgsSplitWithLinesAlgorithm::group() const
41{
42 return QObject::tr( "Vector overlay" );
43}
44
45QString QgsSplitWithLinesAlgorithm::groupId() const
46{
47 return QStringLiteral( "vectoroverlay" );
48}
49
50void QgsSplitWithLinesAlgorithm::initAlgorithm( const QVariantMap & )
51{
52 addParameter( new QgsProcessingParameterFeatureSource( QStringLiteral( "INPUT" ),
53 QObject::tr( "Input layer" ), QList< int >() << static_cast< int >( Qgis::ProcessingSourceType::VectorLine ) << static_cast< int >( Qgis::ProcessingSourceType::VectorPolygon ) ) );
54 addParameter( new QgsProcessingParameterFeatureSource( QStringLiteral( "LINES" ),
55 QObject::tr( "Split layer" ), QList< int >() << static_cast< int >( Qgis::ProcessingSourceType::VectorLine ) << static_cast< int >( Qgis::ProcessingSourceType::VectorPolygon ) ) );
56 addParameter( new QgsProcessingParameterFeatureSink( QStringLiteral( "OUTPUT" ), QObject::tr( "Split" ) ) );
57}
58
59QString QgsSplitWithLinesAlgorithm::shortHelpString() const
60{
61 return QObject::tr( "This algorithm splits the lines or polygons in one layer using the lines or polygon rings in another layer to define the breaking points. "
62 "Intersection between geometries in both layers are considered as split points." );
63}
64
65Qgis::ProcessingAlgorithmDocumentationFlags QgsSplitWithLinesAlgorithm::documentationFlags() const
66{
68}
69
70QgsSplitWithLinesAlgorithm *QgsSplitWithLinesAlgorithm::createInstance() const
71{
72 return new QgsSplitWithLinesAlgorithm();
73}
74
75Qgis::ProcessingAlgorithmFlags QgsSplitWithLinesAlgorithm::flags() const
76{
79 return f;
80}
81
82bool QgsSplitWithLinesAlgorithm::supportInPlaceEdit( const QgsMapLayer *l ) const
83{
84 const QgsVectorLayer *layer = qobject_cast< const QgsVectorLayer * >( l );
85 if ( !layer )
86 return false;
87
89 return false;
90
91 return true;
92}
93
94QVariantMap QgsSplitWithLinesAlgorithm::processAlgorithm( const QVariantMap &parameters, QgsProcessingContext &context, QgsProcessingFeedback *feedback )
95{
96 std::unique_ptr< QgsFeatureSource > source( parameterAsSource( parameters, QStringLiteral( "INPUT" ), context ) );
97 if ( !source )
98 throw QgsProcessingException( invalidSourceError( parameters, QStringLiteral( "INPUT" ) ) );
99
100 std::unique_ptr< QgsFeatureSource > linesSource( parameterAsSource( parameters, QStringLiteral( "LINES" ), context ) );
101 if ( !linesSource )
102 throw QgsProcessingException( invalidSourceError( parameters, QStringLiteral( "LINES" ) ) );
103
104 bool sameLayer = parameters.value( QStringLiteral( "INPUT" ) ) == parameters.value( QStringLiteral( "LINES" ) );
105
106 QString dest;
107 std::unique_ptr< QgsFeatureSink > sink( parameterAsSink( parameters, QStringLiteral( "OUTPUT" ), context, dest, source->fields(),
108 QgsWkbTypes::multiType( source->wkbType() ), source->sourceCrs(), QgsFeatureSink::RegeneratePrimaryKey ) );
109 if ( !sink )
110 throw QgsProcessingException( invalidSinkError( parameters, QStringLiteral( "OUTPUT" ) ) );
111
112 QgsFeatureRequest request;
113 request.setNoAttributes();
114 request.setDestinationCrs( source->sourceCrs(), context.transformContext() );
115
116 QgsFeatureIterator splitFeatures = linesSource->getFeatures( request );
117 QgsFeature aSplitFeature;
118
119 const QgsSpatialIndex splitFeaturesIndex( splitFeatures, feedback, QgsSpatialIndex::FlagStoreFeatureGeometries );
120
121 QgsFeature outFeat;
122 QgsFeatureIterator features = source->getFeatures();
123
124 double step = source->featureCount() > 0 ? 100.0 / source->featureCount() : 1;
125 int i = 0;
126 QgsFeature inFeatureA;
127 while ( features.nextFeature( inFeatureA ) )
128 {
129 i++;
130 if ( feedback->isCanceled() )
131 {
132 break;
133 }
134
135 if ( !inFeatureA.hasGeometry() )
136 {
137 if ( !sink->addFeature( inFeatureA, QgsFeatureSink::FastInsert ) )
138 throw QgsProcessingException( writeFeatureError( sink.get(), parameters, QStringLiteral( "OUTPUT" ) ) );
139 continue;
140 }
141
142 const QgsGeometry originalGeometry = inFeatureA.geometry();
143 outFeat.setAttributes( inFeatureA.attributes() );
144
145 QVector< QgsGeometry > inGeoms = originalGeometry.asGeometryCollection();
146
147 const QgsFeatureIds splitFeatureCandidates = qgis::listToSet( splitFeaturesIndex.intersects( originalGeometry.boundingBox() ) );
148 if ( !splitFeatureCandidates.empty() ) // has intersection of bounding boxes
149 {
150 QVector< QgsGeometry > splittingLines;
151
152 // use prepared geometries for faster intersection tests
153 std::unique_ptr< QgsGeometryEngine > originalGeometryEngine;
154
155 for ( QgsFeatureId splitFeatureCandidateId : splitFeatureCandidates )
156 {
157 // check if trying to self-intersect
158 if ( sameLayer && inFeatureA.id() == splitFeatureCandidateId )
159 continue;
160
161 const QgsGeometry splitFeatureCandidate = splitFeaturesIndex.geometry( splitFeatureCandidateId );
162 if ( !originalGeometryEngine )
163 {
164 originalGeometryEngine.reset( QgsGeometry::createGeometryEngine( originalGeometry.constGet() ) );
165 originalGeometryEngine->prepareGeometry();
166 }
167
168 if ( originalGeometryEngine->intersects( splitFeatureCandidate.constGet() ) )
169 {
170
171 QVector< QgsGeometry > splitGeomParts = splitFeatureCandidate.convertToType( Qgis::GeometryType::Line, true ).asGeometryCollection();
172 splittingLines.append( splitGeomParts );
173 }
174 }
175
176 if ( !splittingLines.empty() )
177 {
178 for ( const QgsGeometry &splitGeom : std::as_const( splittingLines ) )
179 {
180 QgsPointSequence splitterPList;
181 QVector< QgsGeometry > outGeoms;
182
183 // use prepared geometries for faster intersection tests
184 std::unique_ptr< QgsGeometryEngine > splitGeomEngine( QgsGeometry::createGeometryEngine( splitGeom.constGet() ) );
185 splitGeomEngine->prepareGeometry();
186 while ( !inGeoms.empty() )
187 {
188 if ( feedback->isCanceled() )
189 {
190 break;
191 }
192
193 QgsGeometry inGeom = inGeoms.takeFirst();
194 if ( inGeom.isNull() )
195 continue;
196
197 if ( splitGeomEngine->intersects( inGeom.constGet() ) )
198 {
199 QgsGeometry before = inGeom;
200 if ( splitterPList.empty() )
201 {
202 const QgsCoordinateSequence sequence = splitGeom.constGet()->coordinateSequence();
203 for ( const QgsRingSequence &part : sequence )
204 {
205 for ( const QgsPointSequence &ring : part )
206 {
207 for ( const QgsPoint &pt : ring )
208 {
209 splitterPList << pt;
210 }
211 }
212 }
213 }
214
215 QVector< QgsGeometry > newGeometries;
216 QgsPointSequence topologyTestPoints;
217 Qgis::GeometryOperationResult result = inGeom.splitGeometry( splitterPList, newGeometries, false, topologyTestPoints, true );
218
219 // splitGeometry: If there are several intersections
220 // between geometry and splitLine, only the first one is considered.
222 {
223 // sometimes the resultant geometry has changed from the input, but only because of numerical precision issues.
224 // and is effectively indistinguishable from the input. By testing the Hausdorff distance is less than this threshold
225 // we are checking that the maximum "change" between the result and the input is actually significant enough to be meaningful...
226 if ( inGeom.hausdorffDistance( before ) < 1e-12 )
227 {
228 // effectively no change!!
229 outGeoms.append( inGeom );
230 }
231 else
232 {
233 outGeoms.append( inGeom );
234 outGeoms.append( newGeometries );
235 }
236 }
237 else
238 {
239 outGeoms.append( inGeom );
240 }
241 }
242 else
243 {
244 outGeoms.append( inGeom );
245 }
246
247 }
248 inGeoms = outGeoms;
249 }
250 }
251 }
252
253 QVector< QgsGeometry > parts;
254 for ( const QgsGeometry &aGeom : std::as_const( inGeoms ) )
255 {
256 if ( feedback->isCanceled() )
257 {
258 break;
259 }
260
261 bool passed = true;
262 if ( QgsWkbTypes::geometryType( aGeom.wkbType() ) == Qgis::GeometryType::Line )
263 {
264 int numPoints = aGeom.constGet()->nCoordinates();
265
266 if ( numPoints <= 2 )
267 {
268 if ( numPoints == 2 )
269 passed = !static_cast< const QgsCurve * >( aGeom.constGet() )->isClosed(); // tests if vertex 0 = vertex 1
270 else
271 passed = false; // sometimes splitting results in lines of zero length
272 }
273 }
274
275 if ( passed )
276 parts.append( aGeom );
277 }
278
279 for ( const QgsGeometry &g : parts )
280 {
281 outFeat.setGeometry( g );
282 if ( !sink->addFeature( outFeat, QgsFeatureSink::FastInsert ) )
283 throw QgsProcessingException( writeFeatureError( sink.get(), parameters, QStringLiteral( "OUTPUT" ) ) );
284 }
285
286 feedback->setProgress( i * step );
287 }
288
289 QVariantMap outputs;
290 outputs.insert( QStringLiteral( "OUTPUT" ), dest );
291 return outputs;
292}
293
294
295
297
298
@ VectorPolygon
Vector polygon layers.
@ VectorLine
Vector line layers.
GeometryOperationResult
Success or failure of a geometry operation.
Definition qgis.h:1889
@ Success
Operation succeeded.
@ Polygon
Polygons.
@ RegeneratesPrimaryKey
Algorithm always drops any existing primary keys or FID values and regenerates them in outputs.
QFlags< ProcessingAlgorithmFlag > ProcessingAlgorithmFlags
Flags indicating how and when an algorithm operates and should be exposed to users.
Definition qgis.h:3339
QFlags< ProcessingAlgorithmDocumentationFlag > ProcessingAlgorithmDocumentationFlags
Flags describing algorithm behavior for documentation purposes.
Definition qgis.h:3359
@ SupportsInPlaceEdits
Algorithm supports in-place editing.
virtual QgsCoordinateSequence coordinateSequence() const =0
Retrieves the sequence of geometries, rings and nodes.
Abstract base class for curved geometry type.
Definition qgscurve.h:35
Wrapper for iterator of features from vector data provider or vector layer.
bool nextFeature(QgsFeature &f)
Fetch next feature and stores in f, returns true on success.
This class wraps a request for features to a vector layer (or directly its vector data provider).
QgsFeatureRequest & setDestinationCrs(const QgsCoordinateReferenceSystem &crs, const QgsCoordinateTransformContext &context)
Sets the destination crs for feature's geometries.
QgsFeatureRequest & setNoAttributes()
Set that no attributes will be fetched.
@ FastInsert
Use faster inserts, at the cost of updating the passed features to reflect changes made at the provid...
@ RegeneratePrimaryKey
This flag indicates, that a primary key field cannot be guaranteed to be unique and the sink should i...
The feature class encapsulates a single feature including its unique ID, geometry and a list of field...
Definition qgsfeature.h:58
QgsAttributes attributes
Definition qgsfeature.h:67
QgsFeatureId id
Definition qgsfeature.h:66
void setAttributes(const QgsAttributes &attrs)
Sets the feature's attributes.
QgsGeometry geometry
Definition qgsfeature.h:69
bool hasGeometry() const
Returns true if the feature has an associated geometry.
void setGeometry(const QgsGeometry &geometry)
Set the feature's geometry.
bool isCanceled() const
Tells whether the operation has been canceled already.
Definition qgsfeedback.h:53
void setProgress(double progress)
Sets the current progress for the feedback object.
Definition qgsfeedback.h:61
A geometry is the spatial representation of a feature.
QgsGeometry convertToType(Qgis::GeometryType destType, bool destMultipart=false) const
Try to convert the geometry to the requested type.
QVector< QgsGeometry > asGeometryCollection() const
Returns contents of the geometry as a list of geometries.
const QgsAbstractGeometry * constGet() const
Returns a non-modifiable (const) reference to the underlying abstract geometry primitive.
double hausdorffDistance(const QgsGeometry &geom) const
Returns the Hausdorff distance between this geometry and geom.
QgsRectangle boundingBox() const
Returns the bounding box of the geometry.
static QgsGeometryEngine * createGeometryEngine(const QgsAbstractGeometry *geometry, double precision=0.0)
Creates and returns a new geometry engine representing the specified geometry using precision on a gr...
Q_DECL_DEPRECATED Qgis::GeometryOperationResult splitGeometry(const QVector< QgsPointXY > &splitLine, QVector< QgsGeometry > &newGeometries, bool topological, QVector< QgsPointXY > &topologyTestPoints, bool splitFeature=true)
Splits this geometry according to a given line.
Base class for all map layer types.
Definition qgsmaplayer.h:76
Point geometry type, with support for z-dimension and m-values.
Definition qgspoint.h:49
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.
QgsCoordinateTransformContext transformContext() const
Returns the coordinate transform context.
Custom exception class for processing related exceptions.
Base class for providing feedback from a processing algorithm.
A feature sink output for processing algorithms.
An input feature source (such as vector layers) parameter for processing algorithms.
A spatial index for QgsFeature objects.
@ FlagStoreFeatureGeometries
Indicates that the spatial index should also store feature geometries. This requires more memory,...
Represents a vector layer which manages a vector based data sets.
Q_INVOKABLE Qgis::GeometryType geometryType() const
Returns point, line or polygon.
static Qgis::GeometryType geometryType(Qgis::WkbType type)
Returns the geometry type for a WKB type, e.g., both MultiPolygon and CurvePolygon would have a Polyg...
static Qgis::WkbType multiType(Qgis::WkbType type)
Returns the multi type for a WKB type.
QVector< QgsRingSequence > QgsCoordinateSequence
QVector< QgsPointSequence > QgsRingSequence
QVector< QgsPoint > QgsPointSequence
QSet< QgsFeatureId > QgsFeatureIds
qint64 QgsFeatureId
64 bit feature ids negative numbers are used for uncommitted/newly added features