[geos-devel] [GEOS] #1057: GEOSClipByRect might return mixed dimension geometries

[GEOS]

#1057: GEOSClipByRect might return mixed dimension geometries
------------------------+--------------------------
 Reporter:  Algunenano  |      Owner:  geos-devel@…
     Type:  defect      |     Status:  new
 Priority:  major       |  Milestone:  3.8.2
Component:  Default     |    Version:  3.8.0
 Severity:  Unassigned  |   Keywords:
------------------------+--------------------------
 Comes from Postgis' https://trac.osgeo.org/postgis/ticket/4690

 Reproducer:
 {{{
 diff --git a/tests/unit/capi/GEOSClipByRectTest.cpp
 b/tests/unit/capi/GEOSClipByRectTest.cpp
 index e5358436..f2b0ad12 100644
 --- a/tests/unit/capi/GEOSClipByRectTest.cpp
 +++ b/tests/unit/capi/GEOSClipByRectTest.cpp
 @@ -200,5 +200,24 @@ template<> template<> void object::test<13>
      isEqual(geom2_, "POLYGON ((5 5, 5 15, 10 15, 10 10, 15 10, 15 5, 5
 5))");
  }

 -} // namespace tut
 +/// Comes from https://trac.osgeo.org/geos/ticket/1056
 +template<> template<> void object::test<15>
 +()
 +{
 +    const char* wkt = "LINESTRING(2542 4287 5,2540 4299
 5.361999999965,2537 4305 5.910000000033,2530 4311 6.832999999984,2519 4314
 8.670000000042,2511 4318 10.351999999955,2504 4325 12.388000000035,2499
 4331 14.146999999997,2497 4337 15.354999999981,2497 4363
 18.652000000002,2484 4382 19.369999999995,2481 4390 19.572999999975,2480
 4394 19.675999999978,2485 4406 20,2484 4410 20,2480 4415
 20.123999999952,2466 4423 20.464999999967,2464 4425 20.503999999957,2458
 4438 21.265000000014,2453 4443 21.55700000003,2451 4444
 21.589000000036,2446 4446 21.898000000045,2441 4450 22.373000000021,2438
 4455 22.876000000047,2437 4460 23.35699999996,2435 4475
 24.327999999979,2432 4479 24.689999999944,2428 4483 24.974999999977,2423
 4484 25.236000000033,2419 4483 25.587000000058,2414 4481
 25.886000000057,2408 4480 26.34600000002,2404 4483 26.780999999959,2401
 4488 27.089999999967,2398 4491 27.447000000044,2392 4493
 28.158999999985,2385 4494 28.90399999998,2379 4492 29.648000000045,2375
 4489 30.082999999984,2370 4481 30.947000000044,2368 4478
 31.459999999963,2340 4479 34.521999999997,2320 4483 38.736999999965,2294
 4491 42,2289 4493 42,2283 4498 42.108000000008,2277 4504
 42.302000000025,2269 4517 43.079000000027,2263 4522 43.53899999999,2261
 4523 43.706000000006,2255 4523 44.201000000001,2243 4521 45.5,2239 4517
 45.886000000057,2234 4502 47.400999999954,2213 4462 49.908999999985,2211
 4454 50.160999999964,2210 4443 50.631999999983,2211 4429
 51.402000000002,2214 4418 52.089000000036,2212 4414 52.373000000021,2199
 4404 53.283999999985,2188 4397 54,2184 4397 54.295000000042,2179 4398
 54.626999999979,2166 4405 55,2159 4405 55,2155 4402 55.070999999996,2153
 4390 55.351000000024,2149 4381 55.976000000024,2138 4334 60.5,2136 4330
 60.895000000019,2131 4324 61.405999999959,2121 4315 62.069999999949,2103
 4290 63.123999999952,2097 4284 63.302999999956,2090 4279
 63.515999999945,2085 4273 63.805000000051,2076 4266 64.10699999996,2073
 4262 64.185000000056,2072 4260 64.197999999975,2072 4253
 64.396999999997,2072 4242 64.701000000001,2074 4228 65,2075 4208 65,2071
 4202 65,2064 4197 65,2061 4193 65,2061 4186 65,2063 4179 65,2064 4169
 65,2060 4154 65,2056 4145 65.222999999998,2052 4140 65.363999999943,2047
 4139 65.493000000017,2039 4140 65.812000000035,2023 4146 66,2021 4146
 66,2013 4140 66,2011 4139 66,2009 4140 66,2005 4143 66,2003 4143 66,2001
 4142 66,1990 4130 66,1986 4124 66,1984 4118 66)";
 +    geom1_ = GEOSGeomFromWKT(wkt);
 +    geom2_ = GEOSClipByRect(geom1_, -256,256,4352,4352);
 +
 +    for (int i = 0; i < GEOSGetNumGeometries(geom2_); i++)
 +    {
 +        const GEOSGeometry* g = GEOSGetGeometryN(geom2_, i);
 +        bool hasz = GEOSHasZ(g);
 +        if (!hasz) {
 +            std::printf("Expected subgeom %d to have z-coordinate and it
 doesn't", i);
 +        }
 +        ensure(hasz);
 +    }
 +    ensure(GEOSHasZ(geom2_) == true);
 +}

 +} // namespace tut
 \ No newline at end of file
 }}}

 The result of this GEOSClipByRect is a multilinestring where the first
 (index 0) linestring has 3 coordinates (xyz) but the second one only has 2
 (xy). I would expect to either have everything 2D or everything 3D, but
 AFAIK having mixed dimensions in a multigeometry is invalid (at least in
 Postgis it is).

-- 
Ticket URL: <https://trac.osgeo.org/geos/ticket/1057>
GEOS <http://trac.osgeo.org/geos>
GEOS (Geometry Engine - Open Source) is a C++ port of the Java Topology Suite (JTS).