1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
|
/*
* # GEOS C++ example 1
*
* Reads two WKT representations and calculates the
* intersection, prints it out, and cleans up.
*
* In general, to avoid API changes, stick to operations
* on Geometry. The more esoteric APIs are more likely
* to change between versions.
*/
#include <iostream>
/* For geometry operations */
#include <geos/geom/GeometryFactory.h>
#include <geos/geom/Geometry.h>
/* For WKT read/write */
#include <geos/io/WKTReader.h>
#include <geos/io/WKTWriter.h>
/* Geometry/GeometryFactory */
using geos::geom::Geometry;
using geos::geom::GeometryFactory;
/* WKTReader/WKTWriter */
using geos::io::WKTReader;
using geos::io::WKTWriter;
int main()
{
/* New factory with default (float) precision model */
GeometryFactory::Ptr factory = GeometryFactory::create();
/*
* Reader requires a factory to bind the geometry to
* for shared resources like the PrecisionModel
*/
WKTReader reader(*factory);
/* Input WKT strings */
std::string wkt_a("POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))");
std::string wkt_b("POLYGON((5 5, 15 5, 15 15, 5 15, 5 5))");
/* Convert WKT to Geometry */
std::unique_ptr<Geometry> geom_a(reader.read(wkt_a));
std::unique_ptr<Geometry> geom_b(reader.read(wkt_b));
/* Calculate intersection */
std::unique_ptr<Geometry> inter = geom_a->intersection(geom_b.get());
/* Convert Geometry to WKT */
WKTWriter writer;
writer.setTrim(true);
std::string inter_wkt = writer.write(inter.get());
/* Print out results */
std::cout << "Geometry A: " << wkt_a << std::endl;
std::cout << "Geometry B: " << wkt_b << std::endl;
std::cout << "Intersection(A, B): " << inter_wkt << std::endl;
}
|
