7. 矢量图层数据查询选择和保存

于 2023-10-23 11:26:29 发布
阅读量305 收藏 2
点赞数 2
分类专栏： QGis二次开发 文章标签： qgis二次开发 矢量图层数据
本文为博主原创文章，未经博主允许不得转载。
本文链接：https://blog.csdn.net/richie12/article/details/133984015
版权

文章目录

前言

为了获取地理空间数据的详细信息及统计结果，需要使用数据选择、筛选、查询与统计等方法。
数据查询可以获得矢量数据与栅格数据的属性信息；
数据选择（仅针对矢量图层）是指通过空间位置、属性信息等特征选取的部分地理要素，被选择的地理要素可以进一步被其他的编辑工具或分析工具处理；
数据筛选是指通过表达式等方式仅加载矢量数据中的部分要素，以便进一步制图或处理；
数据统计可以针对这些属性信息进行统计操作。

数据查询

QGis

  • 以myplaces.shp为例，添加图层后，右键菜单选择Open Attribute Table, 如下图
    在这里插入图片描述
  • 可以看到一个二维表格，显示矢量图层的数据，如下图
    在这里插入图片描述
  • 其数据存储和关系型数据库相同
  • 表头NAME、PLACE、ADDR_CITY等等相当于数据库中的字段，在QGis中以QgsField类表示
  • 每一行是一个Feature，feature除表中的数据外还有空间属性，在QGis中以QgsFeature类表示
  • 每一行的每一列是一个具体的数据，在QGis中以QgsAttributes类表示
  • 可以从表格中筛选数据，点击下图所示按钮
    在这里插入图片描述
  • 输入表达式，如下图
    在这里插入图片描述
  • 选择Filter Current Selection后，点击按钮后，效果如下图所示
    在这里插入图片描述
    在这里插入图片描述

代码实现

获取图层数据的字段QgsField

  • 类QgsField代表字段，详情见文档
  • 类QgsVectorLayer中的函数QgsFields fields()用于获取字段
void MainWindow::getFieldsSlot()
{
    //添加测试图层
    QString filename = QStringLiteral("maps/shapefile/myplaces.shp");
    QFileInfo ff(filename);
    QgsVectorLayer* layer = (QgsVectorLayer*)mApp->addVectorLayer(filename,ff.baseName());
    zoomToFirstLayer<QgsVectorLayer*>();

    mVectorDataDockWidget->show();
    mVectorDataDockWidget->plainTextEdit()->clear();
    for(auto field : layer->fields())
    {
        mVectorDataDockWidget->plainTextEdit()->appendPlainText(QString("%1:%2").arg(field.name()).arg(field.typeName()));
    }
}

在这里插入图片描述
在这里插入图片描述

获取图层数据QgsFeature以及QgsAttribute

  • 类QgsFeature代表图层中一个特征，该特征包含一个featureId，空间特征(点、线、面以及其坐标)和字段/值属性列表。详情见文档
  • 类QgsAttributes是feature的字段/值属性，详情见文档
  • 类QgsVectorLayer的函数QgsFeatureIterator getFeatures( const QgsFeatureRequest &request = QgsFeatureRequest() )用户获取features，其参数QgsFeatureRequest用于过滤数据，默认参数会获取所有feature
  • 类QgsFeature的函数attributes()用于获取feature的属性
  • 首先演示获取所有feature，代码如下
void MainWindow::getFeaturesSlot()
{
    //添加测试图层
    QString filename = QStringLiteral("maps/shapefile/myplaces.shp");
    QFileInfo ff(filename);
    QgsVectorLayer* layer = (QgsVectorLayer*)mApp->addVectorLayer(filename,ff.baseName());
    zoomToFirstLayer<QgsVectorLayer*>();

    mVectorDataDockWidget->show();
    mVectorDataDockWidget->plainTextEdit()->clear();

    QgsFeatureIterator it = layer->getFeatures();
    QgsFeature f;
    while(it.nextFeature(f))
    {
        QString str;
        str.append(QString("%1").arg(f.id()));
        str.append(" ");
        QgsAttributes attrs = f.attributes();
        for(int i = 0;i < attrs.size();++i)
        {
            str.append(attrs[i].toString());
            str.append(" ");
        }
        mVectorDataDockWidget->plainTextEdit()->appendPlainText(str);
    }
}

在这里插入图片描述
在这里插入图片描述

  • 筛选数据的代码如下，将字段RAINFALL的值大于200.0的保留
    QString str = QString("\"RAINFALL\" > 200.0");
    auto request = QgsFeatureRequest().setFilterExpression(str);
    QgsFeatureIterator it = layer->getFeatures(request);
  • 完整测试代码如下
void MainWindow::getFeaturesSlot()
{
    //添加测试图层
    QString filename = QStringLiteral("maps/shapefile/myplaces.shp");
    QFileInfo ff(filename);
    QgsVectorLayer* layer = (QgsVectorLayer*)mApp->addVectorLayer(filename,ff.baseName());
    zoomToFirstLayer<QgsVectorLayer*>();

    mVectorDataDockWidget->show();
    mVectorDataDockWidget->plainTextEdit()->clear();

    //带feature的expression
    QString str = QString("\"RAINFALL\" > 200.0");
    auto request = QgsFeatureRequest().setFilterExpression(str);
    QgsFeatureIterator it = layer->getFeatures(request);
    QgsFeature f;
    while(it.nextFeature(f))
    {
        QString str;
        str.append(QString("%1").arg(f.id()));
        str.append(" ");
        QgsAttributes attrs = f.attributes();
        for(int i = 0;i < attrs.size();++i)
        {
            str.append(attrs[i].toString());
            str.append(" ");
        }
        mVectorDataDockWidget->plainTextEdit()->appendPlainText(str);
    }
}

在这里插入图片描述

数据选择

QGis

  • 在“Attributes Toolbar”工具栏中，单击按钮右侧的下拉按钮，弹出的菜单中展示了交互式选择的四种方式，如下图所示
    在这里插入图片描述
  • 选择要素之后，颜色变为黄色，如下图
    在这里插入图片描述
  • 在默认情况下，选择的要素以黄色的点、线或黄色面填充符号突出显示在地图画布上。在QGIS的“Options”对话框的“Canvas & Legend”选项卡中，通过“Selection color”选项即可更改选择要素的颜色
    在这里插入图片描述

代码实现

  • 类QgsVectorLayer的一组选择函数实现了多种选择方式，代码如下
void select( QgsFeatureId featureId );
void select( const QgsFeatureIds &featureIds );
void selectAll();
void selectByRect( QgsRectangle &rect, Qgis::SelectBehavior behavior = Qgis::SelectBehavior::SetSelection );
void selectByExpression( const QString &expression, Qgis::SelectBehavior behavior = Qgis::SelectBehavior::SetSelection, QgsExpressionContext *context = nullptr );
void selectByIds( const QgsFeatureIds &ids, Qgis::SelectBehavior behavior = Qgis::SelectBehavior::SetSelection );
void deselect( QgsFeatureId featureId );
void deselect( const QgsFeatureIds &featureIds );
void removeSelection();
void reselect();
void modifySelection( const QgsFeatureIds &selectIds, const QgsFeatureIds &deselectIds );
  • 完整测试代码如下
void MainWindow::selectFeaturesSlot()
{
    //添加测试图层
    QString filename = QStringLiteral("maps/shapefile/myplaces.shp");
    QFileInfo ff(filename);
    QgsVectorLayer* layer = (QgsVectorLayer*)mApp->addVectorLayer(filename,ff.baseName());
    zoomToFirstLayer<QgsVectorLayer*>();

    QString expression = QString("\"RAINFALL\" > 200.0");
    //选择所有
    //    layer->selectAll();
    //按照id选择
    //    QgsFeatureId id = 0;
    //    layer->select(id);
    //按照多个id选择
    //    QgsFeatureIds ids;
    //    ids << 0 << 1 << 2 << 3;
    //    layer->select(ids);
    //按照表达式选择
    //    layer->selectByExpression(expression);
    //按照区域选择,并且使用rubberband显示区域
    QgsRectangle rect(1006585,6222254,1010253,6219118);
    layer->selectByRect(rect);
    QgsPointXY point1(1006585,6222254);
    QgsPointXY point2(1010253,6222254);
    QgsPointXY point3(1006585,6219118);
    QgsPointXY point4(1010253,6219118);
    QgsRubberBand *rubberBand = new QgsRubberBand(mApp->mapCanvas(),QgsWkbTypes::PolygonGeometry);
    rubberBand->addPoint(point1);
    rubberBand->addPoint(point2);
    rubberBand->addPoint(point4);
    rubberBand->addPoint(point3);
    rubberBand->show();
    //选择feature并且闪烁
    //    QgsMapCanvasUtils::flashMatchingFeatures(mApp->mapCanvas(),layer,expression);
    //选择Feature并zoom
    //    QgsMapCanvasUtils::zoomToMatchingFeatures(mApp->mapCanvas(),layer,expression);
}

在这里插入图片描述
在这里插入图片描述

创建虚拟图层并保存为shape文件

虚拟图层是在不改变矢量数据源的情况下，改变某些空间或属性信息，或者组合多个矢量图层，并以虚拟图层的方式展现数据。虚拟图层并不直接指向数据源，而是通过SQL查询语句等引用一个或多个矢量图层。

QGis

  • 此处介绍如何通过虚拟图层将“2017.xlsx”文件连接到吉林省地级行政区划文件（jilin_dist.shp）
  • 分别加载这两个图层
    在这里插入图片描述
  • 打开图层jilin_dist的Attribute Table，如下图
    在这里插入图片描述
  • 打开图层2017–Sheet1的Attribute Table，如下图
    在这里插入图片描述
  • 添加一个虚拟图层，将上述两个图层的数据连接到一起
    在这里插入图片描述
  • 在“Layer name”文本框中输入新生成的虚拟图层的名称；在“Embedded layers”组合框中单击“Import”按钮，将上述两个图层加入列表中，并在“Local name”选项中设置图层名称，同时将其作为SQL查询语句的表名称。单击“Add”按钮可以输入图层位置；单击“Remove”按钮可以删除图层。在“Query”文本框中输入以下查询语句：
SELECT * FROM jilin_dist left outer join prec on jilin_dist.NAME = prec.city;
  • 该语句可以通过jilin_dist图层的“NAME”字段和prec图层的“city”字段将prec图层的属性连接到jilin_dist图层中。
  • 单击“Add”按钮即可在图层列表中看到该虚拟图层。打开virtual_layer图层的Attribute Table可以看到连接后的数据表，如下图
    在这里插入图片描述

代码实现

  • 类QgsVirtualLayerDefinition定义虚拟图层连接信息，详情见文档
  • 其成员函数void addSource( const QString &name, const QString &source, const QString &provider, const QString &encoding = "" );用于添加添加图层。代码如下
    //QgsVirtualLayerDefinition
    QgsVirtualLayerDefinition def;
    //add embedded layers
    def.addSource("prec","maps/virtuallayer/2017.xlsx|layername=Sheet1","ogr");
    def.addSource("jilin_dist","maps/virtuallayer/jilin_dist.shp","ogr");
  • 成员函数void setQuery( const QString &query )用于设置查询语句
    //查询语句
    def.setQuery("SELECT * FROM jilin_dist left outer join prec on jilin_dist.NAME = prec.city;");
  • 新建矢量图层QgsVectorLayer，数据源provider传入virtual，代码如下
const QgsVectorLayer::LayerOptions options { QgsProject::instance()->transformContext() };
QgsVectorLayer *vl = new QgsVectorLayer( def.toString(), QStringLiteral( "test" ), QStringLiteral( "virtual" ), options );
  • 类QgsVectorFileWriter用于将图层保存至文件，详情见文档
  • 其成员函数writeAsVectorFormatV3用于保存，虚拟图层保存至文件的代码如下
        QgsVectorFileWriter::SaveVectorOptions opt;
        QgsVectorFileWriter::writeAsVectorFormatV3(vl,"maps/virtuallayer/myvurtual.shp",QgsProject::instance()->transformContext(),opt);

在这里插入图片描述
在这里插入图片描述

  • 保存的图层文件如下
    在这里插入图片描述

总结

  • 分别演示了：
    • QGis软件矢量图层数据的查询、选择和保存
    • 二次开发c++ api的图层数据的查询、选择和保存
雷动软件工作室
关注
  • 2
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论
专栏目录
python3 qgis开发_QGIS 笔记三：利用pyQGIS进行shapefile坐标转换
weixin_39662594的博客
12-22 636
QGIS 笔记三：利用pyQGIS进行shapefile坐标转换前言不断有人问我关于Arcpy的问题，但是我却真的一点都不了解Arcpy。最近有个task是将一个shapefile中的所有点的坐标从wgs84转换成为gcj。我决心试用一下Arcpy，结果发现非常难使用。接口非常不人性化，而且，文档也不够详细。使用基本上是靠模仿内置的python文件去写自己的脚本，路径是path\to\ArcGIS...
学习QGIS Python API（应用程序编程接口）
Tiny的博客
05-17 1万+
目录 前言 说明 API (Application Programming Interface) 摘要 关于QGIS Python API 解密C ++文档 组织QGIS Python库 qgis.core包 地图和地图图层 协调参考系统 矢量图层 显示矢量数据 访问矢量数据 空间索引 栅格图层 显示栅格数据 访问栅格数据 其他有用的qgis.core类 ...
Qgis 绘制要素并保存为shp文件
josiecici的博客
10-26 7878
Qgis 绘制要素
QGIS查看属性和选择要素
weixin_46713695的博客
01-05 6609
QGIS查看属性和选择要素
QGIS保存图层样式
zhanghongxia1111的博客
12-28 955
geoserver发布mongodb矢量数据地图服务.docx
02-15
将mongodb矢量数据在geoserver中发布；将mongodb矢量数据在geoserver中发布
在ArcGIS中栅格数据矢量化.doc
10-06
在 ArcGIS 中栅格数据矢量化 在 ArcGIS 中栅格数据矢量化是指将栅格图像转换为矢量数据的过程。在这个过程中，我们可以使用 ArcGIS ...通过这种方法，我们可以将栅格图像转换为矢量数据，用于后续的空间分析和处理。
矢量数据配准具体步骤（ARCGIS软件）.doc
09-30
在开始矢量数据配准之前，需要准备两个要素类：已经有坐标系的要素类（基准要素）和需要校正的要素类（被校正要素）。这两个要素类需要加进 ArcMap 中，并使需要校正的图层处于编辑状态。 二、 spatial adjustment ...
最短路径和导航-源码和文档.zip_C#_arcengine _arcengine 位置_图层软件源码_导航
07-15
能够查看矢量图层的属性表，并进行条件选择，在地图上高亮显示。能够显示鼠标当前所指的要素属性信息，鼠标滑过教学楼可以显示其名字和类型（通过配置实现，本软件只配置了建筑物图层）。能够进行最短路径分析，根据...
ArcGIS Engine投影变换/坐标转换/矢量地图配准/七参数计算和应用功能源码
05-04
自己开发的基于DevExpress 16.1.1的ribbon主流界面（类似于平时所用的word扁平化界面），主要实现基本的地图功能和投影变换/坐标转换/矢量地图配准/七参数计算和应用等功能。对于2000坐标的转换，软件自带了国家2000...
QGIS安装以及使用教程
人无远虑，必有近忧
07-08 2万+
Quantum GIS（QGIS）是相当友好的开源地理信息系统，使用GNU（General Public License）授权，属于 Open Source eospatial Foundation(OSGeo)的官方计划。在 GNU 这个授权下，开发者可以自行检阅与调整程序代码，并保障让所有使用者可以免费且自由地修改程序。 QGIS 的目标是成为一个使用简单的 GIS，提供了常见的功能与图征。借着核心图征提供基本功能与附加组件（Python 或 C++）持续地扩充，让使用者可以浏览、管理、编辑、分析数据
[QGIS]常用操作--矢量图层加载与创建
Olio
07-19 9300
文章中使用的QGIS版本为:3.14.16-Pi 本文主要介绍如何使用QGIS通过使用界面操作和Python代码打开矢量数据(如shp文件或csv文件),以及如何添加矢量数据中要素. 打开矢量图层与获取要素矢量图层加载加载CSV文件加载SHP文件重投影与格式互转新建图层 矢量图层加载 加载CSV文件 我在平时调试程序时,会使用csv文件将临时结果输出,然后加载到QGIS中查看.这里简单介绍下,怎么加载csv文件. 菜单栏点击"图层"->“添加图层”->“添加文本数据图层” 在弹出的"数据源.
qgis 软件下载和图层管理和操作
xiaojiejie123s的博客
08-26 752
Quantum GIS (QGIS)是一个用户界面友好、开源代码的GIS客户端程序，支持数据的可视化、管理、编辑与分析和印刷地图的制作。通过集成GRASS，支持强大的分析功能，可运行于Linux、Unix、Mac OSX和Windows系统，并支持多种矢量、栅格与数据库格式及功能。 1.友好的图形用户界面 地物的选择和指定 属性的编辑、显示和检索 动态投影 打印定制 地物标记 矢量和栅格符
从零开始gis——(3)使用qgis绘制图层,并将shp导入到postgres数据库
热门推荐
lqx_sunhan的博客
08-11 2万+
开发环境搭建好了,我们现在需要的是图层，我们可以自己绘制一个。 QGIS版本: 3.2.1 PostgreSQL版本:9.5.0 PostGIS版本: 2.3.7 1.打开qgis软件，点击上方菜单的project –&gt; new 2. 上方菜单中的 layer –&gt; Create layer –&gt; New Shapefile Layer 3. 点击上方...
QGIS中怎样设置图层样式并导出地图样式
BADAO_LIUMANG_QIZHI的博客
04-12 3366
场景 QGIS怎样设置简体中文以及新建可编辑的多边形的图层： QGIS怎样设置简体中文以及新建可编辑的多边形的图层_霸道流氓气质的博客-CSDN博客 在上面新增图层之后怎样设置图层的样式并且将样式导出为文件。 注： 博客：霸道流氓气质的博客_CSDN博客-C#,架构之路,SpringBoot领域博主 关注公众号 霸道的程序猿 获取编程相关电子书、教程推送与免费下载。 实现 1、选中图层右键选择属性，可以在符号化中设置颜色、不透明度、边框等。 2、然后点击左下角样式选择保存样式
QGIS专题制图及输出
小小的天空
04-10 2万+
最近因项目需要，研究了下开源GIS软件QGIS，学GIS做GIS多年，ArcGIS用的很熟悉了，第一次用QGIS，界面框架跟ArcGIS类似，原理相通，上手摸索了一周，发现在专题制图和输出方面跟ArcGIS还是不太一样。在专题制图上，QGIS功能只支持几种主要的方式： single symbol :单一符号专题题； categorized: 分类符号专题图； graduated: 分级符号
QGis（三）查询矢量图层的要素属性字段值
tszhangjunqiao的专栏
08-13 1万+
当加载一个矢量图层后，如果要查看每个要素的属性
chromedriver-win64_116.0.5840.0.zip
最新发布
06-15
chromedriver-win64_116.0.5840.0.zip
基于Java Servlet实现的灾情控制系统.zip
06-15
该资源内项目源码是个人的课程设计、毕业设计，代码都测试ok，都是运行成功后才上传资源，答辩评审平均分达到96分，放心下载使用！ ## 项目备注 1、该资源内项目代码都经过测试运行成功，功能ok的情况下才上传的，请放心下载使用！ 2、本项目适合计算机相关专业(如计科、人工智能、通信工程、自动化、电子信息等)的在校学生、老师或者企业员工下载学习，也适合小白学习进阶，当然也可作为毕设项目、课程设计、作业、项目初期立项演示等。 3、如果基础还行，也可在此代码基础上进行修改，以实现其他功能，也可用于毕设、课设、作业等。 下载后请首先打开README.md文件（如有），仅供学习参考, 切勿用于商业用途。 该资源内项目源码是个人的课程设计，代码都测试ok，都是运行成功后才上传资源，答辩评审平均分达到96分，放心下载使用！ ## 项目备注 1、该资源内项目代码都经过测试运行成功，功能ok的情况下才上传的，请放心下载使用！ 2、本项目适合计算机相关专业(如计科、人工智能、通信工程、自动化、电子信息等)的在校学生、老师或者企业员工下载学习，也适合小白学习进阶，当然也可作为毕设项目、课程设计、作业、项目初期立项演示等。 3、如果基础还行，也可在此代码基础上进行修改，以实现其他功能，也可用于毕设、课设、作业等。 下载后请首先打开README.md文件（如有），仅供学习参考, 切勿用于商业用途。
最新2022年山西省道路矢量图层shp格式数据下载
11-21
4. 在数据下载页面中，查找和选择与您需求相关的2022年山西省道路矢量图层shp数据。可能会提供多个下载选项，您可以根据自己的需要选择最合适的数据版本和分辨率。 5. 确认选择后，点击下载按钮或链接，开始下载...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

雷动软件工作室

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值