proj4rb 0.3.0

data/rakefile.rb

@@ -0,0 +1,130 @@
+require 'rubygems'
+require 'rake'
+require 'rake/testtask'
+require 'rake/rdoctask'
+require 'rake/gempackagetask'
+require 'rake/clean'
+require 'date'
+require 'ftools'
+
+CLOBBER.include('pkg/*', 'proj4rb-doc/**/*', 'lib/*.so', 'lib/*.bundle', 'lib/*.dll', 'ext/*.o', 'ext/*.so', 'ext/*.bundle', 'ext/*.dll', 'ext/Makefile', 'ext/mkmf.log')
+
+desc "Create Makefile"
+file 'ext/Makefile' => ['ext/extconf.rb'] do
+    sh 'cd ext; ruby extconf.rb'
+end
+
+desc "Build from C library"
+task :build => ['ext/Makefile', 'ext/projrb.c'] do
+    sh 'cd ext; make'
+end
+
+# ------- Default Package ----------
+FILES = FileList[
+  'rakefile.rb',
+  'README',
+  'MIT-LICENSE',
+  'data/**/*',
+  'doc/**/*',
+  'example/**/*',
+  'ext/*',
+  'ext/mingw/rakefile.rb',
+  'ext/vc/*.sln',
+  'ext/vc/*.vcproj',
+  'lib/**/*.rb'
+]
+
+# Default GEM Specification
+default_spec = Gem::Specification::new do |s|
+    s.name = 'proj4rb'
+    s.version = "0.3.0"
+    s.summary = "Ruby bindings for the Proj.4 Carthographic Projection library"
+    s.description = <<-EOF
+      Proj4rb is a ruby binding for the Proj.4 Carthographic Projection library, that supports conversions between a very large number of geographic coordinate systems and datums.
+    EOF
+    s.author = 'Guilhem Vellut'
+    s.email = 'guilhem.vellut@gmail.com'
+    s.homepage = 'http://proj4rb.rubyforge.org/'
+    s.rubyforge_project = 'proj4rb'
+    s.required_ruby_version = '>= 1.8.4'
+    s.date = DateTime.now
+    
+    s.platform = Gem::Platform::RUBY
+    s.requirements << 'Proj.4 C library'
+    s.require_path = 'lib'
+    s.extensions = ["ext/extconf.rb"]
+    s.files = FILES.to_a
+    s.test_files = FileList['test/test*.rb']
+    
+    s.has_rdoc = true
+    s.extra_rdoc_files = ["README"]
+    s.rdoc_options.concat ['--main',  'README']
+end
+
+desc "Package the library as a gem"
+Rake::GemPackageTask.new(default_spec) do |pkg|
+    pkg.need_zip = true
+    pkg.need_tar = true
+end
+
+
+# ------- Windows Package ----------
+binaries = (FileList['ext/mingw/*.so',
+                     'ext/mingw/*.dll*'])
+
+win_spec = default_spec.clone
+win_spec.extensions = []
+win_spec.platform = Gem::Platform::CURRENT
+win_spec.files += binaries.map {|binaryname| "lib/#{File.basename(binaryname)}"}
+
+desc "Create Windows Gem"
+task :create_win32_gem do
+  # Copy the win32 extension built by MingW - easier to install
+  # since there are no dependencies of msvcr80.dll
+  current_dir = File.expand_path(File.dirname(__FILE__))
+
+  binaries.each do |binaryname|
+    target = File.join(current_dir, 'lib', File.basename(binaryname))
+    cp(binaryname, target)
+  end
+  
+  # Create the gem, then move it to admin/pkg
+  Gem::Builder.new(win_spec).build
+  gem_file = "#{win_spec.name}-#{win_spec.version}-#{win_spec.platform}.gem"
+  mv(gem_file, "pkg/#{gem_file}")
+
+  # Remove win extension from top level directory  
+  binaries.each do |binaryname|
+    target = File.join(current_dir, 'lib', File.basename(binaryname))
+    rm(target)
+  end
+end
+
+
+# ---------  Test Task ---------
+Rake::TestTask.new do |t|
+  t.libs << "test"
+  t.libs << "lib"
+end
+
+# ---------  RDoc Documentation ---------
+desc "Generate rdoc documentation"
+Rake::RDocTask.new("rdoc") do |rdoc|
+  rdoc.rdoc_dir = 'doc/rdoc'
+  rdoc.title    = "Proj4rb Documentation"
+  # Show source inline with line numbers
+  rdoc.options << "--inline-source" << "--line-numbers"
+  # Make the readme file the start page for the generated html
+  rdoc.options << '--main' << 'README'
+  rdoc.rdoc_files.include('doc/*.rdoc',
+                          'ext/**/*.c',
+                          'lib/**/*.rb',
+                          'README',
+                          'MIT-LICENSE')
+end
+
+task :default => :package
+
+if RUBY_PLATFORM.match(/win32/)
+  task :package => :create_win32_gem
+end

data/test/test_constants.rb

@@ -0,0 +1,20 @@
+$: << 'lib' << 'ext'
+require File.join(File.dirname(__FILE__), '..', 'lib', 'proj4')
+require 'test/unit'
+
+class ConstantsTest < Test::Unit::TestCase
+
+    def test_version
+        assert 440 < Proj4::LIBVERSION
+    end
+
+    def test_deg
+        assert_equal Math::PI/180, Proj4::DEG_TO_RAD
+    end
+
+    def test_rad
+        assert_equal 180/Math::PI, Proj4::RAD_TO_DEG
+    end
+
+end
+

data/test/test_create_projection.rb

@@ -0,0 +1,64 @@
+$: << 'lib' << 'ext'
+require File.join(File.dirname(__FILE__), '..', 'lib', 'proj4')
+require 'test/unit'
+
+class CreateProjectionTest < Test::Unit::TestCase
+
+    def setup
+        @proj_wgs84   = Proj4::Projection.new(["init=epsg:4326"])           # WGS84
+        @proj_gk      = Proj4::Projection.new(["init=epsg:31467"])          # Gauss-Kruger Zone 3
+        @proj_conakry = Proj4::Projection.new(["init=epsg:31528"])          # Conakry 1905 / UTM zone 28N
+        @proj_ortel   = Proj4::Projection.new(["proj=ortel", "lon_0=90w"])  # Ortelius Oval Projection
+    end
+
+    def test_has_inverse
+        assert   @proj_wgs84.hasInverse?
+        assert   @proj_gk.hasInverse?
+        assert   @proj_conakry.hasInverse?
+        assert ! @proj_ortel.hasInverse?
+    end
+
+    def test_is_latlong
+        assert   @proj_wgs84.isLatLong?
+        assert ! @proj_gk.isLatLong?
+        assert ! @proj_conakry.isLatLong?
+        assert ! @proj_ortel.isLatLong?
+    end
+
+    def test_is_geocent
+        assert_equal @proj_gk.isGeocent?, @proj_gk.isGeocentric?    # two names for same method
+        assert ! @proj_wgs84.isGeocent?
+        assert ! @proj_gk.isGeocent?
+        assert ! @proj_conakry.isGeocent?
+        assert ! @proj_ortel.isGeocent?
+    end
+
+    def test_get_def
+        assert_equal '+init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0', @proj_wgs84.getDef.strip
+        d = @proj_gk.getDef.strip
+        assert ('+init=epsg:31467 +proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs +towgs84=606.0,23.0,413.0' == d || '+init=epsg:31467 +proj=tmerc +lat_0=0 +lon_0=9 +k=1.000000 +x_0=3500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m +no_defs +towgs84=606.0,23.0,413.0' == d)
+        assert_equal '+init=epsg:31528 +proj=utm +zone=28 +a=6378249.2 +b=6356515 +towgs84=-23,259,-9,0,0,0,0 +units=m +no_defs', @proj_conakry.getDef.strip
+        assert_equal '+proj=ortel +lon_0=90w +ellps=WGS84', @proj_ortel.getDef.strip
+    end
+
+    def test_inspect
+        assert_equal '#<Proj4::Projection +init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0>', @proj_wgs84.to_s
+        assert_equal '#<Proj4::Projection +init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0>', @proj_wgs84.inspect
+    end
+
+    def test_projection
+        assert_equal 'longlat', @proj_wgs84.projection
+        assert_equal 'tmerc',   @proj_gk.projection
+        assert_equal 'utm',     @proj_conakry.projection
+        assert_equal 'ortel',   @proj_ortel.projection
+    end
+
+    def test_datum
+        assert_equal 'WGS84',   @proj_wgs84.datum
+        assert_equal 'potsdam', @proj_gk.datum
+        assert_nil              @proj_conakry.datum
+        assert_nil              @proj_ortel.datum
+    end
+
+end
+

data/test/test_datums.rb

@@ -0,0 +1,44 @@
+$: << 'lib' << 'ext'
+require File.join(File.dirname(__FILE__), '..', 'lib', 'proj4')
+require 'test/unit'
+
+if Proj4::LIBVERSION >= 449
+    class DatumsTest < Test::Unit::TestCase
+
+        def test_get_all
+            datums = Proj4::Datum.list.sort.collect{ |u| u.id }
+            assert datums.index('WGS84')
+            assert datums.index('potsdam')
+            assert datums.index('ire65')
+        end
+
+        def test_one
+            datum = Proj4::Datum.get('potsdam')
+            assert_kind_of Proj4::Datum, datum
+            assert_equal 'potsdam', datum.id
+            assert_equal 'potsdam', datum.to_s
+            assert_equal 'bessel', datum.ellipse_id
+            assert_equal 'towgs84=606.0,23.0,413.0', datum.defn
+            assert_equal 'Potsdam Rauenberg 1950 DHDN', datum.comments
+            assert_equal '#<Proj4::Datum id="potsdam", ellipse_id="bessel", defn="towgs84=606.0,23.0,413.0", comments="Potsdam Rauenberg 1950 DHDN">', datum.inspect
+        end
+
+        def test_compare
+            u1 = Proj4::Datum.get('potsdam')
+            u2 = Proj4::Datum.get('potsdam')
+            assert u1 == u2
+        end
+
+        def test_failed_get
+            datum = Proj4::Datum.get('foo')
+            assert_nil datum
+        end
+
+        def test_new
+            assert_raise TypeError do
+                Proj4::Datum.new
+            end
+        end
+
+    end
+end

data/test/test_ellipsoids.rb

@@ -0,0 +1,45 @@
+$: << 'lib' << 'ext'
+require File.join(File.dirname(__FILE__), '..', 'lib', 'proj4')
+require 'test/unit'
+
+if Proj4::LIBVERSION >= 449
+    class EllipsoidsTest < Test::Unit::TestCase
+
+        def test_get_all
+            ellipsoids = Proj4::Ellipsoid.list.sort.collect{ |u| u.id }
+            assert ellipsoids.index('WGS84')
+            assert ellipsoids.index('bessel')
+            assert ellipsoids.index('lerch')
+        end
+
+        def test_one
+            ellipsoid = Proj4::Ellipsoid.get('bessel')
+            assert_kind_of Proj4::Ellipsoid, ellipsoid
+            assert_equal 'bessel', ellipsoid.id
+            assert_equal 'a=6377397.155', ellipsoid.major
+            assert_equal 'rf=299.1528128', ellipsoid.ell
+            assert_equal 'Bessel 1841', ellipsoid.name
+            assert_equal 'bessel', ellipsoid.to_s
+            assert_equal '#<Proj4::Ellipsoid id="bessel", major="a=6377397.155", ell="rf=299.1528128", name="Bessel 1841">', ellipsoid.inspect
+        end
+
+        def test_compare
+            e1 = Proj4::Ellipsoid.get('bessel')
+            e2 = Proj4::Ellipsoid.get('bessel')
+            assert e1 == e2
+        end
+
+        def test_failed_get
+            ellipsoid = Proj4::Ellipsoid.get('foo')
+            assert_nil ellipsoid
+        end
+
+        def test_new
+            assert_raise TypeError do
+                Proj4::Ellipsoid.new
+            end
+        end
+
+    end
+end
+

data/test/test_errors.rb

@@ -0,0 +1,70 @@
+$: << 'lib' << 'ext'
+require File.join(File.dirname(__FILE__), '..', 'lib', 'proj4')
+require 'test/unit'
+
+class ErrorsTest < Test::Unit::TestCase
+
+    def test_list
+        assert_equal "Unknown",             Proj4::Error.error(0)
+        assert_equal "NoOptionsInInitFile", Proj4::Error.error(2)
+        assert_equal "NoOptionsInInitFile", Proj4::Error.error(-2)
+        assert_equal "Unknown",             Proj4::Error.error(-2000)
+    end
+
+    def test_parenting
+        assert_kind_of Proj4::Error,  Proj4::UnknownError.new
+        assert_kind_of Proj4::Error,  Proj4::ToleranceConditionError.new
+        assert_kind_of StandardError, Proj4::UnknownError.new
+    end
+
+    def test_num
+        assert 0, Proj4::UnknownError.errnum
+        assert 1, Proj4::NoArgsInInitListError.errnum
+    end
+
+    def test_raise
+        assert_raise Proj4::UnknownError do
+            Proj4::Error.raise_error(0)
+        end
+        assert_raise Proj4::NoOptionsInInitFileError do
+            Proj4::Error.raise_error(2)
+        end
+        assert_raise Proj4::ProjectionNotNamedError do
+            Proj4::Error.raise_error(-4)
+        end
+        assert_raise Proj4::UnknownError do
+            Proj4::Error.raise_error(2000)
+        end
+    end
+
+    def test_strerrno
+        assert_equal 'no arguments in initialization list', Proj4::Error.message(-1)
+        assert_equal 'reciprocal flattening (1/f) = 0',     Proj4::Error.message(-10)
+        assert_equal 'unknown error',                       Proj4::Error.message(0)
+        assert_match /^invalid projection system error/,    Proj4::Error.message(-2000)
+    end
+
+    def test_raise_err0
+        begin
+            Proj4::Error.raise_error(0)
+        rescue => exception
+            assert_equal Proj4::UnknownError, exception.class
+            assert_equal "unknown error", exception.message
+            assert_equal 0, exception.errnum
+            assert_match %r{test/test_errors.rb:[0-9]+:in .test_raise_err0.$} , exception.backtrace[0]
+        end
+    end
+
+    def test_raise_err1
+        begin
+            Proj4::Error.raise_error(1)
+        rescue => exception
+            assert_equal Proj4::NoArgsInInitListError, exception.class
+            assert_equal 'no arguments in initialization list', exception.message
+            assert_equal 1, exception.errnum
+            assert_match %r{test/test_errors.rb:[0-9]+:in .test_raise_err1.$} , exception.backtrace[0]
+        end
+    end
+
+end
+

data/test/test_init_projection.rb

@@ -0,0 +1,108 @@
+$: << 'lib' << 'ext'
+require File.join(File.dirname(__FILE__), '..', 'lib', 'proj4')
+require 'test/unit'
+
+class InitProjectionTest < Test::Unit::TestCase
+
+    def setup
+        # NAD27(76) / UTM zone 17N
+        # <2029> +proj=utm +zone=17 +ellps=clrk66 +units=m +no_defs  no_defs <>
+        @epsg2029i = ['init=epsg:2029']
+        @epsg2029x = ['proj=utm', 'zone=17', 'ellps=clrk66', 'units=m', 'no_defs'].sort
+
+        @nad = ' +init=epsg:2029 +proj=utm +zone=17 +ellps=clrk66 +units=m +no_defs'
+    end
+
+    def test_arg_fail
+        assert_raise ArgumentError do
+            Proj4::Projection._parse_init_parameters()
+        end
+        assert_raise ArgumentError do
+            Proj4::Projection._parse_init_parameters(nil)
+        end
+        assert_raise ArgumentError do
+            Proj4::Projection._parse_init_parameters(1)
+        end
+    end
+
+    def test_arg_string
+        args = Proj4::Projection._parse_init_parameters('init=epsg:2029')
+        assert_equal @epsg2029i, args
+        args = Proj4::Projection._parse_init_parameters('  proj=utm zone=17 ellps=clrk66 units=m no_defs  ')
+        assert_equal @epsg2029x, args.sort
+    end
+
+    def test_arg_string_with_plus
+        args = Proj4::Projection._parse_init_parameters('+init=epsg:2029')
+        assert_equal @epsg2029i, args
+        args = Proj4::Projection._parse_init_parameters('+proj=utm +zone=17 +ellps=clrk66 +units=m +no_defs')
+        assert_equal @epsg2029x, args.sort
+    end
+
+    def test_arg_array
+        args = Proj4::Projection._parse_init_parameters(['init=epsg:2029'])
+        assert_equal @epsg2029i, args
+        args = Proj4::Projection._parse_init_parameters(['proj=utm', 'zone=17', 'ellps=clrk66', 'units=m', 'no_defs'])
+        assert_equal @epsg2029x, args.sort
+    end
+
+    def test_arg_array_with_plus
+        args = Proj4::Projection._parse_init_parameters(['+init=epsg:2029'])
+        assert_equal @epsg2029i, args
+        args = Proj4::Projection._parse_init_parameters(['+proj=utm', '+zone=17', '+ellps=clrk66', '+units=m', '+no_defs'])
+        assert_equal @epsg2029x, args.sort
+    end
+
+    def test_arg_hash_with_string
+        args = Proj4::Projection._parse_init_parameters( {'init' => 'epsg:2029'} )
+        assert_equal @epsg2029i, args
+        args = Proj4::Projection._parse_init_parameters( {'proj' => 'utm', 'zone' => '17', 'ellps' => 'clrk66', 'units' => 'm', 'no_defs' => nil} )
+        assert_equal @epsg2029x, args.sort
+    end
+
+    def test_arg_hash_with_symbol
+        args = Proj4::Projection._parse_init_parameters( {:init => 'epsg:2029'} )
+        assert_equal @epsg2029i, args
+        args = Proj4::Projection._parse_init_parameters( {:proj => 'utm', :zone => '17', :ellps => 'clrk66', :units => 'm', :no_defs => nil} )
+        assert_equal @epsg2029x, args.sort
+    end
+
+    def test_arg_hash_with_symbol_simple
+        args = Proj4::Projection._parse_init_parameters( :init => 'epsg:2029' )
+        assert_equal @epsg2029i, args
+        args = Proj4::Projection._parse_init_parameters( :proj => 'utm', :zone => '17', :ellps => 'clrk66', :units => 'm', :no_defs => nil )
+        assert_equal @epsg2029x, args.sort
+    end
+
+    def test_arg_projection
+        proj = Proj4::Projection.new(['init=epsg:2029'])
+        args = Proj4::Projection._parse_init_parameters(proj)
+        assert_equal [@epsg2029i, @epsg2029x].flatten.sort, args.sort
+    end
+
+    def test_init_arg_string
+        proj = Proj4::Projection.new('+init=epsg:2029')
+        assert_equal @nad, proj.getDef
+    end
+
+    def test_init_arg_array
+        proj = Proj4::Projection.new(['init=epsg:2029'])
+        assert_equal @nad, proj.getDef
+    end
+
+    def test_init_arg_hash
+        proj = Proj4::Projection.new( :proj => 'utm', 'zone' => '17', '+ellps' => 'clrk66', :units => 'm', :no_defs => nil )
+        assert_equal @epsg2029x, proj.getDef.strip.split(' ').collect{ |a| a.sub(/^\+/, '') }.sort
+    end
+
+    def test_init_arg_fail
+        assert_raise Proj4::UnknownProjectionIdError do
+            Proj4::Projection.new( :proj => 'xxxx' )
+        end
+        assert_raise Proj4::ProjectionNotNamedError do
+            Proj4::Projection.new( :foo => 'xxxx' )
+        end
+    end
+
+end
+

data/test/test_prime_meridians.rb

@@ -0,0 +1,44 @@
+$: << 'lib' << 'ext'
+require File.join(File.dirname(__FILE__), '..', 'lib', 'proj4')
+require 'test/unit'
+
+if Proj4::LIBVERSION >= 449
+    class PrimeMeridiansTest < Test::Unit::TestCase
+
+        def test_get_all
+            prime_meridians = Proj4::PrimeMeridian.list.sort.collect{ |u| u.id}
+            assert prime_meridians.index('greenwich')
+            assert prime_meridians.index('athens')
+            assert prime_meridians.index('lisbon')
+            assert prime_meridians.index('rome')
+        end
+
+        def test_one
+            prime_meridian = Proj4::PrimeMeridian.get('lisbon')
+            assert_kind_of Proj4::PrimeMeridian, prime_meridian
+            assert_equal 'lisbon', prime_meridian.id
+            assert_equal 'lisbon', prime_meridian.to_s
+            assert_equal '9d07\'54.862"W', prime_meridian.defn
+            assert_equal '#<Proj4::PrimeMeridian id="lisbon", defn="9d07\'54.862"W">', prime_meridian.inspect
+        end
+
+        def test_compare
+            u1 = Proj4::PrimeMeridian.get('lisbon')
+            u2 = Proj4::PrimeMeridian.get('lisbon')
+            assert u1 == u2
+        end
+
+        def test_failed_get
+            prime_meridian = Proj4::PrimeMeridian.get('foo')
+            assert_nil prime_meridian
+        end
+
+        def test_new
+            assert_raise TypeError do
+                Proj4::PrimeMeridian.new
+            end
+        end
+
+    end
+end
+

data/test/test_projection_type.rb

@@ -0,0 +1,43 @@
+$: << 'lib' << 'ext'
+require File.join(File.dirname(__FILE__), '..', 'lib', 'proj4')
+require 'test/unit'
+
+if Proj4::LIBVERSION >= 449
+    class ProjectionTypesTest < Test::Unit::TestCase
+
+        def test_get_all
+            pt = Proj4::ProjectionType.list.sort.collect{ |u| u.id }
+            assert pt.index('merc')
+            assert pt.index('aea')
+            assert pt.index('bipc')
+        end
+
+        def test_one
+            pt = Proj4::ProjectionType.get('merc')
+            assert_kind_of Proj4::ProjectionType, pt
+            assert_equal 'merc', pt.id
+            assert_equal 'merc', pt.to_s
+            assert_equal 'Mercator', pt.name
+            assert_equal "Mercator\n\tCyl, Sph&Ell\n\tlat_ts=", pt.descr
+            assert_equal '#<Proj4::ProjectionType id="merc", name="Mercator">', pt.inspect
+        end
+
+        def test_compare
+            pt1 = Proj4::ProjectionType.get('merc')
+            pt2 = Proj4::ProjectionType.get('merc')
+            assert pt1 == pt2
+        end
+
+        def test_failed_get
+            assert_nil Proj4::ProjectionType.get('foo')
+        end
+
+        def test_new
+            assert_raise TypeError do
+                Proj4::ProjectionType.new
+            end
+        end
+
+    end
+end
+

data/test/test_simple_projection.rb

@@ -0,0 +1,57 @@
+$: << 'lib' << 'ext'
+require File.join(File.dirname(__FILE__), '..', 'lib', 'proj4')
+require 'test/unit'
+
+class SimpleProjectionTest < Test::Unit::TestCase
+
+    def setup
+        @proj_gk = Proj4::Projection.new(["init=epsg:31467"])
+        @lon =  8.4302123334
+        @lat = 48.9906726079
+        @rw = 3458305
+        @hw = 5428192
+    end
+
+    def rad2deg(rad)
+        rad * Proj4::RAD_TO_DEG
+    end
+
+    def deg2rad(deg)
+        deg * Proj4::DEG_TO_RAD
+    end
+
+    # echo "8.4302123334 48.9906726079" | proj +init=epsg:31467 -
+    def test_forward_gk
+        result = @proj_gk.forward( Proj4::Point.new( deg2rad(@lon), deg2rad(@lat) ) ) 
+        assert_in_delta @rw, result.x, 0.1
+        assert_in_delta @hw, result.y, 0.1
+    end
+
+    def test_forward_gk_degrees
+        result = @proj_gk.forwardDeg( Proj4::Point.new( @lon, @lat ) ) 
+        assert_in_delta @rw, result.x, 0.1
+        assert_in_delta @hw, result.y, 0.1
+    end
+
+    # echo "3458305 5428192" | invproj -f '%.10f' +init=epsg:31467 -
+    def test_inverse_gk
+        result = @proj_gk.inverse( Proj4::Point.new(@rw, @hw) ) 
+        assert_in_delta @lon, rad2deg(result.x), 0.000000001
+        assert_in_delta @lat, rad2deg(result.y), 0.000000001
+    end
+
+    def test_inverse_gk_degrees
+        result = @proj_gk.inverseDeg( Proj4::Point.new(@rw, @hw) ) 
+        assert_in_delta @lon, result.x, 0.000000001
+        assert_in_delta @lat, result.y, 0.000000001
+    end
+
+    # echo "190 92" | proj +init=epsg:31467 -
+    def test_out_of_bounds
+        assert_raise Proj4::LatitudeOrLongitudeExceededLimitsError do
+            @proj_gk.forward( Proj4::Point.new( deg2rad(190), deg2rad(92) ) ) 
+        end
+    end
+
+end
+