[geos-commits] [SCM] GEOS branch 3.9 updated. 7c2f400cb870e33f01a84da5308546dd965749bb
git at osgeo.org
git at osgeo.org
Wed Jan 27 13:52:20 PST 2021
This is an automated email from the git hooks/post-receive script. It was
generated because a ref change was pushed to the repository containing
the project "GEOS".
The branch, 3.9 has been updated
via 7c2f400cb870e33f01a84da5308546dd965749bb (commit)
from a5890d8f9ee7737d697e23b19100a86333fcb665 (commit)
Those revisions listed above that are new to this repository have
not appeared on any other notification email; so we list those
revisions in full, below.
- Log -----------------------------------------------------------------
commit 7c2f400cb870e33f01a84da5308546dd965749bb
Author: Martin Davis <mtnclimb at gmail.com>
Date: Mon Jan 25 13:19:56 2021 -0800
Fix OverlayNG line ordering
diff --git a/include/geos/operation/overlayng/EdgeMerger.h b/include/geos/operation/overlayng/EdgeMerger.h
index b67233a..c02eda7 100644
--- a/include/geos/operation/overlayng/EdgeMerger.h
+++ b/include/geos/operation/overlayng/EdgeMerger.h
@@ -59,32 +59,23 @@ namespace overlayng { // geos.operation.overlayng
* This ensures that the overlay output line direction will be as consistent
* as possible with input lines.
*
+ * The merger also preserves the order of the edges in the input.
+ * This means that for polygon-line overlay
+ * the result lines will be in the same order as in the input
+ * (possibly with multiple result lines for a single input line).
+ *
* @author mdavis
*
*/
class GEOS_DLL EdgeMerger {
-private:
-
-
- // Members
- std::vector<Edge*>& edges;
- std::map<EdgeKey, Edge*> edgeMap;
-
public:
- // Methods
- EdgeMerger(std::vector<Edge*>& p_edges);
-
static std::vector<Edge*> merge(std::vector<Edge*>& edges);
- std::vector<Edge*> merge();
-
-
};
} // namespace geos.operation.overlayng
} // namespace geos.operation
} // namespace geos
-
diff --git a/src/operation/overlayng/EdgeMerger.cpp b/src/operation/overlayng/EdgeMerger.cpp
index 605180c..205dcdb 100644
--- a/src/operation/overlayng/EdgeMerger.cpp
+++ b/src/operation/overlayng/EdgeMerger.cpp
@@ -23,23 +23,12 @@ namespace geos { // geos
namespace operation { // geos.operation
namespace overlayng { // geos.operation.overlayng
-EdgeMerger::EdgeMerger(std::vector<Edge*>& p_edges)
- : edges(p_edges) {}
-
/*public static */
std::vector<Edge*>
EdgeMerger::merge(std::vector<Edge*>& edges)
{
- EdgeMerger merger(edges);
- return merger.merge();
-}
-
-
-/*public static */
-std::vector<Edge*>
-EdgeMerger::merge()
-{
std::vector<Edge*> mergedEdges;
+ std::map<EdgeKey, Edge*> edgeMap;
for (Edge* edge : edges) {
EdgeKey edgeKey(edge);
@@ -47,6 +36,7 @@ EdgeMerger::merge()
if (it == edgeMap.end()) {
// this is the first (and maybe only) edge for this line
edgeMap[edgeKey] = edge;
+ mergedEdges.push_back(edge);
//Debug.println("edge added: " + edge);
//Debug.println(edge.toLineString());
}
@@ -69,11 +59,6 @@ EdgeMerger::merge()
//Debug.println(edge.toLineString());
}
}
-
- // copy map values into return vector
- for (auto it: edgeMap) {
- mergedEdges.push_back(it.second);
- }
return mergedEdges;
}
diff --git a/tests/unit/index/chain/MonotoneChainBuilderTest.cpp b/tests/unit/index/chain/MonotoneChainBuilderTest.cpp
index 08e5015..75995d7 100644
--- a/tests/unit/index/chain/MonotoneChainBuilderTest.cpp
+++ b/tests/unit/index/chain/MonotoneChainBuilderTest.cpp
@@ -2,6 +2,7 @@
// Test Suite for geos::index::chain::MonotoneChainBuilder class.
#include <tut/tut.hpp>
+#include <utility.h>
// geos
#include <geos/io/WKTReader.h>
#include <geos/io/WKTWriter.h>
@@ -46,28 +47,14 @@ void object::test<1>
std::string wkt2("POLYGON((0.1 0.1, 4.0 0.1, 4.0 1.9, 0.1 1.9, 0.1 0.1))");
std::unique_ptr<geos::geom::Geometry> g2(wktreader.read(wkt2));
- // POLYGON ((2.0 1.9, 4.0 1.9, 4.0 0.1, 2.0 0.1, 2.0 1.9))
- double x[5] = {2.0, 4.0, 4.0, 2.0, 2.0};
- double y[5] = {1.9, 1.9, 0.1, 0.1, 1.9};
+ auto result = g2->difference(g1.get());
+ result->normalize();
- auto g3 = g2->difference(g1.get());
- auto cs = g3->getCoordinates();
-
- // std::string wkt3 = wktwriter.write(g3.get());
- // std::cout << wkt3 << std::endl;
-
- for (int i = 0; i < 5; ++i)
- {
- auto csx = cs->getOrdinate(i, 0);
- auto csy = cs->getOrdinate(i, 1);
- // std::cout << csx << ", " << csy << std::endl;
- // std::cout << x[i] << ", " << y[i] << std::endl;
- ensure_equals("x", csx, x[i], 0.01);
- ensure_equals("y", csy, y[i], 0.01);
- }
+ std::string wktExpected("POLYGON ((2 0.1, 2 1.9, 4 1.9, 4 0.1, 2 0.1))");
+ std::unique_ptr<geos::geom::Geometry> expected(wktreader.read(wktExpected));
+ ensure_equals_geometry(expected.get(), result.get());
}
} // namespace tut
-
diff --git a/tests/unit/operation/overlayng/OverlayNGTest.cpp b/tests/unit/operation/overlayng/OverlayNGTest.cpp
index 479597c..b753159 100644
--- a/tests/unit/operation/overlayng/OverlayNGTest.cpp
+++ b/tests/unit/operation/overlayng/OverlayNGTest.cpp
@@ -45,6 +45,24 @@ struct test_overlayng_data {
}
void
+ testOverlayExact(const std::string& a, const std::string& b, const std::string& expected, int opCode, double scaleFactor)
+ {
+ std::unique_ptr<PrecisionModel> pm;
+ if (scaleFactor > 0)
+ pm.reset(new PrecisionModel(scaleFactor));
+ else
+ pm.reset(new PrecisionModel());
+
+ std::unique_ptr<Geometry> geom_a = r.read(a);
+ std::unique_ptr<Geometry> geom_b = r.read(b);
+ std::unique_ptr<Geometry> geom_expected = r.read(expected);
+ std::unique_ptr<Geometry> geom_result = OverlayNG::overlay(geom_a.get(), geom_b.get(), opCode, pm.get());
+ // std::string wkt_result = w.write(geom_result.get());
+ // std::cout << std::endl << wkt_result << std::endl;
+ ensure_equals_exact_geometry(geom_expected.get(), geom_result.get(), 0);
+ }
+
+ void
testOverlayNoOpt(const std::string& a, const std::string& b, const std::string& expected, int opCode, double scaleFactor)
{
PrecisionModel pm(scaleFactor);
@@ -531,4 +549,14 @@ void object::test<42> ()
testOverlay(a, b, exp, OverlayNG::INTERSECTION, 0);
}
+template<>
+template<>
+void object::test<43> ()
+{
+ std::string a = "POLYGON ((1 1, 1 9, 9 9, 9 7, 3 7, 3 3, 9 3, 9 1, 1 1))";
+ std::string b = "MULTILINESTRING ((2 10, 2 0), (4 10, 4 0))";
+ std::string exp = "MULTILINESTRING ((2 9, 2 1), (4 9, 4 7), (4 3, 4 1))";
+ testOverlay(a, b, exp, OverlayNG::INTERSECTION, 0);
+}
+
} // namespace tut
diff --git a/tests/unit/utility.h b/tests/unit/utility.h
index 08f14e1..65781e8 100644
--- a/tests/unit/utility.h
+++ b/tests/unit/utility.h
@@ -327,6 +327,49 @@ ensure_equals_geometry_xyz(const T *lhs_in,
ensure_equals_exact_geometry_xyz(g1.get(), g2.get(), tolerance);
}
+/*
+ * Checks for geometries exactly equal in XY only
+ */
+
+template <typename T> inline void ensure_equals_exact_geometry(const T *lhs_in, const T *rhs_in, double tolerance = 0.0);
+
+template <>
+inline void
+ensure_equals_exact_geometry(const geos::geom::Geometry *lhs_in,
+ const geos::geom::Geometry *rhs_in,
+ double tolerance)
+{
+ assert(nullptr != lhs_in);
+ assert(nullptr != rhs_in);
+
+ using geos::geom::Point;
+ using geos::geom::LineString;
+ using geos::geom::Polygon;
+ using geos::geom::CoordinateSequence;
+ using geos::geom::GeometryCollection;
+
+ ensure_equals("type id do not match",
+ lhs_in->getGeometryTypeId(), rhs_in->getGeometryTypeId());
+
+ if (const Point* gpt1 = dynamic_cast<const Point *>(lhs_in)) {
+ const Point *gpt2 = static_cast<const Point *>(rhs_in);
+ return ensure_equals_dims( gpt1->getCoordinatesRO(), gpt2->getCoordinatesRO(), 2, tolerance);
+ }
+ else if (const LineString* gln1 = dynamic_cast<const LineString *>(lhs_in)) {
+ const LineString *gln2 = static_cast<const LineString *>(rhs_in);
+ return ensure_equals_dims( gln1->getCoordinatesRO(), gln2->getCoordinatesRO(), 2, tolerance);
+ }
+ else if (dynamic_cast<const Polygon *>(lhs_in)) {
+ assert("Not implemented yet" == 0);
+ }
+ else if (const GeometryCollection* gc1 = dynamic_cast<const GeometryCollection *>(lhs_in)) {
+ const GeometryCollection *gc2 = static_cast<const GeometryCollection *>(rhs_in);
+ for (unsigned int i = 0; i < gc1->getNumGeometries(); i++) {
+ ensure_equals_exact_geometry(gc1->getGeometryN(i), gc2->getGeometryN(i), tolerance);
+ }
+ }
+}
+
//
// Utility functions
//
@@ -353,4 +396,3 @@ struct wkb_hex_decoder {
} // namespace tut
#endif // #ifndef GEOS_TUT_UTILITY_H_INCLUDED
-
diff --git a/tests/xmltester/XMLTester.cpp b/tests/xmltester/XMLTester.cpp
index 5a747fc..f4d1ac3 100644
--- a/tests/xmltester/XMLTester.cpp
+++ b/tests/xmltester/XMLTester.cpp
@@ -1312,6 +1312,7 @@ XMLTester::parseTest(const tinyxml2::XMLNode* node)
else if(opName == "union") {
GeomPtr gRes(parseGeometry(opRes, "expected"));
+ gRes->normalize();
profile.start();
@@ -1324,6 +1325,7 @@ XMLTester::parseTest(const tinyxml2::XMLNode* node)
}
profile.stop();
+ gRealRes->normalize();
success = checkOverlaySuccess(*gRes, *gRealRes);
@@ -2414,4 +2416,3 @@ main(int argC, char* argV[])
* Revision 1.29 2006/03/17 14:56:39 strk
* Fixed filename normalizer for sql output
**********************************************************************/
-
-----------------------------------------------------------------------
Summary of changes:
include/geos/operation/overlayng/EdgeMerger.h | 19 +++-------
src/operation/overlayng/EdgeMerger.cpp | 19 +---------
.../unit/index/chain/MonotoneChainBuilderTest.cpp | 25 +++---------
tests/unit/operation/overlayng/OverlayNGTest.cpp | 28 ++++++++++++++
tests/unit/utility.h | 44 +++++++++++++++++++++-
tests/xmltester/XMLTester.cpp | 3 +-
6 files changed, 86 insertions(+), 52 deletions(-)
hooks/post-receive
--
GEOS
More information about the geos-commits
mailing list