sharpetronics/QGIS
1
0
Fork 0
mirror of https://github.com/qgis/QGIS.git synced 2025-03-12 00:02:25 -04:00
Code Issues Packages Projects Releases Wiki Activity

optimize nearest neighbour analysis by storing geometries in the spatial

Browse Source 
index and using index to retrieve them instead of making single-feature
requests
This commit is contained in:
Alexander Bruy 2019-12-10 08:57:17 +02:00
parent 363d8d053c
commit bd72e02e55
1 changed files with 2 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/analysis/processing/qgsalgorithmnearestneighbouranalysis.cpp
View File

@ -86,7 +86,7 @@ QVariantMap QgsNearestNeighbourAnalysisAlgorithm::processAlgorithm( const QVaria
QString outputFile = parameterAsFileOutput( parameters, QStringLiteral( "OUTPUT_HTML_FILE" ), context );
QgsSpatialIndex spatialIndex( *source, feedback );
QgsSpatialIndex spatialIndex( *source, feedback, QgsSpatialIndex::FlagStoreFeatureGeometries );
QgsDistanceArea da;
da.setSourceCrs( source->sourceCrs(), context.transformContext() );
da.setEllipsoid( context.project()->ellipsoid() );
@ -109,9 +109,7 @@ QVariantMap QgsNearestNeighbourAnalysisAlgorithm::processAlgorithm( const QVaria
}
QgsFeatureId neighbourId = spatialIndex.nearestNeighbor( f.geometry().asPoint(), 2 ).at( 1 );
request.setFilterFid( neighbourId ).setSubsetOfAttributes( QList< int >() );
source->getFeatures( request ).nextFeature( neighbour );
sumDist += da.measureLine( neighbour.geometry().asPoint(), f.geometry().asPoint() );
sumDist += da.measureLine( spatialIndex.geometry( neighbourId ).asPoint(), f.geometry().asPoint() );
i++;
feedback->setProgress( i * step );