OGR根据FID删除矢量中的某个要素（多边形）

最新推荐文章于 2025-10-25 15:31:26 发布
原创 最新推荐文章于 2025-10-25 15:31:26 发布 · 3.2k 阅读 · 3 ·
CC 4.0 BY-SA版权
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
文章标签：

#null #layer #存储

本文介绍了一个函数，该函数通过OGR库实现复制矢量文件并删除指定FID的要素。输入参数包括源矢量文件 DataSource、目标输出文件名、待删除FID的数组以及结果文件格式。

这是一个函数，实现复制原来的矢量文件，然后删除掉FID与对应数组pFID中相等的要素，参数1为原来的矢量文件DataSource指针，参数2为输出的矢量文件名，参数3为存储要删除的FID的数组，参数4为结果矢量的格式。

bool CreateResultSHP (OGRDataSource  *poSrcDs ,std ::string outputFileName ,vector < int >  *pFID , const  char * pszFormat )
{
  //创建输出矢量文件
 OGRSFDriver  *poDriver ;
  if  (pszFormat  ==  NULL )
  {
  pszFormat  =  "ESRI Shapefile" ;
  }

 poDriver  = OGRSFDriverRegistrar ::GetRegistrar ()->GetDriverByName ( pszFormat  );
  if  (poDriver  ==  NULL )&
最低0.47元/天 解锁文章
使用GDAL删除矢量文件中的要素的方法
CyberBladeX的博客
09-07 357
使用GDAL删除矢量文件中的要素的方法GDAL（地理空间数据抽象库）是一个功能强大的开源库，用于处理和操作地理空间数据。它提供了许多功能，包括读取、写入和编辑矢量文件。在本文中，我们将介绍如何使用GDAL库删除矢量文件中的要素。在开始之前，请确保已经安装了GDAL库，并且熟悉Python编程语言。首先，我们需要导入必要的模块和函数。在Python中，我们可以使用gdal模块和ogr模块来操作矢量数据。接下来，我们需要打开要进行操作的矢量文件。我们可以使用ogr.Open()
ogr2ogr矢量数据要素格式转换
山路十八弯，走过多少遍！
07-06 1689
在文件格式之间转换简单要素数据。 概要 ogr2ogr [--help-general] [-skipfailures] [-append] [-update] [-select field_list] [-where restricted_where|\@filename] [-progress] [-sql <sql statement>|\@filename] [-dialect dialect] [-preserve_fi...
ogrfeature必须销毁，否则会内存泄漏
coder
04-19 405
上周末，发现多个shp文件操作后，竟然崩溃， 今天查了一天，原以为是std::vector没有释放掉，最后发现原来是 OGRFeature * pFeature = poLayer->GetFeature(featureID);后，会载入内存，但是不销毁的话，关掉shp文件也会驻留在内存里。 因此，需要再加上这句 OGRFeature::DestroyFeature(pFeature); 一切平静了 ...
OGRFeature的DestroyFeature方法
weixin_33933118的博客
07-31 409
 Ogr的销毁DestroyFeature方法： 1 void OGRFeature::DestroyFeature( OGRFeature *poFeature ) 2 { 3 delete poFeature; 4 } 析构函数中删除了几何对象和空间参考定义， ogrfeature.cpp文件中查看： 1 OGRFeature::~OGRFeature()...
Python地理数据处理 七：使用OGR处理几何要素
amyniez的博客，欢迎交流讨论
02-04 3180
目录1. 点集1.1 创建编辑一个点1.2 创建编辑多个点 1. 点集 1.1 创建编辑一个点   OGR中带有额外z坐标的几何要素被认为是2.5D，而不是3D。再OGR进行空间操作时，不考虑z值。   OGR常量表示的不同几何类型： firepit = ogr.Geometry(ogr.wkbPoint) firepit.AddPoint(10, 21) x, y = firepit.GetX(), firepit.GetY() print('{}, {}'.format(x, y)) print(
利用OGR处理几何要素
Mikowoo007的博客
10-05 1220
OGR常量表示的不同几何类型 虽然许多几何要素只存在于二维(2D)的直角坐标系平面(X和Y坐标)，也有三维(3D)带有Z值得几何对象。这些Z值通常用来表示高程，但也可以用于表示其它数据，如每年温度得最大值。从技术上讲，这些几何要素在OGR中被认为是2.5D，而不是3D，因为OGR进行空间操作时，不考虑Z值。注意，虽然可以将Z值添加到二维几何结构中，但将数据写入文件时，Z值将被忽略。 from osgeo import ogr ogr.Geometry(ogr.wkbPoint) 几何类型 2D常量
GDAL中矢量的读写（多边形）
06-21 2071
GDAL中矢量线读写主要思路是创建图层，再创建对象 下面粘贴上矢量数据的读写 #include <iostream> #include "ogrsf_frmts.h" using namespace std; bool ReadDXF(string filePath, vector<vector<OGRPoint>>& vertexPoint) { GDALDataset *poDS = (GDALDataset*)GDALOpenEx(filePa
4、矢量数据处理：ogrinfo 与 ogr2ogr 工具详解
最新发布
sky77的博客
10-25 16
本文详细介绍了地理空间数据处理中两个核心命令行工具——ogrinfo 和 ogr2ogr 的功能与使用方法。涵盖数据格式简介、基础与高级选项解析、实际应用场景（如格式转换、重投影、数据裁剪与过滤）、以及使用注意事项。通过丰富的示例和流程图，帮助用户快速掌握如何利用这些工具高效处理矢量数据，适用于OpenStreetMap等大规模数据集的分析与转换。
6、矢量数据编辑全攻略
hill8的博客
10-19 22
本文详细介绍了使用PyQGIS进行矢量数据编辑的完整操作流程，涵盖添加线和多边形要素、编辑属性与字段、连接CSV属性表、移动几何、删除要素与属性以及图层重投影等核心功能。通过代码示例和原理解析，帮助用户掌握GIS中常见的矢量数据处理技术，适用于地理信息系统的开发与数据分析工作。
【Python&GIS】Python处理矢量数据的基本操作（查询、修改、删除、新建）
Remote Sensing迷途小书童
08-10 1926
本文章分享了如何使用Python对矢量数据进行查询、修改、删除、新建等操作
GDAL OGR C++ API 学习之路 （4）OGRFeature 篇 代码示例
qq_69574549的博客
07-04 580
对属性字段值的列表获取
GDAl 查询并删除要素及其属性（python）
feixiegis的博客
05-22 4057
case：选择出矢量图层中符合给定条件的要素，并将其删除。这其中涉及到要素的空间查询和要素的删除等知识点。接下来将分为打开矢量数据、查询空间要素和删除空间要素三个部分做实例。注意：矢量数据需以可读写的模式打开，即driver.Open(strFilePath,1)。此外删除要素后需要利用REPARK命令实现图形和属性删除的同步。第一步：打开矢量数据#注册所有驱动 gdal.AllRegister(...
GDAL OGR C++ API 学习之路 （1）OGRFeature 篇 代码示例
qq_69574549的博客
04-29 529
类 OGRFeatureOGRFeature是OGR库中用于描述矢量数据要素的类。要素是GIS数据的基本单元，它由一组属性值和空间几何体组成。OGRFeature类可以描述各种类型的要素，例如点、线、多边形、面等OGRFeature类包含两个主要成员变量：一个是，它描述了要素的属性结构；另一个是，它描述了要素的空间几何形状。
GDAL OGR C++ API 学习之路 （1）OGRFieldDefn篇
qq_69574549的博客
04-30 1166
OGRFieldDefn是一个表示矢量数据属性字段元信息的类，它用于描述属性字段的数据类型、名称、宽度等信息需要注意的是，OGRFieldDefn只能处理属性字段，不能处理几何字段。如果要处理几何字段，需要使用OGRGeomFieldDefn类字段描述如下：一个名字。参见 SetName（） /备用名称（可选）：字段的替代描述性名称（有时称为“别名”）。请参阅 SetAlternativeName（） /
Python——使用OGR操作矢量数据
qq_45751615的博客
10-07 1136
OGR库
OGR API 手册
dayong576的专栏
11-14 3170
OGR API 手册    这篇文档是为了说明怎样从一个文件里用OGR的C++类读取和写入数据。强烈建议在读此文档之前首先阅读介绍OGR体系结构介绍文档，里面介绍了OGR里主要的类以及他们所完成的功能。从OGR 读取数据    为了描述如何通过OGR读取数据，我们编写一个小的示例，从OGR文件里读出数据，并按照一定格式输出。    首先我们需要注册我们渴望读取的所有格式。这个通过调
Shapefile属性操作之删
给永远比拿愉快
06-13 2022
文章目录开篇案例介绍代码展示方法总结 版权声明：本文为博主原创文章，转载请注明原文出处！ 作者：阿振 写作时间：2020-06-13 又是一个周末 开篇 延续上篇的《Shapfile属性操作之增》，这篇我们来聊聊如何进行属性删除。使用的工具依旧是我们的GDAL库。 案例介绍 这里我们要处理的数据是中国地图分省的矢量Shapefile，是一个面状数据。在上篇中我们给该数据添加了一个属性字段Abbr用以表示省的简称。这篇我们再把该字段给删除掉。 删除的方法包括两部：首先，从属性表中找到该字段，然后删除该字段
DeleteFeature
kone666
03-09 1766
一、  几种删除方法代码 1.  查询结果中删除   private void Delete1(IFeatureClass PFeatureclass)         {             IQueryFilter pQueryFilter = new QueryFilterClass();             pQueryFilter.WhereClause = "obje
c++ qgis二次开发 矢量编辑功能，根据图层类型，在mapcanvas上对点线面矢量进行编辑（新增，删除，移动，复制 ，剪切等功能），新增矢量要素在mapcanvas画布上勾勒，不用QgsMapToolAddFeature实现
06-11
<think>我们被要求在不使用QgsMapToolAddFeature类的情况下实现矢量要素的新增绘制功能。这意味着我们需要手动处理地图画布上的鼠标事件，并在画布上实时勾勒用户绘制的几何图形。解决方案思路：1.自定义地图工具（继承自QgsMapTool）：我们将创建一个自定义地图工具类，用于捕获鼠标事件并绘制临时几何图形。2.根据图层类型（点、线、面）实现不同的绘制逻辑：-点：单击即可创建点要素-线：通过多个点形成线，双击结束绘制-面：通过多个点形成面，双击结束绘制（首点自动闭合）3.在绘制过程中，使用RubberBand（橡皮筋）在画布上实时显示正在绘制的几何图形。步骤：1.创建自定义地图工具类（如CustomMapTool），并重写以下事件：-canvasPressEvent：鼠标按下事件-canvasMoveEvent：鼠标移动事件（用于显示橡皮筋的临时点）-canvasReleaseEvent：鼠标释放事件（用于添加点）-canvasDoubleClickEvent：双击事件（结束绘制线或面）2.在工具中根据当前图层的几何类型执行不同的绘制逻辑。3.完成绘制后，将橡皮筋的几何对象转换为要素并添加到图层。注意：要支持复制和剪切等操作，我们需要维护一个剪贴板（可以使用QList<QgsFeature>作为内存中的剪贴板）。下面我们逐步实现这些功能。一、自定义地图工具类（CustomMapTool）：我们创建一个继承自QgsMapTool的自定义工具类，用于处理绘制。二、实现其他编辑功能（删除、移动、复制、剪切）：这些功能可以通过选择要素后调用相应的函数来实现，不依赖于特定的地图工具。三、移动功能：我们也需要自定义一个移动工具，同样使用橡皮筋来显示移动过程。由于代码较长，我们将分部分展示关键代码：首先，我们实现自定义绘制工具（CustomMapTool）的框架：</think> ### C++ QGIS 二次开发：手动实现矢量编辑功能（无 `QgsMapToolAddFeature`） 以下是基于 QGIS C++ API 实现矢量编辑功能的完整方案，通过手动处理鼠标事件实现绘制逻辑，不依赖 `QgsMapToolAddFeature` 类。 --- ### 一、基础环境设置 初始化 QGIS 环境并加载矢量图层： ```cpp #include <qgsapplication.h> #include <qgsmapcanvas.h> #include <qgsvectorlayer.h> #include <qgsproject.h> #include <qgsrubberband.h> int main(int argc, char *argv[]) { QgsApplication app(argc, argv, true); QgsApplication::setPrefixPath("/path/to/qgis/installation", true); QgsApplication::initQgis(); // 创建主窗口和地图画布 QMainWindow mainWindow; QgsMapCanvas* canvas = new QgsMapCanvas(&mainWindow); canvas->setCanvasColor(Qt::white); // 加载矢量图层 [^1][^3] QgsVectorLayer* layer = new QgsVectorLayer("path/to/shapefile.shp", "Editable Layer", "ogr"); if (layer->isValid()) { QgsProject::instance()->addMapLayer(layer); canvas->setLayers({layer}); canvas->setExtent(layer->extent()); } mainWindow.setCentralWidget(canvas); mainWindow.show(); return app.exec(); } ``` --- ### 二、自定义矢量绘制工具 实现鼠标事件处理逻辑，支持点/线/面绘制： ```cpp class CustomDrawingTool : public QgsMapTool { Q_OBJECT public: CustomDrawingTool(QgsMapCanvas* canvas, QgsVectorLayer* layer) : QgsMapTool(canvas), mCanvas(canvas), mLayer(layer) { // 创建橡皮筋用于实时预览 mRubberBand = new QgsRubberBand(mCanvas, QgsWkbTypes::LineGeometry); mRubberBand->setColor(Qt::red); mRubberBand->setWidth(2); } ~CustomDrawingTool() { delete mRubberBand; } protected: // 鼠标按下事件 void canvasPressEvent(QgsMapMouseEvent* e) override { QgsPointXY mapPoint = toMapCoordinates(e->pos()); if (mLayer->geometryType() == QgsWkbTypes::PointGeometry) { // 点图层：直接创建要素 createPointFeature(mapPoint); } else { // 线/面：添加到临时点集 mPoints.append(mapPoint); if (mPoints.size() > 1) { // 绘制线段预览 QgsPointXY lastPoint = mPoints[mPoints.size() - 2]; mRubberBand->addPoint(mapPoint); } } } // 鼠标移动事件（线/面实时预览） void canvasMoveEvent(QgsMapMouseEvent* e) override { if (mPoints.isEmpty()) return; QgsPointXY currentPoint = toMapCoordinates(e->pos()); QgsPointXY lastPoint = mPoints.last(); // 更新临时预览线 mRubberBand->movePoint(mRubberBand->numberOfVertices() - 1, currentPoint); } // 鼠标释放事件（面/线完成绘制） void canvasReleaseEvent(QgsMapMouseEvent* e) override { if (e->button() == Qt::RightButton) { if (mPoints.size() >= 2) { // 双击完成绘制 createGeometryFeature(); } resetDrawing(); } } private: // 创建点要素 void createPointFeature(const QgsPointXY& point) { QgsFeature feature; feature.setGeometry(QgsGeometry::fromPointXY(point)); mLayer->startEditing(); mLayer->addFeature(feature); mLayer->commitChanges(); mCanvas->refresh(); } // 创建线/面要素 void createGeometryFeature() { if (mPoints.size() < 2) return; QgsGeometry geometry; switch (mLayer->geometryType()) { case QgsWkbTypes::LineGeometry: geometry = QgsGeometry::fromPolylineXY(mPoints); break; case QgsWkbTypes::PolygonGeometry: // 多边形首尾闭合 mPoints.append(mPoints.first()); geometry = QgsGeometry::fromPolygonXY({mPoints}); break; default: return; } mLayer->startEditing(); QgsFeature feature; feature.setGeometry(geometry); mLayer->addFeature(feature); mLayer->commitChanges(); mCanvas->refresh(); } // 重置绘制状态 void resetDrawing() { mPoints.clear(); mRubberBand->reset(); } private: QgsMapCanvas* mCanvas; QgsVectorLayer* mLayer; QgsRubberBand* mRubberBand; QVector<QgsPointXY> mPoints; }; ``` --- ### 三、核心编辑功能实现 #### 1. 删除要素 ```cpp void deleteSelectedFeatures(QgsVectorLayer* layer) { if (!layer) return; layer->startEditing(); QgsFeatureIds ids = layer->selectedFeatureIds(); if (!ids.isEmpty()) { layer->deleteFeatures(ids); layer->commitChanges(); } } ``` #### 2. 移动要素 ```cpp class MoveTool : public QgsMapTool { Q_OBJECT public: MoveTool(QgsMapCanvas* canvas, QgsVectorLayer* layer) : QgsMapTool(canvas), mLayer(layer) {} protected: void canvasPressEvent(QgsMapMouseEvent* e) override { // 获取被点击的要素 QgsPointXY point = toMapCoordinates(e->pos()); QgsFeature feature; if (!mLayer->getFeatures(QgsFeatureRequest() .setFilterRect(QgsRectangle(point.x()-0.1, point.y()-0.1, point.x()+0.1, point.y()+0.1))) .nextFeature(feature)) return; mSelectedFeature = feature; mStartPoint = point; } void canvasMoveEvent(QgsMapMouseEvent* e) override { if (mSelectedFeature.isNull()) return; // 计算位移向量 QgsVector v = toMapCoordinates(e->pos()) - mStartPoint; // 应用位移 QgsGeometry geometry = mSelectedFeature.geometry(); geometry.translate(v.x(), v.y()); mLayer->startEditing(); mSelectedFeature.setGeometry(geometry); mLayer->updateFeature(mSelectedFeature); mLayer->commitChanges(); canvas()->refresh(); } private: QgsVectorLayer* mLayer; QgsFeature mSelectedFeature; QgsPointXY mStartPoint; }; ``` #### 3. 复制与剪切 ```cpp // 复制（到内存剪贴板） QList<QgsFeature> copyBuffer; void copySelectedFeatures(QgsVectorLayer* layer) { QgsFeatureIterator it = layer->getSelectedFeatures(); copyBuffer.clear(); QgsFeature feature; while (it.nextFeature(feature)) { // 克隆要素（去除ID） QgsFeature cloned = feature; cloned.setId(FID_NULL); copyBuffer.append(cloned); } } // 粘贴 void pasteFeatures(QgsVectorLayer* layer) { if (copyBuffer.isEmpty()) return; layer->startEditing(); for (QgsFeature& feature : copyBuffer) { layer->addFeature(feature); } layer->commitChanges(); } // 剪切 = 复制 + 删除 void cutFeatures(QgsVectorLayer* layer) { copySelectedFeatures(layer); deleteSelectedFeatures(layer); } ``` --- ### 四、图层类型适配器 ```cpp enum EditAction { Delete, Move, Copy, Cut }; void executeEditAction(QgsVectorLayer* layer, EditAction action) { if (!layer->isValid()) return; switch (layer->geometryType()) { case QgsWkbTypes::PointGeometry: case QgsWkbTypes::LineGeometry: case QgsWkbTypes::PolygonGeometry: // 所有类型支持的操作 switch (action) { case Delete: deleteSelectedFeatures(layer); break; case Copy: copySelectedFeatures(layer); break; case Cut: cutFeatures(layer); break; } break; } } ``` --- ### 五、使用示例 ```cpp // 在主窗口初始化工具 CustomDrawingTool* drawingTool = new CustomDrawingTool(canvas, layer); MoveTool* moveTool = new MoveTool(canvas, layer); // 切换工具 canvas->setMapTool(drawingTool); // 绘制模式 canvas->setMapTool(moveTool); // 移动模式 // 执行命令按钮 QPushButton* btnDelete = new QPushButton("删除"); connect(btnDelete, &QPushButton::clicked, [layer](){ executeEditAction(layer, EditAction::Delete); }); ``` --- ### 关键技术点 1. **手动几何构造**：使用 `QgsGeometry::fromPointXY()/fromPolylineXY()/fromPolygonXY()` 直接创建几何对象 2. **橡皮筋预览**：通过 `QgsRubberBand` 实现绘制过程中的实时可视化[^3] 3. **位移计算**：`QgsVector` 处理要素移动的几何变换 4. **图层适配**：根据 `geometryType()` 自动匹配几何操作 5. **内存剪贴板**：`QList<QgsFeature>` 作为临时缓冲区实现复制/粘贴功能 这种方法避免了直接依赖 `QgsMapToolAddFeature`，提供了更底层的几何操作控制，同时保持了绘制的实时可视化效果[^1]。
评论 6
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
表情包 代码片
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值