最近，想要在C++下编一个可用GDAL模块的地图管理系统，找来找去，找到了GEOS。GEOS（Geometry Engine-Open Source）开源几何引擎 是一个用于计算几何的JTS库的 C/C++实现，专注于地理信息系统 （GIS） 软件中使用的算法。 它实现了 OGC 简单要素几何模型，并提供了该标准中的所有空间函数以及许多其他标准。 GEOS 是 PostGIS、QGIS、GDAL、Shapely 等的核心依赖项。GEOS库可以和GDAL库一起编译提供OGR库的所有功能。希望了解GEOS更多信息，可以前往https://trac.osgeo.org/geos/。
本文GEOS库所采用的是QT软件平台下的MinGW32编译器进行编译，编译出来的库文件也用于Windows系统下QT（MinGW）软件平台进行开发。编译后的文件主要有两个：libgeos.dll和libgeos_c.dll，编译出来以后经过测试，可以正常使用。
https://libgeos.org/usage/download/二、导入文件：启动QT Creator，新建一工程，选择工程模板为“Non-Qt Project”——“plain C++ Application”，构建系统选择“CMake”，构建套件(kit)选择“Desktop Qt 5.15.2 MinGW 32-bit”，一路点”下一步“，生成工程框架。将解压后的目录”geos-3.12.1“下面的所有子目录和文件导入到工程目录下面，并删除自动生成的”main.cpp“文件。
三、分析CMakeLists.txt文件并编译：经过仔细对工程目录下的所有CMakeLists.txt文件分析，未发现有误之处，便不再修改CMakeLists.txt文件，进行下一步构建工程。经过4分多钟的等待，最终生成了”libgeos.dll、libgeos_c.dll“两个dll文件和3 个测试可执行文件。编译成功！
四、调用GEOS库前设置：新建一工程，选择工程模板为“Application(Qt)”——“Qt Witgets Application”，构建系统选择“CMake”，构建套件(kit)选择“Desktop Qt 5.15.2 MinGW 32-bit”，一路点”下一步“，生成工程框架。在工程目录下新建“lib”和“include”两个子目录，将”libgeos.dll、libgeos_c.dll“两个dll文件复制到“lib”目录，将解压后的目录”geos-3.12.1“下面的学习“include/geos/”下面的所有文件拷贝至工程目录下的“include”中，至此，调用前设置工作完成。
五、调用GEOS库：使用官方网站的C++API帮助实例，代码如下：
#include <iostream>
#include "mygdalwidget.h"
#include "./ui_mygdalwidget.h"
#include "geos/geom/GeometryFactory.h"
#include "geos/geom/Geometry.h"
#include "geos/io/WKTReader.h"
#include "geos/io/WKTWriter.h"
#include <QMessageBox>
using namespace std;
using namespace geos::geom;
using namespace geos::io;
myGdalWidget::myGdalWidget(QWidget *parent)
: QWidget(parent)
, ui(new Ui::myGdalWidget)
{
ui->setupUi(this);
}
myGdalWidget::~myGdalWidget()
{
delete ui;
}
void myGdalWidget::on_btnConnectGDAL_clicked()
{
GeometryFactory::Ptr fact = GeometryFactory::create();
WKTReader reader(*fact);
string wkt_a("POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))");
string wkt_b("POLYGON((5 5, 15 5, 15 15, 5 15, 5 5))");
/* Convert WKT to Geometry */
unique_ptr<Geometry> geom_a(reader.read(wkt_a));
unique_ptr<Geometry> geom_b(reader.read(wkt_b));
/* Calculate intersection */
unique_ptr<Geometry> inter = geom_a->intersection(geom_b.get());
/* Convert Geometry to WKT */
WKTWriter writer;
writer.setTrim(true); /* Only needed before GEOS 3.12 */
QString inter_wkt = QString::fromStdString(writer.write(inter.get()));
/* Print out results */
cout << "Geometry A: " << wkt_a << endl;
cout << "Geometry B: " << wkt_b << endl;
cout << "Intersection(A, B): << inter_wkt.toStdString()<<endl;
QMessageBox::information(this,"显示两个图形相交",inter_wkt);
}
代码输入完毕，点构建运行，应用程序输出：
11:54:07: Starting F:\QtProjects\build-firstGDALPro-Desktop_Qt_5_15_2_MinGW_32_bit-Debug\firstGDALPro.exe...
Geometry A: POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))
Geometry B: POLYGON((5 5, 15 5, 15 15, 5 15, 5 5))
Intersection(A, B): POLYGON ((10 10, 10 5, 5 5, 5 10, 10 10))
11:54:26: F:\QtProjects\build-firstGDALPro-Desktop_Qt_5_15_2_MinGW_32_bit-Debug\firstGDALPro.exe 退出，退出代码: 0
至此，成功运行！
