56 addVerticesFromFeature( feat, valueAttribute, transform, feedback );
84 addVerticesFromFeature( feat, valueAttribute, transform, feedback );
88 addBreakLinesFromFeature( feat, valueAttribute, transform, feedback );
100 return mTriangulation->addPoint( vertex );
105 return mTriangulation->triangulationToMesh( feedback );
130 if ( valueAttribute >= 0 )
131 value = feature.
attribute( valueAttribute ).toDouble();
137 if ( valueAttribute < 0 )
138 mTriangulation->addPoint( *vit );
149 double valueOnVertex = 0;
150 if ( valueAttribute >= 0 )
151 valueOnVertex = feature.
attribute( valueAttribute ).toDouble();
154 std::vector<const QgsCurve *> curves;
169 std::vector< const QgsCurvePolygon * > polygons;
175 polygons.emplace_back( qgsgeometry_cast< const QgsCurvePolygon * >( ms->
geometryN( i ) ) );
180 polygons.emplace_back( qgsgeometry_cast< const QgsCurvePolygon * >( geom.
constGet() ) );
190 if ( polygon->exteriorRing() )
191 curves.emplace_back( polygon->exteriorRing() );
193 for (
int i = 0; i < polygon->numInteriorRings(); ++i )
197 curves.emplace_back( polygon->interiorRing( i ) );
210 curves.emplace_back( qgsgeometry_cast< const QgsCurve * >( mc->
geometryN( i ) ) );
215 curves.emplace_back( qgsgeometry_cast< const QgsCurve * >( geom.
constGet() ) );
219 for (
const QgsCurve *curve : curves )
228 curve->points( linePoints );
229 bool hasZ = curve->is3D();
230 if ( valueAttribute >= 0 )
231 for (
int i = 0; i < linePoints.count(); ++i )
236 linePoints[i].setZ( valueOnVertex );
239 const QgsPoint &point = linePoints.at( i );
240 linePoints[i] =
QgsPoint( point.
x(), point.
y(), valueOnVertex );
251 mDataset = std::make_unique< QgsMeshZValueDataset >( mesh );
261 if ( datasetIndex != 0 )
264 return mDataset->metadata();
274 return mDataset.get();
281 return QDomElement();
288 if ( vertex.z() < mZMinimum )
289 mZMinimum = vertex.z();
290 if ( vertex.z() > mZMaximum )
291 mZMaximum = vertex.z();
297 if ( valueIndex < 0 || valueIndex >= mMesh.
vertexCount() )
307 QVector<double> zValues( count );
308 for (
int i = valueIndex; i < valueIndex + count; ++i )
309 zValues[i - valueIndex] = mMesh.
vertex( i ).
z();
316 Q_UNUSED( faceIndex );
325 return ( faceIndex > 0 && faceIndex < mMesh.
faceCount() );
342 return QObject::tr(
"Delaunay triangulation" );
348 QList<int> vertexIndextoTriangulate;
350 QList<int> removedVerticesFromTriangulation;
355 vertexIndextoTriangulate.append( vertexIndex );
357 removedVerticesFromTriangulation.append( vertexIndex );
360 bool triangulationReady =
false;
364 while ( !triangulationReady )
368 QVector<int> triangulationVertexToMeshVertex( vertexIndextoTriangulate.count() );
371 for (
int i = 0; i < vertexIndextoTriangulate.count(); ++i )
373 triangulationVertexToMeshVertex[i] = vertexIndextoTriangulate.at( i );
374 triangulation.
addVertex( destinationMesh->
vertices.at( vertexIndextoTriangulate.at( i ) ) );
380 QVector<QgsMeshFace> rawDestinationFaces = resultingTriangulation.
faces;
382 for (
QgsMeshFace &destinationFace : rawDestinationFaces )
384 for (
int &vertexIndex : destinationFace )
385 vertexIndex = triangulationVertexToMeshVertex[vertexIndex];
389 QVector<QgsMeshFace> destinationFaces;
390 for (
const QgsMeshFace &face : rawDestinationFaces )
393 destinationFaces.append( face );
396 bool facesReady =
false;
398 while ( !facesReady && !giveUp )
410 triangulationReady =
true;
426 removedVerticesFromTriangulation.append( error.
elementIndex );
427 vertexIndextoTriangulate.removeOne( error.
elementIndex );
438 if ( !removedVerticesFromTriangulation.isEmpty() )
439 mMessage = QObject::tr(
"%n vertices have not been included in the triangulation",
nullptr, removedVerticesFromTriangulation.count() );
443 if ( triangulationReady && !giveUp )
@ TooManyVerticesInFace
A face has more vertices than the maximum number supported per face.
@ InvalidFace
An error occurs due to an invalid face (for example, vertex indexes are unordered)
@ UniqueSharedVertex
A least two faces share only one vertices.
@ ManifoldFace
ManifoldFace.
@ InvalidVertex
An error occurs due to an invalid vertex (for example, vertex index is out of range the available ver...
@ FlatFace
A flat face is present.
GeometryType
The geometry types are used to group Qgis::WkbType in a coarse way.
@ Forward
Forward transform (from source to destination)
The vertex_iterator class provides STL-style iterator for vertices.
This class represents a coordinate reference system (CRS).
Custom exception class for Coordinate Reference System related exceptions.
Curve polygon geometry type.
Abstract base class for curved geometry type.
DualEdgeTriangulation is an implementation of a triangulation class based on the dual edge data struc...
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.
The feature class encapsulates a single feature including its unique ID, geometry and a list of field...
Q_INVOKABLE QVariant attribute(const QString &name) const
Lookup attribute value by attribute name.
Base class for feedback objects to be used for cancellation of something running in a worker thread.
bool isCanceled() const
Tells whether the operation has been canceled already.
void setProgress(double progress)
Sets the current progress for the feedback object.
int numGeometries() const
Returns the number of geometries within the collection.
const QgsAbstractGeometry * geometryN(int n) const
Returns a const reference to a geometry from within the collection.
A geometry is the spatial representation of a feature.
Qgis::GeometryOperationResult transform(const QgsCoordinateTransform &ct, Qgis::TransformDirection direction=Qgis::TransformDirection::Forward, bool transformZ=false)
Transforms this geometry as described by the coordinate transform ct.
const QgsAbstractGeometry * constGet() const
Returns a non-modifiable (const) reference to the underlying abstract geometry primitive.
bool isMultipart() const
Returns true if WKB of the geometry is of WKBMulti* type.
QgsAbstractGeometry::vertex_iterator vertices_begin() const
Returns STL-style iterator pointing to the first vertex of the geometry.
Qgis::WkbType wkbType() const
Returns type of the geometry as a WKB type (point / linestring / polygon etc.)
QgsAbstractGeometry::vertex_iterator vertices_end() const
Returns STL-style iterator pointing to the imaginary vertex after the last vertex of the geometry.
@ SourceBreakLines
Break lines.
QList< int > mInputVertices
QgsMeshDataBlock is a block of integers/doubles that can be used to retrieve: active flags (e....
@ ScalarDouble
Scalar double values.
@ ActiveFlagInteger
Integer boolean flag whether face is active.
void setValues(const QVector< double > &vals)
Sets values.
void setValid(bool valid)
Sets block validity.
Abstract class that represents a dataset group.
void calculateStatistic() const
Calculates the statistics (minimum and maximum)
QgsMeshDatasetValue represents single dataset value.
Abstract class that represents a dataset.
QgsMeshEditingDelaunayTriangulation()
QString text() const override
Returns a short text string describing what this advanced edit does. Default implementation return a ...
Class that represents an error during mesh editing.
Qgis::MeshEditingErrorType errorType
Class that makes edit operation on a mesh.
bool isFaceGeometricallyCompatible(const QgsMeshFace &face) const
Returns true if the face does not intersect or contains any other elements (faces or vertices) The to...
bool isVertexOnBoundary(int vertexIndex) const
Returns whether the vertex with index vertexIndex is on a boundary.
bool isVertexFree(int vertexIndex) const
Returns whether the vertex with index vertexIndex is a free vertex.
QgsTopologicalMesh & topologicalMesh()
Returns a reference to the topological mesh.
Class that handles mesh creation with Delaunay constrained triangulation.
bool addBreakLines(QgsFeatureIterator &lineFeatureIterator, int valueAttribute, const QgsCoordinateTransform &transformContext, QgsFeedback *feedback=nullptr, long featureCount=1)
Adds break lines from a vector layer, return true if successful.
int addVertex(const QgsPoint &vertex)
Adds a new vertex in the triangulation and returns the index of the new vertex.
void setCrs(const QgsCoordinateReferenceSystem &crs)
Sets the coordinate reference system used for the triangulation.
bool addVertices(QgsFeatureIterator &vertexFeatureIterator, int valueAttribute, const QgsCoordinateTransform &transform, QgsFeedback *feedback=nullptr, long featureCount=1)
Adds vertices to the triangulation from a feature iterator, return true if successful.
QgsMesh triangulatedMesh(QgsFeedback *feedback=nullptr) const
Returns the triangulated mesh.
int datasetCount() const override
Returns the count of datasets in the group.
QDomElement writeXml(QDomDocument &doc, const QgsReadWriteContext &context) const override
Write dataset group information in a DOM element.
void initialize() override
Initialize the dataset group.
QgsMeshDataset * dataset(int index) const override
Returns the dataset with index.
QgsMeshDatasetMetadata datasetMetadata(int datasetIndex) const override
Returns the metadata of the dataset with index datasetIndex.
QgsMeshZValueDatasetGroup(const QString &datasetGroupName, const QgsMesh &mesh)
Constructor.
QgsMeshDataBlock areFacesActive(int faceIndex, int count) const override
Returns whether faces are active.
QgsMeshDatasetMetadata metadata() const override
Returns the metadata of the dataset.
QgsMeshZValueDataset(const QgsMesh &mesh)
Constructor with the mesh.
bool isActive(int faceIndex) const override
Returns whether the face is active.
int valuesCount() const override
Returns the values count.
QgsMeshDatasetValue datasetValue(int valueIndex) const override
Returns the value with index valueIndex.
QgsMeshDataBlock datasetValues(bool isScalar, int valueIndex, int count) const override
Returns count values from valueIndex.
Multi curve geometry collection.
Multi surface geometry collection.
Point geometry type, with support for z-dimension and m-values.
The class is used as a container of context for various read/write operations on other objects.
Class that contains topological differences between two states of a topological mesh,...
Class that contains independent faces an topological information about this faces.
QgsMeshEditingError checkConsistency() const
Checks the consistency of the topological mesh and return false if there is a consistency issue.
QgsMeshEditingError facesCanBeAdded(const TopologicalFaces &topologicalFaces) const
Returns whether the faces can be added to the mesh.
Changes addFaces(const TopologicalFaces &topologicFaces)
Adds faces topologicFaces to the topologic mesh.
QgsMesh * mesh() const
Returns a pointer to the wrapped mesh.
static TopologicalFaces createNewTopologicalFaces(const QVector< QgsMeshFace > &faces, bool uniqueSharedVertexAllowed, QgsMeshEditingError &error)
Creates new topological faces that are not yet included in the mesh.
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...
QVector< QgsPoint > QgsPointSequence
#define QgsDebugMsgLevel(str, level)
QVector< int > QgsMeshFace
List of vertex indexes.
const QgsCoordinateReferenceSystem & crs
Mesh - vertices, edges and faces.
int vertexCount() const
Returns number of vertices.
QVector< QgsMeshVertex > vertices
QVector< QgsMeshFace > faces
int faceCount() const
Returns number of faces.
QgsMeshVertex vertex(int index) const
Returns a vertex at the index.