48 return QStringLiteral(
"vector" );
53 QVector<QgsGeometry> res;
57 for (
const Feature &feature : it.value() )
59 res.append( feature.geometry );
71 return asIndividualFeatures(
type, feedback );
86 return snapPointToIndividualFeatures( point, context );
102 visitFeaturesInRange( distanceRange, elevationRange, visitFeature );
106 QVector< QVariantMap> idsList;
107 for (
auto it = ids.constBegin(); it != ids.constEnd(); ++it )
108 idsList.append( QVariantMap( {{QStringLiteral(
"id" ), *it}} ) );
115 QHash< QgsFeatureId, QVariantMap >
features;
116 auto visitFeature = [&
features](
QgsFeatureId featureId,
double delta,
double distance,
double elevation )
118 auto it =
features.find( featureId );
121 features[ featureId ] = QVariantMap( {{QStringLiteral(
"id" ), featureId },
122 {QStringLiteral(
"delta" ), delta },
123 {QStringLiteral(
"distance" ), distance },
124 {QStringLiteral(
"elevation" ), elevation }
129 const double currentDelta = it.value().value( QStringLiteral(
"delta" ) ).toDouble();
130 if ( delta < currentDelta )
132 *it = QVariantMap( {{QStringLiteral(
"id" ), featureId },
133 {QStringLiteral(
"delta" ), delta },
134 {QStringLiteral(
"distance" ), distance },
135 {QStringLiteral(
"elevation" ), elevation }
143 QVector< QVariantMap> attributes;
145 attributes.append( *it );
147 QVector<QgsProfileIdentifyResults> res;
149 if ( !attributes.empty() )
155 res.reserve( surfaceResults.size() );
168 double bestSnapDistance = std::numeric_limits< double >::max();
170 auto visitFeature = [&bestSnapDistance, &res](
QgsFeatureId,
double delta,
double distance,
double elevation )
172 if ( distance < bestSnapDistance )
174 bestSnapDistance = delta;
184void QgsVectorLayerProfileResults::visitFeaturesAtPoint(
const QgsProfilePoint &point,
double maximumPointDistanceDelta,
double maximumPointElevationDelta,
double maximumSurfaceElevationDelta,
const std::function<
void(
QgsFeatureId,
double delta,
double distance,
double elevation ) > &visitor,
bool visitWithin )
192 for (
const Feature &feature : it.value() )
194 const QgsRectangle featureBounds = feature.crossSectionGeometry.boundingBox();
195 if ( ( featureBounds.
xMinimum() - maximumPointDistanceDelta <= point.
distance() ) && ( featureBounds.
xMaximum() + maximumPointDistanceDelta >= point.
distance() ) )
197 switch ( feature.crossSectionGeometry.type() )
201 for (
auto partIt = feature.crossSectionGeometry.const_parts_begin(); partIt != feature.crossSectionGeometry.const_parts_end(); ++partIt )
203 if (
const QgsPoint *candidatePoint = qgsgeometry_cast< const QgsPoint * >( *partIt ) )
205 const double snapDistanceDelta = std::fabs( point.
distance() - candidatePoint->x() );
206 if ( snapDistanceDelta > maximumPointDistanceDelta )
209 const double snapHeightDelta = std::fabs( point.
elevation() - candidatePoint->y() );
210 if ( snapHeightDelta > maximumPointElevationDelta )
213 const double snapDistance = candidatePoint->distance( targetPoint );
214 visitor( feature.featureId, snapDistance, candidatePoint->x(), candidatePoint->y() );
222 for (
auto partIt = feature.crossSectionGeometry.const_parts_begin(); partIt != feature.crossSectionGeometry.const_parts_end(); ++partIt )
224 if (
const QgsCurve *line = qgsgeometry_cast< const QgsCurve * >( *partIt ) )
227 if (
const QgsLineString *lineString = qgsgeometry_cast< const QgsLineString * >( line ) )
229 if ( lineString->numPoints() == 2 &&
qgsDoubleNear( lineString->pointN( 0 ).x(), lineString->pointN( 1 ).x() ) )
231 const double snapDistanceDelta = std::fabs( point.
distance() - lineString->pointN( 0 ).x() );
232 if ( snapDistanceDelta > maximumPointDistanceDelta )
235 const double snapHeightDelta = std::fabs( point.
elevation() - lineString->pointN( 0 ).y() );
236 if ( snapHeightDelta <= maximumPointElevationDelta )
238 const double snapDistanceP1 = lineString->pointN( 0 ).distance( targetPoint );
239 visitor( feature.featureId, snapDistanceP1, lineString->pointN( 0 ).x(), lineString->pointN( 0 ).y() );
242 const double snapHeightDelta2 = std::fabs( point.
elevation() - lineString->pointN( 1 ).y() );
243 if ( snapHeightDelta2 <= maximumPointElevationDelta )
245 const double snapDistanceP2 = lineString->pointN( 1 ).distance( targetPoint );
246 visitor( feature.featureId, snapDistanceP2, lineString->pointN( 1 ).x(), lineString->pointN( 1 ).y() );
251 double elevation1 = lineString->pointN( 0 ).y();
252 double elevation2 = lineString->pointN( 1 ).y();
253 if ( elevation1 > elevation2 )
254 std::swap( elevation1, elevation2 );
258 const double snapDistance = std::fabs( lineString->pointN( 0 ).x() - point.
distance() );
259 visitor( feature.featureId, snapDistance, lineString->pointN( 0 ).x(), point.
elevation() );
266 const QgsRectangle partBounds = ( *partIt )->boundingBox();
276 const QgsGeometry points( cutLineGeos.intersection( line ) );
278 for (
auto vertexIt = points.vertices_begin(); vertexIt != points.vertices_end(); ++vertexIt )
280 const double snapHeightDelta = std::fabs( point.
elevation() - ( *vertexIt ).y() );
281 if ( snapHeightDelta > maximumSurfaceElevationDelta )
284 const double snapDistance = ( *vertexIt ).distance( targetPoint );
285 visitor( feature.featureId, snapDistance, ( *vertexIt ).x(), ( *vertexIt ).y() );
302 for (
auto partIt = feature.crossSectionGeometry.const_parts_begin(); partIt != feature.crossSectionGeometry.const_parts_end(); ++partIt )
304 if (
const QgsCurve *exterior = qgsgeometry_cast< const QgsPolygon * >( *partIt )->exteriorRing() )
306 const QgsRectangle partBounds = ( *partIt )->boundingBox();
316 const QgsGeometry points( cutLineGeos.intersection( exterior ) );
317 for (
auto vertexIt = points.vertices_begin(); vertexIt != points.vertices_end(); ++vertexIt )
319 const double snapHeightDelta = std::fabs( point.
elevation() - ( *vertexIt ).y() );
320 if ( snapHeightDelta > maximumSurfaceElevationDelta )
323 const double snapDistance = ( *vertexIt ).distance( targetPoint );
324 visitor( feature.featureId, snapDistance, ( *vertexIt ).x(), ( *vertexIt ).y() );
345 profileRangeGeos.prepareGeometry();
349 for (
const Feature &feature : it.value() )
351 if ( feature.crossSectionGeometry.boundingBoxIntersects( profileRange ) )
353 switch ( feature.crossSectionGeometry.type() )
357 for (
auto partIt = feature.crossSectionGeometry.const_parts_begin(); partIt != feature.crossSectionGeometry.const_parts_end(); ++partIt )
359 if (
const QgsPoint *candidatePoint = qgsgeometry_cast< const QgsPoint * >( *partIt ) )
361 if ( profileRange.contains( candidatePoint->x(), candidatePoint->y() ) )
363 visitor( feature.featureId );
373 if ( profileRangeGeos.intersects( feature.crossSectionGeometry.constGet() ) )
375 visitor( feature.featureId );
395 renderResultsAsIndividualFeatures( context );
400 renderMarkersOverContinuousSurfacePlot( context );
413 painter->setBrush( Qt::NoBrush );
414 painter->setPen( Qt::NoPen );
421 const QRectF visibleRegion( minDistance,
minZ, maxDistance - minDistance,
maxZ -
minZ );
422 QPainterPath clipPath;
423 clipPath.addPolygon( context.
worldTransform().map( visibleRegion ) );
424 painter->setClipPath( clipPath, Qt::ClipOperation::IntersectClip );
426 const QgsRectangle clipPathRect( clipPath.boundingRect() );
430 if ( profileFeature.crossSectionGeometry.isEmpty() )
433 QgsGeometry transformed = profileFeature.crossSectionGeometry;
440 switch ( transformed.
type() )
444 if (
const QgsPoint *point = qgsgeometry_cast< const QgsPoint * >( transformed.
constGet() ) )
446 markerSymbol->renderPoint( QPointF( point->x(), point->y() ),
nullptr, context.
renderContext() );
448 else if (
const QgsMultiPoint *multipoint = qgsgeometry_cast< const QgsMultiPoint * >( transformed.
constGet() ) )
450 const int numGeometries = multipoint->numGeometries();
451 for (
int i = 0; i < numGeometries; ++i )
453 markerSymbol->renderPoint( QPointF( multipoint->pointN( i )->x(), multipoint->pointN( i )->y() ),
nullptr, context.
renderContext() );
461 if (
const QgsLineString *line = qgsgeometry_cast< const QgsLineString * >( transformed.
constGet() ) )
463 lineSymbol->renderPolyline( line->asQPolygonF(),
nullptr, context.
renderContext() );
465 else if (
const QgsMultiLineString *multiLinestring = qgsgeometry_cast< const QgsMultiLineString * >( transformed.
constGet() ) )
467 const int numGeometries = multiLinestring->numGeometries();
468 for (
int i = 0; i < numGeometries; ++i )
470 lineSymbol->renderPolyline( multiLinestring->lineStringN( i )->asQPolygonF(),
nullptr, context.
renderContext() );
478 if (
const QgsPolygon *polygon = qgsgeometry_cast< const QgsPolygon * >( transformed.
constGet() ) )
480 if (
const QgsCurve *exterior = polygon->exteriorRing() )
481 fillSymbol->renderPolygon( exterior->asQPolygonF(),
nullptr,
nullptr, context.
renderContext() );
483 else if (
const QgsMultiPolygon *multiPolygon = qgsgeometry_cast< const QgsMultiPolygon * >( transformed.
constGet() ) )
485 const int numGeometries = multiPolygon->numGeometries();
486 for (
int i = 0; i < numGeometries; ++i )
488 fillSymbol->renderPolygon( multiPolygon->polygonN( i )->exteriorRing()->asQPolygonF(),
nullptr,
nullptr, context.
renderContext() );
505 std::unique_ptr< QgsFeatureRenderer > renderer(
mLayer->renderer()->clone() );
510 QSet<QString> attributes = renderer->usedAttributes( context.
renderContext() );
512 std::unique_ptr< QgsMarkerSymbol > marker(
mMarkerSymbol->clone() );
513 std::unique_ptr< QgsLineSymbol > line(
mLineSymbol->clone() );
514 std::unique_ptr< QgsFillSymbol > fill(
mFillSymbol->clone() );
515 attributes.unite( marker->usedAttributes( context.
renderContext() ) );
516 attributes.unite( line->usedAttributes( context.
renderContext() ) );
517 attributes.unite( fill->usedAttributes( context.
renderContext() ) );
527 if ( !rendererSymbol )
531 marker->setOpacity( rendererSymbol->
opacity() );
532 line->setColor( rendererSymbol->
color() );
533 line->setOpacity( rendererSymbol->
opacity() );
534 fill->setColor( rendererSymbol->
color() );
535 fill->setOpacity( rendererSymbol->
opacity() );
541 const QVector< Feature > profileFeatures =
features.value( feature.
id() );
542 for (
const Feature &profileFeature : profileFeatures )
544 renderResult( profileFeature,
546 rendererSymbol->
type() ==
Qgis::SymbolType::Line ? qgis::down_cast<
QgsLineSymbol * >( rendererSymbol ) : line.get(),
547 rendererSymbol->
type() ==
Qgis::SymbolType::Fill ? qgis::down_cast<
QgsFillSymbol * >( rendererSymbol ) : fill.get() );
559 QSet<QString> attributes;
574 const QVector< Feature > profileFeatures =
features.value( feature.
id() );
575 for (
const Feature &profileFeature : profileFeatures )
594 painter->setBrush( Qt::NoBrush );
595 painter->setPen( Qt::NoPen );
602 const QRectF visibleRegion( minDistance,
minZ, maxDistance - minDistance,
maxZ -
minZ );
603 QPainterPath clipPath;
604 clipPath.addPolygon( context.
worldTransform().map( visibleRegion ) );
605 painter->setClipPath( clipPath, Qt::ClipOperation::IntersectClip );
611 if ( std::isnan( pointIt.value() ) )
621 QVector<QgsAbstractProfileResults::Feature> res;
628 for (
const Feature &feature : it.value() )
635 outFeature.
attributes = {{QStringLiteral(
"id" ), feature.featureId }};
643 outFeature.
geometry = feature.crossSectionGeometry;
660 mId = vlGenerator->mId;
663 mMarkerSymbol.reset( vlGenerator->mProfileMarkerSymbol->clone() );
675 , mProfileCurve( request.profileCurve() ? request.profileCurve()->clone() : nullptr )
676 , mTerrainProvider( request.terrainProvider() ? request.terrainProvider()->clone() : nullptr )
677 , mTolerance( request.tolerance() )
678 , mSourceCrs( layer->crs3D() )
679 , mTargetCrs( request.
crs() )
680 , mTransformContext( request.transformContext() )
681 , mExtent( layer->extent() )
683 , mOffset( layer->elevationProperties()->zOffset() )
684 , mScale( layer->elevationProperties()->zScale() )
690 , mExpressionContext( request.expressionContext() )
691 , mFields( layer->fields() )
692 , mDataDefinedProperties( layer->elevationProperties()->dataDefinedProperties() )
693 , mWkbType( layer->wkbType() )
699 if ( mTerrainProvider )
700 mTerrainProvider->prepare();
704 mProfileCurve->dropZValue();
722 if ( !mProfileCurve || mFeedback->isCanceled() )
726 qgsgeometry_cast<QgsLineString *>( mProfileCurve.get() ) )
732 std::unique_ptr< QgsCurve > origCurve = std::move( mProfileCurve );
733 std::unique_ptr< QgsVectorLayerProfileResults > totalResults;
734 double distanceProcessed = 0;
736 QVector<QgsLineString *> disjointParts = profileLine->splitToDisjointXYParts();
737 for (
int i = 0; i < disjointParts.size(); i++ )
739 mProfileCurve.reset( disjointParts[i] );
740 if ( !generateProfileInner() )
742 mProfileCurve = std::move( origCurve );
745 for (
int j = i + 1; j < disjointParts.size(); j++ )
746 delete disjointParts[j];
753 totalResults.reset( mResults.release() );
757 totalResults->mRawPoints.append( mResults->mRawPoints );
758 totalResults->minZ = std::min( totalResults->minZ, mResults->minZ );
759 totalResults->maxZ = std::max( totalResults->maxZ, mResults->maxZ );
760 for (
auto it = mResults->mDistanceToHeightMap.constKeyValueBegin();
761 it != mResults->mDistanceToHeightMap.constKeyValueEnd();
764 totalResults->mDistanceToHeightMap[it->first + distanceProcessed] = it->second;
766 for (
auto it = mResults->features.constKeyValueBegin();
767 it != mResults->features.constKeyValueEnd();
772 feature.crossSectionGeometry.translate( distanceProcessed, 0 );
773 totalResults->features[it->first].push_back( feature );
778 distanceProcessed += mProfileCurve->length();
781 mProfileCurve = std::move( origCurve );
782 mResults.reset( totalResults.release() );
786 return generateProfileInner();
792 mTransformedCurve.reset( mProfileCurve->clone() );
794 if ( mTerrainProvider )
795 mTargetToTerrainProviderTransform =
QgsCoordinateTransform( mTargetCrs, mTerrainProvider->crs(), mTransformContext );
803 QgsDebugError( QStringLiteral(
"Error transforming profile line to vector CRS" ) );
807 const QgsRectangle profileCurveBoundingBox = mTransformedCurve->boundingBox();
808 if ( !profileCurveBoundingBox.
intersects( mExtent ) )
811 if ( mFeedback->isCanceled() )
814 mResults = std::make_unique< QgsVectorLayerProfileResults >();
815 mResults->mLayer = mLayer;
816 mResults->copyPropertiesFromGenerator(
this );
818 mProfileCurveEngine.reset(
new QgsGeos( mProfileCurve.get() ) );
819 mProfileCurveEngine->prepareGeometry();
821 if ( mTolerance == 0.0 )
823 mProfileBufferedCurve = std::unique_ptr<QgsAbstractGeometry>( mProfileCurve->clone() );
830 mProfileBufferedCurveEngine.reset(
new QgsGeos( mProfileBufferedCurve.get() ) );
831 mProfileBufferedCurveEngine->prepareGeometry();
833 mDataDefinedProperties.
prepare( mExpressionContext );
835 if ( mFeedback->isCanceled() )
841 if ( !generateProfileForPoints() )
846 if ( !generateProfileForLines() )
851 if ( !generateProfileForPolygons() )
865 return mResults.release();
870 return mFeedback.get();
873bool QgsVectorLayerProfileGenerator::generateProfileForPoints()
888 while ( !mFeedback->isCanceled() && it.
nextFeature( feature ) )
895 if ( mProfileBufferedCurveEngine->intersects( *it ) )
897 processIntersectionPoint( qgsgeometry_cast< const QgsPoint * >( *it ), feature );
901 return !mFeedback->isCanceled();
904void QgsVectorLayerProfileGenerator::processIntersectionPoint(
const QgsPoint *point,
const QgsFeature &feature )
909 const double height = featureZToHeight( point->
x(), point->
y(), point->
z(), offset );
910 mResults->mRawPoints.append(
QgsPoint( point->
x(), point->
y(), height ) );
911 mResults->minZ = std::min( mResults->minZ, height );
912 mResults->maxZ = std::max( mResults->maxZ, height );
914 const double distanceAlongProfileCurve = mProfileCurveEngine->lineLocatePoint( *point, &error );
915 mResults->mDistanceToHeightMap.insert( distanceAlongProfileCurve, height );
919 if ( mExtrusionEnabled )
924 QgsPoint( point->
x(), point->
y(), height + extrusion ) ) );
926 QgsPoint( distanceAlongProfileCurve, height + extrusion ) ) );
927 mResults->minZ = std::min( mResults->minZ, height + extrusion );
928 mResults->maxZ = std::max( mResults->maxZ, height + extrusion );
936 mResults->features[resultFeature.
featureId].append( resultFeature );
939void QgsVectorLayerProfileGenerator::processIntersectionCurve(
const QgsLineString *intersectionCurve,
const QgsFeature &feature )
945 double maxDistanceAlongProfileCurve = std::numeric_limits<double>::lowest();
950 const int numPoints = intersectionCurve->
numPoints();
951 QVector< double > newX( numPoints );
952 QVector< double > newY( numPoints );
953 QVector< double > newZ( numPoints );
954 QVector< double > newDistance( numPoints );
956 const double *inX = intersectionCurve->
xData();
957 const double *inY = intersectionCurve->
yData();
958 const double *inZ = intersectionCurve->
is3D() ? intersectionCurve->
zData() :
nullptr;
959 double *outX = newX.data();
960 double *outY = newY.data();
961 double *outZ = newZ.data();
962 double *outDistance = newDistance.data();
964 QVector< double > extrudedZ;
965 double *extZOut =
nullptr;
966 if ( mExtrusionEnabled )
968 extrudedZ.resize( numPoints );
969 extZOut = extrudedZ.data();
972 for (
int i = 0 ; ! mFeedback->isCanceled() && i < numPoints; ++i )
974 QgsPoint intersectionPoint( *inX, *inY, ( inZ ? *inZ : std::numeric_limits<double>::quiet_NaN() ) );
976 const double height = featureZToHeight( intersectionPoint.
x(), intersectionPoint.
y(), intersectionPoint.
z(), offset );
977 const double distanceAlongProfileCurve = mProfileCurveEngine->lineLocatePoint( intersectionPoint, &error );
979 maxDistanceAlongProfileCurve = std::max( maxDistanceAlongProfileCurve, distanceAlongProfileCurve );
981 mResults->mRawPoints.append(
QgsPoint( intersectionPoint.
x(), intersectionPoint.
y(), height ) );
982 mResults->minZ = std::min( mResults->minZ, height );
983 mResults->maxZ = std::max( mResults->maxZ, height );
985 mResults->mDistanceToHeightMap.insert( distanceAlongProfileCurve, height );
986 *outDistance++ = distanceAlongProfileCurve;
988 *outX++ = intersectionPoint.
x();
989 *outY++ = intersectionPoint.
y();
992 *extZOut++ = height + extrusion;
994 if ( mExtrusionEnabled )
996 mResults->minZ = std::min( mResults->minZ, height + extrusion );
997 mResults->maxZ = std::max( mResults->maxZ, height + extrusion );
1005 mResults->mDistanceToHeightMap.insert( maxDistanceAlongProfileCurve + 0.000001, std::numeric_limits<double>::quiet_NaN() );
1007 if ( mFeedback->isCanceled() )
1011 if ( mExtrusionEnabled )
1013 std::unique_ptr< QgsLineString > ring = std::make_unique< QgsLineString >( newX, newY, newZ );
1014 std::unique_ptr< QgsLineString > extrudedRing = std::make_unique< QgsLineString >( newX, newY, extrudedZ );
1015 std::unique_ptr< QgsLineString > reversedExtrusion( extrudedRing->reversed() );
1016 ring->append( reversedExtrusion.get() );
1020 std::unique_ptr< QgsLineString > distanceVHeightRing = std::make_unique< QgsLineString >( newDistance, newZ );
1021 std::unique_ptr< QgsLineString > extrudedDistanceVHeightRing = std::make_unique< QgsLineString >( newDistance, extrudedZ );
1022 std::unique_ptr< QgsLineString > reversedDistanceVHeightExtrusion( extrudedDistanceVHeightRing->reversed() );
1023 distanceVHeightRing->append( reversedDistanceVHeightExtrusion.get() );
1024 distanceVHeightRing->close();
1033 mResults->features[resultFeature.
featureId].append( resultFeature );
1036bool QgsVectorLayerProfileGenerator::generateProfileForLines()
1041 if ( mTolerance > 0 )
1052 auto processCurve = [
this](
const QgsFeature & feature,
const QgsCurve * featGeomPart )
1055 std::unique_ptr< QgsAbstractGeometry > intersection( mProfileBufferedCurveEngine->intersection( featGeomPart, &error ) );
1056 if ( !intersection )
1059 if ( mFeedback->isCanceled() )
1064 if ( intersection->isEmpty() )
1066 intersection.reset( featGeomPart->clone() );
1069 QgsGeos featGeomPartGeos( featGeomPart );
1070 featGeomPartGeos.prepareGeometry();
1072 for (
auto it = intersection->const_parts_begin();
1073 !mFeedback->isCanceled() && it != intersection->const_parts_end();
1076 if (
const QgsPoint *intersectionPoint = qgsgeometry_cast< const QgsPoint * >( *it ) )
1080 const double distance = featGeomPartGeos.lineLocatePoint( *intersectionPoint, &error );
1081 std::unique_ptr< QgsPoint > interpolatedPoint( featGeomPart->interpolatePoint( distance ) );
1083 processIntersectionPoint( interpolatedPoint.get(), feature );
1085 else if (
const QgsLineString *intersectionCurve = qgsgeometry_cast< const QgsLineString * >( *it ) )
1087 processIntersectionCurve( intersectionCurve, feature );
1094 while ( !mFeedback->isCanceled() && it.
nextFeature( feature ) )
1101 if ( mProfileBufferedCurveEngine->intersects( *it ) )
1103 processCurve( feature, qgsgeometry_cast< const QgsCurve * >( *it ) );
1108 return !mFeedback->isCanceled();
1111QgsPoint QgsVectorLayerProfileGenerator::interpolatePointOnTriangle(
const QgsPolygon *triangle,
double x,
double y )
const
1118 const double z = QgsMeshLayerUtils::interpolateFromVerticesData( p1, p2, p3, p1.
z(), p2.
z(), p3.
z(),
QgsPointXY( x, y ) );
1122void QgsVectorLayerProfileGenerator::processTriangleIntersectForPoint(
const QgsPolygon *triangle,
const QgsPoint *p, QVector< QgsGeometry > &transformedParts, QVector< QgsGeometry > &crossSectionParts )
1124 const QgsPoint interpolatedPoint = interpolatePointOnTriangle( triangle, p->
x(), p->
y() );
1125 mResults->mRawPoints.append( interpolatedPoint );
1126 mResults->minZ = std::min( mResults->minZ, interpolatedPoint.
z() );
1127 mResults->maxZ = std::max( mResults->maxZ, interpolatedPoint.
z() );
1130 const double distance = mProfileCurveEngine->lineLocatePoint( *p, &lastError );
1131 mResults->mDistanceToHeightMap.insert( distance, interpolatedPoint.
z() );
1133 if ( mExtrusionEnabled )
1138 QgsPoint( interpolatedPoint.
x(), interpolatedPoint.
y(), interpolatedPoint.
z() + extrusion ) ) ) );
1140 QgsPoint( distance, interpolatedPoint.
z() + extrusion ) ) ) );
1141 mResults->minZ = std::min( mResults->minZ, interpolatedPoint.
z() + extrusion );
1142 mResults->maxZ = std::max( mResults->maxZ, interpolatedPoint.
z() + extrusion );
1151void QgsVectorLayerProfileGenerator::processTriangleIntersectForLine(
const QgsPolygon *triangle,
const QgsLineString *intersectionLine, QVector< QgsGeometry > &transformedParts, QVector< QgsGeometry > &crossSectionParts )
1156 int numPoints = intersectionLine->
numPoints();
1157 QVector< double > newX( numPoints );
1158 QVector< double > newY( numPoints );
1159 QVector< double > newZ( numPoints );
1160 QVector< double > newDistance( numPoints );
1162 const double *inX = intersectionLine->
xData();
1163 const double *inY = intersectionLine->
yData();
1164 const double *inZ = intersectionLine->
is3D() ? intersectionLine->
zData() :
nullptr;
1165 double *outX = newX.data();
1166 double *outY = newY.data();
1167 double *outZ = newZ.data();
1168 double *outDistance = newDistance.data();
1170 double lastDistanceAlongProfileCurve = 0.0;
1171 QVector< double > extrudedZ;
1172 double *extZOut =
nullptr;
1173 double extrusion = 0;
1175 if ( mExtrusionEnabled )
1177 extrudedZ.resize( numPoints );
1178 extZOut = extrudedZ.data();
1184 for (
int i = 0 ; ! mFeedback->isCanceled() && i < numPoints; ++i )
1188 double z = inZ ? *inZ++ : 0;
1190 QgsPoint interpolatedPoint( x, y, z );
1196 interpolatedPoint = interpolatePointOnTriangle( triangle, x, y );
1198 double tempOutZ = std::isnan( interpolatedPoint.
z() ) ? 0.0 : interpolatedPoint.
z();
1201 if ( mExtrusionEnabled )
1202 *extZOut++ = tempOutZ + extrusion;
1204 mResults->mRawPoints.append( interpolatedPoint );
1205 mResults->minZ = std::min( mResults->minZ, interpolatedPoint.
z() );
1206 mResults->maxZ = std::max( mResults->maxZ, interpolatedPoint.
z() );
1207 if ( mExtrusionEnabled )
1209 mResults->minZ = std::min( mResults->minZ, interpolatedPoint.
z() + extrusion );
1210 mResults->maxZ = std::max( mResults->maxZ, interpolatedPoint.
z() + extrusion );
1213 const double distance = mProfileCurveEngine->lineLocatePoint( interpolatedPoint, &lastError );
1214 *outDistance++ = distance;
1216 mResults->mDistanceToHeightMap.insert( distance, interpolatedPoint.
z() );
1217 lastDistanceAlongProfileCurve = distance;
1221 mResults->mDistanceToHeightMap.insert( lastDistanceAlongProfileCurve + 0.000001, std::numeric_limits<double>::quiet_NaN() );
1223 if ( mFeedback->isCanceled() )
1226 if ( mExtrusionEnabled )
1228 std::unique_ptr< QgsLineString > ring = std::make_unique< QgsLineString >( newX, newY, newZ );
1229 std::unique_ptr< QgsLineString > extrudedRing = std::make_unique< QgsLineString >( newX, newY, extrudedZ );
1230 std::unique_ptr< QgsLineString > reversedExtrusion( extrudedRing->reversed() );
1231 ring->append( reversedExtrusion.get() );
1235 std::unique_ptr< QgsLineString > distanceVHeightRing = std::make_unique< QgsLineString >( newDistance, newZ );
1236 std::unique_ptr< QgsLineString > extrudedDistanceVHeightRing = std::make_unique< QgsLineString >( newDistance, extrudedZ );
1237 std::unique_ptr< QgsLineString > reversedDistanceVHeightExtrusion( extrudedDistanceVHeightRing->reversed() );
1238 distanceVHeightRing->append( reversedDistanceVHeightExtrusion.get() );
1239 distanceVHeightRing->close();
1249void QgsVectorLayerProfileGenerator::processTriangleIntersectForPolygon(
const QgsPolygon *sourcePolygon,
const QgsPolygon *intersectionPolygon, QVector< QgsGeometry > &transformedParts, QVector< QgsGeometry > &crossSectionParts )
1251 bool oldExtrusion = mExtrusionEnabled;
1284 mExtrusionEnabled =
false;
1285 if ( mProfileBufferedCurveEngine->contains( sourcePolygon ) )
1289 QgsLineString *exteriorLine = qgsgeometry_cast<QgsLineString *>( exterior );
1290 processTriangleIntersectForLine( sourcePolygon, exteriorLine, transformedParts, crossSectionParts );
1295 processTriangleIntersectForLine( sourcePolygon, interiorLine, transformedParts, crossSectionParts );
1303 QgsLineString *exteriorLine = qgsgeometry_cast<QgsLineString *>( exterior )->
clone();
1305 processTriangleIntersectForLine( sourcePolygon, exteriorLine, transformedParts, crossSectionParts );
1306 delete exteriorLine;
1311 if ( mProfileBufferedCurveEngine->contains( interiorLine ) )
1313 processTriangleIntersectForLine( sourcePolygon, interiorLine, transformedParts, crossSectionParts );
1317 interiorLine = qgsgeometry_cast<QgsLineString *>( intersectionPolygon->
interiorRing( i ) )->
clone();
1319 processTriangleIntersectForLine( sourcePolygon, interiorLine, transformedParts, crossSectionParts );
1320 delete interiorLine;
1325 mExtrusionEnabled = oldExtrusion;
1328bool QgsVectorLayerProfileGenerator::generateProfileForPolygons()
1333 if ( mTolerance > 0 )
1344 std::function< void(
const QgsPolygon *triangle,
const QgsAbstractGeometry *intersect, QVector< QgsGeometry > &, QVector< QgsGeometry > & ) > processTriangleLineIntersect;
1345 processTriangleLineIntersect = [
this](
const QgsPolygon * triangle,
const QgsAbstractGeometry * intersection, QVector< QgsGeometry > &transformedParts, QVector< QgsGeometry > &crossSectionParts )
1347 for (
auto it = intersection->const_parts_begin();
1348 ! mFeedback->isCanceled() && it != intersection->const_parts_end();
1355 if (
const QgsPoint *p = qgsgeometry_cast< const QgsPoint * >( *it ) )
1357 processTriangleIntersectForPoint( triangle, p, transformedParts, crossSectionParts );
1362 if (
const QgsLineString *intersectionLine = qgsgeometry_cast< const QgsLineString * >( *it ) )
1364 processTriangleIntersectForLine( triangle, intersectionLine, transformedParts, crossSectionParts );
1369 if (
const QgsPolygon *poly = qgsgeometry_cast< const QgsPolygon * >( *it ) )
1371 processTriangleIntersectForPolygon( triangle, poly, transformedParts, crossSectionParts );
1382 auto triangleIsCollinearInXYPlane = [](
const QgsPolygon * polygon )->
bool
1384 const QgsLineString *ring = qgsgeometry_cast< const QgsLineString * >( polygon->exteriorRing() );
1386 ring->
xAt( 1 ), ring->
yAt( 1 ),
1387 ring->
xAt( 2 ), ring->
yAt( 2 ), 0.005 );
1390 auto processPolygon = [
this, &processTriangleLineIntersect, &triangleIsCollinearInXYPlane](
const QgsCurvePolygon * polygon, QVector< QgsGeometry > &transformedParts, QVector< QgsGeometry > &crossSectionParts,
double offset,
bool & wasCollinear )
1392 std::unique_ptr< QgsPolygon > clampedPolygon;
1393 if (
const QgsPolygon *p = qgsgeometry_cast< const QgsPolygon * >( polygon ) )
1395 clampedPolygon.reset( p->
clone() );
1399 clampedPolygon.reset( qgsgeometry_cast< QgsPolygon * >( polygon->segmentize() ) );
1401 clampAltitudes( clampedPolygon.get(), offset );
1403 if ( mFeedback->isCanceled() )
1406 if ( mTolerance > 0.0 )
1409 if ( mProfileBufferedCurveEngine->intersects( clampedPolygon.get(), &error ) )
1411 std::unique_ptr< QgsAbstractGeometry > intersection;
1412 intersection.reset( mProfileBufferedCurveEngine->intersection( clampedPolygon.get(), &error ) );
1413 if ( error.isEmpty() )
1415 processTriangleLineIntersect( clampedPolygon.get(), intersection.get(), transformedParts, crossSectionParts );
1422 QgsLineString *ring = qgsgeometry_cast< QgsLineString * >( clampedPolygon->exteriorRing() );
1424 QVector< double > newX( numPoints );
1425 QVector< double > newY( numPoints );
1426 QVector< double > newZ( numPoints );
1427 double *outX = newX.data();
1428 double *outY = newY.data();
1429 double *outZ = newZ.data();
1431 const double *inX = ring->
xData();
1432 const double *inY = ring->
yData();
1433 const double *inZ = ring->
zData();
1434 for (
int i = 0 ; ! mFeedback->isCanceled() && i < ring->numPoints() - 1; ++i )
1436 *outX++ = inX[i] + i * 1.0e-9;
1437 *outY++ = inY[i] + i * 1.0e-9;
1440 std::unique_ptr< QgsPolygon > shiftedPoly;
1443 intersection.reset( mProfileBufferedCurveEngine->intersection( shiftedPoly.get(), &error ) );
1444 if ( intersection.get() )
1445 processTriangleLineIntersect( clampedPolygon.get(), intersection.get(), transformedParts, crossSectionParts );
1447 QgsDebugMsgLevel( QStringLiteral(
"processPolygon after shift bad geom! error: %1" ).arg( error ), 0 );
1455 if ( clampedPolygon->numInteriorRings() == 0 && clampedPolygon->exteriorRing() && clampedPolygon->exteriorRing()->numPoints() == 4 && clampedPolygon->exteriorRing()->isClosed() )
1458 std::unique_ptr< QgsMultiPolygon > multiPolygon = std::make_unique< QgsMultiPolygon >();
1459 multiPolygon->addGeometry( clampedPolygon.release() );
1460 tessellation =
QgsGeometry( std::move( multiPolygon ) );
1464 const QgsRectangle bounds = clampedPolygon->boundingBox();
1466 t.addPolygon( *clampedPolygon, 0 );
1469 if ( mFeedback->isCanceled() )
1476 const int numTriangles = qgsgeometry_cast< const QgsMultiPolygon * >( tessellation.
constGet() )->numGeometries();
1477 for (
int i = 0; ! mFeedback->isCanceled() && i < numTriangles; ++i )
1479 const QgsPolygon *triangle = qgsgeometry_cast< const QgsPolygon * >( qgsgeometry_cast< const QgsMultiPolygon * >( tessellation.
constGet() )->geometryN( i ) );
1481 if ( triangleIsCollinearInXYPlane( triangle ) )
1483 wasCollinear =
true;
1484 const QgsLineString *ring = qgsgeometry_cast< const QgsLineString * >( polygon->exteriorRing() );
1487 if (
const QgsLineString *ls = qgsgeometry_cast< const QgsLineString * >( mProfileCurve.get() ) )
1489 for (
int curveSegmentIndex = 0; curveSegmentIndex < mProfileCurve->numPoints() - 1; ++curveSegmentIndex )
1491 const QgsPoint p1 = ls->pointN( curveSegmentIndex );
1492 const QgsPoint p2 = ls->pointN( curveSegmentIndex + 1 );
1495 double minZ = std::numeric_limits< double >::max();
1496 double maxZ = std::numeric_limits< double >::lowest();
1498 for (
auto vertexPair : std::array<std::pair<int, int>, 3> {{ { 0, 1}, {1, 2}, {2, 0} }} )
1500 bool isIntersection =
false;
1504 const double intersectionZ = ring->
zAt( vertexPair.first ) + ( ring->
zAt( vertexPair.second ) - ring->
zAt( vertexPair.first ) ) * fraction;
1505 minZ = std::min( minZ, intersectionZ );
1506 maxZ = std::max( maxZ, intersectionZ );
1510 if ( !intersectionPoint.
isEmpty() )
1513 mResults->mRawPoints.append( intersectionPoint );
1514 mResults->minZ = std::min( mResults->minZ, minZ );
1515 mResults->maxZ = std::max( mResults->maxZ, maxZ );
1517 const double distance = mProfileCurveEngine->lineLocatePoint( intersectionPoint, &lastError );
1519 crossSectionParts.append(
QgsGeometry(
new QgsLineString( QVector< double > {distance, distance}, QVector< double > {minZ, maxZ} ) ) );
1521 mResults->mDistanceToHeightMap.insert( distance, minZ );
1522 mResults->mDistanceToHeightMap.insert( distance, maxZ );
1529 QgsDebugError( QStringLiteral(
"Collinear triangles with curved profile lines are not supported yet" ) );
1535 if ( mProfileBufferedCurveEngine->intersects( triangle, &error ) )
1537 std::unique_ptr< QgsAbstractGeometry > intersection( mProfileBufferedCurveEngine->intersection( triangle, &error ) );
1538 processTriangleLineIntersect( triangle, intersection.get(), transformedParts, crossSectionParts );
1548 while ( ! mFeedback->isCanceled() && it.
nextFeature( feature ) )
1550 if ( !mProfileBufferedCurveEngine->intersects( feature.
geometry().
constGet() ) )
1557 QVector< QgsGeometry > transformedParts;
1558 QVector< QgsGeometry > crossSectionParts;
1559 bool wasCollinear =
false;
1564 if ( mProfileBufferedCurveEngine->intersects( *it ) )
1566 if (
const QgsCurvePolygon *curvePolygon = qgsgeometry_cast< const QgsCurvePolygon * >( *it ) )
1568 processPolygon( curvePolygon, transformedParts, crossSectionParts, offset, wasCollinear );
1570 else if (
const QgsPolyhedralSurface *polySurface = qgsgeometry_cast< const QgsPolyhedralSurface * >( *it ) )
1572 for (
int i = 0; i < polySurface->numPatches(); ++i )
1574 const QgsPolygon *polygon = polySurface->patchN( i );
1575 if ( mProfileBufferedCurveEngine->intersects( polygon ) )
1577 processPolygon( polygon, transformedParts, crossSectionParts, offset, wasCollinear );
1583 QgsDebugError( QStringLiteral(
"Unhandled Geometry type: %1" ).arg( ( *it )->wktTypeStr() ) );
1588 if ( mFeedback->isCanceled() )
1595 if ( !crossSectionParts.empty() )
1597 if ( !wasCollinear )
1618 mResults->features[resultFeature.
featureId].append( resultFeature );
1623double QgsVectorLayerProfileGenerator::terrainHeight(
double x,
double y )
1625 if ( !mTerrainProvider )
1626 return std::numeric_limits<double>::quiet_NaN();
1636 return std::numeric_limits<double>::quiet_NaN();
1639 return mTerrainProvider->heightAt( x, y );
1642double QgsVectorLayerProfileGenerator::featureZToHeight(
double x,
double y,
double z,
double offset )
1644 switch ( mClamping )
1652 const double terrainZ = terrainHeight( x, y );
1653 if ( !std::isnan( terrainZ ) )
1655 switch ( mClamping )
1658 if ( std::isnan( z ) )
1676 return ( std::isnan( z ) ? 0 : z ) * mScale + offset;
1679void QgsVectorLayerProfileGenerator::clampAltitudes(
QgsLineString *lineString,
const QgsPoint ¢roid,
double offset )
1683 if ( mFeedback->isCanceled() )
1686 double terrainZ = 0;
1687 switch ( mClamping )
1696 pt.
setX( lineString->
xAt( i ) );
1697 pt.
setY( lineString->
yAt( i ) );
1701 pt.
set( centroid.
x(), centroid.
y() );
1705 terrainZ = terrainHeight( pt.
x(), pt.
y() );
1715 switch ( mClamping )
1719 geomZ = lineString->
zAt( i );
1726 const double z = ( terrainZ + ( std::isnan( geomZ ) ? 0 : geomZ ) ) * mScale + offset;
1727 lineString->
setZAt( i, z );
1731bool QgsVectorLayerProfileGenerator::clampAltitudes(
QgsPolygon *polygon,
double offset )
1733 if ( !polygon->
is3D() )
1748 QgsLineString *lineString = qgsgeometry_cast<QgsLineString *>( curve );
1752 clampAltitudes( lineString, centroid, offset );
1756 if ( mFeedback->isCanceled() )
1760 QgsLineString *lineString = qgsgeometry_cast<QgsLineString *>( curve );
1764 clampAltitudes( lineString, centroid, offset );
The Qgis class provides global constants for use throughout the application.
@ Relative
Elevation is relative to terrain height (final elevation = terrain elevation + feature elevation)
@ Terrain
Elevation is clamped to terrain (final elevation = terrain elevation)
@ Absolute
Elevation is taken directly from feature and is independent of terrain height (final elevation = feat...
VertexType
Types of vertex.
@ Round
Use rounded joins.
@ Centroid
Clamp just centroid of feature.
@ Vertex
Clamp every vertex of feature.
@ Flat
Flat cap (in line with start/end of line)
@ ContinuousSurface
The features should be treated as representing values on a continuous surface (eg contour lines)
@ IndividualFeatures
Treat each feature as an individual object (eg buildings)
ProfileExportType
Types of export for elevation profiles.
@ Profile2D
Export profiles as 2D profile lines, with elevation stored in exported geometry Y dimension and dista...
@ Features3D
Export profiles as 3D features, with elevation values stored in exported geometry Z values.
@ DistanceVsElevationTable
Export profiles as a table of sampled distance vs elevation values.
@ Reverse
Reverse/inverse transform (from destination to source)
Abstract base class for all geometries.
bool is3D() const
Returns true if the geometry is 3D and contains a z-value.
virtual QgsPoint centroid() const
Returns the centroid of the geometry.
Abstract base class for objects which generate elevation profiles.
Abstract base class for storage of elevation profiles.
Abstract base class for objects which generate elevation profiles which represent a continuous surfac...
std::unique_ptr< QgsLineSymbol > mLineSymbol
Qgis::ProfileSurfaceSymbology mSymbology
std::unique_ptr< QgsFillSymbol > mFillSymbol
std::unique_ptr< QgsFillSymbol > mFillSymbol
QMap< double, double > mDistanceToHeightMap
std::unique_ptr< QgsLineSymbol > mLineSymbol
void renderResults(QgsProfileRenderContext &context) override
Renders the results to the specified context.
void copyPropertiesFromGenerator(const QgsAbstractProfileGenerator *generator) override
Copies properties from specified generator to the results object.
QVector< QgsAbstractProfileResults::Feature > asFeatures(Qgis::ProfileExportType type, QgsFeedback *feedback=nullptr) const override
Returns a list of features representing the calculated elevation results.
QVector< QgsProfileIdentifyResults > identify(const QgsProfilePoint &point, const QgsProfileIdentifyContext &context) override
Identify results visible at the specified profile point.
QgsProfileSnapResult snapPoint(const QgsProfilePoint &point, const QgsProfileSnapContext &context) override
Snaps a point to the generated elevation profile.
double valueAsDouble(int key, const QgsExpressionContext &context, double defaultValue=0.0, bool *ok=nullptr) const
Calculates the current value of the property with the specified key and interprets it as a double.
Custom exception class for Coordinate Reference System related exceptions.
Curve polygon geometry type.
int numInteriorRings() const
Returns the number of interior rings contained with the curve polygon.
bool addZValue(double zValue=0) override
Adds a z-dimension to the geometry, initialized to a preset value.
const QgsCurve * exteriorRing() const
Returns the curve polygon's exterior ring.
const QgsCurve * interiorRing(int i) const
Retrieves an interior ring from the curve polygon.
Abstract base class for curved geometry type.
virtual int numPoints() const =0
Returns the number of points in the curve.
virtual bool pointAt(int node, QgsPoint &point, Qgis::VertexType &type) const =0
Returns the point and vertex id of a point within the curve.
QgsRange which stores a range of double values.
RAII class to pop scope from an expression context on destruction.
void setFeature(const QgsFeature &feature)
Convenience function for setting a feature for the context.
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 & setCoordinateTransform(const QgsCoordinateTransform &transform)
Sets the coordinate transform which will be used to transform the feature's geometries.
QgsFeatureRequest & setFilterFids(const QgsFeatureIds &fids)
Sets the feature IDs that should be fetched.
QgsFeatureRequest & setSubsetOfAttributes(const QgsAttributeList &attrs)
Set a subset of attributes that will be fetched.
QgsFeatureRequest & setDestinationCrs(const QgsCoordinateReferenceSystem &crs, const QgsCoordinateTransformContext &context)
Sets the destination crs for feature's geometries.
void setFeedback(QgsFeedback *feedback)
Attach a feedback object that can be queried regularly by the iterator to check if it should be cance...
QgsFeatureRequest & setFilterRect(const QgsRectangle &rectangle)
Sets the rectangle from which features will be taken.
QgsFeatureRequest & setDistanceWithin(const QgsGeometry &geometry, double distance)
Sets a reference geometry and a maximum distance from this geometry to retrieve features within.
The feature class encapsulates a single feature including its unique ID, geometry and a list of field...
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.
A fill symbol type, for rendering Polygon and MultiPolygon geometries.
static double pointFractionAlongLine(double x1, double y1, double x2, double y2, double px, double py)
Given the line (x1, y1) to (x2, y2) and a point (px, py) returns the fraction of the line length at w...
static bool pointsAreCollinear(double x1, double y1, double x2, double y2, double x3, double y3, double epsilon)
Given the points (x1, y1), (x2, y2) and (x3, y3) returns true if these points can be considered colli...
static bool segmentIntersection(const QgsPoint &p1, const QgsPoint &p2, const QgsPoint &q1, const QgsPoint &q2, QgsPoint &intersectionPoint, bool &isIntersection, double tolerance=1e-8, bool acceptImproperIntersection=false)
Compute the intersection between two segments.
A geometry is the spatial representation of a feature.
static QgsGeometry fromRect(const QgsRectangle &rect)
Creates a new geometry from a QgsRectangle.
QgsAbstractGeometry::const_part_iterator const_parts_begin() const
Returns STL-style const iterator pointing to the first part of the geometry.
bool boundingBoxIntersects(const QgsRectangle &rectangle) const
Returns true if the bounding box of this geometry intersects with a rectangle.
static QgsGeometry collectGeometry(const QVector< QgsGeometry > &geometries)
Creates a new multipart geometry from a list of QgsGeometry objects.
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.
static QgsGeometry fromPointXY(const QgsPointXY &point)
Creates a new geometry from a QgsPointXY object.
QgsGeometry mergeLines() const
Merges any connected lines in a LineString/MultiLineString geometry and converts them to single line ...
bool isEmpty() const
Returns true if the geometry is empty (eg a linestring with no vertices, or a collection with no geom...
QgsAbstractGeometry::const_part_iterator const_parts_end() const
Returns STL-style iterator pointing to the imaginary part after the last part of the geometry.
static QgsGeometry unaryUnion(const QVector< QgsGeometry > &geometries, const QgsGeometryParameters ¶meters=QgsGeometryParameters())
Compute the unary union on a list of geometries.
Qgis::GeometryOperationResult translate(double dx, double dy, double dz=0.0, double dm=0.0)
Translates this geometry by dx, dy, dz and dm.
Does vector analysis using the geos library and handles import, export, exception handling*.
Line string geometry type, with support for z-dimension and m-values.
const double * yData() const
Returns a const pointer to the y vertex data.
const double * xData() const
Returns a const pointer to the x vertex data.
const double * zData() const
Returns a const pointer to the z vertex data, or nullptr if the linestring does not have z values.
int numPoints() const override
Returns the number of points in the curve.
QgsPoint pointN(int i) const
Returns the specified point from inside the line string.
int nCoordinates() const override
Returns the number of nodes contained in the geometry.
double yAt(int index) const override
Returns the y-coordinate of the specified node in the line string.
void setZAt(int index, double z)
Sets the z-coordinate of the specified node in the line string.
bool deleteVertex(QgsVertexId position) override
Deletes a vertex within the geometry.
QgsLineString * clone() const override
Clones the geometry by performing a deep copy.
double zAt(int index) const override
Returns the z-coordinate of the specified node in the line string.
double xAt(int index) const override
Returns the x-coordinate of the specified node in the line string.
A line symbol type, for rendering LineString and MultiLineString geometries.
@ ExtrusionHeight
Extrusion height.
A marker symbol type, for rendering Point and MultiPoint geometries.
Multi line string geometry collection.
Multi point geometry collection.
Multi polygon geometry collection.
A class to represent a 2D point.
void setY(double y)
Sets the y value of the point.
void set(double x, double y)
Sets the x and y value of the point.
void setX(double x)
Sets the x value of the point.
Point geometry type, with support for z-dimension and m-values.
QgsPoint * clone() const override
Clones the geometry by performing a deep copy.
bool isEmpty() const override
Returns true if the geometry is empty.
Polyhedral surface geometry type.
Encapsulates the context in which an elevation profile is to be generated.
Encapsulates the context of identifying profile results.
double maximumPointElevationDelta
Maximum allowed snapping delta for the elevation values when identifying a point.
double maximumPointDistanceDelta
Maximum allowed snapping delta for the distance values when identifying a point.
double maximumSurfaceElevationDelta
Maximum allowed snapping delta for the elevation values when identifying a continuous elevation surfa...
Stores identify results generated by a QgsAbstractProfileResults object.
Encapsulates a point on a distance-elevation profile.
double elevation() const
Returns the elevation of the point.
double distance() const
Returns the distance of the point.
Abstract base class for storage of elevation profiles.
const QTransform & worldTransform() const
Returns the transform from world coordinates to painter coordinates.
QgsDoubleRange elevationRange() const
Returns the range of elevations to include in the render.
QgsDoubleRange distanceRange() const
Returns the range of distances to include in the render.
QgsRenderContext & renderContext()
Returns a reference to the component QgsRenderContext.
Encapsulates properties and constraints relating to fetching elevation profiles from different source...
Encapsulates the context of snapping a profile point.
double maximumPointDistanceDelta
Maximum allowed snapping delta for the distance values when snapping to a point.
double maximumSurfaceElevationDelta
Maximum allowed snapping delta for the elevation values when snapping to a continuous elevation surfa...
double maximumPointElevationDelta
Maximum allowed snapping delta for the elevation values when snapping to a point.
Encapsulates results of snapping a profile point.
QgsProfilePoint snappedPoint
Snapped point.
bool prepare(const QgsExpressionContext &context=QgsExpressionContext()) const final
Prepares the collection against a specified expression context.
QSet< QString > referencedFields(const QgsExpressionContext &context=QgsExpressionContext(), bool ignoreContext=false) const final
Returns the set of any fields referenced by the active properties from the collection.
T lower() const
Returns the lower bound of the range.
T upper() const
Returns the upper bound of the range.
A rectangle specified with double values.
double xMinimum() const
Returns the x minimum value (left side of rectangle).
bool intersects(const QgsRectangle &rect) const
Returns true when rectangle intersects with other rectangle.
double yMinimum() const
Returns the y minimum value (bottom side of rectangle).
double xMaximum() const
Returns the x maximum value (right side of rectangle).
QPainter * painter()
Returns the destination QPainter for the render operation.
QgsExpressionContext & expressionContext()
Gets the expression context.
Scoped object for saving and restoring a QPainter object's state.
Abstract base class for all rendered symbols.
void setColor(const QColor &color) const
Sets the color for the symbol.
qreal opacity() const
Returns the opacity for the symbol.
QColor color() const
Returns the symbol's color.
Qgis::SymbolType type() const
Returns the symbol's type.
Class that takes care of tessellation of polygons into triangles.
Vector layer specific subclass of QgsMapLayerElevationProperties.
Partial snapshot of vector layer's state (only the members necessary for access to features)
Implementation of QgsAbstractProfileGenerator for vector layers.
QgsAbstractProfileResults * takeResults() override
Takes results from the generator.
bool generateProfile(const QgsProfileGenerationContext &context=QgsProfileGenerationContext()) override
Generate the profile (based on data stored in the class).
QString sourceId() const override
Returns a unique identifier representing the source of the profile.
~QgsVectorLayerProfileGenerator() override
QgsFeedback * feedback() const override
Access to feedback object of the generator (may be nullptr)
QgsVectorLayerProfileGenerator(QgsVectorLayer *layer, const QgsProfileRequest &request)
Constructor for QgsVectorLayerProfileGenerator.
QVector< QgsProfileIdentifyResults > identify(const QgsProfilePoint &point, const QgsProfileIdentifyContext &context) override
Identify results visible at the specified profile point.
Qgis::VectorProfileType profileType
std::unique_ptr< QgsMarkerSymbol > mMarkerSymbol
bool respectLayerSymbology
bool mShowMarkerSymbolInSurfacePlots
QVector< QgsGeometry > asGeometries() const override
Returns a list of geometries representing the calculated elevation results.
QPointer< QgsVectorLayer > mLayer
void renderResults(QgsProfileRenderContext &context) override
Renders the results to the specified context.
QVector< QgsAbstractProfileResults::Feature > asFeatures(Qgis::ProfileExportType type, QgsFeedback *feedback=nullptr) const override
Returns a list of features representing the calculated elevation results.
QgsProfileSnapResult snapPoint(const QgsProfilePoint &point, const QgsProfileSnapContext &context) override
Snaps a point to the generated elevation profile.
QString type() const override
Returns the unique string identifier for the results type.
QHash< QgsFeatureId, QVector< Feature > > features
void copyPropertiesFromGenerator(const QgsAbstractProfileGenerator *generator) override
Copies properties from specified generator to the results object.
Represents a vector layer which manages a vector based data sets.
QgsMapLayerElevationProperties * elevationProperties() override
Returns the layer's elevation properties.
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...
#define BUILTIN_UNREACHABLE
bool qgsDoubleNear(double a, double b, double epsilon=4 *std::numeric_limits< double >::epsilon())
Compare two doubles (but allow some difference)
QSet< QgsFeatureId > QgsFeatureIds
qint64 QgsFeatureId
64 bit feature ids negative numbers are used for uncommitted/newly added features
#define QgsDebugMsgLevel(str, level)
#define QgsDebugError(str)
const QgsCoordinateReferenceSystem & crs
Encapsulates information about a feature exported from the profile results.
QString layerIdentifier
Identifier for grouping output features.
QVariantMap attributes
Exported attributes.
QgsGeometry geometry
Exported geometry.
QgsGeometry crossSectionGeometry
Cross section distance vs height geometry for feature.
QgsFeatureId featureId
Original feature ID.
QgsGeometry geometry
Feature's geometry with any terrain height adjustment and extrusion applied.
Utility class for identifying a unique vertex within a geometry.