schleyfox-rgeo 0.2.5

data/test/coord_sys/tc_active_record_table.rb

@@ -0,0 +1,102 @@
+# -----------------------------------------------------------------------------
+# 
+# Tests for OGC CS classes
+# 
+# -----------------------------------------------------------------------------
+# Copyright 2010 Daniel Azuma
+# 
+# All rights reserved.
+# 
+# Redistribution and use in source and binary forms, with or without
+# modification, are permitted provided that the following conditions are met:
+# 
+# * Redistributions of source code must retain the above copyright notice,
+#   this list of conditions and the following disclaimer.
+# * Redistributions in binary form must reproduce the above copyright notice,
+#   this list of conditions and the following disclaimer in the documentation
+#   and/or other materials provided with the distribution.
+# * Neither the name of the copyright holder, nor the names of any other
+#   contributors to this software, may be used to endorse or promote products
+#   derived from this software without specific prior written permission.
+# 
+# THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+# AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+# IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+# ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
+# LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+# CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+# SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+# INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+# CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+# ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+# POSSIBILITY OF SUCH DAMAGE.
+# -----------------------------------------------------------------------------
+;
+
+
+require 'test/unit'
+require 'rgeo'
+require 'yaml'
+
+
+module RGeo
+  module Tests  # :nodoc:
+    module CoordSys  # :nodoc:
+      module ActiveRecordTableTests  # :nodoc:
+        
+        database_configs_ = ::YAML.load_file(::File.dirname(__FILE__)+'/database.yml') rescue nil
+        if database_configs_
+          begin
+            require 'active_record'
+          rescue ::LoadError
+            database_configs_ = nil
+          end
+        end
+        
+        if database_configs_
+          
+          PostGIS_CONFIG = database_configs_['postgis'] rescue nil
+          
+          
+          if PostGIS_CONFIG
+            
+            class TestPostGIS < ::Test::Unit::TestCase  # :nodoc:
+              
+              class ARBase < ::ActiveRecord::Base
+                establish_connection(PostGIS_CONFIG)
+              end
+              
+              @@db = ::RGeo::CoordSys::SRSDatabase::ActiveRecordTable.new(:ar_base_class => ARBase, :auth_name_column => 'auth_name', :auth_srid_column => 'auth_srid', :srtext_column => 'srtext', :proj4text_column => 'proj4text')
+              
+              
+              def test_4326
+                entry_ = @@db.get(4326)
+                assert_equal('+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs', entry_.proj4.original_str)
+                assert_kind_of(::RGeo::CoordSys::CS::GeographicCoordinateSystem, entry_.coord_sys)
+                assert_equal('WGS 84', entry_.name)
+              end
+              
+              
+              def test_3785
+                entry_ = @@db.get(3785)
+                assert_equal('+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +units=m +k=1.0 +nadgrids=@null +no_defs', entry_.proj4.original_str)
+                assert_kind_of(::RGeo::CoordSys::CS::ProjectedCoordinateSystem, entry_.coord_sys)
+                assert_equal('Popular Visualisation CRS / Mercator (deprecated)', entry_.name)
+              end
+              
+              
+            end
+            
+          else
+            puts "WARNING: No postgis section in database.yml; skipping PostGIS tests."
+          end
+          
+        
+        else
+          puts "WARNING: Couldn't find database.yml or ActiveRecord gem is not present; skipping ActiveRecord tests."
+        end
+        
+      end
+    end
+  end
+end

data/test/coord_sys/tc_ogc_cs.rb

@@ -0,0 +1,356 @@
+# -----------------------------------------------------------------------------
+# 
+# Tests for OGC CS classes
+# 
+# -----------------------------------------------------------------------------
+# Copyright 2010 Daniel Azuma
+# 
+# All rights reserved.
+# 
+# Redistribution and use in source and binary forms, with or without
+# modification, are permitted provided that the following conditions are met:
+# 
+# * Redistributions of source code must retain the above copyright notice,
+#   this list of conditions and the following disclaimer.
+# * Redistributions in binary form must reproduce the above copyright notice,
+#   this list of conditions and the following disclaimer in the documentation
+#   and/or other materials provided with the distribution.
+# * Neither the name of the copyright holder, nor the names of any other
+#   contributors to this software, may be used to endorse or promote products
+#   derived from this software without specific prior written permission.
+# 
+# THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+# AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+# IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+# ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
+# LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+# CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+# SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+# INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+# CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+# ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+# POSSIBILITY OF SUCH DAMAGE.
+# -----------------------------------------------------------------------------
+;
+
+
+require 'test/unit'
+require 'rgeo'
+
+
+module RGeo
+  module Tests  # :nodoc:
+    module CoordSys  # :nodoc:
+      
+      class TestOgcCs < ::Test::Unit::TestCase  # :nodoc:
+        
+        
+        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', 400000)
+          assert_equal('false_easting', obj_.name)
+          assert_equal(400000, 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', 6378137, 298.257223563, obj1_, 'EPSG', '7030')
+          assert_equal('WGS 84', obj_.name)
+          assert_equal(6378137, obj_.semi_major_axis)
+          assert_in_delta(298.257223563, obj_.inverse_flattening, 0.1)
+          assert_in_delta(6356752.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', 6378137, 0, obj1_, 'EPSG', '7059')
+          assert_equal('Popular Visualisation Sphere', obj_.name)
+          assert_equal(6378137, obj_.semi_major_axis)
+          assert_equal(0, obj_.inverse_flattening)
+          assert_equal(6378137, 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', 6378137, 6356752.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', 6378137, 6378137, 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_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', 6378137, 6378137, 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', 6378137, 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', 6378137, 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_equal('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', 6377563.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', 400000)
+          obj13_ = ::RGeo::CoordSys::CS::ProjectionParameter.create('false_northing', -100000)
+          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', 27700)
+          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(-100000, 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(6378137, obj_.ellipsoid.semi_major_axis)
+          assert_equal(0, obj_.ellipsoid.inverse_flattening)
+          assert_equal(6378137, 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_equal(input_, obj_.to_wkt)
+        end
+        
+        
+      end
+      
+    end
+  end
+end