libgeos · dbaston · Dec 13, 2022 · Dec 12, 2022 · Dec 12, 2022 · Dec 12, 2022
diff --git a/NEWS.md b/NEWS.md
@@ -30,6 +30,7 @@ xxxx-xx-xx
   - Improve ConvexHull radial sort robustness (GH-724, Martin Davis)
   - Use more robust Delaunay Triangulation frame size heuristic (GH-728, Martin Davis)
   - DiscreteFrechetDistance: Fix crash with empty inputs (GH-751, Dan Baston)
+  - PreparedLineStringIntersects: Fix incorrect result with mixed-dim collection (GH-774, Dan Baston)
 
 
 ## Changes in 3.11.0

diff --git a/include/geos/geom/Geometry.h b/include/geos/geom/Geometry.h
@@ -342,6 +342,11 @@ class GEOS_DLL Geometry {
     /// Returns the dimension of this Geometry (0=point, 1=line, 2=surface)
     virtual Dimension::DimensionType getDimension() const = 0; //Abstract
 
+    /// Checks whether any component of this geometry has dimension d
+    virtual bool hasDimension(Dimension::DimensionType d) const {
+        return getDimension() == d;
+    }
+
     /// Checks whether this Geometry consists only of components having dimension d.
     virtual bool isDimensionStrict(Dimension::DimensionType d) const {
         return d == getDimension();

diff --git a/include/geos/geom/GeometryCollection.h b/include/geos/geom/GeometryCollection.h
@@ -108,6 +108,8 @@ class GEOS_DLL GeometryCollection : public Geometry {
      */
     Dimension::DimensionType getDimension() const override;
 
+    bool hasDimension(Dimension::DimensionType d) const override;
+
     bool isDimensionStrict(Dimension::DimensionType d) const override;
 
     /// Returns coordinate dimension.

diff --git a/include/geos/geom/MultiLineString.h b/include/geos/geom/MultiLineString.h
@@ -57,6 +57,10 @@ class GEOS_DLL MultiLineString: public GeometryCollection {
     /// Returns line dimension (1)
     Dimension::DimensionType getDimension() const override;
 
+    bool hasDimension(Dimension::DimensionType d) const override {
+        return d == Dimension::L;
+    }
+
     bool isDimensionStrict(Dimension::DimensionType d) const override {
         return d == Dimension::L;
     }

diff --git a/include/geos/geom/MultiPoint.h b/include/geos/geom/MultiPoint.h
@@ -62,6 +62,10 @@ class GEOS_DLL MultiPoint: public GeometryCollection {
         return d == Dimension::P;
     }
 
+    bool hasDimension(Dimension::DimensionType d) const override {
+        return d == Dimension::P;
+    }
+
     /// Returns Dimension::False (Point has no boundary)
     int getBoundaryDimension() const override;
 

diff --git a/include/geos/geom/MultiPolygon.h b/include/geos/geom/MultiPolygon.h
@@ -65,6 +65,10 @@ class GEOS_DLL MultiPolygon: public GeometryCollection {
     /// Returns surface dimension (2)
     Dimension::DimensionType getDimension() const override;
 
+    bool hasDimension(Dimension::DimensionType d) const override {
+        return d == Dimension::A;
+    }
+
     bool isDimensionStrict(Dimension::DimensionType d) const override {
         return d == Dimension::A;
     }

diff --git a/src/geom/GeometryCollection.cpp b/src/geom/GeometryCollection.cpp
@@ -117,6 +117,15 @@ GeometryCollection::isDimensionStrict(Dimension::DimensionType d) const {
             });
 }
 
+bool
+GeometryCollection::hasDimension(Dimension::DimensionType d) const {
+    return std::any_of(geometries.begin(),
+                       geometries.end(),
+                       [&d](const std::unique_ptr<Geometry>& g) {
+        return g->hasDimension(d);
+    });
+}
+
 int
 GeometryCollection::getBoundaryDimension() const
 {

diff --git a/src/geom/prep/PreparedLineStringIntersects.cpp b/src/geom/prep/PreparedLineStringIntersects.cpp
@@ -69,19 +69,14 @@ PreparedLineStringIntersects::intersects(const geom::Geometry* g) const
         return true;
     }
 
-    // For L/L case we are done
-    if(g->getDimension() == 1) {
-        return false;
-    }
-
     // For L/A case, need to check for proper inclusion of the target in the test
     if(g->getDimension() == 2
             &&	prepLine.isAnyTargetComponentInTest(g)) {
         return true;
     }
 
     // For L/P case, need to check if any points lie on line(s)
-    if(g->getDimension() == 0) {
+    if(g->hasDimension(Dimension::P)) {
         return isAnyTestPointInTarget(g);
     }
 

diff --git a/tests/unit/geom/GeometryCollectionTest.cpp b/tests/unit/geom/GeometryCollectionTest.cpp
@@ -154,4 +154,31 @@ void object::test<7>()
     ensure_equals(components.size(), 2u);
 }
 
+// Test of hasDimension()
+template<>
+template<>
+void object::test<8>
+()
+{
+    auto gc = readWKT("GEOMETRYCOLLECTION(POINT (1 1), LINESTRING (1 1, 2 2))");
+
+    ensure(gc->hasDimension(geos::geom::Dimension::P));
+    ensure(gc->hasDimension(geos::geom::Dimension::L));
+    ensure(!gc->hasDimension(geos::geom::Dimension::A));
+}
+
+
+// Test of hasDimension() for nested collection
+template<>
+template<>
+void object::test<9>
+()
+{
+    auto gc = readWKT("GEOMETRYCOLLECTION(POINT (1 1), GEOMETRYCOLLECTION(LINESTRING (1 1, 2 2), POLYGON((0 0, 0 1, 1 1, 0 0))))");
+
+    ensure(gc->hasDimension(geos::geom::Dimension::P));
+    ensure(gc->hasDimension(geos::geom::Dimension::L));
+    ensure(gc->hasDimension(geos::geom::Dimension::A));
+}
+
 } // namespace tut
diff --git a/tests/unit/geom/LineStringTest.cpp b/tests/unit/geom/LineStringTest.cpp
@@ -37,6 +37,7 @@ struct test_linestring_data {
     geos::io::WKTReader reader_;
 
     std::unique_ptr<geos::geom::LineString> empty_line_;
+    std::unique_ptr<geos::geom::LineString> line_;
 
     test_linestring_data()
         : pm_(1000)
@@ -45,6 +46,11 @@ struct test_linestring_data {
         , empty_line_(factory_->createLineString())
     {
         assert(nullptr != empty_line_);
+
+        auto cs = geos::detail::make_unique<geos::geom::CoordinateSequence>(2u, false, false);
+        cs->setAt(geos::geom::Coordinate{0, 0}, 0);
+        cs->setAt(geos::geom::Coordinate{1, 1}, 1);
+        line_ = factory_->createLineString(std::move(cs));
     }
 
     ~test_linestring_data()
@@ -571,6 +577,16 @@ void object::test<31>
 
 }
 
+// Test of hasDimension()
+template<>
+template<>
+void object::test<32>
+()
+{
+    ensure(!line_->hasDimension(geos::geom::Dimension::P));
+    ensure(line_->hasDimension(geos::geom::Dimension::L));
+    ensure(!line_->hasDimension(geos::geom::Dimension::A));
+}
 
 } // namespace tut
 
diff --git a/tests/unit/geom/MultiLineStringTest.cpp b/tests/unit/geom/MultiLineStringTest.cpp
@@ -63,5 +63,18 @@ void object::test<3>
     ensure(!mls_->isDimensionStrict(geos::geom::Dimension::A));
 }
 
+
+// Test of hasDimension()
+template<>
+template<>
+void object::test<4>
+()
+{
+    ensure(!mls_->hasDimension(geos::geom::Dimension::P));
+    ensure(mls_->hasDimension(geos::geom::Dimension::L));
+    ensure(!mls_->hasDimension(geos::geom::Dimension::A));
+}
+
+
 } // namespace tut
 
diff --git a/tests/unit/geom/MultiPointTest.cpp b/tests/unit/geom/MultiPointTest.cpp
@@ -417,5 +417,16 @@ void object::test<31>
     ensure(!empty_mp_->isDimensionStrict(geos::geom::Dimension::L));
 }
 
+// Test of hasDimension()
+template<>
+template<>
+void object::test<32>
+()
+{
+    ensure(mp_->hasDimension(geos::geom::Dimension::P));
+    ensure(!mp_->hasDimension(geos::geom::Dimension::L));
+    ensure(!mp_->hasDimension(geos::geom::Dimension::A));
+}
+
 } // namespace tut
 
diff --git a/tests/unit/geom/MultiPolygonTest.cpp b/tests/unit/geom/MultiPolygonTest.cpp
@@ -62,4 +62,15 @@ void object::test<3>
     ensure(!empty_mp_->isDimensionStrict(geos::geom::Dimension::L));
 }
 
+// Test of hasDimension()
+template<>
+template<>
+void object::test<4>
+()
+{
+    ensure(!mp_->hasDimension(geos::geom::Dimension::P));
+    ensure(!mp_->hasDimension(geos::geom::Dimension::L));
+    ensure(mp_->hasDimension(geos::geom::Dimension::A));
+}
+
 } // namespace tut
diff --git a/tests/unit/geom/PointTest.cpp b/tests/unit/geom/PointTest.cpp
@@ -18,6 +18,7 @@
 #include <memory>
 #include <string>
 
+constexpr int MAX_TESTS = 100;
 
 namespace tut {
 //
@@ -59,7 +60,7 @@ struct test_point_data {
     }
 };
 
-typedef test_group<test_point_data> group;
+typedef test_group<test_point_data, MAX_TESTS> group;
 typedef group::object object;
 
 group test_point_group("geos::geom::Point");
@@ -622,5 +623,16 @@ void object::test<47>
     ensure(pt->getEnvelopeInternal()->contains(7, 9));
 }
 
+// Test of hasDimension()
+template<>
+template<>
+void object::test<48>
+()
+{
+    ensure(point_->hasDimension(geos::geom::Dimension::P));
+    ensure(!point_->hasDimension(geos::geom::Dimension::L));
+    ensure(!point_->hasDimension(geos::geom::Dimension::A));
+}
+
 } // namespace tut
 
diff --git a/tests/unit/geom/PolygonTest.cpp b/tests/unit/geom/PolygonTest.cpp
@@ -650,4 +650,15 @@ void object::test<44>()
     ensure_equals(holes.size(), 1u);
 }
 
+// Test of hasDimension()
+template<>
+template<>
+void object::test<45>
+()
+{
+    ensure(!poly_->hasDimension(geos::geom::Dimension::P));
+    ensure(!poly_->hasDimension(geos::geom::Dimension::L));
+    ensure(poly_->hasDimension(geos::geom::Dimension::A));
+}
+
 } // namespace tut
diff --git a/tests/xmltester/tests/general/TestPreparedPredicatesWithGeometryCollection.xml b/tests/xmltester/tests/general/TestPreparedPredicatesWithGeometryCollection.xml
@@ -76,6 +76,17 @@
 <test>  <op name="intersects"     arg1="A" arg2="B">   true  </op> </test>
 </case>
 
+<case>
+  <desc>LineString against GC, with point at endpoint
+  </desc>
+  <a>
+  LINESTRING (0 0, 1 1)
+  </a>
+  <b>
+    GEOMETRYCOLLECTION (POINT (1 1), LINESTRING (2 2, 3 3))
+  </b>
+<test>  <op name="intersects"     arg1="A" arg2="B">   true  </op> </test>
+</case>
 
 
 </run>