expose isValidReason maximumInscribedCircle

Author: jw3126

src/geos_functions.jl

@@ -671,6 +671,22 @@ function lineMerge(obj::Geometry, context::GEOSContext = get_context(obj))
     geomFromGEOS(result, context)
 end
 
+function maximumInscribedCircle(obj::Geometry, tolerance, context::GEOSContext = get_context(obj))
+    result = GEOSMaximumInscribedCircle_r(context, obj, tolerance)
+    if result == C_NULL
+        error("LibGEOS: Error in GEOSMaximumInscribedCircle")
+    end
+    geomFromGEOS(result, context)
+end
+
+function isValidReason(obj::Geometry, context::GEOSContext = get_context(obj))
+    result = GEOSisValidReason_r(context, obj)
+    if result == C_NULL
+        error("LibGEOS: Error in GEOSisValidReason")
+    end
+    result
+end
+
 function simplify(obj::Geometry, tol::Real, context::GEOSContext = get_context(obj))
     result = GEOSSimplify_r(context, obj, tol)
     if result == C_NULL

test/test_geos_functions.jl

@@ -918,11 +918,9 @@ end
     lss = readgeom("LINESTRING EMPTY")
     @test lineMerge(lss) == readgeom("GEOMETRYCOLLECTION EMPTY")
 
-    geo_invalid = readgeom("POLYGON((0 0, 0 1, 1 1, 1 0, -1 1, 0 0))")
-    @test !LibGEOS.isValid(geo_invalid)
-    geo_valid = LibGEOS.makeValid(geo_invalid)
-    @test geo_valid isa LibGEOS.MultiPolygon
-    @test LibGEOS.isValid(geo_valid)
-
     @test LibGEOS.reverse(readgeom("LINESTRING(0 0, 1 1)")) == readgeom("LINESTRING(1 1, 0 0)")
+
+    geo = readgeom("POLYGON((0 0, 0 1, 1 1, 1 0, 0 0))")
+    mic = LibGEOS.maximumInscribedCircle(geo, 1e-4)
+    @test mic == readgeom("LINESTRING (0.5 0.5, 0 0.5)")
 end

test/test_invalid_geometry.jl

@@ -1,3 +1,6 @@
+using Test
+using LibGEOS
+
 @testset "LibGEOS invalid geometry" begin
     # LibGEOS shouldn't crash but error out on invalid geometry
 
@@ -11,4 +14,10 @@
     @test !LibGEOS.isValid(polygon)
     @test LibGEOS.GEOSisValidReason_r(LibGEOS.get_global_context(), polygon) ==
           "Hole lies outside shell[15 15]"
+    @test LibGEOS.isValidReason(polygon) == "Hole lies outside shell[15 15]"
+
+    geo_valid = LibGEOS.makeValid(polygon)
+    @test geo_valid isa LibGEOS.MultiPolygon
+    @test LibGEOS.isValid(geo_valid)
+    @test LibGEOS.isValidReason(geo_valid) == "Valid Geometry"
 end