ffi-geos 1.2.0 → 2.2.0

data/test/geometry_collection_tests.rb

@@ -1,7 +1,5 @@
-# encoding: UTF-8
 # frozen_string_literal: true
 
-$: << File.dirname(__FILE__)
 require 'test_helper'
 
 class GeometryCollectionTests < Minitest::Test
@@ -18,7 +16,7 @@ class GeometryCollectionTests < Minitest::Test
     geom = read('GEOMETRYCOLLECTION(POINT(0 0))')
     assert_kind_of(Enumerable, geom.each)
     assert_kind_of(Enumerable, geom.to_enum)
-    assert_equal(geom, geom.each {})
+    assert_equal(geom, geom.each(&EMPTY_BLOCK))
   end
 
   def test_geometry_collection_array
@@ -40,7 +38,7 @@ class GeometryCollectionTests < Minitest::Test
       'POINT (10 20)'
     ], geom[0, 2].collect { |g| write(g) })
 
-    assert_equal(nil, geom[0, -1])
+    assert_nil(geom[0, -1])
     assert_equal([], geom[-1, 0])
     assert_equal([
       'POINT (10 20)',
@@ -85,4 +83,390 @@ class GeometryCollectionTests < Minitest::Test
     geom = read('GEOMETRYCOLLECTION (POINT(1 2), LINESTRING(1 2, 3 4))')
     assert_equal(2, geom.num_geometries)
   end
+
+  def test_x_max
+    geom = read('GEOMETRYCOLLECTION (
+      POLYGON ((0 0, 5 0, 8 9, -10 5, 0 0)),
+      LINESTRING (0 0, 5 0, 8 9, -10 5, 0 0),
+      POINT(3 12)
+    )')
+
+    assert_equal(8, geom.x_max)
+  end
+
+  def test_x_min
+    geom = read('GEOMETRYCOLLECTION (
+      POLYGON ((0 0, 5 0, 8 9, -10 5, 0 0)),
+      LINESTRING (0 0, 5 0, 8 9, -10 5, 0 0),
+      POINT(3 12)
+    )')
+
+    assert_equal(-10, geom.x_min)
+  end
+
+  def test_y_max
+    geom = read('GEOMETRYCOLLECTION (
+      POLYGON ((0 0, 5 0, 8 9, -10 5, 0 0)),
+      LINESTRING (0 0, 5 0, 8 9, -10 5, 0 0),
+      POINT(3 12)
+    )')
+
+    assert_equal(12, geom.y_max)
+  end
+
+  def test_y_min
+    geom = read('GEOMETRYCOLLECTION (
+      POLYGON ((0 0, 5 0, 8 9, -10 5, 0 0)),
+      LINESTRING (0 0, 5 0, 8 9, -10 5, 0 0),
+      POINT(3 12)
+    )')
+
+    assert_equal(0, geom.y_min)
+  end
+
+  def test_z_max
+    geom = read('GEOMETRYCOLLECTION (
+      POLYGON ((0 0, 5 0, 8 9, -10 5, 0 0)),
+      LINESTRING (0 0, 5 0, 8 9, -10 5, 0 0),
+      POINT(3 12)
+    )')
+    assert_equal(0, geom.z_max)
+
+    geom = read('GEOMETRYCOLLECTION Z (
+      POLYGON Z ((0 0 0, 5 0 3, 8 9 4, -10 5 3, 0 0 0)),
+      LINESTRING Z (0 0 0, 5 0 3, 8 9 4, -10 5 3, 0 0 0),
+      POINT Z (3 12 6)
+    )')
+    assert_equal(6, geom.z_max)
+
+    # GEOS lets you mix dimensionality, while PostGIS doesn't.
+    geom = read('GEOMETRYCOLLECTION (
+      POLYGON ((0 0, 5 0, 8 9, -10 5, 0 0)),
+      LINESTRING (0 0, 5 0, 8 9, -10 5, 0 0),
+      POINT(3 12 10)
+    )')
+    assert_equal(10, geom.z_max)
+  end
+
+  def test_z_min
+    geom = read('GEOMETRYCOLLECTION (
+      POLYGON ((0 0, 5 0, 8 9, -10 5, 0 0)),
+      LINESTRING (0 0, 5 0, 8 9, -10 5, 0 0),
+      POINT(3 12)
+    )')
+    assert_equal(0, geom.z_min)
+
+    geom = read('GEOMETRYCOLLECTION Z (
+      POLYGON Z ((0 0 0, 5 0 3, 8 9 4, -10 5 3, 0 0 0)),
+      LINESTRING Z (0 0 0, 5 0 3, 8 9 4, -10 5 3, 0 0 0),
+      POINT Z (3 12 6)
+    )')
+    assert_equal(0, geom.z_min)
+
+    # GEOS lets you mix dimensionality, while PostGIS doesn't.
+    geom = read('GEOMETRYCOLLECTION (
+      POLYGON ((0 0, 5 0, 8 9, -10 5, 0 0)),
+      LINESTRING (0 0, 5 0, 8 9, -10 5, 0 0),
+      POINT(3 12 -10)
+    )')
+    assert_equal(-10, geom.z_min)
+  end
+
+  def test_snap_to_grid
+    wkt = 'GEOMETRYCOLLECTION (LINESTRING (-10.12 0, -10.12 5, -10.12 5, -10.12 6, -10.12 6, -10.12 6, -10.12 7, -10.12 7, -10.12 7, -10.12 8, -10.12 8, -9 8, -9 9, -10.12 0), ' \
+      'POLYGON ((-10.12 0, -10.12 5, -10.12 5, -10.12 6, -10.12 6, -10.12 6, -10.12 7, -10.12 7, -10.12 7, -10.12 8, -10.12 8, -9 8, -9 9, -10.12 0)), POINT (10.12 10.12))'
+
+    expected = 'GEOMETRYCOLLECTION (LINESTRING (-10 0, -10 5, -10 5, -10 6, -10 6, -10 6, -10 7, -10 7, -10 7, -10 8, -10 8, -9 8, -9 9, -10 0), POLYGON ((-10 0, -10 5, -10 5, -10 6, -10 6, -10 6, -10 7, -10 7, -10 7, -10 8, -10 8, -9 8, -9 9, -10 0)), POINT (10 10))'
+
+    simple_bang_tester(:snap_to_grid, expected, wkt, 1)
+  end
+
+  def test_snap_to_grid_empty
+    assert(read('GEOMETRYCOLLECTION EMPTY').snap_to_grid!.empty?, 'Expected an empty GeometryCollection')
+  end
+
+  def test_snap_to_grid_with_srid
+    wkt = 'GEOMETRYCOLLECTION (
+      LINESTRING (-10.12 0, -10.12 5, -10.12 5, -10.12 6, -10.12 6, -10.12 6, -10.12 7, -10.12 7, -10.12 7, -10.12 8, -10.12 8, -9 8, -9 9, -10.12 0),
+      POLYGON ((-10.12 0, -10.12 5, -10.12 5, -10.12 6, -10.12 6, -10.12 6, -10.12 7, -10.12 7, -10.12 7, -10.12 8, -10.12 8, -9 8, -9 9, -10.12 0)),
+      POINT (10.12 10.12)
+    )'
+
+    expected = 'GEOMETRYCOLLECTION (LINESTRING (-10 0, -10 5, -10 5, -10 6, -10 6, -10 6, -10 7, -10 7, -10 7, -10 8, -10 8, -9 8, -9 9, -10 0), POLYGON ((-10 0, -10 5, -10 5, -10 6, -10 6, -10 6, -10 7, -10 7, -10 7, -10 8, -10 8, -9 8, -9 9, -10 0)), POINT (10 10))'
+
+    srid_copy_tester(:snap_to_grid, expected, 0, :zero, wkt, 1)
+    srid_copy_tester(:snap_to_grid, expected, 4326, :lenient, wkt, 1)
+    srid_copy_tester(:snap_to_grid, expected, 4326, :strict, wkt, 1)
+  end
+
+  def test_snap_to_grid_with_illegal_result
+    assert_raises(Geos::InvalidGeometryError) do
+      read('GEOMETRYCOLLECTION (POINT (0 2), LINESTRING (0 1, 0 11), POLYGON ((0 1, 0 1, 0 6, 0 6, 0 1)))')
+        .snap_to_grid(1)
+    end
+  end
+
+  def test_rotate
+    writer.rounding_precision = 3
+
+    wkt = 'GEOMETRYCOLLECTION (POINT (1 1), LINESTRING (0 0, 10 10), POLYGON ((0 0, 5 0, 5 5, 0 5, 0 0)))'
+
+    affine_tester(:rotate,
+      'GEOMETRYCOLLECTION (POINT (29 11), LINESTRING (30 10, 20 20), POLYGON ((30 10, 30 15, 25 15, 25 10, 30 10)))',
+      wkt,
+      Math::PI / 2,
+      [10.0, 20.0])
+
+    affine_tester(:rotate,
+      'GEOMETRYCOLLECTION (POINT (-2 0), LINESTRING (-3 1, 7 -9), POLYGON ((-3 1, -3 -4, 2 -4, 2 1, -3 1)))',
+      wkt,
+      -Math::PI / 2,
+      [-1.0, 2.0])
+
+    affine_tester(:rotate,
+      'GEOMETRYCOLLECTION (POINT (19 1), LINESTRING (20 0, 10 10), POLYGON ((20 0, 20 5, 15 5, 15 0, 20 0)))',
+      wkt,
+      Math::PI / 2,
+      read('POINT(10 10)'))
+
+    affine_tester(:rotate,
+      'GEOMETRYCOLLECTION (POINT (-0.5 0.5), LINESTRING (0.5 -0.5, -9.5 9.5), POLYGON ((0.5 -0.5, 0.5 4.5, -4.5 4.5, -4.5 -0.5, 0.5 -0.5)))',
+      wkt,
+      Math::PI / 2,
+      read('LINESTRING(0 0, 1 0)'))
+  end
+
+  def test_rotate_x
+    writer.rounding_precision = 3
+    writer.output_dimensions = 3
+
+    wkt = 'GEOMETRYCOLLECTION Z (POINT Z (1 1 1), LINESTRING Z (1 1 1, 10 10 10), POLYGON Z ((0 0 0, 5 0 0, 5 5 0, 0 5 0, 0 0 0)))'
+
+    affine_tester(:rotate_x,
+      'GEOMETRYCOLLECTION Z (POINT Z (1 -1 -1), LINESTRING Z (1 -1 -1, 10 -10 -10), POLYGON Z ((0 0 0, 5 0 0, 5 -5 0, 0 -5 0, 0 0 0)))',
+      wkt,
+      Math::PI)
+
+    affine_tester(:rotate_x,
+      'GEOMETRYCOLLECTION Z (POINT Z (1 -1 1), LINESTRING Z (1 -1 1, 10 -10 10), POLYGON Z ((0 0 0, 5 0 0, 5 0 5, 0 0 5, 0 0 0)))',
+      wkt,
+      Math::PI / 2)
+
+    affine_tester(:rotate_x,
+      'GEOMETRYCOLLECTION Z (POINT Z (1 1 -1), LINESTRING Z (1 1 -1, 10 10 -10), POLYGON Z ((0 0 0, 5 0 0, 5 0 -5, 0 0 -5, 0 0 0)))',
+      wkt,
+      Math::PI + Math::PI / 2)
+
+    affine_tester(:rotate_x,
+      wkt,
+      wkt,
+      Math::PI * 2)
+  end
+
+  def test_rotate_y
+    writer.rounding_precision = 6
+    writer.output_dimensions = 3
+
+    wkt = 'GEOMETRYCOLLECTION Z (POINT Z (1 1 1), LINESTRING Z (1 1 1, 10 10 10), POLYGON Z ((0 0 0, 5 0 0, 5 5 0, 0 5 0, 0 0 0)))'
+
+    affine_tester(:rotate_y,
+      'GEOMETRYCOLLECTION Z (POINT Z (-1 1 -1), LINESTRING Z (-1 1 -1, -10 10 -10), POLYGON Z ((0 0 0, -5 0 0, -5 5 0, 0 5 0, 0 0 0)))',
+      wkt,
+      Math::PI)
+
+    affine_tester(:rotate_y,
+      'GEOMETRYCOLLECTION Z (POINT Z (1 1 -1), LINESTRING Z (1 1 -1, 10 10 -10), POLYGON Z ((0 0 0, 0 0 -5, 0 5 -5, 0 5 0, 0 0 0)))',
+      wkt,
+      Math::PI / 2)
+
+    affine_tester(:rotate_y,
+      'GEOMETRYCOLLECTION Z (POINT Z (-1 1 1), LINESTRING Z (-1 1 1, -10 10 10), POLYGON Z ((0 0 0, 0 0 5, 0 5 5, 0 5 0, 0 0 0)))',
+      wkt,
+      Math::PI + Math::PI / 2)
+
+    affine_tester(:rotate_y,
+      wkt,
+      wkt,
+      Math::PI * 2)
+  end
+
+  def test_rotate_z
+    writer.rounding_precision = 3
+
+    wkt = 'GEOMETRYCOLLECTION (POINT (1 1), LINESTRING (0 0, 10 10), POLYGON ((0 0, 5 0, 5 5, 0 5, 0 0)))'
+
+    affine_tester(:rotate_z,
+      'GEOMETRYCOLLECTION (POINT (-1 -1), LINESTRING (0 0, -10 -10), POLYGON ((0 0, -5 0, -5 -5, 0 -5, 0 0)))',
+      wkt,
+      Math::PI)
+
+    affine_tester(:rotate_z,
+      'GEOMETRYCOLLECTION (POINT (-1 1), LINESTRING (0 0, -10 10), POLYGON ((0 0, 0 5, -5 5, -5 0, 0 0)))',
+      wkt,
+      Math::PI / 2)
+
+    affine_tester(:rotate_z,
+      'GEOMETRYCOLLECTION (POINT (1 -1), LINESTRING (0 0, 10 -10), POLYGON ((0 0, 0 -5, 5 -5, 5 0, 0 0)))',
+      wkt,
+      Math::PI + Math::PI / 2)
+
+    affine_tester(:rotate_z,
+      wkt,
+      wkt,
+      Math::PI * 2)
+  end
+
+  def test_scale
+    affine_tester(:scale,
+      'GEOMETRYCOLLECTION (POINT (5 5), LINESTRING (0 0, 50 50), POLYGON ((0 0, 25 0, 25 25, 0 25, 0 0)))',
+      'GEOMETRYCOLLECTION (POINT (1 1), LINESTRING (0 0, 10 10), POLYGON ((0 0, 5 0, 5 5, 0 5, 0 0)))',
+      5,
+      5)
+
+    affine_tester(:scale,
+      'GEOMETRYCOLLECTION (POINT (3 2), LINESTRING (0 0, 30 20), POLYGON ((0 0, 15 0, 15 10, 0 10, 0 0)))',
+      'GEOMETRYCOLLECTION (POINT (1 1), LINESTRING (0 0, 10 10), POLYGON ((0 0, 5 0, 5 5, 0 5, 0 0)))',
+      3,
+      2)
+
+    writer.output_dimensions = 3
+    affine_tester(:scale,
+      'GEOMETRYCOLLECTION Z (POINT Z (4 2 -8), LINESTRING Z (4 2 -8, 40 20 -80), POLYGON Z ((0 0 0, 20 0 0, 20 10 0, 0 10 0, 0 0 0)))',
+      'GEOMETRYCOLLECTION Z (POINT Z (1 1 1), LINESTRING Z (1 1 1, 10 10 10), POLYGON Z ((0 0 0, 5 0 0, 5 5 0, 0 5 0, 0 0 0)))',
+      4,
+      2,
+      -8)
+  end
+
+  def test_scale_hash
+    affine_tester(:scale,
+      'GEOMETRYCOLLECTION (POINT (5 5), LINESTRING (0 0, 50 50), POLYGON ((0 0, 25 0, 25 25, 0 25, 0 0)))',
+      'GEOMETRYCOLLECTION (POINT (1 1), LINESTRING (0 0, 10 10), POLYGON ((0 0, 5 0, 5 5, 0 5, 0 0)))',
+      x: 5,
+      y: 5)
+
+    affine_tester(:scale,
+      'GEOMETRYCOLLECTION (POINT (3 2), LINESTRING (0 0, 30 20), POLYGON ((0 0, 15 0, 15 10, 0 10, 0 0)))',
+      'GEOMETRYCOLLECTION (POINT (1 1), LINESTRING (0 0, 10 10), POLYGON ((0 0, 5 0, 5 5, 0 5, 0 0)))',
+      x: 3,
+      y: 2)
+
+    writer.output_dimensions = 3
+    affine_tester(:scale,
+      'GEOMETRYCOLLECTION Z (POINT Z (4 2 -8), LINESTRING Z (4 2 -8, 40 20 -80), POLYGON Z ((0 0 0, 20 0 0, 20 10 0, 0 10 0, 0 0 0)))',
+      'GEOMETRYCOLLECTION Z (POINT Z (1 1 1), LINESTRING Z (1 1 1, 10 10 10), POLYGON Z ((0 0 0, 5 0 0, 5 5 0, 0 5 0, 0 0 0)))',
+      x: 4,
+      y: 2,
+      z: -8)
+  end
+
+  def test_trans_scale
+    affine_tester(:trans_scale,
+      'GEOMETRYCOLLECTION (POINT (3 3), LINESTRING (2 2, 12 12), POLYGON ((2 2, 7 2, 7 7, 2 7, 2 2)))',
+      'GEOMETRYCOLLECTION (POINT (2 2), LINESTRING (1 1, 11 11), POLYGON ((1 1, 6 1, 6 6, 1 6, 1 1)))',
+      1, 1, 1, 1)
+
+    wkt = 'GEOMETRYCOLLECTION (POINT (1 1), LINESTRING (0 0, 10 10), POLYGON ((0 0, 5 0, 5 5, 0 5, 0 0)))'
+
+    affine_tester(:trans_scale,
+      'GEOMETRYCOLLECTION (POINT (2 2), LINESTRING (1 1, 11 11), POLYGON ((1 1, 6 1, 6 6, 1 6, 1 1)))',
+      wkt,
+      1, 1, 1, 1)
+
+    affine_tester(:trans_scale,
+      'GEOMETRYCOLLECTION (POINT (0 0), LINESTRING (1 1, -9 -9), POLYGON ((1 1, -4 1, -4 -4, 1 -4, 1 1)))',
+      wkt,
+      -1, -1, -1, -1)
+
+    affine_tester(:trans_scale,
+      'GEOMETRYCOLLECTION (POINT (1 2), LINESTRING (0 1, 10 11), POLYGON ((0 1, 5 1, 5 6, 0 6, 0 1)))',
+      wkt,
+      0, 1, 1, 1)
+
+    affine_tester(:trans_scale,
+      'GEOMETRYCOLLECTION (POINT (2 1), LINESTRING (1 0, 11 10), POLYGON ((1 0, 6 0, 6 5, 1 5, 1 0)))',
+      wkt,
+      1, 0, 1, 1)
+
+    affine_tester(:trans_scale,
+      'GEOMETRYCOLLECTION (POINT (3 2), LINESTRING (2 1, 12 11), POLYGON ((2 1, 7 1, 7 6, 2 6, 2 1)))',
+      wkt,
+      2, 1, 1, 1)
+
+    affine_tester(:trans_scale,
+      'GEOMETRYCOLLECTION (POINT (2 3), LINESTRING (1 2, 11 12), POLYGON ((1 2, 6 2, 6 7, 1 7, 1 2)))',
+      wkt,
+      1, 2, 1, 1)
+
+    affine_tester(:trans_scale,
+      'GEOMETRYCOLLECTION (POINT (4 2), LINESTRING (2 1, 22 11), POLYGON ((2 1, 12 1, 12 6, 2 6, 2 1)))',
+      wkt,
+      1, 1, 2, 1)
+
+    affine_tester(:trans_scale,
+      'GEOMETRYCOLLECTION (POINT (2 4), LINESTRING (1 2, 11 22), POLYGON ((1 2, 6 2, 6 12, 1 12, 1 2)))',
+      wkt,
+      1, 1, 1, 2)
+
+    affine_tester(:trans_scale,
+      'GEOMETRYCOLLECTION (POINT (15 28), LINESTRING (10 21, 60 91), POLYGON ((10 21, 35 21, 35 56, 10 56, 10 21)))',
+      wkt,
+      2, 3, 5, 7)
+
+    writer.output_dimensions = 3
+    affine_tester(:trans_scale,
+      'GEOMETRYCOLLECTION Z (POINT Z (15 28 1), LINESTRING Z (15 28 1, 60 91 10), POLYGON Z ((10 21 0, 35 21 0, 35 56 0, 10 56 0, 10 21 0)))',
+      'GEOMETRYCOLLECTION Z (POINT Z (1 1 1), LINESTRING Z (1 1 1, 10 10 10), POLYGON Z ((0 0 0, 5 0 0, 5 5 0, 0 5 0, 0 0 0)))',
+      2, 3, 5, 7)
+  end
+
+  def test_trans_scale_hash
+    affine_tester(:trans_scale,
+      'GEOMETRYCOLLECTION (POINT (2 2), LINESTRING (1 1, 11 11), POLYGON ((1 1, 6 1, 6 6, 1 6, 1 1)))',
+      'GEOMETRYCOLLECTION (POINT (1 1), LINESTRING (0 0, 10 10), POLYGON ((0 0, 5 0, 5 5, 0 5, 0 0)))',
+      delta_x: 1, delta_y: 1, x_factor: 1, y_factor: 1)
+
+    writer.output_dimensions = 3
+    affine_tester(:trans_scale,
+      'GEOMETRYCOLLECTION Z (POINT Z (15 28 1), LINESTRING Z (15 28 1, 60 91 10), POLYGON Z ((10 21 0, 35 21 0, 35 56 0, 10 56 0, 10 21 0)))',
+      'GEOMETRYCOLLECTION Z (POINT Z (1 1 1), LINESTRING Z (1 1 1, 10 10 10), POLYGON Z ((0 0 0, 5 0 0, 5 5 0, 0 5 0, 0 0 0)))',
+      delta_x: 2, delta_y: 3, x_factor: 5, y_factor: 7)
+
+    affine_tester(:trans_scale,
+      'GEOMETRYCOLLECTION Z (POINT Z (3 1 1), LINESTRING Z (3 1 1, 12 10 10), POLYGON Z ((2 0 0, 7 0 0, 7 5 0, 2 5 0, 2 0 0)))',
+      'GEOMETRYCOLLECTION Z (POINT Z (1 1 1), LINESTRING Z (1 1 1, 10 10 10), POLYGON Z ((0 0 0, 5 0 0, 5 5 0, 0 5 0, 0 0 0)))',
+      delta_x: 2, z_factor: 2)
+  end
+
+  def test_translate
+    affine_tester(:translate,
+      'GEOMETRYCOLLECTION (POINT (6 13), LINESTRING (5 12, 15 22), POLYGON ((5 12, 10 12, 10 17, 5 17, 5 12)))',
+      'GEOMETRYCOLLECTION (POINT (1 1), LINESTRING (0 0, 10 10), POLYGON ((0 0, 5 0, 5 5, 0 5, 0 0)))',
+      5,
+      12)
+
+    writer.output_dimensions = 3
+    affine_tester(:translate,
+      'GEOMETRYCOLLECTION Z (POINT Z (-2 -6 4), LINESTRING Z (-2 -6 4, 7 3 13), POLYGON Z ((-3 -7 3, 2 -7 3, 2 -2 3, -3 -2 3, -3 -7 3)))',
+      'GEOMETRYCOLLECTION Z (POINT Z (1 1 1), LINESTRING Z (1 1 1, 10 10 10), POLYGON Z ((0 0 0, 5 0 0, 5 5 0, 0 5 0, 0 0 0)))',
+      -3,
+      -7,
+      3)
+  end
+
+  def test_translate_hash
+    affine_tester(:translate,
+      'GEOMETRYCOLLECTION (POINT (6 13), LINESTRING (5 12, 15 22), POLYGON ((5 12, 10 12, 10 17, 5 17, 5 12)))',
+      'GEOMETRYCOLLECTION (POINT (1 1), LINESTRING (0 0, 10 10), POLYGON ((0 0, 5 0, 5 5, 0 5, 0 0)))',
+      x: 5,
+      y: 12)
+
+    writer.output_dimensions = 3
+    affine_tester(:translate,
+      'GEOMETRYCOLLECTION Z (POINT Z (-2 -6 4), LINESTRING Z (-2 -6 4, 7 3 13), POLYGON Z ((-3 -7 3, 2 -7 3, 2 -2 3, -3 -2 3, -3 -7 3)))',
+      'GEOMETRYCOLLECTION Z (POINT Z (1 1 1), LINESTRING Z (1 1 1, 10 10 10), POLYGON Z ((0 0 0, 5 0 0, 5 5 0, 0 5 0, 0 0 0)))',
+      x: -3,
+      y: -7,
+      z: 3)
+  end
 end

data/test/geometry_tests.rb

@@ -1,7 +1,5 @@
-# encoding: UTF-8
 # frozen_string_literal: true
 
-$: << File.dirname(__FILE__)
 require 'test_helper'
 
 class GeometryTests < Minitest::Test
@@ -15,15 +13,37 @@ class GeometryTests < Minitest::Test
   def test_intersection
     comparison_tester(
       :intersection,
-      'POLYGON ((5 10, 10 10, 10 5, 5 5, 5 10))',
+      if Geos::GEOS_VERSION > '3.9.0'
+        'POLYGON ((10 10, 10 5, 5 5, 5 10, 10 10))'
+      else
+        'POLYGON ((5 10, 10 10, 10 5, 5 5, 5 10))'
+      end,
       'POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0))',
       'POLYGON ((5 5, 15 5, 15 15, 5 15, 5 5))'
     )
   end
 
-  def test_buffer
-    writer.rounding_precision = 2
+  def test_intersection_with_precision
+    skip unless ENV['FORCE_TESTS'] || Geos::FFIGeos.respond_to?(:GEOSIntersectionPrec_r)
+
+    comparison_tester(
+      :intersection,
+      'GEOMETRYCOLLECTION (POLYGON ((1 2, 1 1, 0.5 1, 1 2)), POLYGON ((9.5 1, 2 1, 2 2, 9 2, 9.5 1)), LINESTRING (1 1, 2 1), LINESTRING (2 2, 1 2))',
+      'MULTIPOLYGON(((0 0,5 10,10 0,0 0),(1 1,1 2,2 2,2 1,1 1),(100 100,100 102,102 102,102 100,100 100)))',
+      'POLYGON((0 1,0 2,10 2,10 1,0 1))',
+      precision: 0
+    )
 
+    comparison_tester(
+      :intersection,
+      'GEOMETRYCOLLECTION (LINESTRING (2 0, 4 0), POINT (0 0), POINT (10 0))',
+      'LINESTRING(0 0, 10 0)',
+      'LINESTRING(9 0, 12 0, 12 20, 4 0, 2 0, 2 10, 0 10, 0 -10)',
+      precision: 2
+    )
+  end
+
+  def test_buffer
     simple_tester(
       :buffer,
       'POLYGON EMPTY',
@@ -31,95 +51,95 @@ class GeometryTests < Minitest::Test
       0
     )
 
-    simple_tester(
+    snapped_tester(
       :buffer,
-      'POLYGON ((10 0, 9.8 -2, 9.2 -3.8, 8.3 -5.6, 7.1 -7.1, 5.6 -8.3, 3.8 -9.2, 2 -9.8, 1.6e-14 -10, -2 -9.8, -3.8 -9.2, -5.6 -8.3, -7.1 -7.1, -8.3 -5.6, -9.2 -3.8, -9.8 -2, -10 -3.2e-14, -9.8 2, -9.2 3.8, -8.3 5.6, -7.1 7.1, -5.6 8.3, -3.8 9.2, -2 9.8, -3.7e-14 10, 2 9.8, 3.8 9.2, 5.6 8.3, 7.1 7.1, 8.3 5.6, 9.2 3.8, 9.8 2, 10 0))',
+      'POLYGON ((10 0, 10 -2, 9 -4, 8 -6, 7 -7, 6 -8, 4 -9, 2 -10, 0 -10, -2 -10, -4 -9, -6 -8, -7 -7, -8 -6, -9 -4, -10 -2, -10 0, -10 2, -9 4, -8 6, -7 7, -6 8, -4 9, -2 10, 0 10, 2 10, 4 9, 6 8, 7 7, 8 6, 9 4, 10 2, 10 0))',
       'POINT(0 0)',
       10
     )
 
     # One segment per quadrant
-    simple_tester(
+    snapped_tester(
       :buffer,
-      'POLYGON ((10 0, 1.6e-14 -10, -10 -3.2e-14, -4.6e-14 10, 10 0))',
+      'POLYGON ((10 0, 0 -10, -10 0, 0 10, 10 0))',
       'POINT(0 0)',
       10,
-      { :quad_segs => 1 }
+      quad_segs: 1
     )
 
     # End cap styles
-    simple_tester(
+    snapped_tester(
       :buffer,
-      'POLYGON ((1e+02 10, 1.1e+02 0, 1e+02 -10, 0 -10, -10 1.2e-15, 0 10, 1e+02 10))',
+      'POLYGON ((100 10, 110 0, 100 -10, 0 -10, -10 0, 0 10, 100 10))',
       'LINESTRING(0 0, 100 0)',
       10,
-      { :quad_segs => 1, :endcap => :round }
+      quad_segs: 1, endcap: :round
     )
 
-    simple_tester(
+    snapped_tester(
       :buffer,
-      'POLYGON ((1e+02 10, 1e+02 -10, 0 -10, 0 10, 1e+02 10))',
+      'POLYGON ((100 10, 100 -10, 0 -10, 0 10, 100 10))',
       'LINESTRING(0 0, 100 0)',
       10,
-      { :quad_segs => 1, :endcap => :flat }
+      quad_segs: 1, endcap: :flat
     )
 
-    simple_tester(
+    snapped_tester(
       :buffer,
-      'POLYGON ((1e+02 10, 1.1e+02 10, 1.1e+02 -10, 0 -10, -10 -10, -10 10, 1e+02 10))',
+      'POLYGON ((100 10, 110 10, 110 -10, 0 -10, -10 -10, -10 10, 100 10))',
       'LINESTRING(0 0, 100 0)',
       10,
-      { :quad_segs => 1, :endcap => :square }
+      quad_segs: 1, endcap: :square
     )
 
     # Join styles
-    simple_tester(
+    snapped_tester(
       :buffer,
-      'POLYGON ((90 10, 90 1e+02, 93 1.1e+02, 1e+02 1.1e+02, 1.1e+02 1.1e+02, 1.1e+02 1e+02, 1.1e+02 0, 1.1e+02 -7.1, 1e+02 -10, 0 -10, -7.1 -7.1, -10 1.2e-15, -7.1 7.1, 0 10, 90 10))',
+      'POLYGON ((90 10, 90 100, 93 107, 100 110, 107 107, 110 100, 110 0, 107 -7, 100 -10, 0 -10, -7 -7, -10 0, -7 7, 0 10, 90 10))',
       'LINESTRING(0 0, 100 0, 100 100)',
       10,
-      { :quad_segs => 2, :join => :round }
+      quad_segs: 2, join: :round
     )
 
-    simple_tester(
+    snapped_tester(
       :buffer,
-      'POLYGON ((90 10, 90 1e+02, 93 1.1e+02, 1e+02 1.1e+02, 1.1e+02 1.1e+02, 1.1e+02 1e+02, 1.1e+02 0, 1e+02 -10, 0 -10, -7.1 -7.1, -10 1.2e-15, -7.1 7.1, 0 10, 90 10))',
+      'POLYGON ((90 10, 90 100, 93 107, 100 110, 107 107, 110 100, 110 0, 100 -10, 0 -10, -7 -7, -10 0, -7 7, 0 10, 90 10))',
       'LINESTRING(0 0, 100 0, 100 100)',
       10,
-      { :quad_segs => 2, :join => :bevel }
+      quad_segs: 2, join: :bevel
     )
 
-    simple_tester(
+    snapped_tester(
       :buffer,
-      'POLYGON ((90 10, 90 1e+02, 93 1.1e+02, 1e+02 1.1e+02, 1.1e+02 1.1e+02, 1.1e+02 1e+02, 1.1e+02 -10, 0 -10, -7.1 -7.1, -10 1.2e-15, -7.1 7.1, 0 10, 90 10))',
+      'POLYGON ((90 10, 90 100, 93 107, 100 110, 107 107, 110 100, 110 -10, 0 -10, -7 -7, -10 0, -7 7, 0 10, 90 10))',
       'LINESTRING(0 0, 100 0, 100 100)',
       10,
-      { :quad_segs => 2, :join => :mitre }
+      quad_segs: 2, join: :mitre
     )
 
-    simple_tester(
+    snapped_tester(
       :buffer,
-      'POLYGON ((90 10, 90 1e+02, 93 1.1e+02, 1e+02 1.1e+02, 1.1e+02 1.1e+02, 1.1e+02 1e+02, 1.1e+02 -5, 1e+02 -9.1, 0 -10, -7.1 -7.1, -10 1.2e-15, -7.1 7.1, 0 10, 90 10))',
+      'POLYGON ((90 10, 90 100, 93 107, 100 110, 107 107, 110 100, 109 -5, 105 -9, 0 -10, -7 -7, -10 0, -7 7, 0 10, 90 10))',
       'LINESTRING(0 0, 100 0, 100 100)',
       10,
-      { :quad_segs => 2, :join => :mitre, :mitre_limit => 1.0 }
+      quad_segs: 2, join: :mitre, mitre_limit: 1.0
     )
 
     # Single-sided buffering
-    simple_tester(
+    snapped_tester(
       :buffer,
-      'POLYGON ((1e+02 0, 0 0, 0 10, 1e+02 10, 1e+02 0))',
+      'POLYGON ((100 0, 0 0, 0 10, 100 10, 100 0))',
       'LINESTRING(0 0, 100 0)',
       10,
-      { :single_sided => true }
+      single_sided: true
     )
 
-    simple_tester(
+    snapped_tester(
       :buffer,
-      'POLYGON ((0 0, 1e+02 0, 1e+02 -10, 0 -10, 0 0))',
+      'POLYGON ((0 0, 100 0, 100 -10, 0 -10, 0 0))',
       'LINESTRING(0 0, 100 0)',
       -10,
-      { :single_sided => true }
+      single_sided: true
     )
   end
 
@@ -137,7 +157,7 @@ class GeometryTests < Minitest::Test
   def test_difference
     comparison_tester(
       :difference,
-      'GEOMETRYCOLLECTION EMPTY',
+      EMPTY_GEOMETRY,
       'POINT(0 0)',
       'POINT(0 0)'
     )
@@ -158,7 +178,7 @@ class GeometryTests < Minitest::Test
 
     comparison_tester(
       :difference,
-      'GEOMETRYCOLLECTION EMPTY',
+      EMPTY_GEOMETRY,
       'POINT(5 0)',
       'LINESTRING(0 0, 10 0)'
     )
@@ -186,31 +206,55 @@ class GeometryTests < Minitest::Test
 
     comparison_tester(
       :difference,
-      'POLYGON ((5 0, 0 0, 0 10, 5 10, 10 10, 10 0, 5 0))',
+      if Geos::GEOS_VERSION > '3.9.0'
+        'POLYGON ((0 10, 5 10, 10 10, 10 0, 5 0, 0 0, 0 10))'
+      else
+        'POLYGON ((0 0, 0 10, 5 10, 10 10, 10 0, 5 0, 0 0))'
+      end,
       'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
       'LINESTRING(5 -10, 5 10)'
     )
 
     comparison_tester(
       :difference,
-      'POLYGON ((10 0, 0 0, 0 10, 10 10, 10 0))',
+      if Geos::GEOS_VERSION > '3.9.0'
+        'POLYGON ((0 10, 10 10, 10 0, 0 0, 0 10))'
+      else
+        'POLYGON ((0 0, 0 10, 10 10, 10 0, 0 0))'
+      end,
       'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
       'LINESTRING(10 0, 20 0)'
     )
 
     comparison_tester(
       :difference,
-      'POLYGON ((5 0, 0 0, 0 10, 10 10, 10 5, 5 5, 5 0))',
+      if Geos::GEOS_VERSION > '3.9.0'
+        'POLYGON ((0 10, 10 10, 10 5, 5 5, 5 0, 0 0, 0 10))'
+      else
+        'POLYGON ((0 0, 0 10, 10 10, 10 5, 5 5, 5 0, 0 0))'
+      end,
       'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
       'POLYGON((5 -5, 5 5, 15 5, 15 -5, 5 -5))'
     )
   end
 
+  def test_difference_with_precision
+    skip unless ENV['FORCE_TESTS'] || Geos::FFIGeos.respond_to?(:GEOSDifferencePrec_r)
+
+    comparison_tester(
+      :difference,
+      'MULTILINESTRING ((2 8, 4 8), (6 8, 10 8))',
+      'LINESTRING (2 8, 10 8)',
+      'LINESTRING (3.9 8.1, 6.1 7.9)',
+      precision: 2
+    )
+  end
+
   def test_sym_difference
     %w{ sym_difference symmetric_difference }.each do |method|
       comparison_tester(
         method,
-        'GEOMETRYCOLLECTION EMPTY',
+        EMPTY_GEOMETRY,
         'POINT(0 0)',
         'POINT(0 0)'
       )
@@ -259,27 +303,51 @@ class GeometryTests < Minitest::Test
 
       comparison_tester(
         method,
-        'GEOMETRYCOLLECTION (LINESTRING (5 -10, 5 0), POLYGON ((5 0, 0 0, 0 10, 5 10, 10 10, 10 0, 5 0)))',
+        if Geos::GEOS_VERSION > '3.9.0'
+          'GEOMETRYCOLLECTION (POLYGON ((0 10, 5 10, 10 10, 10 0, 5 0, 0 0, 0 10)), LINESTRING (5 -10, 5 0))'
+        else
+          'GEOMETRYCOLLECTION (LINESTRING (5 -10, 5 0), POLYGON ((0 0, 0 10, 5 10, 10 10, 10 0, 5 0, 0 0)))'
+        end,
         'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
         'LINESTRING(5 -10, 5 10)'
       )
 
       comparison_tester(
         method,
-        'GEOMETRYCOLLECTION (LINESTRING (10 0, 20 0), POLYGON ((10 0, 0 0, 0 10, 10 10, 10 0)))',
+        if Geos::GEOS_VERSION > '3.9.0'
+          'GEOMETRYCOLLECTION (POLYGON ((0 10, 10 10, 10 0, 0 0, 0 10)), LINESTRING (10 0, 20 0))'
+        else
+          'GEOMETRYCOLLECTION (LINESTRING (10 0, 20 0), POLYGON ((0 0, 0 10, 10 10, 10 0, 0 0)))'
+        end,
         'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
         'LINESTRING(10 0, 20 0)'
       )
 
       comparison_tester(
         method,
-        'MULTIPOLYGON (((5 0, 0 0, 0 10, 10 10, 10 5, 5 5, 5 0)), ((5 0, 10 0, 10 5, 15 5, 15 -5, 5 -5, 5 0)))',
+        if Geos::GEOS_VERSION > '3.9.0'
+          'MULTIPOLYGON (((0 10, 10 10, 10 5, 5 5, 5 0, 0 0, 0 10)), ((10 0, 10 5, 15 5, 15 -5, 5 -5, 5 0, 10 0)))'
+        else
+          'MULTIPOLYGON (((0 0, 0 10, 10 10, 10 5, 5 5, 5 0, 0 0)), ((5 0, 10 0, 10 5, 15 5, 15 -5, 5 -5, 5 0)))'
+        end,
         'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
         'POLYGON((5 -5, 5 5, 15 5, 15 -5, 5 -5))'
       )
     end
   end
 
+  def test_sym_difference_with_precision
+    skip unless ENV['FORCE_TESTS'] || Geos::FFIGeos.respond_to?(:GEOSSymDifferencePrec_r)
+
+    comparison_tester(
+      :sym_difference,
+      'GEOMETRYCOLLECTION (POLYGON ((0 10, 6 10, 10 10, 10 0, 6 0, 0 0, 0 10)), LINESTRING (6 -10, 6 0))',
+      'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
+      'LINESTRING(5 -10, 5 10)',
+      precision: 2
+    )
+  end
+
   def test_boundary
     simple_tester(
       :boundary,
@@ -352,32 +420,59 @@ class GeometryTests < Minitest::Test
 
     comparison_tester(
       :union,
-      'GEOMETRYCOLLECTION (LINESTRING (5 -10, 5 0), POLYGON ((5 0, 0 0, 0 10, 5 10, 10 10, 10 0, 5 0)))',
+      if Geos::GEOS_VERSION > '3.9.0'
+        'GEOMETRYCOLLECTION (POLYGON ((0 10, 5 10, 10 10, 10 0, 5 0, 0 0, 0 10)), LINESTRING (5 -10, 5 0))'
+      else
+        'GEOMETRYCOLLECTION (LINESTRING (5 -10, 5 0), POLYGON ((0 0, 0 10, 5 10, 10 10, 10 0, 5 0, 0 0)))'
+      end,
       'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
       'LINESTRING(5 -10, 5 10)'
     )
 
     comparison_tester(
       :union,
-      'GEOMETRYCOLLECTION (LINESTRING (10 0, 20 0), POLYGON ((10 0, 0 0, 0 10, 10 10, 10 0)))',
+      if Geos::GEOS_VERSION > '3.9.0'
+        'GEOMETRYCOLLECTION (POLYGON ((0 10, 10 10, 10 0, 0 0, 0 10)), LINESTRING (10 0, 20 0))'
+      else
+        'GEOMETRYCOLLECTION (LINESTRING (10 0, 20 0), POLYGON ((0 0, 0 10, 10 10, 10 0, 0 0)))'
+      end,
       'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
       'LINESTRING(10 0, 20 0)'
     )
 
     comparison_tester(
       :union,
-      'POLYGON ((5 0, 0 0, 0 10, 10 10, 10 5, 15 5, 15 -5, 5 -5, 5 0))',
+      if Geos::GEOS_VERSION > '3.9.0'
+        'POLYGON ((0 10, 10 10, 10 5, 15 5, 15 -5, 5 -5, 5 0, 0 0, 0 10))'
+      else
+        'POLYGON ((0 0, 0 10, 10 10, 10 5, 15 5, 15 -5, 5 -5, 5 0, 0 0))'
+      end,
       'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))',
       'POLYGON((5 -5, 5 5, 15 5, 15 -5, 5 -5))'
     )
   end
 
+  def test_union_with_precision
+    skip unless ENV['FORCE_TESTS'] || Geos::FFIGeos.respond_to?(:GEOSUnionPrec_r)
+
+    geom_a = read('POINT (1.9 8.2)')
+    geom_b = read('POINT (4.1 9.8)')
+
+    result = geom_a.union(geom_b, precision: 2)
+
+    assert_equal('MULTIPOINT (2 8, 4 10)', write(result))
+  end
+
   def test_union_cascaded
     skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:union_cascaded)
 
     simple_tester(
       :union_cascaded,
-      'POLYGON ((1 0, 0 0, 0 1, 0 11, 10 11, 10 14, 14 14, 14 10, 11 10, 11 0, 1 0), (11 11, 12 11, 12 12, 11 12, 11 11))',
+      if Geos::GEOS_VERSION > '3.9.0'
+        'POLYGON ((0 0, 0 1, 0 11, 10 11, 10 14, 14 14, 14 10, 11 10, 11 0, 1 0, 0 0), (12 12, 11 12, 11 11, 12 11, 12 12))'
+      else
+        'POLYGON ((1 0, 0 0, 0 1, 0 11, 10 11, 10 14, 14 14, 14 10, 11 10, 11 0, 1 0), (11 11, 12 11, 12 12, 11 12, 11 11))'
+      end,
       'MULTIPOLYGON(
         ((0 0, 1 0, 1 1, 0 1, 0 0)),
         ((10 10, 10 14, 14 14, 14 10, 10 10),
@@ -387,12 +482,33 @@ class GeometryTests < Minitest::Test
     )
   end
 
+  def test_coverage_union
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:coverage_union)
+
+    simple_tester(
+      :union_cascaded,
+      if Geos::GEOS_VERSION > '3.9.0'
+        'POLYGON ((0 1, 1 1, 2 1, 2 0, 1 0, 0 0, 0 1))'
+      else
+        'POLYGON ((0 0, 0 1, 1 1, 2 1, 2 0, 1 0, 0 0))'
+      end,
+      'MULTIPOLYGON(
+        ((0 0, 0 1, 1 1, 1 0, 0 0)),
+        ((1 0, 1 1, 2 1, 2 0, 1 0))
+      ))'
+    )
+  end
+
   def test_unary_union
     skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:unary_union)
 
     simple_tester(
       :unary_union,
-      'POLYGON ((1 0, 0 0, 0 1, 0 11, 10 11, 10 14, 14 14, 14 10, 11 10, 11 0, 1 0), (11 11, 12 11, 12 12, 11 12, 11 11))',
+      if Geos::GEOS_VERSION > '3.9.0'
+        'POLYGON ((0 0, 0 1, 0 11, 10 11, 10 14, 14 14, 14 10, 11 10, 11 0, 1 0, 0 0), (12 12, 11 12, 11 11, 12 11, 12 12))'
+      else
+        'POLYGON ((1 0, 0 0, 0 1, 0 11, 10 11, 10 14, 14 14, 14 10, 11 10, 11 0, 1 0), (11 11, 12 11, 12 12, 11 12, 11 11))'
+      end,
       'MULTIPOLYGON(
         ((0 0, 1 0, 1 1, 0 1, 0 0)),
         ((10 10, 10 14, 14 14, 14 10, 10 10),
@@ -402,6 +518,22 @@ class GeometryTests < Minitest::Test
     )
   end
 
+  def test_unary_union_with_precision
+    skip unless ENV['FORCE_TESTS'] || Geos::FFIGeos.respond_to?(:GEOSUnaryUnionPrec_r)
+
+    simple_tester(
+      :unary_union,
+      'POLYGON ((0 0, 0 12, 9 12, 9 15, 15 15, 15 9, 12 9, 12 0, 0 0))',
+      'MULTIPOLYGON(
+        ((0 0, 1 0, 1 1, 0 1, 0 0)),
+        ((10 10, 10 14, 14 14, 14 10, 10 10),
+        (11 11, 11 12, 12 12, 12 11, 11 11)),
+        ((0 0, 11 0, 11 11, 0 11, 0 0))
+      ))',
+      3
+    )
+  end
+
   def test_node
     skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:node)
 
@@ -415,7 +547,11 @@ class GeometryTests < Minitest::Test
   def test_union_without_arguments
     simple_tester(
       :union,
-      'POLYGON ((1 0, 0 0, 0 1, 0 11, 10 11, 10 14, 14 14, 14 10, 11 10, 11 0, 1 0), (11 11, 12 11, 12 12, 11 12, 11 11))',
+      if Geos::GEOS_VERSION > '3.9.0'
+        'POLYGON ((0 0, 0 1, 0 11, 10 11, 10 14, 14 14, 14 10, 11 10, 11 0, 1 0, 0 0), (12 12, 11 12, 11 11, 12 11, 12 12))'
+      else
+        'POLYGON ((1 0, 0 0, 0 1, 0 11, 10 11, 10 14, 14 14, 14 10, 11 10, 11 0, 1 0), (11 11, 12 11, 12 12, 11 12, 11 11))'
+      end,
       'MULTIPOLYGON(
         ((0 0, 1 0, 1 1, 0 1, 0 0)),
         ((10 10, 10 14, 14 14, 14 10, 10 10),
@@ -511,10 +647,7 @@ class GeometryTests < Minitest::Test
         'LINESTRING(0 0, 10 10)'
       )
 
-      writer.trim = true
-      writer.rounding_precision = 0
-
-      simple_tester(
+      snapped_tester(
         method,
         'POINT (5 4)',
         'POLYGON((0 0, 0 10, 5 5, 10 10, 10 0, 0 0))'
@@ -522,6 +655,17 @@ class GeometryTests < Minitest::Test
     end
   end
 
+  def test_minimum_bounding_circle
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:minimum_bounding_circle)
+
+    geom = read('LINESTRING(0 10, 0 20)')
+
+    assert_equal(
+      'POLYGON ((5 15, 5 14, 5 13, 4 12, 4 11, 3 11, 2 10, 1 10, 0 10, -1 10, -2 10, -3 11, -4 11, -4 12, -5 13, -5 14, -5 15, -5 16, -5 17, -4 18, -4 19, -3 19, -2 20, -1 20, 0 20, 1 20, 2 20, 3 19, 4 19, 4 18, 5 17, 5 16, 5 15))',
+      write(geom.minimum_bounding_circle.snap_to_grid(1))
+    )
+  end
+
   def test_envelope
     simple_tester(
       :envelope,
@@ -656,7 +800,7 @@ class GeometryTests < Minitest::Test
       tests.each do |test|
         expected, method, args = test
         if ENV['FORCE_TESTS'] || geom_a.respond_to?(method)
-          value = geom_a.send(method, *([ geom_b ] + Array(args)))
+          value = geom_a.send(method, *([geom_b] + Array(args)))
           assert_equal(expected, value)
         end
       end
@@ -795,10 +939,10 @@ class GeometryTests < Minitest::Test
   def test_valid_reason
     skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:valid_reason)
 
-    assert_equal("Valid Geometry", read('POINT(0 0)').valid_reason)
-    assert_equal("Invalid Coordinate[0 nan]", read('POINT(0 NaN)').valid_reason)
-    assert_equal("Invalid Coordinate[0 nan]", read('POINT(0 nan)').valid_reason)
-    assert_equal("Self-intersection[2.5 5]", read('POLYGON((0 0, 0 5, 5 5, 5 10, 0 0))').valid_reason)
+    assert_equal('Valid Geometry', read('POINT(0 0)').valid_reason)
+    assert_equal('Invalid Coordinate[0 nan]', read('POINT(0 NaN)').valid_reason)
+    assert_equal('Invalid Coordinate[0 nan]', read('POINT(0 nan)').valid_reason)
+    assert_equal('Self-intersection[2.5 5]', read('POLYGON((0 0, 0 5, 5 5, 5 10, 0 0))').valid_reason)
   end
 
   def test_valid_detail
@@ -811,14 +955,16 @@ class GeometryTests < Minitest::Test
     }
 
     assert_nil(read('POINT(0 0)').valid_detail)
-    tester["Invalid Coordinate", 'POINT (0 nan)', 'POINT(0 NaN)', 0]
-    tester["Self-intersection", 'POINT (2.5 5)', 'POLYGON((0 0, 0 5, 5 5, 5 10, 0 0))', 0]
+    tester['Invalid Coordinate', 'POINT (0 nan)', 'POINT(0 NaN)', 0]
+    tester['Self-intersection', 'POINT (2.5 5)', 'POLYGON((0 0, 0 5, 5 5, 5 10, 0 0))', 0]
 
-    tester["Ring Self-intersection", 'POINT (0 0)', 'POLYGON((0 0, -10 10, 10 10, 0 0, 4 5, -4 5, 0 0)))', 0]
+    tester['Ring Self-intersection', 'POINT (0 0)', 'POLYGON((0 0, -10 10, 10 10, 0 0, 4 5, -4 5, 0 0)))', 0]
 
-    assert_nil(read('POLYGON((0 0, -10 10, 10 10, 0 0, 4 5, -4 5, 0 0)))').valid_detail(
-      :allow_selftouching_ring_forming_hole
-    ))
+    assert_nil(
+      read('POLYGON((0 0, -10 10, 10 10, 0 0, 4 5, -4 5, 0 0)))').valid_detail(
+        :allow_selftouching_ring_forming_hole
+      )
+    )
   end
 
   def test_simple
@@ -847,8 +993,7 @@ class GeometryTests < Minitest::Test
     simple_tester(:num_geometries, 2, 'MULTIPOLYGON(
       ((0 0, 1 0, 1 1, 0 1, 0 0)),
       ((10 10, 10 14, 14 14, 14 10, 10 10),
-      (11 11, 11 12, 12 12, 12 11, 11 11)))'
-    )
+      (11 11, 11 12, 12 12, 12 11, 11 11)))')
     simple_tester(:num_geometries, 6, 'GEOMETRYCOLLECTION (
       MULTIPOLYGON (
         ((0 0, 1 0, 1 1, 0 1, 0 0)),
@@ -859,8 +1004,7 @@ class GeometryTests < Minitest::Test
       MULTILINESTRING ((0 0, 2 3), (10 10, 3 4)),
       LINESTRING (0 0, 2 3),
       MULTIPOINT (0 0, 2 3),
-      POINT (9 0))'
-    )
+      POINT (9 0))')
   end
 
   # get_geometry_n is segfaulting in the binary GEOS build
@@ -881,8 +1025,7 @@ class GeometryTests < Minitest::Test
     simple_tester(:num_interior_rings, 2, 'POLYGON (
       (10 10, 10 14, 14 14, 14 10, 10 10),
       (11 11, 11 12, 12 12, 12 11, 11 11),
-      (13 11, 13 12, 13.5 12, 13.5 11, 13 11))'
-    )
+      (13 11, 13 12, 13.5 12, 13.5 11, 13 11))')
 
     assert_raises(NoMethodError) do
       read('POINT (0 0)').num_interior_rings
@@ -959,7 +1102,7 @@ class GeometryTests < Minitest::Test
   def test_interior_rings
     array_tester(
       :interior_rings,
-      [ 'LINEARRING (11 11, 11 12, 12 12, 12 11, 11 11)' ],
+      ['LINEARRING (11 11, 11 12, 12 12, 12 11, 11 11)'],
       'POLYGON(
         (10 10, 10 14, 14 14, 14 10, 10 10),
         (11 11, 11 12, 12 12, 12 11, 11 11)
@@ -1024,12 +1167,12 @@ class GeometryTests < Minitest::Test
 
   def test_dimensions
     types = {
-      :dontcare => -3,
-      :non_empty => -2,
-      :empty => -1,
-      :point => 0,
-      :curve => 1,
-      :surface => 2
+      dontcare: -3,
+      non_empty: -2,
+      empty: -1,
+      point: 0,
+      curve: 1,
+      surface: 2
     }
 
     simple_tester(:dimensions, types[:point], 'POINT(0 0)')
@@ -1162,12 +1305,22 @@ class GeometryTests < Minitest::Test
     simple_tester(:distance, 2.0, geom, read('LINESTRING (3 0 , 10 0)'))
   end
 
+  def test_distance_indexed
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:distance_indexed)
+
+    geom_a = read('POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))')
+    geom_b = read('POLYGON ((20 30, 10 10, 13 14, 7 8, 20 30))')
+
+    assert_in_delta(9.219544457292887, geom_a.distance_indexed(geom_b), TOLERANCE)
+    assert_in_delta(9.219544457292887, geom_b.distance_indexed(geom_a), TOLERANCE)
+  end
+
   def test_hausdorff_distance
     skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:hausdorff_distance)
 
-    tester = lambda { |expected, g1, g2|
-      geom_1 = read(g1)
-      geom_2 = read(g2)
+    tester = lambda { |expected, g_1, g_2|
+      geom_1 = read(g_1)
+      geom_2 = read(g_2)
       assert_in_delta(expected, geom_1.hausdorff_distance(geom_2), TOLERANCE)
     }
 
@@ -1181,9 +1334,9 @@ class GeometryTests < Minitest::Test
   def test_hausdorff_distance_with_densify_fract
     skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:hausdorff_distance)
 
-    tester = lambda { |expected, g1, g2|
-      geom_1 = read(g1)
-      geom_2 = read(g2)
+    tester = lambda { |expected, g_1, g_2|
+      geom_1 = read(g_1)
+      geom_2 = read(g_2)
       assert_in_delta(expected, geom_1.hausdorff_distance(geom_2, 0.001), TOLERANCE)
     }
 
@@ -1197,17 +1350,18 @@ class GeometryTests < Minitest::Test
   def test_nearest_points
     skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:nearest_points)
 
-    tester = lambda { |expected, g1, g2|
-      geom_1 = read(g1)
-      geom_2 = read(g2)
+    tester = lambda { |expected, g_1, g_2|
+      geom_1 = read(g_1)
+      geom_2 = read(g_2)
 
       cs = geom_1.nearest_points(geom_2)
+      result = cs.to_s if cs
 
-      result = if cs
-        cs.to_s
+      if expected.nil?
+        assert_nil(result)
+      else
+        assert_equal(expected, result)
       end
-
-      assert_equal(expected, result)
     }
 
     tester[
@@ -1217,7 +1371,11 @@ class GeometryTests < Minitest::Test
     ]
 
     tester[
-      '5.0 5.0 NaN, 8.0 8.0 NaN',
+      if Geos::GEOS_VERSION > '3.8'
+        '5.0 5.0, 8.0 8.0'
+      else
+        '5.0 5.0 NaN, 8.0 8.0 NaN'
+      end,
       'POLYGON((1 1, 1 5, 5 5, 5 1, 1 1))',
       'POLYGON((8 8, 9 9, 9 10, 8 8))'
     ]
@@ -1249,11 +1407,46 @@ class GeometryTests < Minitest::Test
     assert_equal(2, polygonized.length)
     assert_equal(
       'POLYGON ((185 221, 88 275, 180 316, 292 281, 185 221))',
-      write(polygonized[0])
+      write(polygonized[0].snap_to_grid(0.1))
     )
     assert_equal(
       'POLYGON ((189 98, 83 187, 185 221, 325 168, 189 98))',
-      write(polygonized[1])
+      write(polygonized[1].snap_to_grid(0.1))
+    )
+  end
+
+  def test_polygonize_with_geometry_arguments
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:polygonize)
+
+    geom_a = read('LINESTRING (100 100, 100 300, 300 300, 300 100, 100 100)')
+    geom_b = read('LINESTRING (150 150, 150 250, 250 250, 250 150, 150 150)')
+
+    polygonized = geom_a.polygonize(geom_b)
+    assert_equal(2, polygonized.length)
+    assert_equal(
+      'POLYGON ((100 100, 100 300, 300 300, 300 100, 100 100), (150 150, 250 150, 250 250, 150 250, 150 150))',
+      write(polygonized[0].snap_to_grid(0.1))
+    )
+    assert_equal(
+      'POLYGON ((150 150, 150 250, 250 250, 250 150, 150 150))',
+      write(polygonized[1].snap_to_grid(0.1))
+    )
+  end
+
+  def test_polygonize_valid
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:polygonize_valid)
+
+    geom_a = read(
+      'GEOMETRYCOLLECTION(
+        LINESTRING (100 100, 100 300, 300 300, 300 100, 100 100),
+        LINESTRING (150 150, 150 250, 250 250, 250 150, 150 150)
+      )'
+    )
+
+    polygonized = geom_a.polygonize_valid
+    assert_equal(
+      'POLYGON ((100 100, 100 300, 300 300, 300 100, 100 100), (150 150, 250 150, 250 250, 150 250, 150 150))',
+      write(polygonized.snap_to_grid(0.1))
     )
   end
 
@@ -1288,7 +1481,8 @@ class GeometryTests < Minitest::Test
         LINESTRING(185 221, 292 281, 180 316),
         LINESTRING(189 98, 83 187, 185 221),
         LINESTRING(189 98, 325 168, 185 221)
-      )')
+      )'
+    )
 
     polygonized = geom_a.polygonize_full
 
@@ -1348,6 +1542,29 @@ class GeometryTests < Minitest::Test
     end
   end
 
+  def test_build_area
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:build_area)
+
+    geom = read('GEOMETRYCOLLECTION (LINESTRING(0 0, 0 1, 1 1), LINESTRING (1 1, 1 0, 0 0))')
+
+    assert_equal('POLYGON ((0 0, 0 1, 1 1, 1 0, 0 0))', write(geom.build_area))
+  end
+
+  def test_make_valid
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:make_valid)
+
+    geom = read('POLYGON((0 0, 1 1, 0 1, 1 0, 0 0))')
+
+    assert_equal(
+      if Geos::GEOS_VERSION > '3.9.0'
+        'MULTIPOLYGON (((1 0, 0 0, 0.5 0.5, 1 0)), ((1 1, 0.5 0.5, 0 1, 1 1)))'
+      else
+        'MULTIPOLYGON (((0 0, 0.5 0.5, 1 0, 0 0)), ((0.5 0.5, 0 1, 1 1, 0.5 0.5)))'
+      end,
+      write(geom.make_valid)
+    )
+  end
+
   def test_shared_paths
     skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:shared_paths)
 
@@ -1403,7 +1620,7 @@ class GeometryTests < Minitest::Test
     geom = read('LINESTRING(0 0, 10 10))')
     assert_kind_of(Enumerable, geom.each)
     assert_kind_of(Enumerable, geom.to_enum)
-    assert_equal(geom, geom.each {})
+    assert_equal(geom, geom.each(&EMPTY_BLOCK))
   end
 
   def test_normalize
@@ -1438,9 +1655,9 @@ class GeometryTests < Minitest::Test
     geom_a = read('POINT(1.0 1.0)')
     geom_b = read('POINT(2.0 2.0)')
 
-    assert(geom_a == geom_a, "Expected geoms to be equal using ==")
-    refute(geom_a == geom_b, "Expected geoms to not be equal using ==")
-    refute(geom_a == "test", "Expected geoms to not be equal using ==")
+    assert(geom_a == geom_a, 'Expected geoms to be equal using ==')
+    refute(geom_a == geom_b, 'Expected geoms to not be equal using ==')
+    refute(geom_a == 'test', 'Expected geoms to not be equal using ==')
   end
 
   def test_eql_exact
@@ -1517,7 +1734,7 @@ class GeometryTests < Minitest::Test
     assert_raises(Geos::MixedSRIDsError) do
       Geos.srid_copy_policy = :strict
       geom_c = geom.intersection(geom_b)
-      assert_equal(231231, geom_c.srid)
+      assert_equal(231_231, geom_c.srid)
     end
   ensure
     Geos.srid_copy_policy = :default
@@ -1573,23 +1790,25 @@ class GeometryTests < Minitest::Test
 
     # empty polygon
     tester['GEOMETRYCOLLECTION EMPTY', 'POLYGON EMPTY', 0]
-    tester['MULTILINESTRING EMPTY', 'POLYGON EMPTY', 0, { :only_edges => true }]
+    tester['MULTILINESTRING EMPTY', 'POLYGON EMPTY', 0, only_edges: true]
 
     # single point
     tester['GEOMETRYCOLLECTION EMPTY', 'POINT (0 0)', 0]
-    tester['MULTILINESTRING EMPTY', 'POINT (0 0)', 0, { :only_edges => true }]
+    tester['MULTILINESTRING EMPTY', 'POINT (0 0)', 0, only_edges: true]
 
     # three collinear points
     tester['GEOMETRYCOLLECTION EMPTY', 'MULTIPOINT(0 0, 5 0, 10 0)', 0]
-    tester['MULTILINESTRING ((5 0, 10 0), (0 0, 5 0))', 'MULTIPOINT(0 0, 5 0, 10 0)', 0, { :only_edges => true }]
+    tester['MULTILINESTRING ((5 0, 10 0), (0 0, 5 0))', 'MULTIPOINT(0 0, 5 0, 10 0)', 0, only_edges: true]
 
     # three points
     tester['GEOMETRYCOLLECTION (POLYGON ((0 0, 10 10, 5 0, 0 0)))', 'MULTIPOINT(0 0, 5 0, 10 10)', 0]
-    tester['MULTILINESTRING ((5 0, 10 10), (0 0, 10 10), (0 0, 5 0))', 'MULTIPOINT(0 0, 5 0, 10 10)', 0, { :only_edges => true }]
+    tester['MULTILINESTRING ((5 0, 10 10), (0 0, 10 10), (0 0, 5 0))', 'MULTIPOINT(0 0, 5 0, 10 10)', 0, only_edges: true]
 
     # polygon with a hole
     tester[
-      'GEOMETRYCOLLECTION (POLYGON ((8 2, 10 10, 8.5 1, 8 2)), POLYGON ((7 8, 10 10, 8 2, 7 8)), POLYGON ((3 8, 10 10, 7 8, 3 8)), POLYGON ((2 2, 8 2, 8.5 1, 2 2)), POLYGON ((2 2, 7 8, 8 2, 2 2)), POLYGON ((2 2, 3 8, 7 8, 2 2)), POLYGON ((0.5 9, 10 10, 3 8, 0.5 9)), POLYGON ((0.5 9, 3 8, 2 2, 0.5 9)), POLYGON ((0 0, 2 2, 8.5 1, 0 0)), POLYGON ((0 0, 0.5 9, 2 2, 0 0)))',
+      'GEOMETRYCOLLECTION (POLYGON ((8 2, 10 10, 8.5 1, 8 2)), POLYGON ((7 8, 10 10, 8 2, 7 8)), POLYGON ((3 8, 10 10, 7 8, 3 8)), ' \
+        'POLYGON ((2 2, 8 2, 8.5 1, 2 2)), POLYGON ((2 2, 7 8, 8 2, 2 2)), POLYGON ((2 2, 3 8, 7 8, 2 2)), POLYGON ((0.5 9, 10 10, 3 8, 0.5 9)), ' \
+        'POLYGON ((0.5 9, 3 8, 2 2, 0.5 9)), POLYGON ((0 0, 2 2, 8.5 1, 0 0)), POLYGON ((0 0, 0.5 9, 2 2, 0 0)))',
       'POLYGON((0 0, 8.5 1, 10 10, 0.5 9, 0 0),(2 2, 3 8, 7 8, 8 2, 2 2)))',
       0
     ]
@@ -1597,16 +1816,15 @@ class GeometryTests < Minitest::Test
     tester[
       'MULTILINESTRING ((8.5 1, 10 10), (8 2, 10 10), (8 2, 8.5 1), (7 8, 10 10), (7 8, 8 2), (3 8, 10 10), (3 8, 7 8), (2 2, 8.5 1), (2 2, 8 2), (2 2, 7 8), (2 2, 3 8), (0.5 9, 10 10), (0.5 9, 3 8), (0.5 9, 2 2), (0 0, 8.5 1), (0 0, 2 2), (0 0, 0.5 9))',
       'POLYGON((0 0, 8.5 1, 10 10, 0.5 9, 0 0),(2 2, 3 8, 7 8, 8 2, 2 2)))',
-      0, {
-        :only_edges => true
-      }
+      0,
+      only_edges: true
     ]
 
     # four points with a tolerance making one collapse
-    tester['MULTILINESTRING ((10 0, 10 10), (0 0, 10 10), (0 0, 10 0))', 'MULTIPOINT(0 0, 10 0, 10 10, 11 10)', 2.0, { :only_edges => true }]
+    tester['MULTILINESTRING ((10 0, 10 10), (0 0, 10 10), (0 0, 10 0))', 'MULTIPOINT(0 0, 10 0, 10 10, 11 10)', 2.0, only_edges: true]
 
     # tolerance as an option
-    tester['MULTILINESTRING ((10 0, 10 10), (0 0, 10 10), (0 0, 10 0))', 'MULTIPOINT(0 0, 10 0, 10 10, 11 10)', { :tolerance => 2.0, :only_edges => true }]
+    tester['MULTILINESTRING ((10 0, 10 10), (0 0, 10 10), (0 0, 10 0))', 'MULTIPOINT(0 0, 10 0, 10 10, 11 10)', tolerance: 2.0, only_edges: true]
   end
 
   def test_voronoi_diagram
@@ -1623,18 +1841,33 @@ class GeometryTests < Minitest::Test
 
     geom = 'MULTIPOINT(0 0, 100 0, 100 100, 0 100)'
 
-    tester['GEOMETRYCOLLECTION (POLYGON ((50 200, 200 200, 200 50, 50 50, 50 200)), POLYGON ((-100 50, -100 200, 50 200, 50 50, -100 50)), POLYGON ((50 -100, -100 -100, -100 50, 50 50, 50 -100)), POLYGON ((200 50, 200 -100, 50 -100, 50 50, 200 50)))', geom]
+    tester[
+      if Geos::GEOS_VERSION > '3.9.0'
+        'GEOMETRYCOLLECTION (POLYGON ((200 200, 200 50, 50 50, 50 200, 200 200)), POLYGON ((-100 200, 50 200, 50 50, -100 50, -100 200)), POLYGON ((-100 -100, -100 50, 50 50, 50 -100, -100 -100)), POLYGON ((200 -100, 50 -100, 50 50, 200 50, 200 -100)))'
+      else
+        'GEOMETRYCOLLECTION (POLYGON ((50 200, 200 200, 200 50, 50 50, 50 200)), POLYGON ((-100 50, -100 200, 50 200, 50 50, -100 50)), POLYGON ((50 -100, -100 -100, -100 50, 50 50, 50 -100)), POLYGON ((200 50, 200 -100, 50 -100, 50 50, 200 50)))'
+      end,
+      geom
+    ]
 
-    tester['MULTILINESTRING ((50 50, 50 200), (200 50, 50 50), (50 50, -100 50), (50 50, 50 -100))', geom, { :tolerance => 0, :only_edges => true }]
+    tester['MULTILINESTRING ((50 50, 50 200), (200 50, 50 50), (50 50, -100 50), (50 50, 50 -100))', geom, tolerance: 0, only_edges: true]
 
-    tester['MULTILINESTRING ((50 50, 50 1100), (1100 50, 50 50), (50 50, -1000 50), (50 50, 50 -1000))', geom, {
-      :only_edges => true,
-      :envelope => read(geom).buffer(1000)
-    }]
+    tester['MULTILINESTRING ((50 50, 50 1100), (1100 50, 50 50), (50 50, -1000 50), (50 50, 50 -1000))', geom,
+      only_edges: true,
+      envelope: read(geom).buffer(1000)
+    ]
 
     # Allows a tolerance for the first argument
     @writer.rounding_precision = 3
-    tester['GEOMETRYCOLLECTION (POLYGON ((290 252, 290 140, 185 140, 185 215, 188 235, 290 252)), POLYGON ((80 215, 80 340, 101 340, 188 235, 185 215, 80 215)), POLYGON ((185 140, 80 140, 80 215, 185 215, 185 140)), POLYGON ((101 340, 290 340, 290 252, 188 235, 101 340)))', "MULTIPOINT ((150 210), (210 270), (150 220), (220 210), (215 269))", 10]
+    tester[
+      if Geos::GEOS_VERSION > '3.9.0'
+        'GEOMETRYCOLLECTION (POLYGON ((290 140, 185 140, 185 215, 188 235, 290 252, 290 140)), POLYGON ((80 340, 101 340, 188 235, 185 215, 80 215, 80 340)), POLYGON ((80 140, 80 215, 185 215, 185 140, 80 140)), POLYGON ((290 340, 290 252, 188 235, 101 340, 290 340)))'
+      else
+        'GEOMETRYCOLLECTION (POLYGON ((290 252, 290 140, 185 140, 185 215, 188 235, 290 252)), POLYGON ((80 215, 80 340, 101 340, 188 235, 185 215, 80 215)), POLYGON ((185 140, 80 140, 80 215, 185 215, 185 140)), POLYGON ((101 340, 290 340, 290 252, 188 235, 101 340)))'
+      end,
+      'MULTIPOINT ((150 210), (210 270), (150 220), (220 210), (215 269))',
+      10
+    ]
   end
 
   def test_precision
@@ -1659,7 +1892,7 @@ class GeometryTests < Minitest::Test
     geom_with_precision = geom.with_precision(5.0)
     assert_equal('LINESTRING EMPTY', write(geom_with_precision))
 
-    geom_with_precision = geom.with_precision(5.0, :keep_collapsed => true)
+    geom_with_precision = geom.with_precision(5.0, keep_collapsed: true)
     assert_equal('LINESTRING (0 0, 0 0)', write(geom_with_precision))
   end
 
@@ -1707,11 +1940,37 @@ class GeometryTests < Minitest::Test
     tester['LINESTRING (200 200, 200 100)', 'POLYGON ((100 100, 300 100, 200 200, 100 100))']
     tester[
       'LINESTRING (-112.712119 33.575919, -112.712127 33.575885)',
-      '0106000000010000000103000000010000001a00000035d42824992d5cc01b834e081dca404073b9c150872d5cc03465a71fd4c940400ec00644882d5cc03b8a73d4d1c94040376dc669882d5cc0bf9cd9aed0c940401363997e892d5cc002f4fbfecdc94040ca4e3fa88b2d5cc0a487a1d5c9c940408f1ce90c8c2d5cc0698995d1c8c94040fab836548c2d5cc0bd175fb4c7c940409f1f46088f2d5cc0962023a0c2c940407b15191d902d5cc068041bd7bfc940400397c79a912d5cc0287d21e4bcc940403201bf46922d5cc065e3c116bbc940409d9d0c8e922d5cc0060fd3beb9c940400ef7915b932d5cc09012bbb6b7c940404fe61f7d932d5cc0e4a08499b6c94040fc71fbe5932d5cc0ea9106b7b5c94040eaec6470942d5cc0c2323674b3c94040601dc70f952d5cc043588d25acc94040aea06989952d5cc03ecf9f36aac94040307f85cc952d5cc0e5eb32fca7c94040dd0a6135962d5cc01b615111a7c9404048a7ae7c962d5cc00a2aaa7ea5c94040f4328ae5962d5cc05eb87361a4c94040c49448a2972d5cc04d81cccea2c940407c80eecb992d5cc06745d4449fc9404035d42824992d5cc01b834e081dca4040'
+      '0106000000010000000103000000010000001a00000035d42824992d5cc01b834e081dca404073b9c150872d5cc03465a71fd4c940400ec00644882d5cc03b8a' \
+        '73d4d1c94040376dc669882d5cc0bf9cd9aed0c940401363997e892d5cc002f4fbfecdc94040ca4e3fa88b2d5cc0a487a1d5c9c940408f1ce90c8c2d5cc06989' \
+        '95d1c8c94040fab836548c2d5cc0bd175fb4c7c940409f1f46088f2d5cc0962023a0c2c940407b15191d902d5cc068041bd7bfc940400397c79a912d5cc0287d' \
+        '21e4bcc940403201bf46922d5cc065e3c116bbc940409d9d0c8e922d5cc0060fd3beb9c940400ef7915b932d5cc09012bbb6b7c940404fe61f7d932d5cc0e4a0' \
+        '8499b6c94040fc71fbe5932d5cc0ea9106b7b5c94040eaec6470942d5cc0c2323674b3c94040601dc70f952d5cc043588d25acc94040aea06989952d5cc03ecf' \
+        '9f36aac94040307f85cc952d5cc0e5eb32fca7c94040dd0a6135962d5cc01b615111a7c9404048a7ae7c962d5cc00a2aaa7ea5c94040f4328ae5962d5cc05eb8' \
+        '7361a4c94040c49448a2972d5cc04d81cccea2c940407c80eecb992d5cc06745d4449fc9404035d42824992d5cc01b834e081dca4040'
     ]
     tester['LINESTRING EMPTY', 'POLYGON EMPTY']
   end
 
+  def test_maximum_inscribed_circle
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:maximum_inscribed_circle)
+
+    geom = read('POLYGON ((100 200, 200 200, 200 100, 100 100, 100 200))')
+    output = geom.maximum_inscribed_circle(0.001)
+    assert_equal('LINESTRING (150 150, 150 200)', write(output))
+  end
+
+  def test_largest_empty_circle
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:largest_empty_circle)
+
+    geom = read('MULTIPOINT ((100 100), (100 200), (200 200), (200 100))')
+    output = geom.largest_empty_circle(0.001)
+    assert_equal('LINESTRING (150 150, 100 100)', write(output))
+
+    geom = read('MULTIPOINT ((100 100), (100 200), (200 200), (200 100))')
+    output = geom.largest_empty_circle(0.001, boundary: read('MULTIPOINT ((100 100), (100 200), (200 200), (200 100))'))
+    assert_equal('LINESTRING (100 100, 100 100)', write(output))
+  end
+
   def test_minimum_width
     skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:minimum_width)
 
@@ -1723,4 +1982,90 @@ class GeometryTests < Minitest::Test
     output = geom.minimum_width
     assert_equal('LINESTRING (5 5, 0 10)', write(output))
   end
+
+  def test_dump_points
+    geom = read('GEOMETRYCOLLECTION(
+      MULTILINESTRING((0 0, 10 10, 20 20), (100 100, 200 200, 300 300)),
+
+      POINT(10 10),
+
+      POLYGON((0 0, 5 0, 5 5, 0 5, 0 0), (1 1, 4 1, 4 4, 1 4, 1 1))
+    )')
+
+    assert_equal([
+      [
+        [
+          Geos.create_point(0, 0),
+          Geos.create_point(10, 10),
+          Geos.create_point(20, 20)
+        ],
+
+        [
+          Geos.create_point(100, 100),
+          Geos.create_point(200, 200),
+          Geos.create_point(300, 300)
+        ]
+      ],
+
+      [
+        Geos.create_point(10, 10)
+      ],
+
+      [
+        [
+          Geos.create_point(0, 0),
+          Geos.create_point(5, 0),
+          Geos.create_point(5, 5),
+          Geos.create_point(0, 5),
+          Geos.create_point(0, 0)
+        ],
+
+        [
+          Geos.create_point(1, 1),
+          Geos.create_point(4, 1),
+          Geos.create_point(4, 4),
+          Geos.create_point(1, 4),
+          Geos.create_point(1, 1)
+        ]
+      ]
+    ], geom.dump_points)
+  end
+
+  def test_reverse
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:reverse)
+
+    simple_tester(:reverse, 'POINT (3 5)', 'POINT (3 5)')
+    simple_tester(:reverse, 'MULTIPOINT (100 100, 10 100, 30 100)', 'MULTIPOINT (100 100, 10 100, 30 100)')
+    simple_tester(:reverse, 'LINESTRING (200 200, 200 100)', 'LINESTRING (200 100, 200 200)')
+
+    if Geos::GEOS_VERSION >= '3.8.1'
+      simple_tester(:reverse, 'MULTILINESTRING ((3 3, 4 4), (1 1, 2 2))', 'MULTILINESTRING ((4 4, 3 3), (2 2, 1 1))')
+    else
+      simple_tester(:reverse, 'MULTILINESTRING ((1 1, 2 2), (3 3, 4 4))', 'MULTILINESTRING ((4 4, 3 3), (2 2, 1 1))')
+    end
+
+    simple_tester(:reverse, 'POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0), (1 1, 2 1, 2 2, 1 2, 1 1))', 'POLYGON ((0 0, 0 10, 10 10, 10 0, 0 0), (1 1, 1 2, 2 2, 2 1, 1 1))')
+    simple_tester(:reverse, 'MULTIPOLYGON (((0 0, 10 0, 10 10, 0 10, 0 0), (1 1, 2 1, 2 2, 1 2, 1 1)), ((100 100, 100 200, 200 200, 100 100)))', 'MULTIPOLYGON (((0 0, 0 10, 10 10, 10 0, 0 0), (1 1, 1 2, 2 2, 2 1, 1 1)), ((100 100, 200 200, 100 200, 100 100)))')
+    simple_tester(:reverse, 'GEOMETRYCOLLECTION (LINESTRING (1 1, 2 2), GEOMETRYCOLLECTION (LINESTRING (3 5, 2 9)))', 'GEOMETRYCOLLECTION (LINESTRING (2 2, 1 1), GEOMETRYCOLLECTION(LINESTRING (2 9, 3 5)))')
+    simple_tester(:reverse, 'POINT EMPTY', 'POINT EMPTY')
+    simple_tester(:reverse, 'LINESTRING EMPTY', 'LINESTRING EMPTY')
+    simple_tester(:reverse, 'LINEARRING EMPTY', 'LINEARRING EMPTY')
+    simple_tester(:reverse, 'POLYGON EMPTY', 'POLYGON EMPTY')
+    simple_tester(:reverse, 'MULTIPOINT EMPTY', 'MULTIPOINT EMPTY')
+    simple_tester(:reverse, 'MULTILINESTRING EMPTY', 'MULTILINESTRING EMPTY')
+    simple_tester(:reverse, 'MULTIPOLYGON EMPTY', 'MULTIPOLYGON EMPTY')
+    simple_tester(:reverse, 'GEOMETRYCOLLECTION EMPTY', 'GEOMETRYCOLLECTION EMPTY')
+  end
+
+  def test_frechet_distance
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:frechet_distance)
+
+    assert_in_delta(read('LINESTRING (0 0, 100 0)').frechet_distance(read('LINESTRING (0 0, 50 50, 100 0)')), 70.7106781186548, TOLERANCE)
+  end
+
+  def test_frechet_distance_with_densify
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:frechet_distance)
+
+    assert_in_delta(read('LINESTRING (0 0, 100 0)').frechet_distance(read('LINESTRING (0 0, 50 50, 100 0)'), 0.5), 50.0, TOLERANCE)
+  end
 end