rgeo 0.2.1 → 0.2.2

data/lib/rgeo/coord_sys/srs_database/url_reader.rb

@@ -0,0 +1,90 @@
+# -----------------------------------------------------------------------------
+# 
+# SRS database interface
+# 
+# -----------------------------------------------------------------------------
+# 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 'net/http'
+
+
+module RGeo
+  
+  module CoordSys
+    
+    module SRSDatabase
+      
+      
+      class UrlReader
+        
+        
+        def initialize(opts_={})
+          @cache = opts_[:cache] ? {} : nil
+        end
+        
+        
+        def get(ident_)
+          ident_ = ident_.to_s
+          return @cache[ident_] if @cache && @cache.include?(ident_)
+          uri_ = ::URI.parse(ident_)
+          result_ = nil
+          ::Net::HTTP.start(uri_.host, uri_.port) do |http_|
+            request_ = uri_.path
+            request_ = "#{request_}?#{uri_.query}" if uri_.query
+            response_ = http_.request_get(request_)
+            if response_.kind_of?(::Net::HTTPSuccess)
+              response_ = response_.body.strip
+              if response_[0,1] == '+'
+                result_ = Entry.new(ident_, :proj4 => response_)
+              else
+                result_ = Entry.new(ident_, :coord_sys => response_)
+              end
+            end
+          end
+          @cache[ident_] = result_ if @cache
+          result_
+        end
+        
+        
+        def clear_cache
+          @cache.clear if @cache
+        end
+        
+        
+      end
+      
+      
+    end
+    
+  end
+  
+end

data/lib/rgeo/feature/types.rb

@@ -39,6 +39,14 @@ module RGeo
   module Feature
     
     
+    # All geometry implementations MUST include this submodule.
+    # This serves as a marker that may be used to test an object for
+    # feature-ness.
+    
+    module Instance
+    end
+    
+    
     # These methods are available as module methods (not instance methods)
     # of the various feature types.
     # For example, you may determine whether a feature object is a
@@ -57,12 +65,8 @@ module RGeo
     module Type
       
       
-      # All geometry implementations MUST include this submodule.
-      # This serves as a marker that may be used to test an object for
-      # feature-ness.
-      
-      module Instance
-      end
+      # Deprecated alias for RGeo::Feature::Instance
+      Instance = Feature::Instance
       
       
       # Returns true if the given object is this type or a subtype
@@ -72,7 +76,7 @@ module RGeo
       # Note that feature objects need not actually include this module.
       
       def check_type(rhs_)
-        rhs_ = rhs_.geometry_type if rhs_.kind_of?(Instance)
+        rhs_ = rhs_.geometry_type if rhs_.kind_of?(Feature::Instance)
         rhs_.kind_of?(Type) && (rhs_ == self || rhs_.include?(self))
       end
       alias_method :===, :check_type

data/lib/rgeo/geos/factory.rb

@@ -74,8 +74,13 @@ module RGeo
           else
             proj4_ = nil
           end
+          coord_sys_ = opts_[:coord_sys]
+          if coord_sys_.kind_of?(::String)
+            coord_sys_ = CoordSys::CS.create_from_wkt(coord_sys_) rescue nil
+          end
           result_ = _create(flags_, opts_[:srid].to_i, buffer_resolution_)
           result_.instance_variable_set(:@proj4, proj4_)
+          result_.instance_variable_set(:@coord_sys, coord_sys_)
           result_
         end
         alias_method :new, :create
@@ -233,7 +238,7 @@ module RGeo
       # See ::RGeo::Feature::Factory#coord_sys
       
       def coord_sys
-        nil
+        @coord_sys
       end
       
       

data/lib/rgeo/geos/impl_additions.rb

@@ -40,7 +40,7 @@ module RGeo
     
     class GeometryImpl  # :nodoc:
       
-      include ::RGeo::Feature::Type::Instance
+      include Feature::Instance
       
       def inspect
         "#<#{self.class}:0x#{object_id.to_s(16)} #{as_text.inspect}>"

data/lib/rgeo/geos/zm_impl.rb

@@ -41,7 +41,7 @@ module RGeo
     
     class ZMGeometryImpl  # :nodoc:
       
-      include Feature::Type::Instance
+      include Feature::Instance
       
       
       def initialize(factory_, zgeometry_, mgeometry_)

data/lib/rgeo/impl_helper/basic_geometry_methods.rb

@@ -41,7 +41,7 @@ module RGeo
     
     module BasicGeometryMethods  # :nodoc:
       
-      include ::RGeo::Feature::Type::Instance
+      include Feature::Instance
       
       
       def inspect  # :nodoc:

data/test/coord_sys/tc_active_record_table.rb

@@ -0,0 +1,97 @@
+# -----------------------------------------------------------------------------
+# 
+# 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_
+          
+          require 'active_record'
+          
+          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; 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