ffi-geos 1.2.1 → 1.2.2

data/test/coordinate_sequence_tests.rb

@@ -1,4 +1,3 @@
-# encoding: UTF-8
 # frozen_string_literal: true
 
 require 'test_helper'
@@ -131,8 +130,8 @@ class CoordinateSequenceTests < Minitest::Test
   end
 
   def test_to_point
-    cs = Geos::CoordinateSequence.new([5,7])
-    assert_equal('POINT (5 7)', write(cs.to_point, :trim => true))
+    cs = Geos::CoordinateSequence.new([ 5, 7 ])
+    assert_equal('POINT (5 7)', write(cs.to_point, trim: true))
   end
 
   def test_to_to_linear_ring
@@ -144,21 +143,21 @@ class CoordinateSequenceTests < Minitest::Test
       [ 0, 0 ]
     ])
 
-    assert_equal('LINEARRING (0 0, 0 5, 5 5, 5 0, 0 0)', write(cs.to_linear_ring, :trim => true))
+    assert_equal('LINEARRING (0 0, 0 5, 5 5, 5 0, 0 0)', write(cs.to_linear_ring, trim: true))
   end
 
   def test_empty
     cs = Geos::CoordinateSequence.new
     assert_geom_empty(cs)
 
-    cs = Geos::CoordinateSequence.new([4,1])
+    cs = Geos::CoordinateSequence.new([ 4, 1 ])
     refute_geom_empty(cs)
   end
 
   def test_to_empty_linear_ring
     cs = Geos::CoordinateSequence.new
 
-    assert_equal('LINEARRING EMPTY', write(cs.to_linear_ring, :trim => true))
+    assert_equal('LINEARRING EMPTY', write(cs.to_linear_ring, trim: true))
   end
 
   def test_to_line_string
@@ -169,13 +168,13 @@ class CoordinateSequenceTests < Minitest::Test
       [ 5, 0 ]
     ])
 
-    assert_equal('LINESTRING (0 0, 0 5, 5 5, 5 0)', write(cs.to_line_string, :trim => true))
+    assert_equal('LINESTRING (0 0, 0 5, 5 5, 5 0)', write(cs.to_line_string, trim: true))
   end
 
   def test_to_empty_line_string
     cs = Geos::CoordinateSequence.new
 
-    assert_equal('LINESTRING EMPTY', write(cs.to_line_string, :trim => true))
+    assert_equal('LINESTRING EMPTY', write(cs.to_line_string, trim: true))
   end
 
   def test_to_polygon
@@ -187,21 +186,23 @@ class CoordinateSequenceTests < Minitest::Test
       [ 0, 0 ]
     ])
 
-    assert_equal('POLYGON ((0 0, 0 5, 5 5, 5 0, 0 0))', write(cs.to_polygon, :trim => true))
+    assert_equal('POLYGON ((0 0, 0 5, 5 5, 5 0, 0 0))', write(cs.to_polygon, trim: true))
   end
 
   def test_to_empty_polygon
     cs = Geos::CoordinateSequence.new
 
-    assert_equal('POLYGON EMPTY', write(cs.to_polygon, :trim => true))
+    assert_equal('POLYGON EMPTY', write(cs.to_polygon, trim: true))
   end
 
-  def test_to_s
+  def test_to_s_2d
     cs = Geos::CoordinateSequence.new([[1, 2], [10, 11]])
-    assert_equal("1.0 2.0, 10.0 11.0", cs.to_s)
+    assert_equal('1.0 2.0, 10.0 11.0', cs.to_s)
+  end
 
+  def test_to_s_3d
     cs = Geos::CoordinateSequence.new([[1, 2, 3], [10, 11, 12]])
-    assert_equal("1.0 2.0 3.0, 10.0 11.0 12.0", cs.to_s)
+    assert_equal('1.0 2.0 3.0, 10.0 11.0 12.0', cs.to_s)
   end
 
   def test_get_by_proxy
@@ -265,7 +266,7 @@ class CoordinateSequenceTests < Minitest::Test
   end
 
   def test_options_hash
-    cs = Geos::CoordinateSequence.new(:size => 10, :dimensions => 2)
+    cs = Geos::CoordinateSequence.new(size: 10, dimensions: 2)
 
     assert_equal(10, cs.size)
     assert_equal(2, cs.dimensions)
@@ -332,4 +333,252 @@ class CoordinateSequenceTests < Minitest::Test
     assert_geom_has_z(read('POINT (0 0 0)').coord_seq)
     refute_geom_has_z(read('POINT (0 0)').coord_seq)
   end
+
+  def test_x_max
+    cs = Geos::CoordinateSequence.new([ -10, -15 ], [ 0, 5 ], [ 10, 20 ])
+    assert_equal(10, cs.x_max)
+  end
+
+  def test_x_min
+    cs = Geos::CoordinateSequence.new([ -10, -15 ], [ 0, 5 ], [ 10, 20 ])
+    assert_equal(-10, cs.x_min)
+  end
+
+  def test_y_max
+    cs = Geos::CoordinateSequence.new([ -10, -15 ], [ 0, 5 ], [ 10, 20 ])
+    assert_equal(20, cs.y_max)
+  end
+
+  def test_y_min
+    cs = Geos::CoordinateSequence.new([ -10, -15 ], [ 0, 5 ], [ 10, 20 ])
+    assert_equal(-15, cs.y_min)
+  end
+
+  def test_z_max
+    cs = Geos::CoordinateSequence.new([ -10, -15 ], [ 0, 5 ], [ 10, 20 ])
+    assert(cs.z_max.nan?, " Expected NaN")
+
+    cs = Geos::CoordinateSequence.new([ -10, -15, -20 ], [ 0, 5, 10 ], [ 10, 20, 30 ])
+    assert_equal(30, cs.z_max)
+  end
+
+  def test_z_min
+    cs = Geos::CoordinateSequence.new([ -10, -15 ], [ 0, 5 ], [ 10, 20 ])
+    assert(cs.z_min.nan?, " Expected NaN")
+
+    cs = Geos::CoordinateSequence.new([ -10, -15, -20 ], [ 0, 5, 10 ], [ 10, 20, 30 ])
+    assert_equal(-20, cs.z_min)
+  end
+
+  def test_snap_to_grid_with_size
+    expected = [
+      [[-10.0, -15.0], [0.0, 5.0], [10.0, 20.0]],
+      [[-10.1, -15.1], [0.1, 5.1], [10.1, 20.1]],
+      [[-10.12, -15.12], [0.12, 5.12], [10.12, 20.12]],
+      [[-10.123, -15.123], [0.123, 5.123], [10.123, 20.123]],
+      [[-10.1235, -15.1235], [0.1235, 5.1235], [10.1235, 20.1235]],
+      [[-10.12346, -15.12346], [0.12346, 5.12346], [10.12346, 20.12346]],
+      [[-10.123457, -15.123457], [0.123457, 5.123457], [10.123457, 20.123457]],
+      [[-10.1234568, -15.1234568], [0.1234568, 5.1234568], [10.1234568, 20.1234568]],
+      [[-10.12345679, -15.12345679], [0.12345679, 5.12345679], [10.12345679, 20.12345679]]
+    ]
+
+    coordinates = [
+      [ -10.123456789, -15.123456789 ],
+      [ 0.123456789, 5.123456789 ],
+      [ 10.123456789, 20.123456789 ]
+    ]
+
+    9.times do |i|
+      cs = Geos::CoordinateSequence.new(*coordinates)
+      cs.snap_to_grid!(10 ** -i)
+
+      # XXX - Ruby 1.8.7 sometimes sees the the float values as differing
+      # slightly, but not enough that it would matter for these tests.
+      # Test equality on the inspect Strings instead of the float values.
+      assert_equal(expected[i].inspect, cs.to_a.inspect)
+
+      cs = Geos::CoordinateSequence.new(*coordinates)
+      snapped = cs.snap_to_grid(10 ** -i)
+      assert_equal(coordinates, cs.to_a)
+      assert_equal(expected[i].inspect, snapped.to_a.inspect)
+    end
+  end
+
+  def test_snap_to_grid_with_hash
+    cs = Geos::CoordinateSequence.new(
+      [ 10, 10 ],
+      [ 20, 20 ],
+      [ 30, 30 ]
+    )
+    cs.snap_to_grid!(:size_x => 1, :size_y => 1, :offset_x => 12.5, :offset_y => 12.5)
+
+    assert_equal([
+      [ 9.5, 9.5 ],
+      [ 20.5, 20.5 ],
+      [ 30.5, 30.5 ]
+    ], cs.to_a)
+  end
+
+  def test_snap_to_grid_with_geometry_origin
+    cs = Geos::CoordinateSequence.new(
+      [ 10, 10 ],
+      [ 20, 20 ],
+      [ 30, 30 ]
+    )
+    cs.snap_to_grid!(:size => 1, :offset => read('LINESTRING (0 0, 25 25)'))
+
+    assert_equal([
+      [ 9.5, 9.5 ],
+      [ 20.5, 20.5 ],
+      [ 30.5, 30.5 ]
+    ], cs.to_a)
+  end
+
+  def test_snap_to_grid_with_z
+    cs = Geos::CoordinateSequence.new(
+      [ 10, 10, 10 ],
+      [ 20, 20, 20 ],
+      [ 30, 30, 30 ]
+    )
+    cs.snap_to_grid!(
+      :size_x => 1,
+      :size_y => 1,
+      :size_z => 1,
+
+      :offset_x => 12.5,
+      :offset_y => 12.5,
+      :offset_z => 12.5
+    )
+
+    assert_equal([
+      [ 9.5, 9.5, 9.5 ],
+      [ 20.5, 20.5, 20.5 ],
+      [ 30.5, 30.5, 30.5 ]
+    ], cs.to_a)
+  end
+
+  def test_snap_to_grid_remove_duplicate_points
+    coords = [
+      [-10.0, 0.0],
+      [-10.0, 5.0], [-10.0, 5.0],
+      [-10.0, 6.0], [-10.0, 6.0], [-10.0, 6.0],
+      [-10.0, 7.0], [-10.0, 7.0], [-10.0, 7.0],
+      [-10.0, 8.0], [-10.0, 8.0],
+      [-9.0, 8.0], [-9.0, 9.0],
+      [-10.0, 0.0]
+    ]
+
+    expected = [
+      [-10.0, 0.0],
+      [-10.0, 5.0],
+      [-10.0, 6.0],
+      [-10.0, 7.0],
+      [-10.0, 8.0],
+      [-9.0, 8.0],
+      [-9.0, 9.0],
+      [-10.0, 0.0]
+    ]
+
+    cs = Geos::CoordinateSequence.new(coords)
+    cs.snap_to_grid!
+
+    assert_equal(expected, cs.to_a)
+
+    cs = Geos::CoordinateSequence.new(coords)
+    cs2 = cs.snap_to_grid
+
+    assert_equal(coords, cs.to_a)
+    assert_equal(expected, cs2.to_a)
+  end
+
+  undef :affine_tester
+  def affine_tester(method, expected, coords, *args)
+    cs = Geos::CoordinateSequence.new(coords)
+    cs.send("#{method}!", *args)
+
+    expected.length.times do |i|
+      assert_in_delta(expected[i], cs.get_ordinate(0, i), TOLERANCE)
+    end
+
+    cs = Geos::CoordinateSequence.new(coords)
+    cs2 = cs.send(method, *args)
+
+    expected.length.times do |i|
+      assert_in_delta(coords[i], cs.get_ordinate(0, i), TOLERANCE)
+      assert_in_delta(expected[i], cs2.get_ordinate(0, i), TOLERANCE)
+    end
+  end
+
+  def test_rotate
+    affine_tester(:rotate, [ 29.0, 11.0 ], [ 1, 1 ], Math::PI / 2, [ 10.0, 20.0 ])
+    affine_tester(:rotate, [ -2.0, 0.0 ], [ 1, 1 ], -Math::PI / 2, [ -1.0, 2.0 ])
+    affine_tester(:rotate, [ 19.0, 1.0 ], [ 1, 1 ], Math::PI / 2, read('POINT(10 10)'))
+    affine_tester(:rotate, [ -0.5, 0.5 ], [ 1, 1 ], Math::PI / 2, read('LINESTRING(0 0, 1 0)'))
+  end
+
+  def test_rotate_x
+    affine_tester(:rotate_x, [ 1, -1, -1 ], [ 1, 1, 1 ], Math::PI)
+    affine_tester(:rotate_x, [ 1, -1, 1 ], [ 1, 1, 1 ], Math::PI / 2)
+    affine_tester(:rotate_x, [ 1, 1, -1 ], [ 1, 1, 1 ], Math::PI + Math::PI / 2)
+    affine_tester(:rotate_x, [ 1, 1, 1 ], [ 1, 1, 1 ], Math::PI * 2)
+  end
+
+  def test_rotate_y
+    affine_tester(:rotate_y, [ -1, 1, -1 ], [ 1, 1, 1 ], Math::PI)
+    affine_tester(:rotate_y, [ 1, 1, -1 ], [ 1, 1, 1 ], Math::PI / 2)
+    affine_tester(:rotate_y, [ -1, 1, 1 ], [ 1, 1, 1 ], Math::PI + Math::PI / 2)
+    affine_tester(:rotate_y, [ 1, 1, 1 ], [ 1, 1, 1 ], Math::PI * 2)
+  end
+
+  def test_rotate_z
+    affine_tester(:rotate_z, [ -1, -1 ], [ 1, 1 ], Math::PI)
+    affine_tester(:rotate_z, [ -1, 1 ], [ 1, 1 ], Math::PI / 2)
+    affine_tester(:rotate_z, [ 1, -1 ], [ 1, 1 ], Math::PI + Math::PI / 2)
+    affine_tester(:rotate_z, [ 1, 1 ], [ 1, 1 ], Math::PI * 2)
+  end
+
+  def test_scale
+    affine_tester(:scale, [ 5, 5 ], [ 1, 1 ], 5, 5)
+    affine_tester(:scale, [ 3, 2 ], [ 1, 1 ], 3, 2)
+    affine_tester(:scale, [ 40, 40, 40 ], [ 10, 20, -5 ], 4, 2, -8)
+  end
+
+  def test_scale_hash
+    affine_tester(:scale, [ 5, 5 ], [ 1, 1 ], :x => 5, :y => 5)
+    affine_tester(:scale, [ 3, 2 ], [ 1, 1 ], :x => 3, :y => 2)
+    affine_tester(:scale, [ 40, 40, 40 ], [ 10, 20, -5 ], :x => 4, :y => 2, :z => -8)
+  end
+
+  def test_trans_scale
+    affine_tester(:trans_scale, [ 2, 2 ], [ 1, 1 ], 1, 1, 1, 1)
+    affine_tester(:trans_scale, [ 3, 3 ], [ 2, 2 ], 1, 1, 1, 1)
+    affine_tester(:trans_scale, [ 0, 0 ], [ 1, 1 ], -1, -1, -1, -1)
+    affine_tester(:trans_scale, [ 1, 2 ], [ 1, 1 ], 0, 1, 1, 1)
+    affine_tester(:trans_scale, [ 2, 1 ], [ 1, 1 ], 1, 0, 1, 1)
+    affine_tester(:trans_scale, [ 0, 2 ], [ 1, 1 ], 1, 1, 0, 1)
+    affine_tester(:trans_scale, [ 2, 0 ], [ 1, 1 ], 1, 1, 1, 0)
+    affine_tester(:trans_scale, [ 3, 2 ], [ 1, 1 ], 2, 1, 1, 1)
+    affine_tester(:trans_scale, [ 2, 3 ], [ 1, 1 ], 1, 2, 1, 1)
+    affine_tester(:trans_scale, [ 4, 2 ], [ 1, 1 ], 1, 1, 2, 1)
+    affine_tester(:trans_scale, [ 2, 4 ], [ 1, 1 ], 1, 1, 1, 2)
+    affine_tester(:trans_scale, [ 15, 28 ], [ 1, 1 ], 2, 3, 5, 7)
+    affine_tester(:trans_scale, [ 15, 28, 1 ], [ 1, 1, 1 ], 2, 3, 5, 7)
+  end
+
+  def test_trans_scale_hash
+    affine_tester(:trans_scale, [ 2, 2 ], [ 1, 1 ], :delta_x => 1, :delta_y => 1, :x_factor => 1, :y_factor => 1)
+    affine_tester(:trans_scale, [ 15, 28, 1 ], [ 1, 1, 1 ], :delta_x => 2, :delta_y => 3, :x_factor => 5, :y_factor => 7)
+    affine_tester(:trans_scale, [ 3, 1, 1 ], [ 1, 1, 1 ], :delta_x => 2, :z_factor => 2)
+  end
+
+  def test_translate
+    affine_tester(:translate, [ 5, 12 ], [ 0, 0 ], 5, 12)
+    affine_tester(:translate, [ -3, -7, 3 ], [ 0, 0, 0 ], -3, -7, 3)
+  end
+
+  def test_translate_hash
+    affine_tester(:translate, [ 5, 12 ], [ 0, 0 ], :x => 5, :y => 12)
+    affine_tester(:translate, [ -3, -7, 3 ], [ 0, 0, 0 ], :x => -3, :y => -7, :z => 3)
+  end
 end

data/test/geometry_collection_tests.rb

@@ -1,4 +1,3 @@
-# encoding: UTF-8
 # frozen_string_literal: true
 
 require 'test_helper'
@@ -84,4 +83,416 @@ 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