ffi-geos 2.3.1 → 2.5.0

data/lib/ffi-geos.rb

@@ -15,6 +15,10 @@ module Geos
     File.join(GEOS_BASE, 'wkb_reader')
   autoload :WkbWriter,
     File.join(GEOS_BASE, 'wkb_writer')
+  autoload :GeoJSONReader,
+    File.join(GEOS_BASE, 'geojson_reader')
+  autoload :GeoJSONWriter,
+    File.join(GEOS_BASE, 'geojson_writer')
   autoload :CoordinateSequence,
     File.join(GEOS_BASE, 'coordinate_sequence')
   autoload :Geometry,
@@ -50,7 +54,7 @@ module Geos
 
   module FFIGeos
     def self.search_paths
-      @search_paths ||= \
+      @search_paths ||=
         if ENV['GEOS_LIBRARY_PATH']
           [ENV['GEOS_LIBRARY_PATH']]
         elsif FFI::Platform::IS_WINDOWS
@@ -77,7 +81,7 @@ module Geos
     end
 
     def self.geos_library_path
-      @geos_library_path ||= \
+      @geos_library_path ||=
         # On MingW the libraries have version numbers
         find_lib('{lib,}geos_c{,-?}')
     end
@@ -104,6 +108,11 @@ module Geos
       :ndr, 1 # Little Endian
     ])
 
+    Geos::Flavors = enum(:flavor, [
+      :extended, 1,
+      :iso, 2
+    ])
+
     Geos::BufferCapStyles = enum(:buffer_cap_style, [
       :round, 1,
       :flat, 2,
@@ -144,6 +153,11 @@ module Geos
       :keep_collapsed, 1 << 1
     ])
 
+    Geos::PolygonHullSimplifyModes = enum(:polygon_hull_simplify_mode, [
+      :vertex_ratio, 1,
+      :area_ratio, 2
+    ])
+
     FFI_LAYOUT = {
       #### Utility functions ####
 
@@ -430,6 +444,11 @@ module Geos
         :pointer, :pointer, :pointer, :double, :int, :buffer_cap_style, :buffer_join_style, :double
       ],
 
+      GEOSDensify_r: [
+        # *geom, *handle, *geom, tolerence
+        :pointer, :pointer, :pointer, :double
+      ],
+
       # Deprecated in GEOS 3.3.0.
       GEOSSingleSidedBuffer_r: [
         # *geom, *handle, *geom, width, quad_segs, buffer_join_style, mitre_limit, is_left
@@ -446,6 +465,26 @@ module Geos
         :pointer, :pointer, :pointer
       ],
 
+      GEOSConcaveHull_r: [
+        # *geom, *handle, *geom, ratio, allow_holes
+        :pointer, :pointer, :pointer, :double, :uint
+      ],
+
+      GEOSConcaveHullByLength_r: [
+        # *geom, *handle, *geom, length, allow_holes
+        :pointer, :pointer, :pointer, :double, :uint
+      ],
+
+      GEOSPolygonHullSimplify_r: [
+        # *geom, *handle, *geom, is_outer, vertex_num_fraction
+        :pointer, :pointer, :pointer, :uint, :double
+      ],
+
+      GEOSPolygonHullSimplifyMode_r: [
+        # *geom, *handle, *geom, is_outer, parameter_mode, parameter
+        :pointer, :pointer, :pointer, :uint, :uint, :double
+      ],
+
       GEOSDifference_r: [
         # *geom, *handle, *geom_a, *geom_b
         :pointer, :pointer, :pointer, :pointer
@@ -486,6 +525,11 @@ module Geos
         :pointer, :pointer, :pointer
       ],
 
+      GEOSDisjointSubsetUnion_r: [
+        # *geom, *handle, *geom
+        :pointer, :pointer, :pointer
+      ],
+
       GEOSUnaryUnion_r: [
         # *geom, *handle, *geom
         :pointer, :pointer, :pointer
@@ -522,6 +566,11 @@ module Geos
         :pointer, :pointer, :pointer
       ],
 
+      GEOSHilbertCode_r: [
+        # int, *handler, *geom, *extent, level, *code
+        :int, :pointer, :pointer, :pointer, :uint, :pointer
+      ],
+
       GEOSMinimumBoundingCircle_r: [
         # *geom, *handle, *geom, *double radius, **geom center
         :pointer, :pointer, :pointer, :pointer, :pointer
@@ -537,6 +586,16 @@ module Geos
         :pointer, :pointer, :pointer
       ],
 
+      GEOSLineMergeDirected_r: [
+        # *geom, *handle, *geom
+        :pointer, :pointer, :pointer
+      ],
+
+      GEOSLineSubstring_r: [
+        # *geom, *handle, *geom, start_fraction, end_fraction
+        :pointer, :pointer, :pointer, :double, :double
+      ],
+
       GEOSGeom_getXMin_r: [
         # 0 on exception, *handle, (double *) value
         :int, :pointer, :pointer, :pointer
@@ -587,6 +646,11 @@ module Geos
         :pointer, :pointer, :pointer, :double, :int
       ],
 
+      GEOSConstrainedDelaunayTriangulation_r: [
+        # *geom, *handle, *geom
+        :pointer, :pointer, :pointer
+      ],
+
       GEOSVoronoiDiagram_r: [
         # *geom, *handle, *geom, *envelope, tolerance, only_edges
         :pointer, :pointer, :pointer, :pointer, :double, :int
@@ -667,6 +731,11 @@ module Geos
         :char, :pointer, :pointer, :pointer, :double
       ],
 
+      GEOSEqualsIdentical_r: [
+        # (2 on exception, 1 on true, 2 on false), *handle, *geom_a, *geom_b
+        :char, :pointer, :pointer, :pointer
+      ],
+
       GEOSisEmpty_r: [
         # (2 on exception, 1 on true, 2 on false), *handle, *geom
         :char, :pointer, :pointer
@@ -702,6 +771,11 @@ module Geos
         :char, :pointer, :pointer
       ],
 
+      GEOSHasM_r: [
+        # (2 on exception, 1 on true, 2 on false), *handle, *geom
+        :char, :pointer, :pointer
+      ],
+
       GEOSisClosed_r: [
         # (2 on exception, 1 on true, 2 on false), *handle, *geom
         :char, :pointer, :pointer
@@ -727,6 +801,11 @@ module Geos
         :int, :pointer, :pointer, :pointer, :pointer
       ],
 
+      GEOSDistanceWithin_r: [
+        # (0 on exception, 1 otherwise), *handle, *geom_a, *geom_b, double distance
+        :char, :pointer, :pointer, :pointer, :double
+      ],
+
       GEOSHausdorffDistance_r: [
         # (0 on exception, 1 otherwise), *handle, *geom_a, *geom_b, (double *) distance
         :int, :pointer, :pointer, :pointer, :pointer
@@ -757,6 +836,11 @@ module Geos
         :int, :pointer, :pointer
       ],
 
+      GEOSOrientPolygons_r: [
+        # -1 on exception, *handle, *geom, exterior_cw
+        :int, :pointer, :pointer, :int
+      ],
+
       GEOSGetInteriorRingN_r: [
         # *geom, *handle, *geom, n
         :pointer, :pointer, :pointer, :int
@@ -787,6 +871,11 @@ module Geos
         :int, :pointer, :pointer, :pointer
       ],
 
+      GEOSGeomGetM_r: [
+        # -1 on exception, *handle, *geom, *point
+        :int, :pointer, :pointer, :pointer
+      ],
+
       GEOSGeomGetPointN_r: [
         # *point, *handle, *geom, n
         :pointer, :pointer, :pointer, :int
@@ -812,6 +901,11 @@ module Geos
         :double, :pointer, :pointer
       ],
 
+      GEOSConcaveHullOfPolygons_r: [
+        # *geom, *handle, *geom, length_ratio, tight, allow_holes
+        :pointer, :pointer, :pointer, :double, :uint, :uint
+      ],
+
       GEOSMinimumRotatedRectangle_r: [
         # *geom, *handle, *geom
         :pointer, :pointer, :pointer
@@ -955,6 +1049,21 @@ module Geos
         # (2 on exception, 1 on true, 0 on false), *handle, *prepared, *geom
         :char, :pointer, :pointer, :pointer
       ],
+
+      GEOSPreparedDistance_r: [
+        # (1 on success, 0 on failure), *handle, *prepared, *geom, *distance
+        :int, :pointer, :pointer, :pointer, :pointer
+      ],
+
+      GEOSPreparedDistanceWithin_r: [
+        # (1 on true, 0 on false), *handle, *prepared, *geom, max_distance
+        :char, :pointer, :pointer, :pointer, :double
+      ],
+
+      GEOSPreparedNearestPoints_r: [
+        # *coord_seq, *handle, *prepared, *geom
+        :pointer, :pointer, :pointer, :pointer
+      ],
       #### /PreparedGeometry functions ####
 
       #### WktReader functions ####
@@ -1088,8 +1197,52 @@ module Geos
         # void, *handle, *geom, bool
         :void, :pointer, :pointer, :char
       ],
+
+      GEOSWKBWriter_getFlavor_r: [
+        # flavor, *handle, *geom
+        :flavor, :pointer, :pointer
+      ],
+
+      GEOSWKBWriter_setFlavor_r: [
+        # void, *handle, *geom, flavor
+        :void, :pointer, :pointer, :flavor
+      ],
       #### /WkbWriter functions ####
 
+      #### GeoJSONReader functions ####
+      GEOSGeoJSONReader_create_r: [
+        # *geojson_reader, *handle
+        :pointer, :pointer
+      ],
+
+      GEOSGeoJSONReader_readGeometry_r: [
+        # *geom, *handle, *geojson_reader, string
+        :pointer, :pointer, :pointer, :string
+      ],
+
+      GEOSGeoJSONReader_destroy_r: [
+        # void, *handle, *geojson_reader
+        :void, :pointer, :pointer
+      ],
+      #### /GeoJSONReader functions ###
+
+      #### GeoJSONWriter functions ####
+      GEOSGeoJSONWriter_create_r: [
+        # *geojson_writer, *handle
+        :pointer, :pointer
+      ],
+
+      GEOSGeoJSONWriter_destroy_r: [
+        # void, *handle, *geojson_writer
+        :void, :pointer, :pointer
+      ],
+
+      GEOSGeoJSONWriter_writeGeometry_r: [
+        # string, *handle, *geojson_writer, :geom, :indent
+        :string, :pointer, :pointer, :pointer, :int
+      ],
+      #### /GeoJSONWriter functions ####
+
       #### Linearref functions ####
       GEOSProject_r: [
         # distance, *handle, *geom_a, *geom_b
@@ -1252,7 +1405,7 @@ module Geos
       end
 
       def default_error_handler(*args)
-        raise Geos::GEOSException, args[0] % args[1]
+        raise Geos::GEOSException, sprintf(args[0], *args[1])
       end
   end
 

data/test/coordinate_sequence_tests.rb

@@ -12,16 +12,19 @@ class CoordinateSequenceTests < Minitest::Test
 
   def test_set_and_get_x
     @cs.set_x(0, 10.01)
+
     assert_in_delta(10.01, @cs.get_x(0), TOLERANCE)
   end
 
   def test_set_and_get_y
     @cs.set_y(0, 20.02)
+
     assert_in_delta(20.02, @cs.get_y(0), TOLERANCE)
   end
 
   def test_set_and_get_z
     @cs.set_z(0, 20.02)
+
     assert_in_delta(20.02, @cs.get_z(0), TOLERANCE)
   end
 
@@ -53,7 +56,7 @@ class CoordinateSequenceTests < Minitest::Test
       [0, 0]
     ])
 
-    assert(cs.counter_clockwise?)
+    assert_predicate(cs, :counter_clockwise?)
 
     cs = Geos::CoordinateSequence.new([
       [0, 0],
@@ -62,7 +65,7 @@ class CoordinateSequenceTests < Minitest::Test
       [0, 0]
     ])
 
-    refute(cs.counter_clockwise?)
+    refute_predicate(cs, :counter_clockwise?)
   end
 
   def test_check_bounds
@@ -115,6 +118,7 @@ class CoordinateSequenceTests < Minitest::Test
 
   def test_with_no_arguments
     cs = Geos::CoordinateSequence.new
+
     assert_equal(0, cs.size)
     assert_equal(3, cs.dimensions)
   end
@@ -153,6 +157,7 @@ class CoordinateSequenceTests < Minitest::Test
 
   def test_to_point
     cs = Geos::CoordinateSequence.new([5, 7])
+
     assert_equal('POINT (5 7)', write(cs.to_point, trim: true))
   end
 
@@ -170,9 +175,11 @@ class CoordinateSequenceTests < Minitest::Test
 
   def test_empty
     cs = Geos::CoordinateSequence.new
+
     assert_geom_empty(cs)
 
     cs = Geos::CoordinateSequence.new([4, 1])
+
     refute_geom_empty(cs)
   end
 
@@ -219,11 +226,13 @@ class CoordinateSequenceTests < Minitest::Test
 
   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)
   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)
   end
 
@@ -358,37 +367,45 @@ class CoordinateSequenceTests < Minitest::Test
 
   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')
+
+    assert_predicate(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')
+
+    assert_predicate(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
 
@@ -422,6 +439,7 @@ class CoordinateSequenceTests < Minitest::Test
 
       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

data/test/geojson_reader_tests.rb

@@ -0,0 +1,174 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeoJSONReaderTests < Minitest::Test
+  include TestHelper
+
+  attr_reader :json_reader
+
+  def setup
+    super
+
+    skip unless ENV['FORCE_TESTS'] || Geos::FFIGeos.respond_to?(:GEOSGeoJSONReader_create_r)
+
+    @json_reader = Geos::GeoJSONReader.new
+
+    @writer.rounding_precision = 3
+  end
+
+  def json_read(*args, **options)
+    json_reader.read(*args, **options)
+  end
+
+  def geojson_tester(expected, json, **options)
+    assert_equal(expected, write(json_read(json, **options)))
+  end
+
+  def test_point
+    geojson_tester(
+      'POINT (-117.000 33.000)',
+      '{"type":"Point","coordinates":[-117.0,33.0]}'
+    )
+  end
+
+  def test_line_string
+    geojson_tester(
+      'LINESTRING (102.000 0.000, 103.000 1.000, 104.000 0.000, 105.000 1.000)',
+      '{"type":"LineString","coordinates":[[102.0,0.0],[103.0,1.0],[104.0,0.0],[105.0,1.0]]}'
+    )
+  end
+
+  def test_polygon
+    geojson_tester(
+      'POLYGON ((30.000 10.000, 40.000 40.000, 20.000 40.000, 10.000 20.000, 30.000 10.000))',
+      '{"type":"Polygon","coordinates":[[[30,10],[40,40],[20,40],[10,20],[30,10]]]}'
+    )
+  end
+
+  def test_polygon_with_inner_ring
+    geojson_tester(
+      'POLYGON ((35.000 10.000, 45.000 45.000, 15.000 40.000, 10.000 20.000, 35.000 10.000), (20.000 30.000, 35.000 35.000, 30.000 20.000, 20.000 30.000))',
+      '{"type":"Polygon","coordinates":[[[35,10],[45,45],[15,40],[10,20],[35,10]],[[20,30],[35,35],[30,20],[20,30]]]}'
+    )
+  end
+
+  def test_multi_point
+    geojson_tester(
+      if Geos::GEOS_NICE_VERSION >= '031200'
+        'MULTIPOINT ((10.000 40.000), (40.000 30.000), (20.000 20.000), (30.000 10.000))'
+      else
+        'MULTIPOINT (10.000 40.000, 40.000 30.000, 20.000 20.000, 30.000 10.000)'
+      end,
+      '{"type":"MultiPoint","coordinates":[[10, 40], [40, 30], [20, 20], [30, 10]]}'
+    )
+  end
+
+  def test_multi_line_string
+    geojson_tester(
+      'MULTILINESTRING ((10.000 10.000, 20.000 20.000, 10.000 40.000), (40.000 40.000, 30.000 30.000, 40.000 20.000, 30.000 10.000))',
+      '{"type":"MultiLineString","coordinates":[[[10, 10], [20, 20], [10, 40]],[[40, 40], [30, 30], [40, 20], [30, 10]]]}'
+    )
+  end
+
+  def test_multi_polygon
+    geojson_tester(
+      'MULTIPOLYGON (((40.000 40.000, 20.000 45.000, 45.000 30.000, 40.000 40.000)), ((20.000 35.000, 10.000 30.000, 10.000 10.000, 30.000 5.000, 45.000 20.000, 20.000 35.000), (30.000 20.000, 20.000 15.000, 20.000 25.000, 30.000 20.000)))',
+      '{"type": "MultiPolygon", "coordinates": [[[[40, 40], [20, 45], [45, 30], [40, 40]]], [[[20, 35], [10, 30], [10, 10], [30, 5], [45, 20], [20, 35]], [[30, 20], [20, 15], [20, 25], [30, 20]]]]}'
+    )
+  end
+
+  def test_geometry_collection
+    geojson_tester(
+      'GEOMETRYCOLLECTION (POINT (40.000 10.000), LINESTRING (10.000 10.000, 20.000 20.000, 10.000 40.000), POLYGON ((40.000 40.000, 20.000 45.000, 45.000 30.000, 40.000 40.000)))',
+      '{"type": "GeometryCollection","geometries": [{"type": "Point","coordinates": [40, 10]},{"type": "LineString","coordinates": [[10, 10], [20, 20], [10, 40]]},{"type": "Polygon","coordinates": [[[40, 40], [20, 45], [45, 30], [40, 40]]]}]}'
+    )
+  end
+
+  def test_feature_collection
+    geojson_tester(
+      'GEOMETRYCOLLECTION (POINT (-117.000 33.000), POINT (-122.000 45.000))',
+      '{"type":"FeatureCollection","features":[{"type":"Feature","geometry":{"type":"Point","coordinates":[-117.0,33.0]}},{"type":"Feature","geometry":{"type":"Point","coordinates":[-122.0,45.0]}}]}'
+    )
+  end
+
+  def test_empty_point
+    geojson_tester(
+      'POINT EMPTY',
+      '{"type":"Point","coordinates":[]}'
+    )
+  end
+
+  def test_empty_line_string
+    geojson_tester(
+      'LINESTRING EMPTY',
+      '{"type":"LineString","coordinates":[]}'
+    )
+  end
+
+  def test_empty_polygon
+    geojson_tester(
+      'POLYGON EMPTY',
+      '{"type":"Polygon","coordinates":[]}'
+    )
+  end
+
+  def test_empty_multi_point
+    geojson_tester(
+      'MULTIPOINT EMPTY',
+      '{"type":"MultiPoint","coordinates":[]}'
+    )
+  end
+
+  def test_empty_multi_line_string
+    geojson_tester(
+      'MULTILINESTRING EMPTY',
+      '{"type":"MultiLineString","coordinates":[]}'
+    )
+  end
+
+  def test_empty_multi_polygon
+    geojson_tester(
+      'MULTIPOLYGON EMPTY',
+      '{"type": "MultiPolygon", "coordinates": []}'
+    )
+  end
+
+  def test_empty_geometry_collection
+    geojson_tester(
+      'GEOMETRYCOLLECTION EMPTY',
+      '{"type": "GeometryCollection","geometries": []}'
+    )
+  end
+
+  def test_incomplete_geojson
+    assert_raises(Geos::GeoJSONReader::ParseError) do
+      json_reader.read('{"type":"Point","coordinates":[-117.0]}')
+    end
+
+    assert_raises(Geos::GeoJSONReader::ParseError) do
+      json_reader.read('{"type":"LineString","coordinates":[[1,2],[2]]}')
+    end
+  end
+
+  def test_broken_geojson
+    assert_raises(Geos::GeoJSONReader::ParseError) do
+      json_reader.read('<gml>NOT_GEO_JSON</gml>')
+    end
+  end
+
+  def test_incompatible_type
+    assert_raises(Geos::GeoJSONReader::ParseError) do
+      json_reader.read('{"type":"Line","coordinates":[[1,2],[2,3]]}')
+    end
+  end
+
+  def test_srid_from_options
+    geom = json_reader.read(
+      '{"type":"Point","coordinates":[-117.0,33.0]}',
+      srid: 3857
+    )
+
+    assert_equal(3857, geom.srid)
+  end
+end

data/test/geojson_writer_tests.rb

@@ -0,0 +1,103 @@
+# frozen_string_literal: true
+
+require 'test_helper'
+
+class GeoJSONWriterTests < Minitest::Test
+  include TestHelper
+
+  attr_reader :json_writer
+
+  def setup
+    super
+
+    skip unless ENV['FORCE_TESTS'] || Geos::FFIGeos.respond_to?(:GEOSGeoJSONWriter_create_r)
+
+    @json_writer = Geos::GeoJSONWriter.new
+  end
+
+  def json_write(*args, **options)
+    json_writer.write(*args, **options)
+  end
+
+  def geojson_tester(expected, geom, **options)
+    assert_equal(expected, json_write(read(geom), **options))
+  end
+
+  def test_point
+    geojson_tester('{"type":"Point","coordinates":[-117.0,33.0]}', 'POINT(-117 33)')
+  end
+
+  def test_line_string
+    geojson_tester('{"type":"LineString","coordinates":[[102.0,0.0],[103.0,1.0],[104.0,0.0],[105.0,1.0]]}', 'LINESTRING(102.0 0.0, 103.0 1.0, 104.0 0.0, 105.0 1.0)')
+  end
+
+  def test_polygon
+    geojson_tester('{"type":"Polygon","coordinates":[[[30.0,10.0],[40.0,40.0],[20.0,40.0],[10.0,20.0],[30.0,10.0]]]}', 'POLYGON((30 10, 40 40, 20 40, 10 20, 30 10))')
+  end
+
+  def test_polygon_with_inner_ring
+    geojson_tester('{"type":"Polygon","coordinates":[[[35.0,10.0],[45.0,45.0],[15.0,40.0],[10.0,20.0],[35.0,10.0]],[[20.0,30.0],[35.0,35.0],[30.0,20.0],[20.0,30.0]]]}', 'POLYGON((35 10, 45 45, 15 40, 10 20, 35 10), (20 30, 35 35, 30 20, 20 30))')
+  end
+
+  def test_multi_point
+    geojson_tester('{"type":"MultiPoint","coordinates":[[10.0,40.0],[40.0,30.0],[20.0,20.0],[30.0,10.0]]}', 'MULTIPOINT ((10 40), (40 30), (20 20), (30 10))')
+  end
+
+  def test_multi_line_string
+    geojson_tester('{"type":"MultiLineString","coordinates":[[[10.0,10.0],[20.0,20.0],[10.0,40.0]],[[40.0,40.0],[30.0,30.0],[40.0,20.0],[30.0,10.0]]]}', 'MULTILINESTRING ((10 10, 20 20, 10 40),(40 40, 30 30, 40 20, 30 10))')
+  end
+
+  def test_multi_polygon
+    geojson_tester('{"type":"MultiPolygon","coordinates":[[[[30.0,20.0],[45.0,40.0],[10.0,40.0],[30.0,20.0]]],[[[15.0,5.0],[40.0,10.0],[10.0,20.0],[5.0,10.0],[15.0,5.0]]]]}', 'MULTIPOLYGON (((30 20, 45 40, 10 40, 30 20)),((15 5, 40 10, 10 20, 5 10, 15 5)))')
+  end
+
+  def test_geometry_collection
+    geojson_tester('{"type":"GeometryCollection","geometries":[{"type":"Point","coordinates":[1.0,1.0]},{"type":"Point","coordinates":[2.0,2.0]}]}', 'GEOMETRYCOLLECTION(POINT(1 1),POINT(2 2))')
+  end
+
+  def test_write_with_indentation
+    geojson_tester(<<~JSON.strip, 'LINESTRING(102.0 0.0, 103.0 1.0, 104.0 0.0, 105.0 1.0)', indentation: 2)
+      {
+        "type": "LineString",
+        "coordinates": [
+          [
+            102.0,
+            0.0
+          ],
+          [
+            103.0,
+            1.0
+          ],
+          [
+            104.0,
+            0.0
+          ],
+          [
+            105.0,
+            1.0
+          ]
+        ]
+      }
+    JSON
+  end
+
+  def test_empty_point
+    geojson_tester('{"type":"Point","coordinates":[]}', 'POINT EMPTY')
+  end
+
+  def test_empty_line_string
+    geojson_tester('{"type":"LineString","coordinates":[]}', 'LINESTRING EMPTY')
+  end
+
+  def test_empty_polygon
+    geojson_tester('{"type":"Polygon","coordinates":[[]]}', 'POLYGON EMPTY')
+  end
+
+  def test_empty_geometry_collection
+    geojson_tester('{"type":"GeometryCollection","geometries":[]}', 'GEOMETRYCOLLECTION EMPTY')
+  end
+
+  def test_linear_ring
+    geojson_tester('{"type":"LineString","coordinates":[[0.0,0.0],[1.0,1.0],[1.0,0.0],[0.0,0.0]]}', 'LINEARRING (0 0, 1 1, 1 0, 0 0)')
+  end
+end