ffi-geos 2.4.0 → 2.5.0

data/test/geometry/buffer_tests.rb

@@ -0,0 +1,116 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryBufferTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_buffer
+    simple_tester(
+      :buffer,
+      'POLYGON EMPTY',
+      'POINT(0 0)',
+      0
+    )
+
+    snapped_tester(
+      :buffer,
+      '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
+    snapped_tester(
+      :buffer,
+      'POLYGON ((10 0, 0 -10, -10 0, 0 10, 10 0))',
+      'POINT(0 0)',
+      10,
+      quad_segs: 1
+    )
+
+    # End cap styles
+    snapped_tester(
+      :buffer,
+      '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
+    )
+
+    snapped_tester(
+      :buffer,
+      'POLYGON ((100 10, 100 -10, 0 -10, 0 10, 100 10))',
+      'LINESTRING(0 0, 100 0)',
+      10,
+      quad_segs: 1, endcap: :flat
+    )
+
+    snapped_tester(
+      :buffer,
+      '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
+    )
+
+    # Join styles
+    snapped_tester(
+      :buffer,
+      '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
+    )
+
+    snapped_tester(
+      :buffer,
+      '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
+    )
+
+    snapped_tester(
+      :buffer,
+      '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
+    )
+
+    snapped_tester(
+      :buffer,
+      if Geos::GEOS_NICE_VERSION >= '031100'
+        'POLYGON ((90 10, 90 100, 93 107, 100 110, 107 107, 110 100, 110 -4, 104 -10, 0 -10, -7 -7, -10 0, -7 7, 0 10, 90 10))'
+      else
+        '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))'
+      end,
+      'LINESTRING(0 0, 100 0, 100 100)',
+      10,
+      quad_segs: 2, join: :mitre, mitre_limit: 1.0
+    )
+
+    # Single-sided buffering
+    snapped_tester(
+      :buffer,
+      'POLYGON ((100 0, 0 0, 0 10, 100 10, 100 0))',
+      'LINESTRING(0 0, 100 0)',
+      10,
+      single_sided: true
+    )
+
+    snapped_tester(
+      :buffer,
+      'POLYGON ((0 0, 100 0, 100 -10, 0 -10, 0 0))',
+      'LINESTRING(0 0, 100 0)',
+      -10,
+      single_sided: true
+    )
+  end
+end

data/test/geometry/build_area_tests.rb

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

data/test/geometry/centroid_tests.rb

@@ -0,0 +1,37 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryCentroidTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_centroid_and_center
+    %w{
+      centroid
+      center
+    }.each do |method|
+      simple_tester(
+        method,
+        'POINT (0 0)',
+        'POINT(0 0)'
+      )
+
+      simple_tester(
+        method,
+        'POINT (5 5)',
+        'LINESTRING(0 0, 10 10)'
+      )
+
+      snapped_tester(
+        method,
+        'POINT (5 4)',
+        'POLYGON((0 0, 0 10, 5 5, 10 10, 10 0, 0 0))'
+      )
+    end
+  end
+end

data/test/geometry/clip_by_rect_tests.rb

@@ -0,0 +1,56 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryClipByRectTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_clip_by_rect
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:clip_by_rect)
+
+    %w{
+      clip_by_rect
+      clip_by_rectangle
+    }.each do |method|
+      simple_tester(
+        method,
+        'POINT (0 0)',
+        'POINT (0 0)',
+        -1, -1, 1, 1
+      )
+
+      simple_tester(
+        method,
+        'GEOMETRYCOLLECTION EMPTY',
+        'POINT (0 0)',
+        0, 0, 2, 2
+      )
+
+      simple_tester(
+        method,
+        'LINESTRING (1 0, 2 0)',
+        'LINESTRING (0 0, 10 0)',
+        1, -1, 2, 1
+      )
+
+      simple_tester(
+        method,
+        'POLYGON ((1 1, 1 5, 5 5, 5 1, 1 1))',
+        'POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0))',
+        1, 1, 5, 5
+      )
+
+      simple_tester(
+        method,
+        'POLYGON ((0 0, 0 5, 5 5, 5 0, 0 0))',
+        'POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0))',
+        -1, -1, 5, 5
+      )
+    end
+  end
+end

data/test/geometry/clone_tests.rb

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

data/test/geometry/concave_hull_of_polygons_tests.rb

@@ -0,0 +1,28 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+describe 'Concave Hull Of Polygons Tests' do
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  it 'handles empty polygons' do
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:concave_hull_of_polygons)
+
+    geom = read('POLYGON EMPTY')
+
+    assert(geom.concave_hull_of_polygons(0.7).eql?(read('POLYGON EMPTY')))
+  end
+
+  it 'handles a multipolygon' do
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:concave_hull_of_polygons)
+
+    geom = read('MULTIPOLYGON(((0 0, 0 1, 1 1, 1 0, 0 0)))')
+
+    assert(geom.concave_hull_of_polygons(0.7).eql?(read('MULTIPOLYGON(((0 0, 0 1, 1 1, 1 0, 0 0)))')))
+  end
+end

data/test/geometry/concave_hull_tests.rb

@@ -0,0 +1,38 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+describe 'Geometry Concave Hull Tests' do
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def tester(expected, geom, **options)
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:concave_hull)
+
+    simple_tester(
+      :concave_hull,
+      expected,
+      geom,
+      **options
+    )
+  end
+
+  it 'returns a concave hull' do
+    tester(
+      'POLYGON ((30 70, 10 90, 60 72, 90 90, 90 60, 90 10, 60 30, 10 10, 40 40, 60 50, 47 66, 40 60, 30 70))',
+      'MULTIPOINT ((10 90), (10 10), (90 10), (90 90), (40 40), (60 30), (30 70), (40 60), (60 50), (60 72), (47 66), (90 60))'
+    )
+  end
+
+  it 'handles a concave hull by length' do
+    tester(
+      'POLYGON ((30 70, 10 90, 60 72, 90 90, 90 60, 90 10, 60 30, 10 10, 40 40, 30 70))',
+      'MULTIPOINT ((10 90), (10 10), (90 10), (90 90), (40 40), (60 30), (30 70), (40 60), (60 50), (60 72), (47 66), (90 60))',
+      length: 50
+    )
+  end
+end

data/test/geometry/convex_hull_tests.rb

@@ -0,0 +1,26 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryConvexHullTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_convex_hull
+    geom = read('POINT(0 0)')
+
+    assert_geom_eql_exact(read('POINT(0 0)'), geom.convex_hull)
+
+    geom = read('LINESTRING(0 0, 10 10)')
+
+    assert_geom_eql_exact(read('LINESTRING(0 0, 10 10)'), geom.convex_hull)
+
+    geom = read('POLYGON((0 0, 0 10, 5 5, 10 10, 10 0, 0 0))')
+
+    assert_geom_eql_exact(read('POLYGON((0 0, 0 10, 10 10, 10 0, 0 0))'), geom.convex_hull)
+  end
+end

data/test/geometry/coord_seq_tests.rb

@@ -0,0 +1,27 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryCoordSeqTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_coord_seq
+    tester = lambda { |expected, g|
+      geom = read(g)
+      cs = geom.coord_seq
+      expected.each_with_index do |c, i|
+        assert_equal(c[0], cs.get_x(i))
+        assert_equal(c[1], cs.get_y(i))
+      end
+    }
+
+    tester[[[0, 0]], 'POINT(0 0)']
+    tester[[[0, 0], [2, 3]], 'LINESTRING (0 0, 2 3)']
+    tester[[[0, 0], [0, 5], [5, 5], [5, 0], [0, 0]], 'LINEARRING(0 0, 0 5, 5 5, 5 0, 0 0)']
+  end
+end

data/test/geometry/delaunay_triangulation_tests.rb

@@ -0,0 +1,82 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeometryDelaunayTriangulationTests < Minitest::Test
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def test_delaunay_triangulation
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:delaunay_triangulation)
+
+    tester = lambda { |expected, geom, *args|
+      geom = read(geom)
+      geom_tri = geom.delaunay_triangulation(*args)
+      geom_tri.normalize!
+
+      assert_equal(expected, write(geom_tri))
+    }
+
+    writer.trim = true
+
+    # empty polygon
+    tester['GEOMETRYCOLLECTION EMPTY', 'POLYGON EMPTY', 0]
+    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]
+
+    # 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]
+
+    # 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]
+
+    # 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)))',
+      'POLYGON((0 0, 8.5 1, 10 10, 0.5 9, 0 0),(2 2, 3 8, 7 8, 8 2, 2 2))',
+      0
+    ]
+
+    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
+    ]
+
+    # 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]
+
+    # 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]
+  end
+
+  def test_constrained_delaunay_triangulation
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:constrained_delaunay_triangulation)
+
+    tester = lambda { |expected, geom|
+      geom = read(geom)
+      geom_tri = geom.constrained_delaunay_triangulation
+      geom_tri.normalize!
+
+      assert_equal(write(read(expected).normalize), write(geom_tri))
+    }
+
+    writer.trim = true
+
+    tester['GEOMETRYCOLLECTION EMPTY', 'POLYGON EMPTY']
+    tester['GEOMETRYCOLLECTION EMPTY', 'POINT(0 0)']
+    tester['GEOMETRYCOLLECTION (POLYGON ((10 10, 20 40, 90 10, 10 10)), POLYGON ((90 90, 20 40, 90 10, 90 90)))', 'POLYGON ((10 10, 20 40, 90 90, 90 10, 10 10))']
+  end
+end

data/test/geometry/densify_tests.rb

@@ -0,0 +1,95 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+describe 'Geometry Densify Tests' do
+  include TestHelper
+
+  def setup
+    super
+    writer.trim = true
+  end
+
+  def tester(expected, geom, tolerence)
+    skip unless ENV['FORCE_TESTS'] || Geos::Geometry.method_defined?(:densify)
+
+    simple_tester(
+      :densify,
+      expected,
+      geom,
+      tolerence
+    )
+  end
+
+  it 'can densify with tolerance greater than or equal to length of all edges' do
+    tester(
+      'POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0), (1 1, 1 2, 2 2, 2 1, 1 1))',
+      'POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0), (1 1, 1 2, 2 2, 2 1, 1 1))',
+      10.0
+    )
+  end
+
+  it 'can densify with a tolerance that evenly subdivides all outer and inner edges' do
+    tester(
+      'POLYGON ((0 0, 5 0, 10 0, 10 5, 10 10, 5 10, 0 10, 0 5, 0 0), (1 1, 1 4, 1 7, 4 7, 7 7, 7 4, 7 1, 4 1, 1 1))',
+      'POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0), (1 1, 1 7, 7 7, 7 1, 1 1))',
+      5.0
+    )
+  end
+
+  it 'can densify a LINESTRING' do
+    tester(
+      'LINESTRING (0 0, 0 3, 0 6)',
+      'LINESTRING (0 0, 0 6 )',
+      3.0
+    )
+  end
+
+  it 'can ensure that tolerance results in the right number of subdivisions' do
+    tester(
+      'LINESTRING (0 0, 0 2, 0 4, 0 6)',
+      'LINESTRING (0 0, 0 6 )',
+      2.9999999
+    )
+  end
+
+  it 'can densify a LINEARRING' do
+    tester(
+      'LINEARRING (0 0, 0 3, 0 6, 3 6, 6 6, 4 4, 2 2, 0 0)',
+      'LINEARRING (0 0, 0 6, 6 6, 0 0)',
+      3.0
+    )
+  end
+
+  it 'can densify a POINT' do
+    tester(
+      'POINT (0 0)',
+      'POINT (0 0)',
+      3.0
+    )
+  end
+
+  it 'can densify a MULTIPOINT' do
+    tester(
+      'MULTIPOINT ((0 0), (10 10))',
+      'MULTIPOINT ((0 0), (10 10))',
+      3.0
+    )
+  end
+
+  it 'can densify an empty polygon' do
+    tester(
+      'POLYGON EMPTY',
+      'POLYGON EMPTY',
+      3.0
+    )
+  end
+
+  it 'densify with an invalid tolerances should fail' do
+    geom = read('POLYGON ((0 0, 10 0, 10 10, 0 10, 0 0))')
+
+    assert_raises(Geos::GEOSException, 'IllegalArgumentException: Tolerance must be positive') do
+      geom.densify(-1.0)
+    end
+  end
+end

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