proj4rb 1.0.0 → 2.0.0

data/test/prime_meridians_test.rb

@@ -0,0 +1,33 @@
+# encoding: UTF-8
+
+require_relative './abstract_test'
+
+class PrimeMeridiansTest < AbstractTest
+  def test_get_all
+    prime_meridians = Proj::PrimeMeridian.list.sort.collect {|prime_meridian| prime_meridian.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 = Proj::PrimeMeridian.get('lisbon')
+    assert_kind_of(Proj::PrimeMeridian, prime_meridian)
+    assert_equal('lisbon', prime_meridian.id)
+    assert_equal('9d07\'54.862"W', prime_meridian.defn)
+    assert_equal('#<Proj::PrimeMeridian id="lisbon", defn="9d07\'54.862"W">', prime_meridian.inspect)
+  end
+
+  def test_compare
+    u1 = Proj::PrimeMeridian.get('lisbon')
+    u2 = Proj::PrimeMeridian.get('lisbon')
+    assert u1 == u2
+  end
+
+  def test_failed_get
+    prime_meridian = Proj::PrimeMeridian.get('foo')
+    assert_nil prime_meridian
+  end
+end
+

data/test/proj_test.rb

@@ -0,0 +1,17 @@
+# encoding: UTF-8
+
+require_relative './abstract_test'
+
+class ProjTest < AbstractTest
+  def test_info
+    info = Proj.info
+    assert_equal(6, info[:major])
+    assert_equal(2, info[:minor])
+    assert_equal(1, info[:patch])
+    assert_equal('Rel. 6.2.1, November 1st, 2019', info[:release])
+    assert_equal('6.2.1', info[:version])
+    refute_nil(info[:searchpath])
+    refute(info[:paths].null?)
+    assert_equal(1, info[:path_count])
+  end
+end

data/test/projection_test.rb

@@ -0,0 +1,224 @@
+# encoding: UTF-8
+
+require_relative './abstract_test'
+
+class ProjectionTest < AbstractTest
+  PRECISION = 0.1 ** 4
+  
+  def setup
+    @proj_wgs84 = Proj::Projection.new(["init=epsg:4326"])       # WGS84
+    @proj_gk = Proj::Projection.new(["+init=epsg:31467"])      # Gauss-Kruger Zone 3
+    @proj_merc  = Proj::Projection.new(["proj=merc"])
+    @proj_conakry = Proj::Projection.new(["+init=epsg:31528"])      # Conakry 1905 / UTM zone 28N
+    @proj_ortel = Proj::Projection.new(["+proj=ortel", "+lon_0=90w"])  # Ortelius Oval Projection
+    @epsg2029i = ['+init=epsg:2029']
+    @epsg2029_args = ['+proj=utm', '+zone=17', '+ellps=clrk66', '+units=m', '+no_defs']
+  end
+
+  def test_arg_fail
+    assert_raises ArgumentError do
+      Proj::Projection.parse()
+    end
+    assert_raises ArgumentError do
+      Proj::Projection.parse(nil)
+    end
+    assert_raises ArgumentError do
+      Proj::Projection.parse(1)
+    end
+  end
+
+  def test_arg_string
+    args = Proj::Projection.parse('+init=epsg:2029')
+    assert_equal(@epsg2029i, args)
+    args = Proj::Projection.parse('  +proj=utm +zone=17 +ellps=clrk66 +units=m +no_defs  ')
+    assert_equal(@epsg2029_args, args)
+  end
+
+  def test_arg_string_with_plus
+    args = Proj::Projection.parse('+init=epsg:2029')
+    assert_equal(@epsg2029i, args)
+    args = Proj::Projection.parse('+proj=utm +zone=17 +ellps=clrk66 +units=m +no_defs')
+    assert_equal(@epsg2029_args, args)
+  end
+
+  def test_arg_array
+    args = Proj::Projection.parse(['+init=epsg:2029'])
+    assert_equal(@epsg2029i, args)
+    args = Proj::Projection.parse(['+proj=utm', '+zone=17', '+ellps=clrk66', '+units=m', '+no_defs'])
+    assert_equal(@epsg2029_args, args)
+  end
+
+  def test_arg_array_with_plus
+    args = Proj::Projection.parse(['+init=epsg:2029'])
+    assert_equal(@epsg2029i, args)
+    args = Proj::Projection.parse(['+proj=utm', '+zone=17', '+ellps=clrk66', '+units=m', '+no_defs'])
+    assert_equal(@epsg2029_args, args)
+  end
+
+  def test_arg_hash_with_string
+    args = Proj::Projection.parse('init' => 'epsg:2029')
+    assert_equal(@epsg2029i, args)
+    args = Proj::Projection.parse('proj' => 'utm', 'zone' => '17', 'ellps' => 'clrk66', 'units' => 'm', 'no_defs' => nil)
+    assert_equal(@epsg2029_args, args)
+  end
+
+  def test_arg_hash_with_symbol
+    args = Proj::Projection.parse(:init => 'epsg:2029')
+    assert_equal(@epsg2029i, args)
+    args = Proj::Projection.parse(:proj => 'utm', :zone => '17', :ellps => 'clrk66', :units => 'm', :no_defs => nil)
+    assert_equal(@epsg2029_args, args)
+  end
+
+  def test_arg_hash_with_symbol_simple
+    args = Proj::Projection.parse(:init => 'epsg:2029')
+    assert_equal(@epsg2029i, args)
+    args = Proj::Projection.parse(:proj => 'utm', :zone => '17', :ellps => 'clrk66', :units => 'm', :no_defs => nil)
+    assert_equal(@epsg2029_args, args)
+  end
+
+  def test_arg_projection
+    proj = Proj::Projection.new(['+init=epsg:2029'])
+    args = Proj::Projection.parse(proj)
+    assert_equal(["+init=epsg:2029", "+proj=utm", "+zone=17", "+ellps=clrk66", "+units=m", "+no_defs"], args)
+  end
+
+  def test_init_arg_string
+    proj = Proj::Projection.new('+init=epsg:2029')
+    assert_equal(' +init=epsg:2029 +proj=utm +zone=17 +ellps=clrk66 +units=m +no_defs', proj.getDef)
+  end
+
+  def test_init_arg_array
+    proj = Proj::Projection.new(['+init=epsg:2029'])
+    assert_equal(' +init=epsg:2029 +proj=utm +zone=17 +ellps=clrk66 +units=m +no_defs', proj.getDef)
+  end
+
+  def test_init_arg_hash
+    proj = Proj::Projection.new(:proj => 'utm', 'zone' => '17', 'ellps' => 'clrk66', :units => 'm', :no_defs => nil)
+    assert_equal(' +proj=utm +zone=17 +ellps=clrk66 +units=m +no_defs', proj.getDef)
+  end
+
+  def test_init_arg_fail
+    assert_raises Proj::Error do
+      Proj::Projection.new(:proj => 'xxxx')
+    end
+
+    assert_raises Proj::Error do
+      Proj::Projection.new(:foo => 'xxxx')
+    end
+  end
+
+  def test_is_latlong
+    assert(@proj_wgs84.isLatLong?)
+    refute(@proj_gk.isLatLong?)
+    refute(@proj_conakry.isLatLong?)
+    refute(@proj_ortel.isLatLong?)
+  end
+
+  def test_is_geocent
+    assert_equal(@proj_gk.isGeocent?, @proj_gk.isGeocentric?)  # two names for same method
+    refute(@proj_wgs84.isGeocent?)
+    refute(@proj_gk.isGeocent?)
+    refute(@proj_conakry.isGeocent?)
+    refute(@proj_ortel.isGeocent?)
+  end
+
+  def test_get_def
+    assert_equal(' +init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0', @proj_wgs84.getDef)
+    assert_equal(' +init=epsg:31467 +proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel +units=m +no_defs', @proj_gk.getDef)
+    assert_equal('+proj=ortel +lon_0=90w +ellps=GRS80', @proj_ortel.getDef.strip)
+  end
+
+  def test_to_s
+    assert_equal('#<Proj::Projection +init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0>', @proj_wgs84.to_s)
+  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_nil(@proj_gk.datum)
+    assert_nil(@proj_conakry.datum)
+  end
+
+  # echo "8.4302123334 48.9906726079" | proj +init=epsg:31467 -
+  def test_forward_gk
+    point = Proj::Point.new(8.4302123334, 48.9906726079)
+    result = @proj_gk.forward(point.to_radians)
+    assert_in_delta(3458305.0, result.x, 0.1)
+    assert_in_delta(5428192.0, result.y, 0.1)
+  end
+
+  def test_forward_gk_degrees
+    point = Proj::Point.new(8.4302123334, 48.9906726079)
+    result = @proj_gk.forwardDeg(point)
+    assert_in_delta(3458305.0, result.x, 0.1)
+    assert_in_delta(5428192.0, result.y, 0.1)
+  end
+
+  # echo "3458305 5428192" | invproj -f '%.10f' +init=epsg:31467 -
+  def test_inverse_gk
+    point = Proj::Point.new(3458305.0, 5428192.0)
+    result = @proj_gk.inverse(point).to_degrees
+    assert_in_delta(result.x, 8.4302123334, PRECISION)
+    assert_in_delta(result.y, 48.9906726079, PRECISION)
+  end
+
+  def test_inverse_gk_degrees
+    point = Proj::Point.new(3458305.0, 5428192.0)
+    result = @proj_gk.inverseDeg(point)
+    assert_in_delta(result.x, 8.4302123334, PRECISION)
+    assert_in_delta(result.y, 48.9906726079, PRECISION)
+  end
+
+  # echo "190 92" | proj +init=epsg:31467 -
+  def test_out_of_bounds
+    error = assert_raises(Proj::Error) do
+      point = Proj::Point.new(190, 92).to_radians
+      @proj_gk.forward(point)
+    end
+    assert_equal('latitude or longitude exceeded limits', error.message)
+  end
+
+  # echo "3458305 5428192" | cs2cs -f '%.10f' +init=epsg:31467 +to +init=epsg:4326 -
+  def test_gk_to_wgs84
+    from = Proj::Point.new(3458305.0, 5428192.0)
+    to = @proj_gk.transform(@proj_wgs84, from).to_degrees
+    assert_in_delta(8.4302123334, to.x, PRECISION)
+    assert_in_delta(48.9906726079, to.y, PRECISION)
+  end
+
+  # echo "8.4293092923 48.9896114523" | cs2cs -f '%.10f' +init=epsg:4326 +to +init=epsg:31467 -
+  def test_wgs84_to_gk
+    from = Proj::Point.new(8.4302123334, 48.9906726079)
+    to = @proj_wgs84.transform(@proj_gk, from.to_radians)
+    assert_in_delta(3458305.0, to.x, PRECISION)
+    assert_in_delta(5428192.0, to.y, PRECISION)
+  end
+
+  def test_mercator_at_pole_raise
+    from = Proj::Point.new(0, 90)
+    assert_raises(Proj::Error) do
+      @proj_wgs84.transform(@proj_merc, from.to_radians)
+    end
+  end
+
+  def test_collection
+    from0 = Proj::Point.new(3458305.0, 5428192.0)
+    from1 = Proj::Point.new(0, 0)
+    collection = @proj_gk.transform_all(@proj_wgs84, [from0, from1])
+
+    to0 = collection[0].to_degrees
+    to1 = collection[1].to_degrees
+
+    assert_in_delta(8.4302123334, to0.x, PRECISION)
+    assert_in_delta(48.9906726079, to0.y, PRECISION)
+
+    assert_in_delta(-20.9657785647, to1.x, PRECISION)
+    assert_in_delta(0, to1.y, PRECISION)
+  end
+end

data/test/transformation_test.rb

@@ -0,0 +1,68 @@
+# encoding: UTF-8
+
+require_relative './abstract_test'
+
+class TransformationTest < AbstractTest
+  PRECISION = 0.5
+
+  def setup
+    @crs_wgs84 = Proj::Crs.new('epsg:4326')
+    @crs_gk  = Proj::Crs.new('epsg:31467')
+  end
+
+  def test_create_from_strings
+    transform = Proj::Transformation.new('epsg:31467', 'epsg:4326')
+    assert(transform.info)
+  end
+
+  def test_create_crs
+    transform = Proj::Transformation.new(@crs_wgs84, @crs_gk)
+    assert(transform.info)
+  end
+
+  # echo "3458305 5428192" | cs2cs -f '%.10f' +init=epsg:31467 +to +init=epsg:4326 -
+  def test_gk_to_wgs84_forward
+    transform = Proj::Transformation.new(@crs_gk, @crs_wgs84)
+    from = Proj::Coordinate.new(x: 5428192.0, y: 3458305.0, z: -5.1790915237)
+    to = transform.forward(from)
+
+    assert_in_delta(48.98963932450735, to.x, PRECISION)
+    assert_in_delta(8.429263044355544, to.y, PRECISION)
+    assert_in_delta(-5.1790915237, to.z, PRECISION)
+    assert_in_delta(0, to.t, PRECISION)
+  end
+
+  def test_gk_to_wgs84_inverse
+    transform = Proj::Transformation.new(@crs_gk, @crs_wgs84)
+    from = Proj::Coordinate.new(lam: 48.9906726079, phi: 8.4302123334)
+    to = transform.inverse(from)
+
+    assert_in_delta(5428307, to.x, PRECISION)
+    assert_in_delta(3458375, to.y, PRECISION)
+    assert_in_delta(0, to.z, PRECISION)
+    assert_in_delta(0, to.t, PRECISION)
+  end
+
+  # echo "8.4293092923 48.9896114523" | cs2cs -f '%.10f' +init=epsg:4326 +to +init=epsg:31467 -
+  def test_wgs84_to_gk_forward
+    transform = Proj::Transformation.new(@crs_wgs84, @crs_gk)
+    from = Proj::Coordinate.new(lam: 48.9906726079, phi: 8.4302123334)
+    to = transform.forward(from)
+
+    assert_in_delta(5428307, to.x, PRECISION)
+    assert_in_delta(3458375, to.y, PRECISION)
+    assert_in_delta(0, to.z, PRECISION)
+    assert_in_delta(0, to.t, PRECISION)
+  end
+
+  def test_wgs84_to_gk_forward_inverse
+    transform = Proj::Transformation.new(@crs_wgs84, @crs_gk)
+    from = Proj::Coordinate.new(x: 5428192.0, y: 3458305.0, z: -5.1790915237)
+    to = transform.inverse(from)
+
+    assert_in_delta(48.98963932450735, to.x, PRECISION)
+    assert_in_delta(8.429263044355544, to.y, PRECISION)
+    assert_in_delta(-5.1790915237, to.z, PRECISION)
+    assert_in_delta(0, to.t, PRECISION)
+  end
+end

data/test/unit_test.rb

@@ -0,0 +1,47 @@
+# encoding: UTF-8
+
+require_relative './abstract_test'
+
+class UnitsTest < AbstractTest
+  def test_get_all
+    units = Proj::Unit.list.sort.collect {|unit| unit.id}
+    assert(units.include?('deg'))
+    assert(units.include?('km'))
+    assert(units.include?('m'))
+    assert(units.include?('yd'))
+    assert(units.include?('us-mi'))
+  end
+
+  def test_linear_unit
+    unit = Proj::Unit.get('km')
+    assert_kind_of Proj::Unit, unit
+    assert_equal('km', unit.id)
+    assert_equal('km', unit.to_s)
+    assert_equal('1000', unit.to_meter)
+    assert_equal(1000.0, unit.factor)
+    assert_equal('Kilometer', unit.name)
+    assert_equal('#<Proj::Unit id="km", to_meter="1000", factor="1000.0", name="Kilometer">', unit.inspect)
+  end
+
+  def test_angular_unit
+    unit = Proj::Unit.get('deg')
+    assert_kind_of Proj::Unit, unit
+    assert_equal('deg', unit.id)
+    assert_equal('deg', unit.to_s)
+    assert_equal('0.017453292519943296', unit.to_meter)
+    assert_equal(0.017453292519943295, unit.factor)
+    assert_equal('Degree', unit.name)
+    assert_equal('#<Proj::Unit id="deg", to_meter="0.017453292519943296", factor="0.017453292519943295", name="Degree">', unit.inspect)
+  end
+
+  def test_compare
+    u1 = Proj::Unit.get('km')
+    u2 = Proj::Unit.get('km')
+    assert(u1 == u2)
+  end
+
+  def test_failed_get
+    unit = Proj::Unit.get('foo')
+    assert_nil unit
+  end
+end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: proj4rb
 version: !ruby/object:Gem::Version
-  version: 1.0.0
+  version: 2.0.0
 platform: ruby
 authors:
 - Guilhem Vellut
@@ -10,73 +10,91 @@ authors:
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2014-12-14 00:00:00.000000000 Z
-dependencies: []
-description: |2
-      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 datumspec.
+date: 2019-12-31 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: ffi
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: minitest
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: yard
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+description: "    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 datumspec.\n"
 email: 
 executables: []
-extensions:
-- ext/extconf.rb
+extensions: []
 extra_rdoc_files: []
 files:
+- ChangeLog
+- Gemfile
 - MIT-LICENSE
 - README.rdoc
 - Rakefile
-- data/GL27
-- data/MD
-- data/TN
-- data/WI
-- data/WO
-- data/conus
-- data/epsg
-- data/epsg-deprecated
-- data/esri
-- data/esri.extra
-- data/hawaii
-- data/nad.lst
-- data/nad27
-- data/nad83
-- data/ntv1_can.dat
-- data/null
-- data/other.extra
-- data/proj_def.dat
-- data/prvi
-- data/stgeorge
-- data/stlrnc
-- data/stpaul
-- data/world
-- example/basic.rb
-- example/list-datums.rb
-- example/list-ellipsoids.rb
-- example/list-errors.rb
-- example/list-prime-meridians.rb
-- example/list-projection-types.rb
-- example/list-units.rb
-- example/version.rb
-- ext/Makefile
-- ext/extconf.rb
-- ext/mkmf.log
-- ext/out.log
-- ext/proj4_ruby-x64-mingw32.def
-- ext/proj4_ruby.so
-- ext/projrb.c
-- ext/projrb.o
-- ext/vc/proj4_ruby.sln
-- ext/vc/proj4_ruby.vcproj
+- lib/area.rb
+- lib/config.rb
+- lib/context.rb
+- lib/coordinate.rb
+- lib/crs.rb
+- lib/ellipsoid.rb
+- lib/error.rb
+- lib/operation.rb
+- lib/pj_object.rb
+- lib/point.rb
+- lib/prime_meridian.rb
+- lib/proj.rb
 - lib/proj4.rb
-- test/test_constants.rb
-- test/test_create_projection.rb
-- test/test_datums.rb
-- test/test_ellipsoids.rb
-- test/test_errors.rb
-- test/test_init_projection.rb
-- test/test_prime_meridians.rb
-- test/test_projection_type.rb
-- test/test_simple_projection.rb
-- test/test_suite.rb
-- test/test_transform.rb
-- test/test_units.rb
+- lib/projection.rb
+- lib/transformation.rb
+- lib/unit.rb
+- proj4rb.gemspec
+- test/abstract_test.rb
+- test/context_test.rb
+- test/coordinate_test.rb
+- test/crs_test.rb
+- test/ellipsoid_test.rb
+- test/operation_test.rb
+- test/prime_meridians_test.rb
+- test/proj_test.rb
+- test/projection_test.rb
+- test/transformation_test.rb
+- test/unit_test.rb
 homepage: https://github.com/cfis/proj4rb
 licenses:
 - MIT
@@ -89,29 +107,16 @@ required_ruby_version: !ruby/object:Gem::Requirement
   requirements:
   - - ">="
     - !ruby/object:Gem::Version
-      version: 1.8.7
+      version: 2.4.9
 required_rubygems_version: !ruby/object:Gem::Requirement
   requirements:
   - - ">="
     - !ruby/object:Gem::Version
       version: '0'
 requirements:
-- Proj.4 C library
-rubyforge_project: 
-rubygems_version: 2.4.5
+- Proj (Proj4) Library
+rubygems_version: 3.0.6
 signing_key: 
 specification_version: 4
 summary: Ruby bindings for the Proj.4 Carthographic Projection library
-test_files:
-- test/test_constants.rb
-- test/test_create_projection.rb
-- test/test_datums.rb
-- test/test_ellipsoids.rb
-- test/test_errors.rb
-- test/test_init_projection.rb
-- test/test_prime_meridians.rb
-- test/test_projection_type.rb
-- test/test_simple_projection.rb
-- test/test_suite.rb
-- test/test_transform.rb
-- test/test_units.rb
+test_files: []