rgeo-dschee 0.5.4

data/test/common/polygon_tests.rb

@@ -0,0 +1,261 @@
+# -----------------------------------------------------------------------------
+#
+# Common tests for polygon implementations
+#
+# -----------------------------------------------------------------------------
+
+require "rgeo"
+
+module RGeo
+  module Tests # :nodoc:
+    module Common # :nodoc:
+      module PolygonTests # :nodoc:
+        def test_creation_simple
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 1)
+          point3_ = @factory.point(1, 0)
+          exterior_ = @factory.linear_ring([point1_, point2_, point3_, point1_])
+          polygon_ = @factory.polygon(exterior_)
+          assert_not_nil(polygon_)
+          assert(::RGeo::Feature::Polygon === polygon_)
+          assert_equal(::RGeo::Feature::Polygon, polygon_.geometry_type)
+          assert(exterior_.eql?(polygon_.exterior_ring))
+          assert_equal(0, polygon_.num_interior_rings)
+          assert_nil(polygon_.interior_ring_n(0))
+          assert_nil(polygon_.interior_ring_n(-1))
+        end
+
+        def test_creation_one_hole
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 10)
+          point3_ = @factory.point(10, 10)
+          point4_ = @factory.point(10, 0)
+          point5_ = @factory.point(4, 4)
+          point6_ = @factory.point(5, 6)
+          point7_ = @factory.point(6, 4)
+          exterior_ = @factory.linear_ring([point1_, point2_, point3_, point4_, point1_])
+          interior_ = @factory.linear_ring([point5_, point6_, point7_, point5_])
+          polygon_ = @factory.polygon(exterior_, [interior_])
+          assert_not_nil(polygon_)
+          assert(::RGeo::Feature::Polygon === polygon_)
+          assert_equal(::RGeo::Feature::Polygon, polygon_.geometry_type)
+          assert(exterior_.eql?(polygon_.exterior_ring))
+          assert_equal(1, polygon_.num_interior_rings)
+          assert(interior_.eql?(polygon_.interior_ring_n(0)))
+          assert_nil(polygon_.interior_ring_n(1))
+          assert_nil(polygon_.interior_ring_n(-1))
+        end
+
+        def test_required_equivalences
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 1)
+          point3_ = @factory.point(1, 0)
+          exterior1_ = @factory.linear_ring([point1_, point2_, point3_, point1_])
+          poly1_ = @factory.polygon(exterior1_)
+          point4_ = @factory.point(0, 0)
+          point5_ = @factory.point(0, 1)
+          point6_ = @factory.point(1, 0)
+          exterior2_ = @factory.linear_ring([point4_, point5_, point6_, point4_])
+          poly2_ = @factory.polygon(exterior2_)
+          assert(poly1_.eql?(poly2_))
+          assert(poly1_ == poly2_)
+        end
+
+        def test_fully_equal
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 1)
+          point3_ = @factory.point(1, 0)
+          exterior1_ = @factory.linear_ring([point1_, point2_, point3_, point1_])
+          poly1_ = @factory.polygon(exterior1_)
+          point4_ = @factory.point(0, 0)
+          point5_ = @factory.point(0, 1)
+          point6_ = @factory.point(1, 0)
+          exterior2_ = @factory.linear_ring([point4_, point5_, point6_, point4_])
+          poly2_ = @factory.polygon(exterior2_)
+          assert(poly1_.rep_equals?(poly2_))
+          assert(poly1_.equals?(poly2_))
+        end
+
+        def test_geometrically_equal_but_ordered_different
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 1)
+          point3_ = @factory.point(1, 0)
+          exterior1_ = @factory.linear_ring([point1_, point2_, point3_, point1_])
+          poly1_ = @factory.polygon(exterior1_)
+          exterior2_ = @factory.linear_ring([point2_, point3_, point1_, point2_])
+          poly2_ = @factory.polygon(exterior2_)
+          assert(!poly1_.rep_equals?(poly2_))
+          assert(poly1_.equals?(poly2_))
+        end
+
+        def test_geometrically_equal_but_different_directions
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 1)
+          point3_ = @factory.point(1, 0)
+          exterior1_ = @factory.linear_ring([point1_, point2_, point3_, point1_])
+          poly1_ = @factory.polygon(exterior1_)
+          exterior2_ = @factory.linear_ring([point1_, point3_, point2_, point1_])
+          poly2_ = @factory.polygon(exterior2_)
+          assert(!poly1_.rep_equals?(poly2_))
+          assert(poly1_.equals?(poly2_))
+        end
+
+        def test_hashes_equal_for_representationally_equivalent_objects
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 1)
+          point3_ = @factory.point(1, 0)
+          exterior1_ = @factory.linear_ring([point1_, point2_, point3_, point1_])
+          poly1_ = @factory.polygon(exterior1_)
+          point4_ = @factory.point(0, 0)
+          point5_ = @factory.point(0, 1)
+          point6_ = @factory.point(1, 0)
+          exterior2_ = @factory.linear_ring([point4_, point5_, point6_, point4_])
+          poly2_ = @factory.polygon(exterior2_)
+          assert_equal(poly1_.hash, poly2_.hash)
+        end
+
+        def test_wkt_creation_simple
+          parsed_poly_ = @factory.parse_wkt("POLYGON((0 0, 0 1, 1 0, 0 0))")
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 1)
+          point3_ = @factory.point(1, 0)
+          exterior_ = @factory.linear_ring([point1_, point2_, point3_, point1_])
+          built_poly_ = @factory.polygon(exterior_)
+          assert(built_poly_.eql?(parsed_poly_))
+        end
+
+        def test_wkt_creation_one_hole
+          parsed_poly_ = @factory.parse_wkt("POLYGON((0 0, 0 10, 10 10, 10 0, 0 0), (4 4, 5 6, 6 4, 4 4))")
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 10)
+          point3_ = @factory.point(10, 10)
+          point4_ = @factory.point(10, 0)
+          point5_ = @factory.point(4, 4)
+          point6_ = @factory.point(5, 6)
+          point7_ = @factory.point(6, 4)
+          exterior_ = @factory.linear_ring([point1_, point2_, point3_, point4_, point1_])
+          interior_ = @factory.linear_ring([point5_, point6_, point7_, point5_])
+          built_poly_ = @factory.polygon(exterior_, [interior_])
+          assert(built_poly_.eql?(parsed_poly_))
+        end
+
+        def test_clone
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 1)
+          point3_ = @factory.point(1, 0)
+          exterior_ = @factory.linear_ring([point1_, point2_, point3_, point1_])
+          poly1_ = @factory.polygon(exterior_)
+          poly2_ = poly1_.clone
+          assert(poly1_.eql?(poly2_))
+          assert(exterior_.eql?(poly2_.exterior_ring))
+          assert_equal(0, poly2_.num_interior_rings)
+        end
+
+        def test_type_check
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 1)
+          point3_ = @factory.point(1, 0)
+          exterior_ = @factory.linear_ring([point1_, point2_, point3_, point1_])
+          poly_ = @factory.polygon(exterior_)
+          assert(::RGeo::Feature::Geometry.check_type(poly_))
+          assert(!::RGeo::Feature::Point.check_type(poly_))
+          assert(!::RGeo::Feature::GeometryCollection.check_type(poly_))
+          assert(::RGeo::Feature::Surface.check_type(poly_))
+          assert(::RGeo::Feature::Polygon.check_type(poly_))
+        end
+
+        def test_as_text_wkt_round_trip
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 1)
+          point3_ = @factory.point(1, 0)
+          exterior_ = @factory.linear_ring([point1_, point2_, point3_, point1_])
+          poly1_ = @factory.polygon(exterior_)
+          text_ = poly1_.as_text
+          poly2_ = @factory.parse_wkt(text_)
+          assert(poly1_.eql?(poly2_))
+        end
+
+        def test_as_binary_wkb_round_trip
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 1)
+          point3_ = @factory.point(1, 0)
+          exterior_ = @factory.linear_ring([point1_, point2_, point3_, point1_])
+          poly1_ = @factory.polygon(exterior_)
+          binary_ = poly1_.as_binary
+          poly2_ = @factory.parse_wkb(binary_)
+          assert(poly1_.eql?(poly2_))
+        end
+
+        def test_dimension
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 10)
+          point3_ = @factory.point(10, 10)
+          point4_ = @factory.point(10, 0)
+          point5_ = @factory.point(4, 4)
+          point6_ = @factory.point(5, 6)
+          point7_ = @factory.point(6, 4)
+          exterior_ = @factory.linear_ring([point1_, point2_, point3_, point4_, point1_])
+          interior_ = @factory.linear_ring([point5_, point6_, point7_, point5_])
+          poly_ = @factory.polygon(exterior_, [interior_])
+          assert_equal(2, poly_.dimension)
+        end
+
+        def test_is_empty
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 1)
+          point3_ = @factory.point(1, 0)
+          exterior_ = @factory.linear_ring([point1_, point2_, point3_, point1_])
+          poly1_ = @factory.polygon(exterior_)
+          assert(!poly1_.is_empty?)
+          poly2_ = @factory.polygon(@factory.linear_ring([]))
+          assert(poly2_.is_empty?)
+        end
+
+        def test_polygon_coordinates
+          coordinates = [
+            [[0.0, 0.0], [1.0, 0.0], [1.0, 1.0], [0.0, 1.0], [0.0, 0.0]],
+            [[0.25, 0.25], [0.75, 0.25], [0.75, 0.75], [0.25, 0.75], [0.25, 0.25]]
+          ]
+
+          ring = @factory.line_string(coordinates.first.map { |(x, y)| @factory.point x, y })
+          inner_ring = @factory.line_string(coordinates.last.map { |(x, y)| @factory.point x, y })
+          polygon = @factory.polygon ring, [inner_ring]
+          assert_equal(polygon.coordinates, coordinates)
+        end
+
+        def test_ignores_consecutive_repeated_points
+          point1_ = @factory.point(0, 0)
+          point2_ = @factory.point(0, 0)
+          point3_ = @factory.point(0, 1)
+          point4_ = @factory.point(0, 1)
+          point5_ = @factory.point(0, 1)
+          point6_ = @factory.point(1, 0)
+          point7_ = @factory.point(1, 0)
+          point8_ = @factory.point(1, 0)
+          exterior_ = @factory.linear_ring([point1_, point2_, point3_, point4_, point5_, point6_, point7_, point8_])
+
+          poly1_ = @factory.polygon(exterior_)
+          assert(!poly1_.nil?)
+
+          line_string = poly1_.exterior_ring
+
+            case line_string.class.name
+            when "RGeo::Tests::GeosCAPI::TestPolygon"
+            when "RGeo::Geos::FFILinearRingImpl"
+            when "RGeo::Geos::CAPILinearRingImpl"
+            when "RGeo::Tests::GeosFFI::TestPolygon"
+              assert(line_string.points.count == 9)
+            else
+              assert(line_string.num_points == 4)
+            end
+
+          points = line_string.points
+          assert(points.first.x == points.last.x)
+          assert(points.first.y == points.last.y)
+          assert(points.first.z == points.last.z)
+        end
+
+      end
+    end
+  end
+end

data/test/coord_sys/ogc_cs_test.rb

@@ -0,0 +1,342 @@
+# -----------------------------------------------------------------------------
+#
+# Tests for OGC CS classes
+#
+# -----------------------------------------------------------------------------
+
+require "test/unit"
+require "rgeo"
+
+module RGeo
+  module Tests # :nodoc:
+    module CoordSys # :nodoc:
+      class TestOgcCs < ::Test::Unit::TestCase # :nodoc:
+        # Handle differences in floating-point output.
+
+        def _lenient_regex_for(str_)
+          ::Regexp.new(str_.gsub(/(\d)\.(\d{10,})/) do |_m_|
+            before_ = Regexp.last_match(1)
+            after_ = Regexp.last_match(2)[0, 10]
+            "#{before_}.#{after_}\\d*"
+          end.gsub(/(\.|\[|\]|\(|\)|\$|\^|\||\+)/) { |_m_| "\\#{Regexp.last_match(1)}" })
+        end
+
+        def test_axis_info_by_value
+          obj_ = ::RGeo::CoordSys::CS::AxisInfo.create("N", ::RGeo::CoordSys::CS::AO_NORTH)
+          assert_equal("N", obj_.name)
+          assert_equal(::RGeo::CoordSys::CS::AO_NORTH, obj_.orientation)
+          assert_equal('AXIS["N",NORTH]', obj_.to_wkt)
+        end
+
+        def test_axis_info_by_name
+          obj_ = ::RGeo::CoordSys::CS::AxisInfo.create("S", "SOUTH")
+          assert_equal("S", obj_.name)
+          assert_equal(::RGeo::CoordSys::CS::AO_SOUTH, obj_.orientation)
+          assert_equal('AXIS["S",SOUTH]', obj_.to_wkt)
+          obj2_ = ::RGeo::CoordSys::CS::AxisInfo.create("S", ::RGeo::CoordSys::CS::AO_SOUTH)
+          assert_equal(obj_, obj2_)
+        end
+
+        def test_parameter
+          obj_ = ::RGeo::CoordSys::CS::ProjectionParameter.create("false_easting", 400_000)
+          assert_equal("false_easting", obj_.name)
+          assert_equal(400_000, obj_.value)
+          assert_equal('PARAMETER["false_easting",400000.0]', obj_.to_wkt)
+        end
+
+        def test_towgs84
+          obj_ = ::RGeo::CoordSys::CS::WGS84ConversionInfo.create(1, 2, 3, 4, 5, 6, 7)
+          assert_equal(1, obj_.dx)
+          assert_equal(2, obj_.dy)
+          assert_equal(3, obj_.dz)
+          assert_equal(4, obj_.ex)
+          assert_equal(5, obj_.ey)
+          assert_equal(6, obj_.ez)
+          assert_equal(7, obj_.ppm)
+          assert_equal("TOWGS84[1.0,2.0,3.0,4.0,5.0,6.0,7.0]", obj_.to_wkt)
+        end
+
+        def test_unit
+          obj_ = ::RGeo::CoordSys::CS::Unit.create("metre", 1)
+          assert_equal("metre", obj_.name)
+          assert_equal(1, obj_.conversion_factor)
+          assert_nil(obj_.authority)
+          assert_equal('UNIT["metre",1.0]', obj_.to_wkt)
+        end
+
+        def test_unit_with_authority
+          obj_ = ::RGeo::CoordSys::CS::Unit.create("metre", 1, "EPSG", 9001)
+          assert_equal("metre", obj_.name)
+          assert_equal(1, obj_.conversion_factor)
+          assert_equal("EPSG", obj_.authority)
+          assert_equal("9001", obj_.authority_code)
+          assert_equal('UNIT["metre",1.0,AUTHORITY["EPSG","9001"]]', obj_.to_wkt)
+        end
+
+        def test_linear_unit
+          obj_ = ::RGeo::CoordSys::CS::LinearUnit.create("metre", 1)
+          assert_equal(1, obj_.meters_per_unit)
+          assert_equal('UNIT["metre",1.0]', obj_.to_wkt)
+        end
+
+        def test_angular_unit
+          obj_ = ::RGeo::CoordSys::CS::AngularUnit.create("radian", 1)
+          assert_equal(1, obj_.radians_per_unit)
+          assert_equal('UNIT["radian",1.0]', obj_.to_wkt)
+        end
+
+        def test_prime_meridian
+          obj1_ = ::RGeo::CoordSys::CS::AngularUnit.create("radian", 1)
+          obj_ = ::RGeo::CoordSys::CS::PrimeMeridian.create("Greenwich", obj1_, 0, "EPSG", "8901")
+          assert_equal("Greenwich", obj_.name)
+          assert_equal(0, obj_.longitude)
+          assert_equal('PRIMEM["Greenwich",0.0,AUTHORITY["EPSG","8901"]]', obj_.to_wkt)
+        end
+
+        def test_create_flattened_sphere
+          obj1_ = ::RGeo::CoordSys::CS::LinearUnit.create("metre", 1)
+          obj_ = ::RGeo::CoordSys::CS::Ellipsoid.create_flattened_sphere("WGS 84", 6_378_137, 298.257223563, obj1_, "EPSG", "7030")
+          assert_equal("WGS 84", obj_.name)
+          assert_equal(6_378_137, obj_.semi_major_axis)
+          assert_in_delta(298.257223563, obj_.inverse_flattening, 0.1)
+          assert_in_delta(6_356_752.314245, obj_.semi_minor_axis, 0.1)
+          assert_equal("EPSG", obj_.authority)
+          assert_equal("7030", obj_.authority_code)
+          assert_equal('SPHEROID["WGS 84",6378137.0,298.257223563,AUTHORITY["EPSG","7030"]]', obj_.to_wkt)
+        end
+
+        def test_create_unflattened_sphere
+          obj1_ = ::RGeo::CoordSys::CS::LinearUnit.create("metre", 1)
+          obj_ = ::RGeo::CoordSys::CS::Ellipsoid.create_flattened_sphere("Popular Visualisation Sphere", 6_378_137, 0, obj1_, "EPSG", "7059")
+          assert_equal("Popular Visualisation Sphere", obj_.name)
+          assert_equal(6_378_137, obj_.semi_major_axis)
+          assert_equal(0, obj_.inverse_flattening)
+          assert_equal(6_378_137, obj_.semi_minor_axis)
+          assert_equal("EPSG", obj_.authority)
+          assert_equal("7059", obj_.authority_code)
+          assert_equal('SPHEROID["Popular Visualisation Sphere",6378137.0,0.0,AUTHORITY["EPSG","7059"]]', obj_.to_wkt)
+        end
+
+        def test_create_ellipsoid
+          obj1_ = ::RGeo::CoordSys::CS::LinearUnit.create("metre", 1)
+          obj_ = ::RGeo::CoordSys::CS::Ellipsoid.create_ellipsoid("WGS 84", 6_378_137, 6_356_752.314245, obj1_, "EPSG", "7030")
+          assert_in_delta(298.257223563, obj_.inverse_flattening, 0.1)
+        end
+
+        def test_create_spherical_ellipsoid
+          obj1_ = ::RGeo::CoordSys::CS::LinearUnit.create("metre", 1)
+          obj_ = ::RGeo::CoordSys::CS::Ellipsoid.create_ellipsoid("Popular Visualisation Sphere", 6_378_137, 6_378_137, obj1_, "EPSG", "7059")
+          assert_equal(0, obj_.inverse_flattening)
+        end
+
+        def test_local_datum
+          obj_ = ::RGeo::CoordSys::CS::LocalDatum.create("Random Local Datum", ::RGeo::CoordSys::CS::LD_MIN)
+          assert_equal("Random Local Datum", obj_.name)
+          assert_equal(::RGeo::CoordSys::CS::LD_MIN, obj_.datum_type)
+          assert_equal('LOCAL_DATUM["Random Local Datum",10000]', obj_.to_wkt)
+        end
+
+        def test_local_datum_with_extension
+          obj_ = ::RGeo::CoordSys::CS::LocalDatum.create("Random Local Datum", ::RGeo::CoordSys::CS::LD_MIN, nil, nil, nil, nil, nil, foo: :bar)
+          assert_equal("bar", obj_.extension(:foo))
+          assert_nil(obj_.extension(:bar))
+          assert_equal('LOCAL_DATUM["Random Local Datum",10000,EXTENSION["foo","bar"]]', obj_.to_wkt)
+        end
+
+        def test_vertical_datum
+          obj_ = ::RGeo::CoordSys::CS::VerticalDatum.create("Ordnance Datum Newlyn", ::RGeo::CoordSys::CS::VD_GEOID_MODE_DERIVED, "EPSG", "5101")
+          assert_equal("Ordnance Datum Newlyn", obj_.name)
+          assert_equal(::RGeo::CoordSys::CS::VD_GEOID_MODE_DERIVED, obj_.datum_type)
+          assert_equal("EPSG", obj_.authority)
+          assert_equal("5101", obj_.authority_code)
+          assert_equal('VERT_DATUM["Ordnance Datum Newlyn",2005,AUTHORITY["EPSG","5101"]]', obj_.to_wkt)
+        end
+
+        def test_horizontal_datum
+          obj1_ = ::RGeo::CoordSys::CS::LinearUnit.create("metre", 1)
+          obj2_ = ::RGeo::CoordSys::CS::Ellipsoid.create_ellipsoid("Popular Visualisation Sphere", 6_378_137, 6_378_137, obj1_, "EPSG", "7059")
+          obj3_ = ::RGeo::CoordSys::CS::WGS84ConversionInfo.create(0, 0, 0, 0, 0, 0, 0)
+          obj_ = ::RGeo::CoordSys::CS::HorizontalDatum.create("Popular_Visualisation_Datum", ::RGeo::CoordSys::CS::HD_GEOCENTRIC, obj2_, obj3_, "EPSG", "6055")
+          assert_equal("Popular_Visualisation_Datum", obj_.name)
+          assert_equal(::RGeo::CoordSys::CS::HD_GEOCENTRIC, obj_.datum_type)
+          assert_equal("EPSG", obj_.authority)
+          assert_equal("6055", obj_.authority_code)
+          assert_equal('DATUM["Popular_Visualisation_Datum",SPHEROID["Popular Visualisation Sphere",6378137.0,0.0,AUTHORITY["EPSG","7059"]],TOWGS84[0.0,0.0,0.0,0.0,0.0,0.0,0.0],AUTHORITY["EPSG","6055"]]', obj_.to_wkt)
+        end
+
+        def test_projection
+          obj_ = ::RGeo::CoordSys::CS::Projection.create("Transverse_Mercator", "Transverse_Mercator", [])
+          assert_equal("Transverse_Mercator", obj_.name)
+          assert_equal("Transverse_Mercator", obj_.class_name)
+          assert_equal(0, obj_.num_parameters)
+          assert_equal('PROJECTION["Transverse_Mercator"]', obj_.to_wkt)
+        end
+
+        def test_local_coordinate_system
+          obj1_ = ::RGeo::CoordSys::CS::LocalDatum.create("Random Local Datum", ::RGeo::CoordSys::CS::LD_MIN)
+          obj2_ = ::RGeo::CoordSys::CS::LinearUnit.create("metre", 1, "EPSG", 9001)
+          obj3_ = ::RGeo::CoordSys::CS::AxisInfo.create("N", ::RGeo::CoordSys::CS::AO_NORTH)
+          obj4_ = ::RGeo::CoordSys::CS::AxisInfo.create("E", ::RGeo::CoordSys::CS::AO_EAST)
+          obj_ = ::RGeo::CoordSys::CS::LocalCoordinateSystem.create("My CS", obj1_, obj2_, [obj3_, obj4_])
+          assert_equal("My CS", obj_.name)
+          assert_equal(2, obj_.dimension)
+          assert_equal("Random Local Datum", obj_.local_datum.name)
+          assert_equal("N", obj_.get_axis(0).name)
+          assert_equal("E", obj_.get_axis(1).name)
+          assert_equal("metre", obj_.get_units(0).name)
+          assert_equal("metre", obj_.get_units(1).name)
+          assert_equal('LOCAL_CS["My CS",LOCAL_DATUM["Random Local Datum",10000],UNIT["metre",1.0,AUTHORITY["EPSG","9001"]],AXIS["N",NORTH],AXIS["E",EAST]]', obj_.to_wkt)
+        end
+
+        def test_geocentric_coordinate_system
+          obj1_ = ::RGeo::CoordSys::CS::Ellipsoid.create_flattened_sphere("WGS 84", 6_378_137, 298.257223563, nil, "EPSG", "7030")
+          obj2_ = ::RGeo::CoordSys::CS::HorizontalDatum.create("World Geodetic System 1984", ::RGeo::CoordSys::CS::HD_GEOCENTRIC, obj1_, nil, "EPSG", "6326")
+          obj3_ = ::RGeo::CoordSys::CS::PrimeMeridian.create("Greenwich", nil, 0.0, "EPSG", "8901")
+          obj4_ = ::RGeo::CoordSys::CS::LinearUnit.create("m", 1.0)
+          obj5_ = ::RGeo::CoordSys::CS::AxisInfo.create("Geocentric X", ::RGeo::CoordSys::CS::AO_OTHER)
+          obj6_ = ::RGeo::CoordSys::CS::AxisInfo.create("Geocentric Y", ::RGeo::CoordSys::CS::AO_EAST)
+          obj7_ = ::RGeo::CoordSys::CS::AxisInfo.create("Geocentric Z", ::RGeo::CoordSys::CS::AO_NORTH)
+          obj_ = ::RGeo::CoordSys::CS::GeocentricCoordinateSystem.create("WGS 84 (geocentric)", obj2_, obj3_, obj4_, obj5_, obj6_, obj7_, "EPSG", 4328)
+          assert_equal("WGS 84 (geocentric)", obj_.name)
+          assert_equal(3, obj_.dimension)
+          assert_equal("World Geodetic System 1984", obj_.horizontal_datum.name)
+          assert_equal("Greenwich", obj_.prime_meridian.name)
+          assert_equal("m", obj_.linear_unit.name)
+          assert_equal("Geocentric X", obj_.get_axis(0).name)
+          assert_equal("Geocentric Y", obj_.get_axis(1).name)
+          assert_equal("Geocentric Z", obj_.get_axis(2).name)
+          assert_equal("m", obj_.get_units(0).name)
+          assert_equal("m", obj_.get_units(1).name)
+          assert_equal("m", obj_.get_units(2).name)
+          assert_equal('GEOCCS["WGS 84 (geocentric)",DATUM["World Geodetic System 1984",SPHEROID["WGS 84",6378137.0,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0.0,AUTHORITY["EPSG","8901"]],UNIT["m",1.0],AXIS["Geocentric X",OTHER],AXIS["Geocentric Y",EAST],AXIS["Geocentric Z",NORTH],AUTHORITY["EPSG","4328"]]', obj_.to_wkt)
+        end
+
+        def test_vertical_coordinate_system
+          obj1_ = ::RGeo::CoordSys::CS::VerticalDatum.create("Ordnance Datum Newlyn", ::RGeo::CoordSys::CS::VD_GEOID_MODE_DERIVED, "EPSG", 5101)
+          obj2_ = ::RGeo::CoordSys::CS::LinearUnit.create("metre", 1, "EPSG", 9001)
+          obj3_ = ::RGeo::CoordSys::CS::AxisInfo.create("Up", ::RGeo::CoordSys::CS::AO_UP)
+          obj_ = ::RGeo::CoordSys::CS::VerticalCoordinateSystem.create("Newlyn", obj1_, obj2_, obj3_, "EPSG", 5701)
+          assert_equal("Newlyn", obj_.name)
+          assert_equal(1, obj_.dimension)
+          assert_equal("Ordnance Datum Newlyn", obj_.vertical_datum.name)
+          assert_equal("metre", obj_.vertical_unit.name)
+          assert_equal("Up", obj_.get_axis(0).name)
+          assert_equal("metre", obj_.get_units(0).name)
+          assert_equal('VERT_CS["Newlyn",VERT_DATUM["Ordnance Datum Newlyn",2005,AUTHORITY["EPSG","5101"]],UNIT["metre",1.0,AUTHORITY["EPSG","9001"]],AXIS["Up",UP],AUTHORITY["EPSG","5701"]]', obj_.to_wkt)
+        end
+
+        def test_geographic_coordinate_system
+          obj1_ = ::RGeo::CoordSys::CS::Ellipsoid.create_flattened_sphere("WGS 84", 6_378_137, 298.257223563, nil, "EPSG", "7030")
+          obj2_ = ::RGeo::CoordSys::CS::AngularUnit.create("degree", 0.01745329251994328, "EPSG", 9122)
+          obj3_ = ::RGeo::CoordSys::CS::HorizontalDatum.create("WGS_1984", ::RGeo::CoordSys::CS::HD_GEOCENTRIC, obj1_, nil, "EPSG", "6326")
+          obj4_ = ::RGeo::CoordSys::CS::PrimeMeridian.create("Greenwich", nil, 0, "EPSG", "8901")
+          obj_ = ::RGeo::CoordSys::CS::GeographicCoordinateSystem.create("WGS 84", obj2_, obj3_, obj4_, nil, nil, "EPSG", 4326)
+          assert_equal("WGS 84", obj_.name)
+          assert_equal(2, obj_.dimension)
+          assert_equal("WGS_1984", obj_.horizontal_datum.name)
+          assert_equal("Greenwich", obj_.prime_meridian.name)
+          assert_equal("degree", obj_.angular_unit.name)
+          assert_nil(obj_.get_axis(0))
+          assert_nil(obj_.get_axis(1))
+          assert_equal("degree", obj_.get_units(0).name)
+          assert_equal("degree", obj_.get_units(1).name)
+          assert_match(_lenient_regex_for('GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137.0,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0.0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]'), obj_.to_wkt)
+        end
+
+        def test_projected_coordinate_system
+          obj1_ = ::RGeo::CoordSys::CS::Ellipsoid.create_flattened_sphere("Airy 1830", 6_377_563.396, 299.3249646, nil, "EPSG", "7001")
+          obj2_ = ::RGeo::CoordSys::CS::WGS84ConversionInfo.create(375, -111, 431, 0, 0, 0, 0)
+          obj3_ = ::RGeo::CoordSys::CS::AngularUnit.create("DMSH", 0.0174532925199433, "EPSG", 9108)
+          obj4_ = ::RGeo::CoordSys::CS::HorizontalDatum.create("OSGB_1936", ::RGeo::CoordSys::CS::HD_CLASSIC, obj1_, obj2_, "EPSG", "6277")
+          obj5_ = ::RGeo::CoordSys::CS::PrimeMeridian.create("Greenwich", nil, 0, "EPSG", "8901")
+          obj6_ = ::RGeo::CoordSys::CS::AxisInfo.create("Lat", ::RGeo::CoordSys::CS::AO_NORTH)
+          obj7_ = ::RGeo::CoordSys::CS::AxisInfo.create("Long", ::RGeo::CoordSys::CS::AO_EAST)
+          obj8_ = ::RGeo::CoordSys::CS::GeographicCoordinateSystem.create("OSGB 1936", obj3_, obj4_, obj5_, obj6_, obj7_, "EPSG", 4277)
+          obj9_ = ::RGeo::CoordSys::CS::ProjectionParameter.create("latitude_of_origin", 49)
+          obj10_ = ::RGeo::CoordSys::CS::ProjectionParameter.create("central_meridian", -2)
+          obj11_ = ::RGeo::CoordSys::CS::ProjectionParameter.create("scale_factor", 0.999601272)
+          obj12_ = ::RGeo::CoordSys::CS::ProjectionParameter.create("false_easting", 400_000)
+          obj13_ = ::RGeo::CoordSys::CS::ProjectionParameter.create("false_northing", -100_000)
+          obj14_ = ::RGeo::CoordSys::CS::Projection.create("Transverse_Mercator", "Transverse_Mercator", [obj9_, obj10_, obj11_, obj12_, obj13_])
+          obj15_ = ::RGeo::CoordSys::CS::LinearUnit.create("metre", 1, "EPSG", 9001)
+          obj16_ = ::RGeo::CoordSys::CS::AxisInfo.create("E", ::RGeo::CoordSys::CS::AO_EAST)
+          obj17_ = ::RGeo::CoordSys::CS::AxisInfo.create("N", ::RGeo::CoordSys::CS::AO_NORTH)
+          obj_ = ::RGeo::CoordSys::CS::ProjectedCoordinateSystem.create("OSGB 1936 / British National Grid", obj8_, obj14_, obj15_, obj16_, obj17_, "EPSG", 27_700)
+          assert_equal("OSGB 1936 / British National Grid", obj_.name)
+          assert_equal(2, obj_.dimension)
+          assert_equal("OSGB_1936", obj_.horizontal_datum.name)
+          assert_equal("OSGB 1936", obj_.geographic_coordinate_system.name)
+          assert_equal("Transverse_Mercator", obj_.projection.name)
+          assert_equal(5, obj_.projection.num_parameters)
+          assert_equal("latitude_of_origin", obj_.projection.get_parameter(0).name)
+          assert_equal(49, obj_.projection.get_parameter(0).value)
+          assert_equal("false_northing", obj_.projection.get_parameter(4).name)
+          assert_equal(-100_000, obj_.projection.get_parameter(4).value)
+          assert_equal("metre", obj_.linear_unit.name)
+          assert_equal("E", obj_.get_axis(0).name)
+          assert_equal("N", obj_.get_axis(1).name)
+          assert_equal("metre", obj_.get_units(0).name)
+          assert_equal("metre", obj_.get_units(1).name)
+          assert_equal('PROJCS["OSGB 1936 / British National Grid",GEOGCS["OSGB 1936",DATUM["OSGB_1936",SPHEROID["Airy 1830",6377563.396,299.3249646,AUTHORITY["EPSG","7001"]],TOWGS84[375.0,-111.0,431.0,0.0,0.0,0.0,0.0],AUTHORITY["EPSG","6277"]],PRIMEM["Greenwich",0.0,AUTHORITY["EPSG","8901"]],UNIT["DMSH",0.0174532925199433,AUTHORITY["EPSG","9108"]],AXIS["Lat",NORTH],AXIS["Long",EAST],AUTHORITY["EPSG","4277"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",49.0],PARAMETER["central_meridian",-2.0],PARAMETER["scale_factor",0.999601272],PARAMETER["false_easting",400000.0],PARAMETER["false_northing",-100000.0],UNIT["metre",1.0,AUTHORITY["EPSG","9001"]],AXIS["E",EAST],AXIS["N",NORTH],AUTHORITY["EPSG","27700"]]', obj_.to_wkt)
+        end
+
+        def test_parse_epsg_6055
+          input_ = 'DATUM["Popular_Visualisation_Datum",SPHEROID["Popular Visualisation Sphere",6378137.0,0.0,AUTHORITY["EPSG","7059"]],TOWGS84[0.0,0.0,0.0,0.0,0.0,0.0,0.0],AUTHORITY["EPSG","6055"]]'
+          obj_ = ::RGeo::CoordSys::CS.create_from_wkt(input_)
+          assert_kind_of(::RGeo::CoordSys::CS::HorizontalDatum, obj_)
+          assert_equal("Popular_Visualisation_Datum", obj_.name)
+          assert_equal(::RGeo::CoordSys::CS::HD_GEOCENTRIC, obj_.datum_type)
+          assert_equal("EPSG", obj_.authority)
+          assert_equal("6055", obj_.authority_code)
+          assert_equal("Popular Visualisation Sphere", obj_.ellipsoid.name)
+          assert_equal(6_378_137, obj_.ellipsoid.semi_major_axis)
+          assert_equal(0, obj_.ellipsoid.inverse_flattening)
+          assert_equal(6_378_137, obj_.ellipsoid.semi_minor_axis)
+          assert_equal("EPSG", obj_.ellipsoid.authority)
+          assert_equal("7059", obj_.ellipsoid.authority_code)
+          assert_equal(0, obj_.wgs84_parameters.dx)
+          assert_equal(input_, obj_.to_wkt)
+        end
+
+        def test_parse_epsg_7405
+          input_ = 'COMPD_CS["OSGB36 / British National Grid + ODN",PROJCS["OSGB 1936 / British National Grid",GEOGCS["OSGB 1936",DATUM["OSGB 1936",SPHEROID["Airy 1830",6377563.396,299.3249646,AUTHORITY["EPSG","7001"]],TOWGS84[446.448,-125.157,542.06,0.15,0.247,0.842,-4.2261596151967575],AUTHORITY["EPSG","6277"]],PRIMEM["Greenwich",0.0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.017453292519943295],AXIS["Geodetic latitude",NORTH],AXIS["Geodetic longitude",EAST],AUTHORITY["EPSG","4277"]],PROJECTION["Transverse Mercator",AUTHORITY["EPSG","9807"]],PARAMETER["central_meridian",-2.0],PARAMETER["latitude_of_origin",49.0],PARAMETER["scale_factor",0.9996012717],PARAMETER["false_easting",400000.0],PARAMETER["false_northing",-100000.0],UNIT["m",1.0],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","27700"]],VERT_CS["Newlyn",VERT_DATUM["Ordnance Datum Newlyn",2005,AUTHORITY["EPSG","5101"]],UNIT["m",1.0],AXIS["Gravity-related height",UP],AUTHORITY["EPSG","5701"]],AUTHORITY["EPSG","7405"]]'
+          obj_ = ::RGeo::CoordSys::CS.create_from_wkt(input_)
+          assert_kind_of(::RGeo::CoordSys::CS::CompoundCoordinateSystem, obj_)
+          assert_kind_of(::RGeo::CoordSys::CS::ProjectedCoordinateSystem, obj_.head)
+          assert_kind_of(::RGeo::CoordSys::CS::VerticalCoordinateSystem, obj_.tail)
+          assert_equal(3, obj_.dimension)
+          assert_match(_lenient_regex_for(input_), obj_.to_wkt)
+        end
+
+        def test_parse_local_datum_with_extension
+          input_ = 'LOCAL_DATUM["Random Local Datum",10000,EXTENSION["foo","bar"]]'
+          obj_ = ::RGeo::CoordSys::CS.create_from_wkt(input_)
+          assert_kind_of(::RGeo::CoordSys::CS::LocalDatum, obj_)
+          assert_equal("bar", obj_.extension(:foo))
+          assert_nil(obj_.extension(:bar))
+          assert_equal(input_, obj_.to_wkt)
+        end
+
+        def test_marshal_roundtrip
+          input_ = 'COMPD_CS["OSGB36 / British National Grid + ODN",PROJCS["OSGB 1936 / British National Grid",GEOGCS["OSGB 1936",DATUM["OSGB 1936",SPHEROID["Airy 1830",6377563.396,299.3249646,AUTHORITY["EPSG","7001"]],TOWGS84[446.448,-125.157,542.06,0.15,0.247,0.842,-4.2261596151967575],AUTHORITY["EPSG","6277"]],PRIMEM["Greenwich",0.0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.017453292519943295],AXIS["Geodetic latitude",NORTH],AXIS["Geodetic longitude",EAST],AUTHORITY["EPSG","4277"]],PROJECTION["Transverse Mercator",AUTHORITY["EPSG","9807"]],PARAMETER["central_meridian",-2.0],PARAMETER["latitude_of_origin",49.0],PARAMETER["scale_factor",0.9996012717],PARAMETER["false_easting",400000.0],PARAMETER["false_northing",-100000.0],UNIT["m",1.0],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","27700"]],VERT_CS["Newlyn",VERT_DATUM["Ordnance Datum Newlyn",2005,AUTHORITY["EPSG","5101"]],UNIT["m",1.0],AXIS["Gravity-related height",UP],AUTHORITY["EPSG","5701"]],AUTHORITY["EPSG","7405"]]'
+          obj1_ = ::RGeo::CoordSys::CS.create_from_wkt(input_)
+          dump_ = ::Marshal.dump(obj1_)
+          obj2_ = ::Marshal.load(dump_)
+          assert_equal(obj1_, obj2_)
+        end
+
+        if ::RGeo.yaml_supported?
+
+          def test_yaml_roundtrip
+            input_ = 'COMPD_CS["OSGB36 / British National Grid + ODN",PROJCS["OSGB 1936 / British National Grid",GEOGCS["OSGB 1936",DATUM["OSGB 1936",SPHEROID["Airy 1830",6377563.396,299.3249646,AUTHORITY["EPSG","7001"]],TOWGS84[446.448,-125.157,542.06,0.15,0.247,0.842,-4.2261596151967575],AUTHORITY["EPSG","6277"]],PRIMEM["Greenwich",0.0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.017453292519943295],AXIS["Geodetic latitude",NORTH],AXIS["Geodetic longitude",EAST],AUTHORITY["EPSG","4277"]],PROJECTION["Transverse Mercator",AUTHORITY["EPSG","9807"]],PARAMETER["central_meridian",-2.0],PARAMETER["latitude_of_origin",49.0],PARAMETER["scale_factor",0.9996012717],PARAMETER["false_easting",400000.0],PARAMETER["false_northing",-100000.0],UNIT["m",1.0],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","27700"]],VERT_CS["Newlyn",VERT_DATUM["Ordnance Datum Newlyn",2005,AUTHORITY["EPSG","5101"]],UNIT["m",1.0],AXIS["Gravity-related height",UP],AUTHORITY["EPSG","5701"]],AUTHORITY["EPSG","7405"]]'
+            obj1_ = ::RGeo::CoordSys::CS.create_from_wkt(input_)
+            dump_ = ::Psych.dump(obj1_)
+            obj2_ = ::Psych.load(dump_)
+            assert_equal(obj1_, obj2_)
+          end
+
+        end
+      end
+    end
+  end
+end