QGIS API Documentation 3.38.0-Grenoble (exported)
Loading...
Searching...
No Matches
qgsvectorlayereditutils.cpp
Go to the documentation of this file.
1/***************************************************************************
2 qgsvectorlayereditutils.cpp
3 ---------------------
4 begin : Dezember 2012
5 copyright : (C) 2012 by Martin Dobias
6 email : wonder dot sk at gmail 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 ***************************************************************************/
16
19#include "qgsfeatureiterator.h"
21#include "qgslinestring.h"
22#include "qgslogger.h"
23#include "qgspoint.h"
24#include "qgis.h"
25#include "qgswkbtypes.h"
26#include "qgsvectorlayerutils.h"
27#include "qgsvectorlayer.h"
28#include "qgsgeometryoptions.h"
29#include "qgsabstractgeometry.h"
32
33#include <limits>
34
35
40
41bool QgsVectorLayerEditUtils::insertVertex( double x, double y, QgsFeatureId atFeatureId, int beforeVertex )
42{
43 if ( !mLayer->isSpatial() )
44 return false;
45
46 QgsFeature f;
47 if ( !mLayer->getFeatures( QgsFeatureRequest().setFilterFid( atFeatureId ).setNoAttributes() ).nextFeature( f ) || !f.hasGeometry() )
48 return false; // geometry not found
49
50 QgsGeometry geometry = f.geometry();
51
52 geometry.insertVertex( x, y, beforeVertex );
53
54 mLayer->changeGeometry( atFeatureId, geometry );
55 return true;
56}
57
58bool QgsVectorLayerEditUtils::insertVertex( const QgsPoint &point, QgsFeatureId atFeatureId, int beforeVertex )
59{
60 if ( !mLayer->isSpatial() )
61 return false;
62
63 QgsFeature f;
64 if ( !mLayer->getFeatures( QgsFeatureRequest().setFilterFid( atFeatureId ).setNoAttributes() ).nextFeature( f ) || !f.hasGeometry() )
65 return false; // geometry not found
66
67 QgsGeometry geometry = f.geometry();
68
69 geometry.insertVertex( point, beforeVertex );
70
71 mLayer->changeGeometry( atFeatureId, geometry );
72 return true;
73}
74
75bool QgsVectorLayerEditUtils::moveVertex( double x, double y, QgsFeatureId atFeatureId, int atVertex )
76{
77 QgsPoint p( x, y );
78 return moveVertex( p, atFeatureId, atVertex );
79}
80
81bool QgsVectorLayerEditUtils::moveVertex( const QgsPoint &p, QgsFeatureId atFeatureId, int atVertex )
82{
83 if ( !mLayer->isSpatial() )
84 return false;
85
86 QgsFeature f;
87 if ( !mLayer->getFeatures( QgsFeatureRequest().setFilterFid( atFeatureId ).setNoAttributes() ).nextFeature( f ) || !f.hasGeometry() )
88 return false; // geometry not found
89
90 QgsGeometry geometry = f.geometry();
91
92 // If original point is not 3D but destination yes, check if it can be promoted
93 if ( p.is3D() && !geometry.constGet()->is3D() && QgsWkbTypes::hasZ( mLayer->wkbType() ) )
94 {
96 return false;
97 }
98
99 // If original point has not M-value but destination yes, check if it can be promoted
100 if ( p.isMeasure() && !geometry.constGet()->isMeasure() && QgsWkbTypes::hasM( mLayer->wkbType() ) )
101 {
103 return false;
104 }
105
106 if ( !geometry.moveVertex( p, atVertex ) )
107 return false;
108
109 return mLayer->changeGeometry( atFeatureId, geometry );
110}
111
113{
114 if ( !mLayer->isSpatial() )
116
117 QgsFeature f;
118 if ( !mLayer->getFeatures( QgsFeatureRequest().setFilterFid( featureId ).setNoAttributes() ).nextFeature( f ) || !f.hasGeometry() )
119 return Qgis::VectorEditResult::FetchFeatureFailed; // geometry not found
120
121 QgsGeometry geometry = f.geometry();
122
123 if ( !geometry.deleteVertex( vertex ) )
125
126 if ( geometry.constGet() && geometry.constGet()->nCoordinates() == 0 )
127 {
128 //last vertex deleted, set geometry to null
129 geometry.set( nullptr );
130 }
131
132 mLayer->changeGeometry( featureId, geometry );
134}
135
136
137static
138Qgis::GeometryOperationResult staticAddRing( QgsVectorLayer *layer, std::unique_ptr< QgsCurve > &ring, const QgsFeatureIds &targetFeatureIds, QgsFeatureIds *modifiedFeatureIds, bool firstOne = true )
139{
140
141 if ( !layer || !layer->isSpatial() )
142 {
144 }
145
146 if ( !ring )
147 {
149 }
150
151 if ( !ring->isClosed() )
152 {
154 }
155
156 if ( !layer->isValid() || !layer->editBuffer() || !layer->dataProvider() )
157 {
159 }
160
161 Qgis::GeometryOperationResult addRingReturnCode = Qgis::GeometryOperationResult::AddRingNotInExistingFeature; //default: return code for 'ring not inserted'
162 QgsFeature f;
163
165 if ( !targetFeatureIds.isEmpty() )
166 {
167 //check only specified features
168 fit = layer->getFeatures( QgsFeatureRequest().setFilterFids( targetFeatureIds ) );
169 }
170 else
171 {
172 //check all intersecting features
173 QgsRectangle bBox = ring->boundingBox();
174 fit = layer->getFeatures( QgsFeatureRequest().setFilterRect( bBox ).setFlags( Qgis::FeatureRequestFlag::ExactIntersect ) );
175 }
176
177 //find first valid feature we can add the ring to
178 while ( fit.nextFeature( f ) )
179 {
180 if ( !f.hasGeometry() )
181 continue;
182
183 //add ring takes ownership of ring, and deletes it if there's an error
184 QgsGeometry g = f.geometry();
185
186 if ( ring->orientation() != g.polygonOrientation() )
187 {
188 addRingReturnCode = g.addRing( static_cast< QgsCurve * >( ring->clone() ) );
189 }
190 else
191 {
192 addRingReturnCode = g.addRing( static_cast< QgsCurve * >( ring->reversed() ) );
193 }
194 if ( addRingReturnCode == Qgis::GeometryOperationResult::Success )
195 {
196 layer->changeGeometry( f.id(), g );
197 if ( modifiedFeatureIds )
198 {
199 modifiedFeatureIds->insert( f.id() );
200 if ( firstOne )
201 {
202 break;
203 }
204 }
205
206 }
207 }
208
209 return addRingReturnCode;
210}
211
212Qgis::GeometryOperationResult QgsVectorLayerEditUtils::addRing( const QVector<QgsPointXY> &ring, const QgsFeatureIds &targetFeatureIds, QgsFeatureId *modifiedFeatureId )
213{
215 for ( QVector<QgsPointXY>::const_iterator it = ring.constBegin(); it != ring.constEnd(); ++it )
216 {
217 l << QgsPoint( *it );
218 }
219 return addRing( l, targetFeatureIds, modifiedFeatureId );
220}
221
223{
224 QgsLineString *ringLine = new QgsLineString( ring );
225 return addRing( ringLine, targetFeatureIds, modifiedFeatureId );
226}
227
229{
230 std::unique_ptr<QgsCurve> uniquePtrRing( ring );
231 if ( modifiedFeatureId )
232 {
233 QgsFeatureIds *modifiedFeatureIds = new QgsFeatureIds;
234 Qgis::GeometryOperationResult result = staticAddRing( mLayer, uniquePtrRing, targetFeatureIds, modifiedFeatureIds, true );
235 *modifiedFeatureId = *modifiedFeatureIds->begin();
236 return result;
237 }
238 return staticAddRing( mLayer, uniquePtrRing, targetFeatureIds, nullptr, true );
239}
240
242{
243
244 std::unique_ptr<QgsCurve> uniquePtrRing( ring );
245 return staticAddRing( mLayer, uniquePtrRing, targetFeatureIds, modifiedFeatureIds, false );
246}
247
248
249
251{
253 for ( QVector<QgsPointXY>::const_iterator it = points.constBegin(); it != points.constEnd(); ++it )
254 {
255 l << QgsPoint( *it );
256 }
257 return addPart( l, featureId );
258}
259
261{
262 if ( !mLayer->isSpatial() )
264
265 QgsGeometry geometry;
266 bool firstPart = false;
267 QgsFeature f;
268 if ( !mLayer->getFeatures( QgsFeatureRequest().setFilterFid( featureId ).setNoAttributes() ).nextFeature( f ) )
270
271 if ( !f.hasGeometry() )
272 {
273 //no existing geometry, so adding first part to null geometry
274 firstPart = true;
275 }
276 else
277 {
278 geometry = f.geometry();
279 }
280
281 Qgis::GeometryOperationResult errorCode = geometry.addPartV2( points, mLayer->wkbType() );
283 {
284 if ( firstPart && QgsWkbTypes::isSingleType( mLayer->wkbType() )
286 {
287 //convert back to single part if required by layer
288 geometry.convertToSingleType();
289 }
290 mLayer->changeGeometry( featureId, geometry );
291 }
292 return errorCode;
293}
294
296{
297
298 if ( !mLayer->isSpatial() )
300
301 QgsGeometry geometry;
302 bool firstPart = false;
303 QgsFeature f;
304 if ( !mLayer->getFeatures( QgsFeatureRequest().setFilterFid( featureId ).setNoAttributes() ).nextFeature( f ) )
306
307 if ( !f.hasGeometry() )
308 {
309 //no existing geometry, so adding first part to null geometry
310 firstPart = true;
311 }
312 else
313 {
314 geometry = f.geometry();
315 if ( ring->orientation() != geometry.polygonOrientation() )
316 {
317 ring = ring->reversed();
318 }
319 }
320 Qgis::GeometryOperationResult errorCode = geometry.addPartV2( ring, mLayer->wkbType() );
321
323 {
324 if ( firstPart && QgsWkbTypes::isSingleType( mLayer->wkbType() )
326 {
327 //convert back to single part if required by layer
328 geometry.convertToSingleType();
329 }
330 mLayer->changeGeometry( featureId, geometry );
331 }
332 return errorCode;
333}
334
335// TODO QGIS 4.0 -- this should return Qgis::GeometryOperationResult
336int QgsVectorLayerEditUtils::translateFeature( QgsFeatureId featureId, double dx, double dy )
337{
338 if ( !mLayer->isSpatial() )
339 return 1;
340
341 QgsFeature f;
342 if ( !mLayer->getFeatures( QgsFeatureRequest().setFilterFid( featureId ).setNoAttributes() ).nextFeature( f ) || !f.hasGeometry() )
343 return 1; //geometry not found
344
345 QgsGeometry geometry = f.geometry();
346
347 Qgis::GeometryOperationResult errorCode = geometry.translate( dx, dy );
349 {
350 mLayer->changeGeometry( featureId, geometry );
351 }
352 return errorCode == Qgis::GeometryOperationResult::Success ? 0 : 1;
353}
354
355Qgis::GeometryOperationResult QgsVectorLayerEditUtils::splitFeatures( const QVector<QgsPointXY> &splitLine, bool topologicalEditing )
356{
357
359 for ( QVector<QgsPointXY>::const_iterator it = splitLine.constBegin(); it != splitLine.constEnd(); ++it )
360 {
361 l << QgsPoint( *it );
362 }
363 return splitFeatures( l, topologicalEditing );
364}
365
367{
368 QgsLineString lineString( splitLine );
369 QgsPointSequence topologyTestPoints;
370 bool preserveCircular = false;
371 return splitFeatures( &lineString, topologyTestPoints, preserveCircular, topologicalEditing );
372}
373
374Qgis::GeometryOperationResult QgsVectorLayerEditUtils::splitFeatures( const QgsCurve *curve, QgsPointSequence &topologyTestPoints, bool preserveCircular, bool topologicalEditing )
375{
376 if ( !mLayer->isSpatial() )
378
379 QgsRectangle bBox; //bounding box of the split line
381 Qgis::GeometryOperationResult splitFunctionReturn; //return code of QgsGeometry::splitGeometry
382 int numberOfSplitFeatures = 0;
383
384 QgsFeatureIterator features;
385 const QgsFeatureIds selectedIds = mLayer->selectedFeatureIds();
386
387 // deactivate preserving circular if the curve contains only straight segments to avoid transforming Polygon to CurvePolygon
388 preserveCircular &= curve->hasCurvedSegments();
389
390 if ( !selectedIds.isEmpty() ) //consider only the selected features if there is a selection
391 {
392 features = mLayer->getSelectedFeatures();
393 }
394 else //else consider all the feature that intersect the bounding box of the split line
395 {
396
397 bBox = curve->boundingBox();
398
399 if ( bBox.isEmpty() )
400 {
401 //if the bbox is a line, try to make a square out of it
402 if ( bBox.width() == 0.0 && bBox.height() > 0 )
403 {
404 bBox.setXMinimum( bBox.xMinimum() - bBox.height() / 2 );
405 bBox.setXMaximum( bBox.xMaximum() + bBox.height() / 2 );
406 }
407 else if ( bBox.height() == 0.0 && bBox.width() > 0 )
408 {
409 bBox.setYMinimum( bBox.yMinimum() - bBox.width() / 2 );
410 bBox.setYMaximum( bBox.yMaximum() + bBox.width() / 2 );
411 }
412 else
413 {
414 //If we have a single point, we still create a non-null box
415 double bufferDistance = 0.000001;
416 if ( mLayer->crs().isGeographic() )
417 bufferDistance = 0.00000001;
418 bBox.setXMinimum( bBox.xMinimum() - bufferDistance );
419 bBox.setXMaximum( bBox.xMaximum() + bufferDistance );
420 bBox.setYMinimum( bBox.yMinimum() - bufferDistance );
421 bBox.setYMaximum( bBox.yMaximum() + bufferDistance );
422 }
423 }
424
425 features = mLayer->getFeatures( QgsFeatureRequest().setFilterRect( bBox ).setFlags( Qgis::FeatureRequestFlag::ExactIntersect ) );
426 }
427
429
430 const int fieldCount = mLayer->fields().count();
431
432 QgsFeature feat;
433 while ( features.nextFeature( feat ) )
434 {
435 if ( !feat.hasGeometry() )
436 {
437 continue;
438 }
439 QVector<QgsGeometry> newGeometries;
440 QgsPointSequence featureTopologyTestPoints;
441 const QgsGeometry originalGeom = feat.geometry();
442 QgsGeometry featureGeom = originalGeom;
443 splitFunctionReturn = featureGeom.splitGeometry( curve, newGeometries, preserveCircular, topologicalEditing, featureTopologyTestPoints );
444 topologyTestPoints.append( featureTopologyTestPoints );
445 if ( splitFunctionReturn == Qgis::GeometryOperationResult::Success )
446 {
447 //change this geometry
448 mLayer->changeGeometry( feat.id(), featureGeom );
449
450 //update any attributes for original feature which are set to GeometryRatio split policy
451 QgsAttributeMap attributeMap;
452 for ( int fieldIdx = 0; fieldIdx < fieldCount; ++fieldIdx )
453 {
454 const QgsField field = mLayer->fields().at( fieldIdx );
455 switch ( field.splitPolicy() )
456 {
460 break;
461
463 {
464 if ( field.isNumeric() )
465 {
466 const double originalValue = feat.attribute( fieldIdx ).toDouble();
467
468 double originalSize = 0;
469
470 switch ( originalGeom.type() )
471 {
475 originalSize = 0;
476 break;
478 originalSize = originalGeom.length();
479 break;
481 originalSize = originalGeom.area();
482 break;
483 }
484
485 double newSize = 0;
486 switch ( featureGeom.type() )
487 {
491 newSize = 0;
492 break;
494 newSize = featureGeom.length();
495 break;
497 newSize = featureGeom.area();
498 break;
499 }
500
501 attributeMap.insert( fieldIdx, originalSize > 0 ? ( originalValue * newSize / originalSize ) : originalValue );
502 }
503 break;
504 }
505 }
506 }
507
508 if ( !attributeMap.isEmpty() )
509 {
510 mLayer->changeAttributeValues( feat.id(), attributeMap );
511 }
512
513 //insert new features
514 for ( const QgsGeometry &geom : std::as_const( newGeometries ) )
515 {
516 QgsAttributeMap attributeMap;
517 for ( int fieldIdx = 0; fieldIdx < fieldCount; ++fieldIdx )
518 {
519 const QgsField field = mLayer->fields().at( fieldIdx );
520 // respect field split policy
521 switch ( field.splitPolicy() )
522 {
524 //do nothing - default values ​​are determined
525 break;
526
528 attributeMap.insert( fieldIdx, feat.attribute( fieldIdx ) );
529 break;
530
532 {
533 if ( !field.isNumeric() )
534 {
535 attributeMap.insert( fieldIdx, feat.attribute( fieldIdx ) );
536 }
537 else
538 {
539 const double originalValue = feat.attribute( fieldIdx ).toDouble();
540
541 double originalSize = 0;
542
543 switch ( originalGeom.type() )
544 {
548 originalSize = 0;
549 break;
551 originalSize = originalGeom.length();
552 break;
554 originalSize = originalGeom.area();
555 break;
556 }
557
558 double newSize = 0;
559 switch ( geom.type() )
560 {
564 newSize = 0;
565 break;
567 newSize = geom.length();
568 break;
570 newSize = geom.area();
571 break;
572 }
573
574 attributeMap.insert( fieldIdx, originalSize > 0 ? ( originalValue * newSize / originalSize ) : originalValue );
575 }
576 break;
577 }
578
580 attributeMap.insert( fieldIdx, QgsUnsetAttributeValue() );
581 break;
582 }
583 }
584
585 featuresDataToAdd << QgsVectorLayerUtils::QgsFeatureData( geom, attributeMap );
586 }
587
588 if ( topologicalEditing )
589 {
590 QgsPointSequence::const_iterator topol_it = featureTopologyTestPoints.constBegin();
591 for ( ; topol_it != featureTopologyTestPoints.constEnd(); ++topol_it )
592 {
593 addTopologicalPoints( *topol_it );
594 }
595 }
596 ++numberOfSplitFeatures;
597 }
598 else if ( splitFunctionReturn != Qgis::GeometryOperationResult::Success && splitFunctionReturn != Qgis::GeometryOperationResult::NothingHappened ) // i.e. no split but no error occurred
599 {
600 returnCode = splitFunctionReturn;
601 }
602 }
603
604 if ( !featuresDataToAdd.isEmpty() )
605 {
606 // finally create and add all bits of geometries cut off the original geometries
607 // (this is much faster than creating features one by one)
608 QgsFeatureList featuresListToAdd = QgsVectorLayerUtils::createFeatures( mLayer, featuresDataToAdd );
609 mLayer->addFeatures( featuresListToAdd );
610 }
611
612 if ( numberOfSplitFeatures == 0 )
613 {
615 }
616
617 return returnCode;
618}
619
620Qgis::GeometryOperationResult QgsVectorLayerEditUtils::splitParts( const QVector<QgsPointXY> &splitLine, bool topologicalEditing )
621{
623 for ( QVector<QgsPointXY>::const_iterator it = splitLine.constBegin(); it != splitLine.constEnd(); ++it )
624 {
625 l << QgsPoint( *it );
626 }
627 return splitParts( l, topologicalEditing );
628}
629
631{
632 if ( !mLayer->isSpatial() )
634
635 double xMin, yMin, xMax, yMax;
636 QgsRectangle bBox; //bounding box of the split line
637 int numberOfSplitParts = 0;
638
640
641 if ( mLayer->selectedFeatureCount() > 0 ) //consider only the selected features if there is a selection
642 {
643 fit = mLayer->getSelectedFeatures();
644 }
645 else //else consider all the feature that intersect the bounding box of the split line
646 {
647 if ( boundingBoxFromPointList( splitLine, xMin, yMin, xMax, yMax ) )
648 {
649 bBox.setXMinimum( xMin );
650 bBox.setYMinimum( yMin );
651 bBox.setXMaximum( xMax );
652 bBox.setYMaximum( yMax );
653 }
654 else
655 {
657 }
658
659 if ( bBox.isEmpty() )
660 {
661 //if the bbox is a line, try to make a square out of it
662 if ( bBox.width() == 0.0 && bBox.height() > 0 )
663 {
664 bBox.setXMinimum( bBox.xMinimum() - bBox.height() / 2 );
665 bBox.setXMaximum( bBox.xMaximum() + bBox.height() / 2 );
666 }
667 else if ( bBox.height() == 0.0 && bBox.width() > 0 )
668 {
669 bBox.setYMinimum( bBox.yMinimum() - bBox.width() / 2 );
670 bBox.setYMaximum( bBox.yMaximum() + bBox.width() / 2 );
671 }
672 else
673 {
674 //If we have a single point, we still create a non-null box
675 double bufferDistance = 0.000001;
676 if ( mLayer->crs().isGeographic() )
677 bufferDistance = 0.00000001;
678 bBox.setXMinimum( bBox.xMinimum() - bufferDistance );
679 bBox.setXMaximum( bBox.xMaximum() + bufferDistance );
680 bBox.setYMinimum( bBox.yMinimum() - bufferDistance );
681 bBox.setYMaximum( bBox.yMaximum() + bufferDistance );
682 }
683 }
684
685 fit = mLayer->getFeatures( QgsFeatureRequest().setFilterRect( bBox ).setFlags( Qgis::FeatureRequestFlag::ExactIntersect ) );
686 }
687
688 QgsFeature feat;
689 while ( fit.nextFeature( feat ) )
690 {
691 QgsGeometry featureGeom = feat.geometry();
692
693 const QVector<QgsGeometry> geomCollection = featureGeom.asGeometryCollection();
694 QVector<QgsGeometry> resultCollection;
695 QgsPointSequence topologyTestPoints;
696 for ( QgsGeometry part : geomCollection )
697 {
698 QVector<QgsGeometry> newGeometries;
699 QgsPointSequence partTopologyTestPoints;
700
701 const Qgis::GeometryOperationResult splitFunctionReturn = part.splitGeometry( splitLine, newGeometries, topologicalEditing, partTopologyTestPoints, false );
702
703 if ( splitFunctionReturn == Qgis::GeometryOperationResult::Success && !newGeometries.isEmpty() )
704 {
705 for ( int i = 0; i < newGeometries.size(); ++i )
706 {
707 resultCollection.append( newGeometries.at( i ).asGeometryCollection() );
708 }
709
710 topologyTestPoints.append( partTopologyTestPoints );
711
712 ++numberOfSplitParts;
713 }
714 // Note: For multilinestring layers, when the split line does not intersect the feature part,
715 // QgsGeometry::splitGeometry returns InvalidBaseGeometry instead of NothingHappened
716 else if ( splitFunctionReturn == Qgis::GeometryOperationResult::NothingHappened ||
718 {
719 // Add part as is
720 resultCollection.append( part );
721 }
722 else if ( splitFunctionReturn != Qgis::GeometryOperationResult::Success )
723 {
724 return splitFunctionReturn;
725 }
726 }
727
728 QgsGeometry newGeom = QgsGeometry::collectGeometry( resultCollection );
729 mLayer->changeGeometry( feat.id(), newGeom ) ;
730
731 if ( topologicalEditing )
732 {
733 QgsPointSequence::const_iterator topol_it = topologyTestPoints.constBegin();
734 for ( ; topol_it != topologyTestPoints.constEnd(); ++topol_it )
735 {
736 addTopologicalPoints( *topol_it );
737 }
738 }
739
740 }
741 if ( numberOfSplitParts == 0 && mLayer->selectedFeatureCount() > 0 )
742 {
743 //There is a selection but no feature has been split.
744 //Maybe user forgot that only the selected features are split
746 }
747
749}
750
751
753{
754 if ( !mLayer->isSpatial() )
755 return 1;
756
757 if ( geom.isNull() )
758 {
759 return 1;
760 }
761
762 bool pointsAdded = false;
763
765 while ( it != geom.vertices_end() )
766 {
767 if ( addTopologicalPoints( *it ) == 0 )
768 {
769 pointsAdded = true;
770 }
771 ++it;
772 }
773
774 return pointsAdded ? 0 : 2;
775}
776
778{
779 if ( !mLayer->isSpatial() )
780 return 1;
781
782 double segmentSearchEpsilon = mLayer->crs().isGeographic() ? 1e-12 : 1e-8;
783
784 //work with a tolerance because coordinate projection may introduce some rounding
785 double threshold = mLayer->geometryOptions()->geometryPrecision();
786
787 if ( qgsDoubleNear( threshold, 0.0 ) )
788 {
789 threshold = 1e-8;
790
791 if ( mLayer->crs().mapUnits() == Qgis::DistanceUnit::Meters )
792 {
793 threshold = 0.001;
794 }
795 else if ( mLayer->crs().mapUnits() == Qgis::DistanceUnit::Feet )
796 {
797 threshold = 0.0001;
798 }
799 }
800
801 QgsRectangle searchRect( p, p, false );
802 searchRect.grow( threshold );
803
804 QgsFeature f;
806 .setFilterRect( searchRect )
808 .setNoAttributes() );
809
810 bool pointsAdded = false;
811 while ( fit.nextFeature( f ) )
812 {
813 QgsGeometry geom = f.geometry();
814 if ( geom.addTopologicalPoint( p, threshold, segmentSearchEpsilon ) )
815 {
816 pointsAdded = true;
817 mLayer->changeGeometry( f.id(), geom );
818 }
819 }
820
821 return pointsAdded ? 0 : 2;
822}
823
825{
826 if ( !mLayer->isSpatial() )
827 return 1;
828
829 if ( ps.isEmpty() )
830 {
831 return 1;
832 }
833
834 bool pointsAdded = false;
835
836 QgsPointSequence::const_iterator it = ps.constBegin();
837 while ( it != ps.constEnd() )
838 {
839 if ( addTopologicalPoints( *it ) == 0 )
840 {
841 pointsAdded = true;
842 }
843 ++it;
844 }
845
846 return pointsAdded ? 0 : 2;
847}
848
853
854bool QgsVectorLayerEditUtils::mergeFeatures( const QgsFeatureId &targetFeatureId, const QgsFeatureIds &mergeFeatureIds, const QgsAttributes &mergeAttributes, const QgsGeometry &unionGeometry, QString &errorMessage )
855{
856 errorMessage.clear();
857
858 if ( mergeFeatureIds.isEmpty() )
859 {
860 errorMessage = QObject::tr( "List of features to merge is empty" );
861 return false;
862 }
863
864 QgsAttributeMap newAttributes;
865 for ( int i = 0; i < mergeAttributes.count(); ++i )
866 {
867 QVariant val = mergeAttributes.at( i );
868
869 bool isDefaultValue = mLayer->fields().fieldOrigin( i ) == Qgis::FieldOrigin::Provider &&
870 mLayer->dataProvider() &&
871 mLayer->dataProvider()->defaultValueClause( mLayer->fields().fieldOriginIndex( i ) ) == val;
872
873 // convert to destination data type
874 QString errorMessageConvertCompatible;
875 if ( !isDefaultValue && !mLayer->fields().at( i ).convertCompatible( val, &errorMessageConvertCompatible ) )
876 {
877 if ( errorMessage.isEmpty() )
878 errorMessage = QObject::tr( "Could not store value '%1' in field of type %2: %3" ).arg( mergeAttributes.at( i ).toString(), mLayer->fields().at( i ).typeName(), errorMessageConvertCompatible );
879 }
880 newAttributes[ i ] = val;
881 }
882
883 mLayer->beginEditCommand( QObject::tr( "Merged features" ) );
884
885 // Delete other features but the target feature
886 QgsFeatureIds::const_iterator feature_it = mergeFeatureIds.constBegin();
887 for ( ; feature_it != mergeFeatureIds.constEnd(); ++feature_it )
888 {
889 if ( *feature_it != targetFeatureId )
890 mLayer->deleteFeature( *feature_it );
891 }
892
893 // Modify target feature or create a new one if invalid
894 QgsGeometry mergeGeometry = unionGeometry;
895 if ( targetFeatureId == FID_NULL )
896 {
897 QgsFeature mergeFeature = QgsVectorLayerUtils::createFeature( mLayer, mergeGeometry, newAttributes );
898 mLayer->addFeature( mergeFeature );
899 }
900 else
901 {
902 mLayer->changeGeometry( targetFeatureId, mergeGeometry );
903 mLayer->changeAttributeValues( targetFeatureId, newAttributes );
904 }
905
906 mLayer->endEditCommand();
907
908 mLayer->triggerRepaint();
909
910 return true;
911}
912
913bool QgsVectorLayerEditUtils::boundingBoxFromPointList( const QgsPointSequence &list, double &xmin, double &ymin, double &xmax, double &ymax ) const
914{
915 if ( list.empty() )
916 {
917 return false;
918 }
919
920 xmin = std::numeric_limits<double>::max();
921 xmax = -std::numeric_limits<double>::max();
922 ymin = std::numeric_limits<double>::max();
923 ymax = -std::numeric_limits<double>::max();
924
925 for ( QgsPointSequence::const_iterator it = list.constBegin(); it != list.constEnd(); ++it )
926 {
927 if ( it->x() < xmin )
928 {
929 xmin = it->x();
930 }
931 if ( it->x() > xmax )
932 {
933 xmax = it->x();
934 }
935 if ( it->y() < ymin )
936 {
937 ymin = it->y();
938 }
939 if ( it->y() > ymax )
940 {
941 ymax = it->y();
942 }
943 }
944
945 return true;
946}
GeometryOperationResult
Success or failure of a geometry operation.
Definition qgis.h:1709
@ AddPartSelectedGeometryNotFound
The selected geometry cannot be found.
@ InvalidInputGeometryType
The input geometry (ring, part, split line, etc.) has not the correct geometry type.
@ Success
Operation succeeded.
@ AddRingNotInExistingFeature
The input ring doesn't have any existing ring to fit into.
@ AddRingNotClosed
The input ring is not closed.
@ NothingHappened
Nothing happened, without any error.
@ InvalidBaseGeometry
The base geometry on which the operation is done is invalid or empty.
@ LayerNotEditable
Cannot edit layer.
@ Feet
Imperial feet.
@ ExactIntersect
Use exact geometry intersection (slower) instead of bounding boxes.
@ GeometryRatio
New values are computed by the ratio of their area/length compared to the area/length of the original...
@ UnsetField
Clears the field value so that the data provider backend will populate using any backend triggers or ...
@ DefaultValue
Use default field value.
@ Duplicate
Duplicate original value.
@ Polygon
Polygons.
@ Unknown
Unknown types.
@ Null
No geometry.
@ Provider
Field originates from the underlying data provider of the vector layer.
VectorEditResult
Specifies the result of a vector layer edit operation.
Definition qgis.h:1485
@ EmptyGeometry
Edit operation resulted in an empty geometry.
@ Success
Edit operation was successful.
@ FetchFeatureFailed
Unable to fetch requested feature.
@ EditFailed
Edit operation failed.
@ InvalidLayer
Edit failed due to invalid layer.
The vertex_iterator class provides STL-style iterator for vertices.
virtual bool addZValue(double zValue=0)=0
Adds a z-dimension to the geometry, initialized to a preset value.
bool isMeasure() const
Returns true if the geometry contains m values.
virtual QgsRectangle boundingBox() const
Returns the minimal bounding box for the geometry.
bool is3D() const
Returns true if the geometry is 3D and contains a z-value.
virtual int nCoordinates() const
Returns the number of nodes contained in the geometry.
virtual bool addMValue(double mValue=0)=0
Adds a measure to the geometry, initialized to a preset value.
virtual bool hasCurvedSegments() const
Returns true if the geometry contains curved segments.
A vector of attributes.
Abstract base class for curved geometry type.
Definition qgscurve.h:35
Qgis::AngularDirection orientation() const
Returns the curve's orientation, e.g.
Definition qgscurve.cpp:286
virtual QgsCurve * reversed() const =0
Returns a reversed copy of the curve, where the direction of the curve has been flipped.
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).
The feature class encapsulates a single feature including its unique ID, geometry and a list of field...
Definition qgsfeature.h:58
QgsFeatureId id
Definition qgsfeature.h:66
QgsGeometry geometry
Definition qgsfeature.h:69
bool hasGeometry() const
Returns true if the feature has an associated geometry.
Q_INVOKABLE QVariant attribute(const QString &name) const
Lookup attribute value by attribute name.
Encapsulate a field in an attribute table or data source.
Definition qgsfield.h:53
QString typeName() const
Gets the field type.
Definition qgsfield.cpp:160
bool convertCompatible(QVariant &v, QString *errorMessage=nullptr) const
Converts the provided variant to a compatible format.
Definition qgsfield.cpp:472
Qgis::FieldDomainSplitPolicy splitPolicy() const
Returns the field's split policy, which indicates how field values should be handled during a split o...
Definition qgsfield.cpp:754
bool isNumeric
Definition qgsfield.h:56
int count
Definition qgsfields.h:50
Qgis::FieldOrigin fieldOrigin(int fieldIdx) const
Returns the field's origin (value from an enumeration).
QgsField at(int i) const
Returns the field at particular index (must be in range 0..N-1).
int fieldOriginIndex(int fieldIdx) const
Returns the field's origin index (its meaning is specific to each type of origin).
double geometryPrecision() const
The precision in which geometries on this layer should be saved.
A geometry is the spatial representation of a feature.
bool deleteVertex(int atVertex)
Deletes the vertex at the given position number and item (first number is index 0)
double length() const
Returns the planar, 2-dimensional length of geometry.
bool addTopologicalPoint(const QgsPoint &point, double snappingTolerance=1e-8, double segmentSearchEpsilon=1e-12)
Adds a vertex to the segment which intersect point but don't already have a vertex there.
static QgsGeometry collectGeometry(const QVector< QgsGeometry > &geometries)
Creates a new multipart geometry from a list of QgsGeometry objects.
QVector< QgsGeometry > asGeometryCollection() const
Returns contents of the geometry as a list of geometries.
QgsAbstractGeometry * get()
Returns a modifiable (non-const) reference to the underlying abstract geometry primitive.
const QgsAbstractGeometry * constGet() const
Returns a non-modifiable (const) reference to the underlying abstract geometry primitive.
bool insertVertex(double x, double y, int beforeVertex)
Insert a new vertex before the given vertex index, ring and item (first number is index 0) If the req...
bool convertToSingleType()
Converts multi type geometry into single type geometry e.g.
Qgis::GeometryOperationResult addRing(const QVector< QgsPointXY > &ring)
Adds a new ring to this geometry.
Qgis::GeometryType type
double area() const
Returns the planar, 2-dimensional area of the geometry.
Qgis::AngularDirection polygonOrientation() const
Returns the orientation of the polygon.
void set(QgsAbstractGeometry *geometry)
Sets the underlying geometry store.
QgsAbstractGeometry::vertex_iterator vertices_begin() const
Returns STL-style iterator pointing to the first vertex of the geometry.
Qgis::GeometryOperationResult addPartV2(const QVector< QgsPointXY > &points, Qgis::WkbType wkbType=Qgis::WkbType::Unknown)
Adds a new part to a the geometry.
Qgis::GeometryOperationResult translate(double dx, double dy, double dz=0.0, double dm=0.0)
Translates this geometry by dx, dy, dz and dm.
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.
bool moveVertex(double x, double y, int atVertex)
Moves the vertex at the given position number and item (first number is index 0) to the given coordin...
QgsAbstractGeometry::vertex_iterator vertices_end() const
Returns STL-style iterator pointing to the imaginary vertex after the last vertex of the geometry.
Line string geometry type, with support for z-dimension and m-values.
QgsCoordinateReferenceSystem crs
Definition qgsmaplayer.h:82
void triggerRepaint(bool deferredUpdate=false)
Will advise the map canvas (and any other interested party) that this layer requires to be repainted.
A class to represent a 2D point.
Definition qgspointxy.h:60
Point geometry type, with support for z-dimension and m-values.
Definition qgspoint.h:49
A rectangle specified with double values.
double xMinimum() const
Returns the x minimum value (left side of rectangle).
void setYMinimum(double y)
Set the minimum y value.
double yMinimum() const
Returns the y minimum value (bottom side of rectangle).
void setXMinimum(double x)
Set the minimum x value.
double width() const
Returns the width of the rectangle.
double xMaximum() const
Returns the x maximum value (right side of rectangle).
double yMaximum() const
Returns the y maximum value (top side of rectangle).
void setYMaximum(double y)
Set the maximum y value.
void setXMaximum(double x)
Set the maximum x value.
void grow(double delta)
Grows the rectangle in place by the specified amount.
bool isEmpty() const
Returns true if the rectangle has no area.
double height() const
Returns the height of the rectangle.
T value(const QString &dynamicKeyPart=QString()) const
Returns settings value.
static const QgsSettingsEntryDouble * settingsDigitizingDefaultMValue
Settings entry digitizing default m value.
static const QgsSettingsEntryDouble * settingsDigitizingDefaultZValue
Settings entry digitizing default z value.
Represents a default, "not-specified" value for a feature attribute.
virtual bool doesStrictFeatureTypeCheck() const
Returns true if the provider is strict about the type of inserted features (e.g.
virtual QString defaultValueClause(int fieldIndex) const
Returns any default value clauses which are present at the provider for a specified field index.
int translateFeature(QgsFeatureId featureId, double dx, double dy)
Translates feature by dx, dy.
bool mergeFeatures(const QgsFeatureId &targetFeatureId, const QgsFeatureIds &mergeFeatureIds, const QgsAttributes &mergeAttributes, const QgsGeometry &unionGeometry, QString &errorMessage)
Merge features into a single one.
QgsVectorLayerEditUtils(QgsVectorLayer *layer)
bool insertVertex(double x, double y, QgsFeatureId atFeatureId, int beforeVertex)
Insert a new vertex before the given vertex number, in the given ring, item (first number is index 0)...
Q_DECL_DEPRECATED Qgis::GeometryOperationResult addPart(const QVector< QgsPointXY > &ring, QgsFeatureId featureId)
Adds a new part polygon to a multipart feature.
Qgis::VectorEditResult deleteVertex(QgsFeatureId featureId, int vertex)
Deletes a vertex from a feature.
Qgis::GeometryOperationResult addRingV2(QgsCurve *ring, const QgsFeatureIds &targetFeatureIds=QgsFeatureIds(), QgsFeatureIds *modifiedFeatureIds=nullptr)
Adds a ring to polygon/multipolygon features.
int addTopologicalPoints(const QgsGeometry &geom)
Adds topological points for every vertex of the geometry.
Q_DECL_DEPRECATED Qgis::GeometryOperationResult splitParts(const QVector< QgsPointXY > &splitLine, bool topologicalEditing=false)
Splits parts cut by the given line.
Q_DECL_DEPRECATED Qgis::GeometryOperationResult splitFeatures(const QVector< QgsPointXY > &splitLine, bool topologicalEditing=false)
Splits features cut by the given line.
bool moveVertex(double x, double y, QgsFeatureId atFeatureId, int atVertex)
Moves the vertex at the given position number, ring and item (first number is index 0),...
Q_DECL_DEPRECATED Qgis::GeometryOperationResult addRing(const QVector< QgsPointXY > &ring, const QgsFeatureIds &targetFeatureIds=QgsFeatureIds(), QgsFeatureId *modifiedFeatureId=nullptr)
Adds a ring to polygon/multipolygon features.
Encapsulate geometry and attributes for new features, to be passed to createFeatures.
QList< QgsVectorLayerUtils::QgsFeatureData > QgsFeaturesDataList
Alias for list of QgsFeatureData.
static QgsFeature createFeature(const QgsVectorLayer *layer, const QgsGeometry &geometry=QgsGeometry(), const QgsAttributeMap &attributes=QgsAttributeMap(), QgsExpressionContext *context=nullptr)
Creates a new feature ready for insertion into a layer.
static QgsFeatureList createFeatures(const QgsVectorLayer *layer, const QgsFeaturesDataList &featuresData, QgsExpressionContext *context=nullptr)
Creates a set of new features ready for insertion into a layer.
Represents a vector layer which manages a vector based data sets.
bool isSpatial() const FINAL
Returns true if this is a geometry layer and false in case of NoGeometry (table only) or UnknownGeome...
QgsFeatureIterator getFeatures(const QgsFeatureRequest &request=QgsFeatureRequest()) const FINAL
Queries the layer for features specified in request.
bool addFeatures(QgsFeatureList &features, QgsFeatureSink::Flags flags=QgsFeatureSink::Flags()) FINAL
Adds a list of features to the sink.
Q_INVOKABLE bool deleteFeature(QgsFeatureId fid, QgsVectorLayer::DeleteContext *context=nullptr)
Deletes a feature from the layer (but does not commit it).
int selectedFeatureCount() const
Returns the number of features that are selected in this layer.
void endEditCommand()
Finish edit command and add it to undo/redo stack.
Q_INVOKABLE const QgsFeatureIds & selectedFeatureIds() const
Returns a list of the selected features IDs in this layer.
QgsGeometryOptions * geometryOptions() const
Configuration and logic to apply automatically on any edit happening on this layer.
Q_INVOKABLE Qgis::WkbType wkbType() const FINAL
Returns the WKBType or WKBUnknown in case of error.
Q_INVOKABLE QgsVectorLayerEditBuffer * editBuffer()
Buffer with uncommitted editing operations. Only valid after editing has been turned on.
QgsFeatureIterator getSelectedFeatures(QgsFeatureRequest request=QgsFeatureRequest()) const
Returns an iterator of the selected features.
QgsVectorDataProvider * dataProvider() FINAL
Returns the layer's data provider, it may be nullptr.
bool addFeature(QgsFeature &feature, QgsFeatureSink::Flags flags=QgsFeatureSink::Flags()) FINAL
Adds a single feature to the sink.
void beginEditCommand(const QString &text)
Create edit command for undo/redo operations.
Q_INVOKABLE bool changeAttributeValues(QgsFeatureId fid, const QgsAttributeMap &newValues, const QgsAttributeMap &oldValues=QgsAttributeMap(), bool skipDefaultValues=false, QgsVectorLayerToolsContext *context=nullptr)
Changes attributes' values for a feature (but does not immediately commit the changes).
bool changeGeometry(QgsFeatureId fid, QgsGeometry &geometry, bool skipDefaultValue=false)
Changes a feature's geometry within the layer's edit buffer (but does not immediately commit the chan...
static bool hasZ(Qgis::WkbType type)
Tests whether a WKB type contains the z-dimension.
static bool hasM(Qgis::WkbType type)
Tests whether a WKB type contains m values.
static bool isSingleType(Qgis::WkbType type)
Returns true if the WKB type is a single type.
bool qgsDoubleNear(double a, double b, double epsilon=4 *std::numeric_limits< double >::epsilon())
Compare two doubles (but allow some difference)
Definition qgis.h:5445
QVector< QgsPoint > QgsPointSequence
QMap< int, QVariant > QgsAttributeMap
QList< QgsFeature > QgsFeatureList
#define FID_NULL
QSet< QgsFeatureId > QgsFeatureIds
qint64 QgsFeatureId
64 bit feature ids negative numbers are used for uncommitted/newly added features