Skip to content

Commit

Permalink
Rename variables for clarity
Browse files Browse the repository at this point in the history
  • Loading branch information
nyalldawson committed Oct 15, 2020
1 parent 05645b3 commit f3a5257
Showing 1 changed file with 14 additions and 14 deletions.
28 changes: 14 additions & 14 deletions src/analysis/processing/qgsalgorithmsplitwithlines.cpp
Expand Up @@ -109,7 +109,7 @@ QVariantMap QgsSplitWithLinesAlgorithm::processAlgorithm( const QVariantMap &par
QgsFeatureIterator splitLines = linesSource->getFeatures( request );
QgsFeature aSplitFeature;

QgsSpatialIndex spatialIndex( splitLines, feedback, QgsSpatialIndex::FlagStoreFeatureGeometries );
const QgsSpatialIndex splitLinesIndex( splitLines, feedback, QgsSpatialIndex::FlagStoreFeatureGeometries );

QgsFeature outFeat;
QgsFeatureIterator features = source->getFeatures();
Expand All @@ -131,35 +131,35 @@ QVariantMap QgsSplitWithLinesAlgorithm::processAlgorithm( const QVariantMap &par
continue;
}

QgsGeometry inGeom = inFeatureA.geometry();
const QgsGeometry originalGeometry = inFeatureA.geometry();
outFeat.setAttributes( inFeatureA.attributes() );

QVector< QgsGeometry > inGeoms = inGeom.asGeometryCollection();
QVector< QgsGeometry > inGeoms = originalGeometry.asGeometryCollection();

const QgsFeatureIds lines = qgis::listToSet( spatialIndex.intersects( inGeom.boundingBox() ) );
if ( !lines.empty() ) // has intersection of bounding boxes
const QgsFeatureIds splitLineCandidates = qgis::listToSet( splitLinesIndex.intersects( originalGeometry.boundingBox() ) );
if ( !splitLineCandidates.empty() ) // has intersection of bounding boxes
{
QVector< QgsGeometry > splittingLines;

// use prepared geometries for faster intersection tests
std::unique_ptr< QgsGeometryEngine > engine;
std::unique_ptr< QgsGeometryEngine > originalGeometryEngine;

for ( QgsFeatureId line : lines )
for ( QgsFeatureId splitLineCandidateId : splitLineCandidates )
{
// check if trying to self-intersect
if ( sameLayer && inFeatureA.id() == line )
if ( sameLayer && inFeatureA.id() == splitLineCandidateId )
continue;

QgsGeometry splitGeom = spatialIndex.geometry( line );
if ( !engine )
const QgsGeometry splitLineCandidate = splitLinesIndex.geometry( splitLineCandidateId );
if ( !originalGeometryEngine )
{
engine.reset( QgsGeometry::createGeometryEngine( inGeom.constGet() ) );
engine->prepareGeometry();
originalGeometryEngine.reset( QgsGeometry::createGeometryEngine( originalGeometry.constGet() ) );
originalGeometryEngine->prepareGeometry();
}

if ( engine->intersects( splitGeom.constGet() ) )
if ( originalGeometryEngine->intersects( splitLineCandidate.constGet() ) )
{
QVector< QgsGeometry > splitGeomParts = splitGeom.asGeometryCollection();
QVector< QgsGeometry > splitGeomParts = splitLineCandidate.asGeometryCollection();
splittingLines.append( splitGeomParts );
}
}
Expand Down

0 comments on commit f3a5257

Please sign in to comment.