Qgis 开发初级 《ToolBox》

Qgis 有个ToolBox 的,在Processing->ToolBox 菜单里面,界面如下。

理论上Qgis这里面的工具都是可以用脚本或者C++ 代码调用的。界面以Vector overlay 为例子简单介绍下使用方式。Vector overlay 的意思是矢量叠置分析,和arcgis软件类似的。点开之后会看到具体的功能。和之前的Geometry的叠置分析的功能是一样的,不同是这个图层级别的叠置分析,效率自然也比遍历Geometry的效率要高。

接下面介绍下使用这些的方法

1、Difference

cpp 复制代码
QgsVectorLayer* Diff(QgsVectorLayer* layer1, QgsVectorLayer* layer2) {
    const QString id = "native:difference";
    const QgsProcessingAlgorithm* algorithm = QgsApplication::processingRegistry()->algorithmById(id);
    QVariantMap conf;
    conf.insert(QStringLiteral("INPUT"), QVariant::fromValue(layer1));
    conf.insert(QStringLiteral("OVERLAY"), QVariant::fromValue(layer2));
    QgsProcessingOutputLayerDefinition value("TEMPORARY_OUTPUT");
    conf.insert(QStringLiteral("OUTPUT"), value);
    QgsProcessingContext context;

    QgsProcessingFeedback feedback;
    QVariantMap results = algorithm->run(conf, context, &feedback);

    QgsVectorLayer* outputLayer(dynamic_cast<QgsVectorLayer*>(context.getMapLayer(results.value(QStringLiteral("OUTPUT")).toString())));
    //暂时先加载,实际用的时候需要屏蔽
    QgsProject::instance()->addMapLayer(outputLayer);

    return outputLayer;
}

2、Clip

cpp 复制代码
QgsVectorLayer* Clip(QgsVectorLayer* inputLayer, QgsVectorLayer* overlayLayer) {
    //差集
    const QString id = "native:clip";
    QVariantMap conf;
    conf.insert(QStringLiteral("INPUT"), QVariant::fromValue(inputLayer));
    conf.insert(QStringLiteral("OVERLAY"), QVariant::fromValue(overlayLayer));
    QgsProcessingOutputLayerDefinition value("TEMPORARY_OUTPUT");
    conf.insert(QStringLiteral("OUTPUT"), value);
    conf.insert(QStringLiteral("LAYER_OPTIONS"), QStringLiteral("COORDINATE_PRECISION=1"));
    QgsProcessingContext context;

    QgsProcessingFeedback feedback;


    const QgsProcessingAlgorithm* algorithm = QgsApplication::processingRegistry()->algorithmById(id);


    QVariantMap results = algorithm->run(conf, context, &feedback);

    QgsVectorLayer* resultlayer(dynamic_cast<QgsVectorLayer*>(context.getMapLayer(results.value(QStringLiteral("OUTPUT")).toString())));

    //暂时先加载,实际用的时候需要屏蔽
    QgsProject::instance()->addMapLayer(resultlayer);

    return resultlayer;
}

3、Union

cpp 复制代码
QgsVectorLayer* Union(QgsVectorLayer* layer1, QgsVectorLayer* layer2) {
    const QString id = "native:union";
    const QgsProcessingAlgorithm* algorithm = QgsApplication::processingRegistry()->algorithmById(id);

    // 检查图层类型
    if (layer1->geometryType() != layer2->geometryType()) {
        qWarning() << "Input layers must have the same geometry type.";
        return nullptr;
    }

    QVariantMap conf;
    conf.insert(QStringLiteral("INPUT"), QVariant::fromValue(layer1));
    conf.insert(QStringLiteral("OVERLAY"), QVariant::fromValue(layer2));
    QgsProcessingOutputLayerDefinition value("TEMPORARY_OUTPUT");
    conf.insert(QStringLiteral("OUTPUT"), value);
    conf.insert(QStringLiteral("LAYER_OPTIONS"), QStringLiteral("COORDINATE_PRECISION=1"));
    QgsProcessingContext context;

    QgsProcessingFeedback feedback;
    QVariantMap results = algorithm->run(conf, context, &feedback);

    QgsVectorLayer* outputLayer(dynamic_cast<QgsVectorLayer*>(context.getMapLayer(results.value(QStringLiteral("OUTPUT")).toString())));
    //暂时先加载,实际用的时候需要屏蔽
    QgsProject::instance()->addMapLayer(outputLayer);

    return outputLayer;
}

4、小节

使用这些工具的代码基本大同小异,最关键的就是id,id基本就是原来的名字的小写。读者可以自己去试试里面其他工具的使用。使用过程需要注意一点的是如果有无效的Geometey,可能会报错,可以在执行之前删除无效的Geometry。还有一点就是面图层在执行过程中因为误差问题可能会出现线的情况,常常需要删除线,一般需要在Geometry的 parts 里面做判断的。还有一些qgis的常用功能,像投影,三维,server和pyphon脚本的使用。这些功能后面可能会继续介绍,也可能不在做介绍了。

相关推荐
罗西的思考3 分钟前
[2W字长文] 探秘Transformer系列之(23)--- 长度外推
人工智能·算法
codelang1 小时前
Cline + MCP 开发实战
前端·后端
风象南2 小时前
SpringBoot中6种自定义starter开发方法
java·spring boot·后端
好_快2 小时前
Lodash源码阅读-memoizeCapped
前端·javascript·源码阅读
好_快2 小时前
Lodash源码阅读-toString
前端·javascript·源码阅读
好_快2 小时前
Lodash源码阅读-memoize
前端·javascript·源码阅读
excel3 小时前
webpack 核心编译器 十四 节
前端
excel3 小时前
webpack 核心编译器 十三 节
前端
腾讯TNTWeb前端团队10 小时前
helux v5 发布了,像pinia一样优雅地管理你的react状态吧
前端·javascript·react.js