QGIS【实例 01】筛选范围内的数据相交intersection工具使用及Feature (347897) has invalid geometry问题处理

在C++和QGIS（一个开源的地理信息系统软件）环境下，计算两个线图层的交叉点可以通过以下步骤实现： 1. **加载线图层**：使用QGIS的API加载两个线图层。 2. **创建交叉点图层**：创建一个新的点图层用于存储交叉点。 3. **遍历线要素**：遍历两个线图层中的每一个线要素，计算它们的交叉点。 4. **添加交叉点到新图层**：将计算得到的交叉点添加到新创建的点图层中。 以下是一个简单的示例代码，展示了如何在C++中使用QGIS API计算两个线图层的交叉点： ```cpp #include <QgsVectorLayer.h> #include <QgsFeature.h> #include <QgsGeometry.h> #include <QgsFeatureIterator.h> #include <QgsProject.h> #include <QgsVectorFileWriter.h> bool calculateIntersections(const QString& layer1Path, const QString& layer2Path, const QString& outputPath) { // 加载线图层 QgsVectorLayer* layer1 = new QgsVectorLayer(layer1Path, "Layer 1", "ogr"); QgsVectorLayer* layer2 = new QgsVectorLayer(layer2Path, "Layer 2", "ogr"); if (!layer1->isValid() || !layer2->isValid()) { qDebug("One of the layers is invalid"); return false; } // 创建交叉点图层 QgsVectorLayer* intersectionLayer = new QgsVectorLayer("Point", "Intersections", "memory"); QgsFeatureList features; // 遍历第一个图层的要素 QgsFeatureIterator it1 = layer1->getFeatures(); QgsFeature f1; while (it1.nextFeature(f1)) { QgsGeometry geom1 = f1.geometry(); // 遍历第二个图层的要素 QgsFeatureIterator it2 = layer2->getFeatures(); QgsFeature f2; while (it2.nextFeature(f2)) { QgsGeometry geom2 = f2.geometry(); // 计算交叉点 QgsGeometry intersectGeom = geom1.intersection(geom2); if (intersectGeom.isNull()) { continue; } QgsFeature newFeature; newFeature.setGeometry(intersectGeom); features.append(newFeature); } } intersectionLayer->dataProvider()->addFeatures(features); // 保存交叉点图层 QgsVectorFileWriter::writeAsVectorFormat(intersectionLayer, outputPath, "UTF-8", intersectionLayer->crs(), "ESRI Shapefile"); return true; } int main(int argc, char *argv[]) { QString layer1Path = "path/to/layer1.shp"; QString layer2Path = "path/to/layer2.shp"; QString outputPath = "path/to/output.shp"; if (calculateIntersections(layer1Path, layer2Path, outputPath)) { qDebug("Intersections calculated successfully"); } else { qDebug("Failed to calculate intersections"); } return 0; } ``` ### 代码说明： 1. **加载线图层**：使用`QgsVectorLayer`加载两个线图层。 2. **创建交叉点图层**：创建一个新的内存中的点图层。 3. **遍历线要素**：使用`QgsFeatureIterator`遍历两个图层的要素，并计算它们的交叉点。 4. **添加交叉点到新图层**：将计算得到的交叉点添加到新图层中。 5. **保存交叉点图层**：使用`QgsVectorFileWriter`将交叉点图层保存为Shapefile文件。