Branch data     Line data    Source code

       1                 :            : /***************************************************************************
       2                 :            :     qgsgeometrylineintersectioncheck.cpp
       3                 :            :     ---------------------
       4                 :            :     begin                : June 2017
       5                 :            :     copyright            : (C) 2017 by Sandro Mani / Sourcepole AG
       6                 :            :     email                : smani at sourcepole dot ch
       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                 :            :  ***************************************************************************/
      15                 :            : 
      16                 :            : #include "qgsgeometrycheckcontext.h"
      17                 :            : #include "qgsgeometrylineintersectioncheck.h"
      18                 :            : #include "qgslinestring.h"
      19                 :            : #include "qgsvectorlayer.h"
      20                 :            : #include "qgsgeometrycheckerror.h"
      21                 :            : 
      22                 :          1 : void QgsGeometryLineIntersectionCheck::collectErrors( const QMap<QString, QgsFeaturePool *> &featurePools, QList<QgsGeometryCheckError *> &errors, QStringList &messages, QgsFeedback *feedback, const LayerFeatureIds &ids ) const
      23                 :            : {
      24                 :          1 :   Q_UNUSED( messages )
      25                 :            : 
      26                 :          1 :   QMap<QString, QgsFeatureIds> featureIds = ids.isEmpty() ? allLayerFeatureIds( featurePools ) : ids.toMap();
      27                 :          1 :   QgsGeometryCheckerUtils::LayerFeatures layerFeaturesA( featurePools, featureIds, compatibleGeometryTypes(), feedback, mContext, true );
      28                 :          1 :   QList<QString> layerIds = featureIds.keys();
      29                 :         10 :   for ( const QgsGeometryCheckerUtils::LayerFeature &layerFeatureA : layerFeaturesA )
      30                 :            :   {
      31                 :            :     // Ensure each pair of layers only gets compared once: remove the current layer from the layerIds, but add it to the layerList for layerFeaturesB
      32                 :          9 :     layerIds.removeOne( layerFeatureA.layer()->id() );
      33                 :            : 
      34                 :          9 :     const QgsAbstractGeometry *geom = layerFeatureA.geometry().constGet();
      35                 :         19 :     for ( int iPart = 0, nParts = geom->partCount(); iPart < nParts; ++iPart )
      36                 :            :     {
      37                 :         10 :       const QgsLineString *line = dynamic_cast<const QgsLineString *>( QgsGeometryCheckerUtils::getGeomPart( geom, iPart ) );
      38                 :         10 :       if ( !line )
      39                 :            :       {
      40                 :            :         // Should not happen
      41                 :          0 :         continue;
      42                 :            :       }
      43                 :            : 
      44                 :            :       // Check whether the line intersects with any other lines
      45                 :         10 :       QgsGeometryCheckerUtils::LayerFeatures layerFeaturesB( featurePools, QList<QString>() << layerFeatureA.layer()->id() << layerIds, line->boundingBox(), {QgsWkbTypes::LineGeometry}, mContext );
      46                 :         28 :       for ( const QgsGeometryCheckerUtils::LayerFeature &layerFeatureB : layerFeaturesB )
      47                 :            :       {
      48                 :            :         // > : only report intersections within same layer once
      49                 :         18 :         if ( layerFeatureA.layer()->id() == layerFeatureB.layer()->id() && layerFeatureB.feature().id() > layerFeatureA.feature().id() )
      50                 :            :         {
      51                 :          3 :           continue;
      52                 :            :         }
      53                 :            : 
      54                 :         15 :         const QgsAbstractGeometry *testGeom = layerFeatureB.geometry().constGet();
      55                 :         35 :         for ( int jPart = 0, mParts = testGeom->partCount(); jPart < mParts; ++jPart )
      56                 :            :         {
      57                 :            :           // Skip current feature part, only report intersections within same part once
      58                 :         20 :           if ( layerFeatureB.feature().id() == layerFeatureA.feature().id() && iPart >= jPart )
      59                 :            :           {
      60                 :         11 :             continue;
      61                 :            :           }
      62                 :          9 :           const QgsLineString *testLine = dynamic_cast<const QgsLineString *>( QgsGeometryCheckerUtils::getGeomPart( testGeom, jPart ) );
      63                 :          9 :           if ( !testLine )
      64                 :            :           {
      65                 :          0 :             continue;
      66                 :            :           }
      67                 :          9 :           const QList< QgsPoint > intersections = QgsGeometryCheckerUtils::lineIntersections( line, testLine, mContext->tolerance );
      68                 :         10 :           for ( const QgsPoint &inter : intersections )
      69                 :            :           {
      70                 :          1 :             errors.append( new QgsGeometryCheckError( this, layerFeatureA, inter, QgsVertexId( iPart ), layerFeatureB.id() ) );
      71                 :            :           }
      72                 :          9 :         }
      73                 :            :       }
      74                 :         10 :     }
      75                 :            :   }
      76                 :          1 : }
      77                 :            : 
      78                 :          0 : void QgsGeometryLineIntersectionCheck::fixError( const QMap<QString, QgsFeaturePool *> &featurePools, QgsGeometryCheckError *error, int method, const QMap<QString, int> & /*mergeAttributeIndices*/, Changes & /*changes*/ ) const
      79                 :            : {
      80                 :          0 :   Q_UNUSED( featurePools )
      81                 :            : 
      82                 :          0 :   if ( method == NoChange )
      83                 :            :   {
      84                 :          0 :     error->setFixed( method );
      85                 :          0 :   }
      86                 :            :   else
      87                 :            :   {
      88                 :          0 :     error->setFixFailed( tr( "Unknown method" ) );
      89                 :            :   }
      90                 :          0 : }
      91                 :            : 
      92                 :          0 : QStringList QgsGeometryLineIntersectionCheck::resolutionMethods() const
      93                 :            : {
      94                 :          0 :   static QStringList methods = QStringList() << tr( "No action" );
      95                 :          0 :   return methods;
      96                 :          0 : }
      97                 :            : 
      98                 :          0 : QgsGeometryCheck::CheckType QgsGeometryLineIntersectionCheck::factoryCheckType()
      99                 :            : {
     100                 :          0 :   return QgsGeometryCheck::FeatureNodeCheck;
     101                 :            : }