GDAL C++ API 学习之路 (5) OGRGeometry 几何类 OGRGeometry

最新推荐文章于 2024-05-25 01:54:49 发布
最新推荐文章于 2024-05-25 01:54:49 发布
阅读量892 收藏 2
点赞数
分类专栏： OGR 文章标签： c++ 学习 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_69574549/article/details/131923388
版权

Distance

virtual double Distance(const OGRGeometry*) const

计算两个几何之间的距离        返回两个几何图形之间的最短距离。距离表示为与几何坐标相同的单位

参数:

poOtherGeom -- 另一个要比较的几何图形。

返回: 几何图形之间的距离，如果发生错误，则为 -1

// 创建两个点
OGRPoint point1(0, 0);
OGRPoint point2(3, 4);

// 计算两个点之间的距离
double distance = point1.Distance(&point2);

cout << "Distance between point1 and point2: " << distance << endl;

ConvexHull

Virtual OGRGeometry *ConvexHull（） const

计算凸包        将创建并返回一个新的几何对象，其中包含调用该方法的几何的凸包

返回: 现在由调用方拥有的新分配的几何图形，或失败时为 NULL

    // 创建一个线环（LinearRing）
    OGRLinearRing ring;
    ring.addPoint(0, 0);
    ring.addPoint(1, 0);
    ring.addPoint(1, 1);
    ring.addPoint(0, 1);
    ring.addPoint(0, 0);

    // 创建一个多边形（Polygon）
    OGRPolygon polygon;
    polygon.addRing(&ring);

    // 计算凸包
    OGRGeometry* convexHull = polygon.ConvexHull();

    // 输出凸包的WKT表示
    char* wkt;
    convexHull->exportToWkt(&wkt);
    std::cout << "Convex Hull: " << wkt << std::endl;

    // 释放内存
    CPLFree(wkt);
    delete convexHull;

ConcaveHull

virtual OGRGeometry *ConcaveHull(double dfRatio, bool bAllowHoles) const

计算几何图形的“凹壳”        计算几何图形的“凹壳”

凹壳完全包含在凸包内，并且还包含输入的所有点，但面积较小。面积比是凸包和凹壳面积的比值。常用于将多点转换为多边形区域。包含输入几何图形中的所有点

参数:

  • dfRatio -- 凸包和凹包面积的比率。

  • bAllowHole -- 是否允许孔。

返回: 现在由调用方拥有的新分配的几何图形，或失败时为 NULL

    // 创建一个线环（LinearRing）
    OGRLinearRing ring;
    ring.addPoint(0, 0);
    ring.addPoint(1, 0);
    ring.addPoint(1, 1);
    ring.addPoint(0, 1);
    ring.addPoint(0, 0);

    // 创建一个多边形（Polygon）
    OGRPolygon polygon;
    polygon.addRing(&ring);

    // 计算凹多边形凸包
    double concaveHullRatio = 0.8; // 指定一个凹多边形凸包的压缩比例
    bool allowHoles = true;        // 允许生成含有洞的凹多边形凸包
    OGRGeometry* concaveHull = polygon.ConcaveHull(concaveHullRatio, allowHoles);

    // 输出凹多边形凸包的WKT表示
    char* wkt;
    concaveHull->exportToWkt(&wkt);
    std::cout << "Concave Hull: " << wkt << std::endl;

    // 释放内存
    CPLFree(wkt);
    delete concaveHull;

Buffer

virtual OGRGeometry *Buffer（double dfDist， int nQuadSegs = 30） const

计算几何的缓冲区

生成一个新几何图形，其中包含调用该几何图形的缓冲区。缓冲区是一个面，其中包含原始几何的缓冲距离内的区域

某些缓冲区部分被正确描述为曲线，但会转换为近似多边形。nQuadSegs 参数可用于控制应使用多少条线来定义 90 度曲线 - 圆的象限。值 30 是合理的默认值。较大的值会导致生成的缓冲区几何中出现大量顶点，而较小的数字会降低结果的准确性

参数:

  • dfDist -- 要应用的缓冲距离。应表示为与几何坐标相同的单位。

  • nQuadSegs -- 用于近似 90 度（象限）曲率的段数。

返回: 新创建的几何图形，如果发生错误，则为 NULL

 // 创建一个点（Point）
 OGRPoint point(0, 0);

// 对点进行缓冲区分析，生成一个缓冲区几何对象
double bufferDistance = 1.0; // 缓冲区的距离
int numQuadSegments = 30;    // 四边形段数，用于控制缓冲区的平滑程度
OGRGeometry* buffer = point.Buffer(bufferDistance, numQuadSegments);

Intersection

Virtual OGRGeometry *Intersection（const OGRGeometry*） const

计算交集

生成一个新的几何图形，它是操作的两个几何图形的交点区域。Intersects（） 方法可用于测试两个几何图形是否相交

不检查几何有效性。如果不确定输入几何的有效性，在之前调用 IsValid（），否则结果可能是错误的

参数:

poOtherGeom -- 另一个几何与“这个”几何相交。

返回: 表示交点的新几何图形，如果没有交集或发生错误，则为 NULL

    // 创建两个多边形对象
    OGRLinearRing ring1, ring2;
    ring1.addPoint(0, 0);
    ring1.addPoint(0, 1);
    ring1.addPoint(1, 1);
    ring1.addPoint(1, 0);
    ring1.addPoint(0, 0);

    ring2.addPoint(0.5, 0.5);
    ring2.addPoint(0.5, 1.5);
    ring2.addPoint(1.5, 1.5);
    ring2.addPoint(1.5, 0.5);
    ring2.addPoint(0.5, 0.5);

    OGRPolygon polygon1, polygon2;
    polygon1.addRing(&ring1);
    polygon2.addRing(&ring2);

    // 计算两个多边形的交集
    OGRGeometry* intersection = polygon1.Intersection(&polygon2);

Union

Virtual OGRGeometry *Union（const OGRGeometry*） const 

计算联合

生成一个新的几何图形，该几何图形是操作的两个几何图形的并集区域

不检查几何有效性。如果不确定输入几何的有效性，在之前调用 IsValid（），否则结果可能是错误的

参数:

poOtherGeom -- 与“这个”几何结合的其他几何。

返回:  表示并集或 NULL 的新几何图形（如果发生错误）

// 创建两个多边形对象
    OGRLinearRing ring1, ring2;
    ring1.addPoint(0, 0);
    ring1.addPoint(0, 1);
    ring1.addPoint(1, 1);
    ring1.addPoint(1, 0);
    ring1.addPoint(0, 0);

    ring2.addPoint(0.5, 0.5);
    ring2.addPoint(0.5, 1.5);
    ring2.addPoint(1.5, 1.5);
    ring2.addPoint(1.5, 0.5);
    ring2.addPoint(0.5, 0.5);

    OGRPolygon polygon1, polygon2;
    polygon1.addRing(&ring1);
    polygon2.addRing(&ring2);

    // 计算两个多边形的并集
    OGRGeometry* unionGeometry = polygon1.Union(&polygon2);

UnionCascaded

Virtual OGRGeometry *UnionCascaded（） const

使用级联计算联合        输入几何必须是多面

不检查几何有效性。如果不确定输入几何的有效性，在之前调用 IsValid（），否则结果可能是错误的

返回: 表示并集或 NULL 的新几何图形（如果发生错误）

    // 创建一个多边形对象
    OGRLinearRing ring;
    ring.addPoint(0, 0);
    ring.addPoint(0, 1);
    ring.addPoint(1, 1);
    ring.addPoint(1, 0);
    ring.addPoint(0, 0);

    OGRPolygon polygon;
    polygon.addRing(&ring);

    // 计算多边形与其自身的并集
    OGRGeometry* unionGeometry = polygon.UnionCascaded();

  • Virtual OGRGeometry *UnionCascaded() const: 这个函数用于返回当前几何对象与其自身的联合（Union）。也就是说，它将当前几何对象中的所有几何组件进行联合操作，得到一个包含所有组件的新几何对象。例如，如果当前几何对象是一个多边形，且其中包含了多个不相交的环，那么该函数将返回一个新的多边形，其中包含了所有环的联合结果。

  • Virtual OGRGeometry *Union(const OGRGeometry*) const: 这个函数用于返回当前几何对象与另一个几何对象的联合（Union）。也就是说，它将当前几何对象与另一个几何对象进行联合操作，得到一个包含两个几何对象的所有组件的新几何对象。例如，如果当前几何对象是一个多边形，而传入的参数是另一个多边形，那么该函数将返回一个新的多边形，其中包含了两个多边形的联合结果

  • 总结：

  • UnionCascaded() 函数用于计算当前几何对象内部的联合结果。
  • Union(const OGRGeometry*) 函数用于计算当前几何对象与另一个几何对象的联合结果

UnaryUnion

OGRGeometry *UnaryUnion（） const

返回单个几何图形的所有组件的并集        通常用于将集合转换为覆盖相同区域的最小多边形集

返回: 表示并集或 NULL 的新几何图形（如果发生错误）

    // 创建两个多边形对象
    OGRLinearRing ring1;
    ring1.addPoint(0, 0);
    ring1.addPoint(0, 1);
    ring1.addPoint(1, 1);
    ring1.addPoint(1, 0);
    ring1.addPoint(0, 0);

    OGRPolygon polygon1;
    polygon1.addRing(&ring1);

    OGRLinearRing ring2;
    ring2.addPoint(1, 1);
    ring2.addPoint(1, 2);
    ring2.addPoint(2, 2);
    ring2.addPoint(2, 1);
    ring2.addPoint(1, 1);

    OGRPolygon polygon2;
    polygon2.addRing(&ring2);

    // 计算两个多边形对象的一元并集
    OGRGeometry* unaryUnionGeometry = polygon1.Union(&polygon2);

Difference

virtual OGRGeometry *Difference(const OGRGeometry*) const

计算差异

生成一个新几何图形，该几何图形是该几何图形的区域，并删除了第二个几何图形的区域

不检查几何有效性。如果不确定输入几何的有效性，在之前调用 IsValid（），否则结果可能是错误的

参数:

poOtherGeom -- 从“this”几何中删除的其他几何图形。

返回:  表示差值的新几何图形，如果差值为空或发生错误，则为 NULL

    // 创建两个多边形对象
    OGRLinearRing ring1;
    ring1.addPoint(0, 0);
    ring1.addPoint(0, 2);
    ring1.addPoint(2, 2);
    ring1.addPoint(2, 0);
    ring1.addPoint(0, 0);

    OGRPolygon polygon1;
    polygon1.addRing(&ring1);

    OGRLinearRing ring2;
    ring2.addPoint(1, 1);
    ring2.addPoint(1, 3);
    ring2.addPoint(3, 3);
    ring2.addPoint(3, 1);
    ring2.addPoint(1, 1);

    OGRPolygon polygon2;
    polygon2.addRing(&ring2);

    // 计算两个多边形对象的差异
    OGRGeometry* differenceGeometry = polygon1.Difference(&polygon2);

SymDifference

virtual OGRGeometry *SymDifference（const OGRGeometry*） const

计算对称差分

生成一个新的几何图形，该几何图形是此几何图形与传递到方法中的第二个几何图形的对称差值

不检查几何有效性。如果不确定输入几何的有效性，在之前调用 IsValid（），否则结果可能是错误的

参数:

poOtherGeom -- 其他几何形状。

返回:  表示对称差值的新几何图形，如果差值为空或发生错误，则为 NULL

    // 创建两个多边形对象
    OGRLinearRing ring1;
    ring1.addPoint(0, 0);
    ring1.addPoint(0, 2);
    ring1.addPoint(2, 2);
    ring1.addPoint(2, 0);
    ring1.addPoint(0, 0);

    OGRPolygon polygon1;
    polygon1.addRing(&ring1);

    OGRLinearRing ring2;
    ring2.addPoint(1, 1);
    ring2.addPoint(1, 3);
    ring2.addPoint(3, 3);
    ring2.addPoint(3, 1);
    ring2.addPoint(1, 1);

    OGRPolygon polygon2;
    polygon2.addRing(&ring2);

    // 计算两个多边形对象的对称差异
    OGRGeometry* symDifferenceGeometry = polygon1.SymDifference(&polygon2);

Centroid

virtual OGRErr Centroid（OGRPoint *poPoint） const

计算几何质心

质心位置应用于传入的 OGRPoint 对象。质心不一定在几何图形内

返回:  OGRERR_NONE成功或OGRERR_FAILURE错误

    // 创建一个多边形对象
    OGRLinearRing ring;
    ring.addPoint(0, 0);
    ring.addPoint(0, 4);
    ring.addPoint(4, 4);
    ring.addPoint(4, 0);
    ring.addPoint(0, 0);

    OGRPolygon polygon;
    polygon.addRing(&ring);

    // 创建一个空的 OGRPoint 对象用于存储质心
    OGRPoint centroidPoint;

    // 计算多边形对象的质心
    OGRErr err = polygon.Centroid(&centroidPoint);

    if (err == OGRERR_NONE) {
        // 输出质心的坐标
        double x = centroidPoint.getX();
        double y = centroidPoint.getY();
        std::cout << "Centroid coordinates: (" << x << ", " << y << ")" << std::endl;
    } else {
        std::cout << "Error calculating centroid." << std::endl;

场主不吃鍋巴
关注
  • 0
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
开源ＧＩＳ简介
shaojieli的专栏
02-04 1870
C++开源GIS中间件类库： 　　GDAL（栅格）/OGR（矢量）提供了类型丰富的读写支持 　　GEOS(Geometry Engine Open Source)是基于Ｃ＋＋的空间拓扑分析实现类库，遵循ＬＧＰＬ协议发布。GEOS类库提供了丰富的空间拓扑操作函数，用以判断几何对象间的相互关系，以及空间分析操作之后形成新的几何对象。点、线、面要素的两两相互关系，包括相合、分离、相交、重合、包含、相
GDAL C++ API 学习之路 (1) OGRGeometry 几何类 OGRGeometry
qq_69574549的博客
07-09 597
OGRGeometry 类是 GDAL/OGR 库中用于处理地理空间数据的核心类之一。它是一个抽象基类，表示各种几何对象，如点、线、面等。OGRGeometry 类提供了许多方法和操作，用于创建、编辑、分析和操作几何对象
【小沐学GIS】GDAL库安装和使用（C++、Python）
最新发布
爱看书的小沐
05-25 2646
GDAL 是用于栅格和矢量地理空间数据格式的转换器库，由开源地理空间基金会在 MIT 风格的开源许可下发布。作为一个库，它向调用应用程序提供单个栅格抽象数据模型和单个矢量抽象数据模型，用于所有支持的格式。它还带有各种有用的命令行实用程序，用于数据转换和处理。╮(￣▽￣)╭如果您感觉方法或代码不咋地//(ㄒoㄒ)//，就在评论处留言，作者继续改进；o_O???如果您需要相关功能的代码定制化开发，可以留言私信作者；(✿◡‿◡)！！
GDAL C++ API 学习之路 OGRGeometry 多边形类 OGRPolygon
qq_69574549的博客
07-30 738
OGRPolygon是 OGR 几何图形库中的一个类，用于表示多边形的几何图形。它是一种封闭的多边形，由一个或多个外环（）和零个或多个内环（）组成。多边形是平面上的一个封闭区域，可以包含其他几何图形，如点、线、其他多边形等创建一个空多边形OGRPolygonconst&other复制构造函数operatorconst&other赋值运算符inline**begin返回迭代器的开始时间inline**end迭代器的返回结束。
GDAL C++ API 学习之路 OGRGeometry 线类 OGRLineString
qq_69574549的博客
07-30 547
类是 OGR 库中的一个几何对象类，用于表示线段或折线。它由多个坐标点组成，并且在坐标点之间形成线段。可以包含 2D、3D 或 3D+M 坐标点，其中 M 表示额外的度量值，例如时间或速度创建一个空线串const&other复制构造函数operatorconst&other赋值运算符。
GDAL C++ API 学习之路 （1）几何工厂类 OGRGeometryFactory
qq_69574549的博客
07-07 607
OGRGeometryFactory 是 GDAL/OGR 库中的一个工厂类，用于创建和管理 OGR 几何对象。它提供了一组静态方法，用于创建不同类型的几何对象，如点、线、多边形等 通过使用 OGRGeometryFactory，可以方便地创建、复制和销毁 OGR 几何对象，而不需要直接操作对象的内存分配和释放 注：OGRGeometryFactory 是一个静态类，无需实例化即可直接调用其方法。它提供了一种统一的方式来创建和操作 OGR 几何对象，简化了几何对象的创建和管理过程
C++ GDAL/OGR 图层求交集Intersection
Xeon_CC的博客
07-11 4056
GDAL
GDAL.rar_GDAL_c++GDAL模块
09-24
`OGRFeature`和`OGRFeatureDefn`类用于描述特征和其属性，而`OGRGeometry`类则处理几何对象。 7. **栅格处理**：GDAL支持多种栅格处理算法，如重采样、裁剪、镶嵌、统计数据计算等。`GDALRasterBand`类提供了许多...
gdal开发文档,gdal库的使用,C,C++
09-10
在学习GDAL时，配合中文版的开发文档是非常有帮助的，它能指导你理解GDAL的原理，掌握API的使用方法，并解决在实际项目中遇到的问题。阅读`gdal开发文档.pdf`，你将逐步熟悉GDAL库的各个部分，提升你在地理信息系统...
gdal-1.9.2基于vs2010已编译好的c++
04-22
- 调用GDAL的矢量API，如OGRGeometry类，处理矢量数据。 在VS2010环境下，首先需要确保安装了C++的开发环境，并将GDAL的头文件目录和库文件目录添加到项目的编译和链接设置中。这样，就可以在C++代码中包含GDAL的...
gdal224.zip_GIS编程_C/C++_
08-09
OGRGeometry类提供了创建、修改和分析几何对象的方法。 5. **地理编码**：GDAL支持地址到坐标（反向地理编码）和坐标到地址（地理编码）的转换。 6. **数据处理**：GDAL提供了大量的内建函数和算法，用于图像处理...
编译好的gdal库，C++加载调用
01-26
安装过程中，确保包含了C++ API，因为我们需要在C++项目中使用GDAL。 在C++项目中引入GDAL库，你需要包含必要的头文件，如`#include "gdal_priv.h"`和`#include "ogr_api.h"`。GDAL库的核心类是GDALDataset和...
c++ 编译GDAL库
langtaosha_a的专栏
01-19 3511
一、GDAL简介 GDAL(Geospatial Data Abstraction Library)是一个在X/MIT许可协议下的开源栅格空间数据转换库。它利用抽象数据模型来表达所支持的各种文件格式。它还有一系列命令行工具来进行数据转换和处理。 OGR是GDAL项目的一个分支，功能与GDAL类似，只不过它提供对矢量数据的支持。 有很多著名的GIS类产品都使用了GDAL/OGR
c++ 多边形求交集代码（凸多边形与凹多边形的交集）
mingo_敏
05-01 5194
多边形求交集cpp代码 凸多边形与凹多边形的交集： 凸多边形与凸多边形的交集：c++ 多边形求交集代码（凸多边形与凸多边形交集） #include "stdafx.h" #include <iostream> #include <stdlib.h> #include <opencv2\opencv.hpp> using namespace std; usi...
【C++】任意两个多边形相交问题
qq_35504602的博客
07-23 4031
标题
QGIS二次开发07：buffer实现
IamChen的专栏
11-09 818
QGIS 缓冲 QT
GDAL库（c++版）栅格API教程
qq_41342446的博客
03-31 4113
** 栅格API教程 ** 1. 打开文件 在打开GDAL支持的栅格数据存储之前，必须注册驱动程序。每个支持的格式都有一个驱动程序。通常这是通过 GDALAllRegister() 函数，该函数尝试注册所有已知的驱动程序，包括使用 GDALDriverManager::AutoLoadDrivers() . 注册，申请应调用独立 GDALOpen() 函数打开数据集，传递数据集的名称和所需的访问权限（GA_ReadOnly只读或GA_Update）。 #include "gdal_priv.h" #incl
【C++】Windows配置GDAL（附编译结果）
Z_W_H_的博客
02-20 1529
SQLite是一款轻量级的嵌入式关系型数据库管理系统，它遵循ACID（原子性、一致性、隔离性和持久性）原则，并提供了完整的事务处理能力。SQLite由D.Richard Hipp创建，并且是一个开源的项目，其源代码大约有3万行，大小约为250KB。它被设计为紧凑且独立，不依赖于外部库或服务。支持ACID保证数据的一致性和正确性；数据存储在一个单一的磁盘文件中，最大可扩展到2TB；能够在不同的字节顺序之间自由共享数据库文件；占用资源较少，处理速度快；
gdal平面几何空间关系
flyingshineangel的专栏
01-06 543
关于平面几何的空间关系判定，gdal提供了8个函数，分别是：Intersects（相交），Equals（相等），Disjoint（不相交），Touches（接触），Crosses（交叉），Within（被包含），Contains（包含），Overlaps（重叠）。关于几种空间关系的概念，网上有许多权威解释，此处不做解释说明。上图给出了一个多边形poly和一组点F，C，G，通过gdal计算poly和F，C，F的空间关系，得出下面两张表（OP表示poly.OP(point*)，!
gdal c++ 提取ndvi
10-09
GDAL是一个功能强大的地理数据抽象库，用于读取、写入和处理各种栅格数据格式。在使用GDAL提取NDVI（归一化植被指数）时，我们可以按照以下步骤进行操作。 首先，我们需要引入GDAL库，并声明需要使用的数据集和波段。 ``` #include "gdal_priv.h" #include "cpl_conv.h" int main() { //声明需要使用的数据集和波段 GDALDatasetH hDataset; GDALRasterBandH hRedBand, hNirBand; //... return 0; } ``` 接下来，我们需要打开影像文件，并获取红光和近红外波段。 ``` hDataset = GDALOpen("input_image.tif", GA_ReadOnly); hRedBand = GDALGetRasterBand(hDataset, 1); hNirBand = GDALGetRasterBand(hDataset, 2); ``` 然后，我们需要读取红光和近红外波段的数据，并计算NDVI。 ``` int width = GDALGetRasterBandXSize(hRedBand); int height = GDALGetRasterBandYSize(hRedBand); float* redData = (float*) CPLMalloc(sizeof(float)*width*height); float* nirData = (float*) CPLMalloc(sizeof(float)*width*height); GDALRasterIO(hRedBand, GF_Read, 0, 0, width, height, redData, width, height, GDT_Float32, 0, 0); GDALRasterIO(hNirBand, GF_Read, 0, 0, width, height, nirData, width, height, GDT_Float32, 0, 0); float* ndviData = (float*) CPLMalloc(sizeof(float)*width*height); for (int i = 0; i < width*height; i++) { if (redData[i] + nirData[i] == 0) { ndviData[i] = 0; } else{ ndviData[i] = (nirData[i] - redData[i]) / (nirData[i] + redData[i]); } } ``` 最后，我们可以将计算得到的NDVI数据保存到新的影像文件中。 ``` GDALDriverH hDriver = GDALGetDriverByName("GTiff"); GDALDatasetH hOutputDataset = GDALCreateCopy(hDriver, "output_ndvi.tif", hDataset, 0, NULL, NULL, NULL); GDALRasterBandH hOutputBand = GDALGetRasterBand(hOutputDataset, 1); GDALRasterIO(hOutputBand, GF_Write, 0, 0, width, height, ndviData, width, height, GDT_Float32, 0, 0); ``` 上述代码片段展示了如何使用GDAL C API提取NDVI。通过逐像素计算红光和近红外波段之间的比值，我们可以得到NDVI并将其保存为新的栅格影像文件output_ndvi.tif。

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

打赏作者

场主不吃鍋巴

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

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

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

打赏作者

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

抵扣说明：

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

余额充值