ffi-geos 2.3.1 → 2.5.0

data/test/geometry/difference_tests.rb

@@ -0,0 +1,108 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryDifferenceTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_difference
+    comparison_tester(
+      :difference,
+      EMPTY_GEOMETRY,
+      'POINT(0 0)',
+      'POINT(0 0)'
+    )
+
+    comparison_tester(
+      :difference,
+      'POINT (0 0)',
+      'POINT(0 0)',
+      'POINT(1 0)'
+    )
+
+    comparison_tester(
+      :difference,
+      'LINESTRING (0 0, 10 0)',
+      'LINESTRING(0 0, 10 0)',
+      'POINT(5 0)'
+    )
+
+    comparison_tester(
+      :difference,
+      EMPTY_GEOMETRY,
+      'POINT(5 0)',
+      'LINESTRING(0 0, 10 0)'
+    )
+
+    comparison_tester(
+      :difference,
+      'POINT (5 0)',
+      'POINT(5 0)',
+      'LINESTRING(0 1, 10 1)'
+    )
+
+    comparison_tester(
+      :difference,
+      'MULTILINESTRING ((0 0, 5 0), (5 0, 10 0))',
+      'LINESTRING(0 0, 10 0)',
+      'LINESTRING(5 -10, 5 10)'
+    )
+
+    comparison_tester(
+      :difference,
+      'LINESTRING (0 0, 5 0)',
+      'LINESTRING(0 0, 10 0)',
+      'LINESTRING(5 0, 20 0)'
+    )
+
+    comparison_tester(
+      :difference,
+      if Geos::GEOS_NICE_VERSION > '030900'
+        '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,
+      if Geos::GEOS_NICE_VERSION > '030900'
+        '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,
+      if Geos::GEOS_NICE_VERSION > '030900'
+        '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
+end

data/test/geometry/dimensions_tests.rb

@@ -0,0 +1,46 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryDimensionsTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_dimensions
+    types = {
+      dontcare: -3,
+      non_empty: -2,
+      empty: -1,
+      point: 0,
+      curve: 1,
+      surface: 2
+    }
+
+    simple_tester(:dimensions, types[:point], 'POINT(0 0)')
+    simple_tester(:dimensions, types[:point], 'MULTIPOINT (0 1, 2 3)')
+    simple_tester(:dimensions, types[:curve], 'LINESTRING (0 0, 2 3)')
+    simple_tester(:dimensions, types[:curve], 'MULTILINESTRING ((0 1, 2 3), (10 10, 3 4))')
+    simple_tester(:dimensions, types[:surface], 'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))')
+    simple_tester(:dimensions, types[:surface], '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))
+    )')
+    simple_tester(:dimensions, types[:surface], 'GEOMETRYCOLLECTION (
+      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))
+      ),
+      POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0)),
+      MULTILINESTRING ((0 0, 2 3), (10 10, 3 4)),
+      LINESTRING (0 0, 2 3),
+      MULTIPOINT (0 0, 2 3),
+      POINT (9 0)
+    )')
+  end
+end

data/test/geometry/distance_tests.rb

@@ -0,0 +1,29 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryDistanceTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_distance
+    geom = 'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))'
+    simple_tester(:distance, 0.0, geom, read('POINT(0.5 0.5)'))
+    simple_tester(:distance, 1.0, geom, read('POINT (-1 0)'))
+    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
+end

data/test/geometry/dump_points_tests.rb

@@ -0,0 +1,60 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryDumpPointsTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  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
+end

data/test/geometry/dup_tests.rb

@@ -0,0 +1,29 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryDupTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_dup
+    geom_a = read('POINT(0 0)')
+    geom_b = geom_a.dup
+
+    assert_equal(geom_a, geom_b)
+  end
+
+  def test_dup_srid
+    srid = 4326
+    geom_a = read('POINT(0 0)')
+    geom_a.srid = srid
+    geom_b = geom_a.dup
+
+    assert_equal(geom_a, geom_b)
+    assert_equal(srid, geom_b.srid)
+  end
+end

data/test/geometry/empty_tests.rb

@@ -0,0 +1,23 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryEmptyTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_empty
+    refute_geom_empty(read('POINT(0 0)'))
+    assert_geom_empty(read('POINT EMPTY'))
+    refute_geom_empty(read('LINESTRING(0 0, 10 0)'))
+    assert_geom_empty(read('LINESTRING EMPTY'))
+    refute_geom_empty(read('POLYGON((0 0, 10 0, 10 10, 0 0))'))
+    assert_geom_empty(read('POLYGON EMPTY'))
+    refute_geom_empty(read('GEOMETRYCOLLECTION(POINT(0 0))'))
+    assert_geom_empty(read('GEOMETRYCOLLECTION EMPTY'))
+  end
+end

data/test/geometry/envelope_tests.rb

@@ -0,0 +1,26 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryEnvelopeTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_envelope
+    simple_tester(
+      :envelope,
+      'POINT (0 0)',
+      'POINT(0 0)'
+    )
+
+    simple_tester(
+      :envelope,
+      'POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0))',
+      'LINESTRING(0 0, 10 10)'
+    )
+  end
+end

data/test/geometry/equal_identical_tests.rb

@@ -0,0 +1,78 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryEqualIdenticalTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def tester(expected, geom_a, geom_b)
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:eql_identical?)
+
+    assert_equal(expected, read(geom_a).eql_identical?(read(geom_b)))
+  end
+
+  def test_different_types
+    tester(false, 'POINT EMPTY', 'LINESTRING EMPTY')
+  end
+
+  def test_different_dimensions_empty
+    tester(false, 'POINT EMPTY', 'POINT Z EMPTY')
+  end
+
+  def test_different_dimensions
+    tester(false, 'POINT Z (1 2 3)', 'POINT M (1 2 3)')
+  end
+
+  def test_different_dimensions_zm_versus_m
+    tester(false, 'POINT ZM (1 2 3 4)', 'POINT Z (1 2 3)')
+  end
+
+  def test_different_structure
+    tester(false, 'LINESTRING (1 1, 2 2)', 'MULTILINESTRING ((1 1, 2 2))')
+  end
+
+  def test_different_types_geometry_collection
+    tester(false, 'GEOMETRYCOLLECTION (LINESTRING (1 1, 2 2))', 'MULTILINESTRING ((1 1, 2 2))')
+  end
+
+  def test_non_finite_values
+    tester(true, 'POINT(inf inf)', 'POINT(inf inf)')
+  end
+
+  def test_equal_lines
+    tester(true, 'LINESTRING M (1 1 0, 2 2 1)', 'LINESTRING M (1 1 0, 2 2 1)')
+  end
+
+  def test_different_lines
+    tester(false, 'LINESTRING M (1 1 0, 2 2 1)', 'LINESTRING M (1 1 1, 2 2 1)')
+  end
+
+  def test_equal_polygons
+    tester(true, 'POLYGON ((0 0, 1 0, 1 1, 0 0))', 'POLYGON ((0 0, 1 0, 1 1, 0 0))')
+  end
+
+  def test_different_polygons
+    tester(false, 'POLYGON ((0 0, 1 0, 1 1, 0 0))', 'POLYGON ((1 0, 1 1, 0 0, 1 0))')
+  end
+
+  def test_equal_polygons_with_holes
+    tester(false, 'POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0), (1 1, 2 1, 2 2, 1 1))', 'POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0), (1 1, 2 1, 2 2, 1 1), (3 3, 4 3, 4 4, 3 3))')
+  end
+
+  def test_identical_collections
+    tester(true, 'MULTILINESTRING ((1 1, 2 2), (2 2, 3 3))', 'MULTILINESTRING ((1 1, 2 2), (2 2, 3 3))')
+  end
+
+  def test_different_collection_structure
+    tester(false, 'MULTILINESTRING ((1 1, 2 2), (2 2, 3 3))', 'MULTILINESTRING ((2 2, 3 3), (1 1, 2 2))')
+  end
+
+  def test_negative_zeo_and_positive_zero
+    tester(true, 'POINT(1 0)', 'POINT(1 -0)')
+  end
+end

data/test/geometry/equal_tests.rb

@@ -0,0 +1,62 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryEqualTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_eql
+    geom_a = read('POINT(1.0 1.0)')
+    geom_b = read('POINT(2.0 2.0)')
+
+    %w{ eql? equals? }.each do |method|
+      assert(geom_a.send(method, geom_a), "Expected geoms to be equal using #{method}")
+      refute(geom_a.send(method, geom_b), "Expected geoms to not be equal using #{method}")
+    end
+  end
+
+  def test_equals_operator
+    geom_a = read('POINT(1.0 1.0)')
+    geom_b = read('POINT(2.0 2.0)')
+
+    refute_equal(geom_a, geom_b, 'Expected geoms to not be equal using ==')
+    refute_equal(geom_a, 'test', 'Expected geoms to not be equal using ==')
+  end
+
+  def test_eql_exact
+    geom_a = read('POINT(1.0 1.0)')
+    geom_b = read('POINT(2.0 2.0)')
+
+    %w{ eql_exact? equals_exact? exactly_equals? }.each do |method|
+      refute(geom_a.send(method, geom_b, 0.001), "Expected geoms to not be equal using #{method}")
+    end
+  end
+
+  def test_eql_almost_default
+    geom = read('POINT (1 1)')
+    geom_a = read('POINT (1.0000001 1.0000001)')
+    geom_b = read('POINT (1.000001 1.000001)')
+
+    %w{ eql_almost? equals_almost? almost_equals? }.each do |method|
+      assert(geom.send(method, geom_a), "Expected geoms to be equal using #{method}")
+      refute(geom.send(method, geom_b), "Expected geoms to not be equal using #{method}")
+    end
+  end
+
+  def test_eql_almost
+    geom_a = read('POINT(1.0 1.0)')
+    geom_b = read('POINT(1.1 1.1)')
+
+    refute_equal(geom_a, geom_b)
+
+    %w{ eql_almost? equals_almost? almost_equals? }.each do |method|
+      assert(geom_a.send(method, geom_b, 0), "Expected geoms to be equal using #{method}")
+      refute(geom_a.send(method, geom_b, 1), "Expected geoms to not be equal using #{method}")
+    end
+  end
+end

data/test/geometry/exterior_ring_tests.rb

@@ -0,0 +1,27 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryInteriorRingTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_exterior_ring
+    simple_tester(
+      :exterior_ring,
+      'LINEARRING (10 10, 10 14, 14 14, 14 10, 10 10)',
+      'POLYGON (
+        (10 10, 10 14, 14 14, 14 10, 10 10),
+        (11 11, 11 12, 12 12, 12 11, 11 11)
+      )'
+    )
+
+    assert_raises(NoMethodError) do
+      read('POINT (0 0)').exterior_ring
+    end
+  end
+end

data/test/geometry/extract_unique_points_tests.rb

@@ -0,0 +1,41 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryExtractUniquePointsTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_extract_unique_points
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:extract_unique_points)
+
+    geom = read('GEOMETRYCOLLECTION (
+      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))
+      ),
+      POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0)),
+      MULTILINESTRING ((0 0, 2 3), (10 10, 3 4)),
+      LINESTRING (0 0, 2 3),
+      MULTIPOINT (0 0, 2 3),
+      POINT (9 0),
+      POINT (1 0),
+      LINESTRING EMPTY
+    )')
+
+    simple_tester(
+      :extract_unique_points,
+      if Geos::GEOS_NICE_VERSION >= '031200'
+        'MULTIPOINT ((0 0), (1 0), (1 1), (0 1), (10 10), (10 14), (14 14), (14 10), (11 11), (11 12), (12 12), (12 11), (2 3), (3 4), (9 0))'
+      else
+        'MULTIPOINT (0 0, 1 0, 1 1, 0 1, 10 10, 10 14, 14 14, 14 10, 11 11, 11 12, 12 12, 12 11, 2 3, 3 4, 9 0)'
+      end,
+      geom.extract_unique_points
+    )
+  end
+end

data/test/geometry/frecet_distance_tests.rb

@@ -0,0 +1,24 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryFrechetDistanceTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  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

data/test/geometry/get_geometry_n_tests.rb

@@ -0,0 +1,21 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryGetGeometryNTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  # get_geometry_n is segfaulting in the binary GEOS build
+  def test_get_geometry_n
+    skip unless defined?(Geos::FFIGeos)
+
+    simple_tester(:get_geometry_n, 'POINT (0 1)', 'MULTIPOINT (0 1, 2 3)', 0)
+    simple_tester(:get_geometry_n, 'POINT (2 3)', 'MULTIPOINT (0 1, 2 3)', 1)
+    simple_tester(:get_geometry_n, nil, 'MULTIPOINT (0 1, 2 3)', 2)
+  end
+end

data/test/geometry/hausdorff_distance_tests.rb

@@ -0,0 +1,46 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryHausdorffDistanceTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_hausdorff_distance
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:hausdorff_distance)
+
+    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)
+    }
+
+    geom_a = 'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))'
+
+    tester[10.0498756211209, geom_a, 'POINT(0 10)']
+    tester[2.23606797749979, geom_a, 'POINT(-1 0)']
+    tester[9.0, geom_a, 'LINESTRING (3 0 , 10 0)']
+  end
+
+  def test_hausdorff_distance_with_densify_fract
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:hausdorff_distance)
+
+    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)
+    }
+
+    geom_a = 'POLYGON ((0 0, 1 0, 1 1, 0 1, 0 0))'
+
+    tester[10.0498756211209, geom_a, 'POINT(0 10)']
+    tester[2.23606797749979, geom_a, 'POINT(-1 0)']
+    tester[9.0, geom_a, 'LINESTRING (3 0 , 10 0)']
+  end
+end

data/test/geometry/hilbert_code_tests.rb

@@ -0,0 +1,45 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+describe '#hilbert_code' do
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  it 'handles various points on the Hilbert curve' do
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:hilbert_code)
+
+    geom_a = read('POINT (0 0)')
+    geom_b = read('POINT (1 1)')
+    extent = read('POLYGON ((0 0, 0 1, 1 1, 1 0, 0 0))')
+
+    assert_equal(0, geom_a.hilbert_code(extent, 1))
+    assert_equal(0, geom_a.hilbert_code(extent, 16))
+    assert_equal(10, geom_b.hilbert_code(extent, 2))
+    assert_equal(43690, geom_b.hilbert_code(extent, 8))
+    assert_equal(2863311530, geom_b.hilbert_code(extent, 16))
+  end
+
+  it 'can calculate the midpoint' do
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:hilbert_code)
+
+    extent = read('POLYGON ((0 0, 0 1, 1 1, 1 0, 0 0))')
+
+    assert_equal(2, extent.hilbert_code(extent, 2))
+  end
+
+  it 'handles out of bounds' do
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:hilbert_code)
+
+    geom = read('POINT (0 0)')
+    extent = read('POLYGON ((0 0, 0 1, 1 1, 1 0, 0 0))')
+
+    assert_raises(Geos::GEOSException) do
+      geom.hilbert_code(extent, 17)
+    end
+  end
+end