dr-ffi-proj4 0.0.1

data/lib/ffi-proj4/tools.rb

@@ -0,0 +1,57 @@
+
+module Proj4
+  module Tools
+    def bool_result(r)
+      case r
+        when 1
+          true
+        when 0
+          false
+        else
+          raise RuntimeError.new("Unexpected boolean result: #{r}")
+      end
+    end
+
+    def rad_to_deg!(rad)
+      unless rad.nil?
+        case rad
+          when Proj4::Point, Proj4::ProjXY
+            rad.to_deg!
+          else
+            rad * Proj4::RAD_TO_DEG
+        end
+      end
+    end
+
+    def rad_to_deg(rad)
+      rad_to_deg!(
+        !rad.is_a?(Numeric) && rad.respond_to?(:dup) ?
+          rad.dup :
+          rad
+      )
+    end
+
+    def deg_to_rad!(deg)
+      unless deg.nil?
+        case deg
+          when Proj4::Point, Proj4::ProjXY
+            deg.to_rad!
+          else
+            deg * Proj4::DEG_TO_RAD
+        end
+      end
+    end
+
+    def deg_to_rad(deg)
+      deg_to_rad!(
+        !deg.is_a?(Numeric) && deg.respond_to?(:dup) ?
+          deg.dup :
+          deg
+      )
+    end
+
+    class << self
+      include Proj4::Tools
+    end
+  end
+end

data/lib/ffi-proj4/version.rb

@@ -0,0 +1,3 @@
+module Proj4
+  VERSION = "0.0.1"
+end

data/test/point_tests.rb

@@ -0,0 +1,69 @@
+
+$: << File.dirname(__FILE__)
+require 'test_helper'
+
+class PointTests < MiniTest::Unit::TestCase
+  include TestHelper
+
+  def rad_deg_tester(method, expected_x, expected_y, expected_z, x, y, z)
+    point = Proj4::Point.new(x, y, z)
+    point.send("#{method}!")
+
+    assert_in_delta(expected_x, point.x, TOLERANCE)
+    assert_in_delta(expected_y, point.y, TOLERANCE)
+    assert_in_delta(expected_z, point.z, TOLERANCE)
+
+    point = Proj4::Point.new(x, y, z)
+    point2 = point.send(method)
+
+    assert_in_delta(x, point.x, TOLERANCE)
+    assert_in_delta(y, point.y, TOLERANCE)
+    assert_in_delta(z, point.z, TOLERANCE)
+
+    assert_in_delta(expected_x, point2.x, TOLERANCE)
+    assert_in_delta(expected_y, point2.y, TOLERANCE)
+    assert_in_delta(expected_z, point2.z, TOLERANCE)
+  end
+
+  def test_to_rad
+    rad_deg_tester(:to_rad,
+      1.0471975511965976,
+      1.5707963267948966,
+      0.7853981633974483,
+      60,
+      90,
+      45
+    )
+  end
+
+  def test_to_deg
+    rad_deg_tester(:to_deg,
+      60,
+      90,
+      45,
+      1.0471975511965976,
+      1.5707963267948966,
+      0.7853981633974483
+    )
+  end
+
+  def test_aliases
+    point = Proj4::Point.new(10, 20, 30)
+    assert_equal(10, point.x)
+    assert_equal(20, point.y)
+    assert_equal(10, point.lon)
+    assert_equal(20, point.lat)
+
+    point.lon = 15
+    point.lat = 25
+
+    assert_equal(15, point.x)
+    assert_equal(25, point.y)
+
+    point.x = 80
+    point.y = 85
+
+    assert_equal(80, point.lon)
+    assert_equal(85, point.lat)
+  end
+end

data/test/projection_tests.rb

@@ -0,0 +1,197 @@
+
+$: << File.dirname(__FILE__)
+require 'test_helper'
+
+class ProjectionTests < MiniTest::Unit::TestCase
+  include TestHelper
+
+  def definition_sorter(definition)
+    definition.split(/\s+/).sort.join(' ')
+  end
+
+  def definition_tester(expected, proj)
+    assert_equal(
+      definition_sorter(expected),
+      definition_sorter(Proj4::Projection.new(proj).definition)
+    )
+  end
+
+  def hash_tester(expected, proj)
+    hash = proj.to_hash
+    assert_equal(expected.size, hash.size)
+
+    expected.each do |k, v|
+      assert_equal(v, hash[k], " Expected #{v} for #{k}, got #{hash[k]}")
+    end
+  end
+
+  def test_read_strings
+    definition_tester(
+      '+init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0',
+      'init=epsg:4326'
+    )
+
+    definition_tester(
+      '+datum=potsdam +ellps=bessel +init=epsg:31467 +k=1 +lat_0=0 +lon_0=9 +no_defs +proj=tmerc +towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 +units=m +x_0=3500000 +y_0=0',
+      'init=epsg:31467'
+    )
+
+    definition_tester(
+      '+init=epsg:31528 +proj=utm +zone=28 +a=6378249.2 +b=6356515 +towgs84=-23,259,-9,0,0,0,0 +units=m +no_defs',
+      'init=epsg:31528'
+    )
+
+    definition_tester(
+      '+proj=ortel +lon_0=90w +ellps=WGS84',
+      'proj=ortel lon_0=90w'
+    )
+
+    assert_raises(Proj4::ProjectionParseError) do
+      definition_tester(
+        '',
+        'gibberish'
+      )
+    end
+  end
+
+  def test_read_arrays
+    definition_tester(
+      '+init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0',
+      [ 'init=epsg:4326' ]
+    )
+
+    definition_tester(
+      '+datum=potsdam +ellps=bessel +init=epsg:31467 +k=1 +lat_0=0 +lon_0=9 +no_defs +proj=tmerc +towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 +units=m +x_0=3500000 +y_0=0',
+      [ 'init=epsg:31467' ]
+    )
+
+    definition_tester(
+      '+init=epsg:31528 +proj=utm +zone=28 +a=6378249.2 +b=6356515 +towgs84=-23,259,-9,0,0,0,0 +units=m +no_defs',
+      [ 'init=epsg:31528' ]
+    )
+
+    definition_tester(
+      '+proj=ortel +lon_0=90w +ellps=WGS84',
+      [ 'proj=ortel', 'lon_0=90w' ]
+    )
+
+    assert_raises(Proj4::ProjectionParseError) do
+      definition_tester(
+        '',
+        [ 'gibberish' ]
+      )
+    end
+  end
+
+  def test_read_hashes
+    definition_tester(
+      '+init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0',
+      { :init => 'epsg:4326' }
+    )
+
+    definition_tester(
+      '+datum=potsdam +ellps=bessel +init=epsg:31467 +k=1 +lat_0=0 +lon_0=9 +no_defs +proj=tmerc +towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 +units=m +x_0=3500000 +y_0=0',
+      { :init => 'epsg:31467' }
+    )
+
+    definition_tester(
+      '+init=epsg:31528 +proj=utm +zone=28 +a=6378249.2 +b=6356515 +towgs84=-23,259,-9,0,0,0,0 +units=m +no_defs',
+      { :init => 'epsg:31528' }
+    )
+
+    definition_tester(
+      '+lon_0=90w +proj=ortel +ellps=WGS84',
+      { :proj => 'ortel', :lon_0 => '90w' }
+    )
+
+    assert_raises(Proj4::ProjectionParseError) do
+      definition_tester(
+        '',
+        { :proj => 'gibberish' }
+      )
+    end
+  end
+
+  def test_has_inverse
+    # XXX - checking if a projection has an inverse isn't supported by the
+    # PROJ.4 public API.
+    # assert(PROJ_WGS84.hasInverse?)
+    # assert(PROJ_GK.hasInverse?)
+    # assert(PROJ_CONAKRY.hasInverse?)
+    # assert(!PROJ_ORTEL.hasInverse?)
+  end
+
+  def test_is_latlong
+    %w{ isLatLong? lat_long? }.each do |method|
+      assert(PROJ_WGS84.send(method))
+      assert(!PROJ_GK.send(method))
+      assert(!PROJ_CONAKRY.send(method))
+      assert(!PROJ_ORTEL.send(method))
+    end
+  end
+
+  def test_is_geocent
+    %w{ isGeocent? isGeocentric? geocentric? }.each do |method|
+      assert(!PROJ_WGS84.send(method))
+      assert(!PROJ_GK.send(method))
+      assert(!PROJ_CONAKRY.send(method))
+      assert(!PROJ_ORTEL.send(method))
+    end
+  end
+
+  def test_get_def
+    %w{ getDef definition }.each do |method|
+      assert_equal(
+        '+init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0',
+        PROJ_WGS84.send(method).strip
+      )
+
+      assert_equal(
+        '+init=epsg:31467 +proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +datum=potsdam +units=m +no_defs +ellps=bessel +towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7',
+        PROJ_GK.send(method).strip
+      )
+
+      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.send(method).strip
+      )
+
+      assert_equal(
+        '+proj=ortel +lon_0=90w +ellps=WGS84',
+        PROJ_ORTEL.send(method).strip
+      )
+    end
+  end
+
+  def test_inspect
+    assert_equal('#<Proj4::Projection +init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0>', PROJ_WGS84.to_s)
+    assert_equal('#<Proj4::Projection +init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0>', PROJ_WGS84.inspect)
+  end
+
+  def test_to_hash
+    hash_tester(PROJ_WGS84_HASH, PROJ_WGS84)
+    hash_tester(PROJ_GK_HASH, PROJ_GK)
+    hash_tester(PROJ_CONAKRY_HASH, PROJ_CONAKRY)
+    hash_tester(PROJ_ORTEL_HASH, PROJ_ORTEL)
+  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
+
+  def test_shortcut_create
+    proj = Proj4::Projection.new("epsg:4326")
+
+    hash_tester(PROJ_WGS84_HASH, proj)
+  end
+end

data/test/simple_transformation_tests.rb

@@ -0,0 +1,82 @@
+
+$: << File.dirname(__FILE__)
+require 'test_helper'
+
+class SimpleTransformationTests < MiniTest::Unit::TestCase
+  include TestHelper
+
+  LONG = 8.4302123334
+  LAT = 48.9906726079
+  RW = 3458305
+  HW = 5428192
+
+  def simple_transformation_tester(method, expected_x, expected_y, proj, x, y, tolerance = TOLERANCE)
+    point = Proj4::Point.new(x, y)
+    result = proj.send("#{method}!", point)
+
+    assert_equal(result.object_id, point.object_id)
+    assert_in_delta(expected_x, result.x, tolerance)
+    assert_in_delta(expected_y, result.y, tolerance)
+
+    point = Proj4::Point.new(x, y)
+    result = proj.send(method, point)
+
+    refute_equal(result.object_id, point.object_id)
+    assert_in_delta(x, point.x, tolerance)
+    assert_in_delta(y, point.y, tolerance)
+    assert_in_delta(expected_x, result.x, tolerance)
+    assert_in_delta(expected_y, result.y, tolerance)
+  end
+
+  def test_forward
+    simple_transformation_tester(
+      :forward,
+      RW,
+      HW,
+      PROJ_GK,
+      Proj4::Tools.deg_to_rad(LONG),
+      Proj4::Tools.deg_to_rad(LAT),
+      0.1
+    )
+  end
+
+  def test_inverse
+    simple_transformation_tester(
+      :inverse,
+      Proj4::Tools.deg_to_rad(LONG),
+      Proj4::Tools.deg_to_rad(LAT),
+      PROJ_GK,
+      RW,
+      HW
+    )
+  end
+
+  def test_forward_deg
+    simple_transformation_tester(
+      :forward_deg,
+      RW,
+      HW,
+      PROJ_GK,
+      LONG,
+      LAT,
+      0.1
+    )
+  end
+
+  def test_inverse_deg
+    simple_transformation_tester(
+      :inverse_deg,
+      LONG,
+      LAT,
+      PROJ_GK,
+      RW,
+      HW
+    )
+  end
+
+  def test_out_of_bounds
+    assert_raises Proj4::LatitudeOrLongitudeExceededLimitsError do
+      PROJ_GK.forward_deg(190, 92)
+    end
+  end
+end

data/test/test_helper.rb

@@ -0,0 +1,78 @@
+
+require 'rubygems'
+require 'minitest/autorun'
+
+if RUBY_VERSION >= '1.9'
+  require 'minitest/reporters'
+end
+
+if ENV['USE_BINARY_PROJ4']
+  require 'proj4_ruby'
+else
+  require File.join(File.dirname(__FILE__), %w{ .. lib ffi-proj4 })
+end
+
+Proj4.proj_lib = File.join(File.dirname(__FILE__), %w{ .. data }) if Proj4.respond_to?(:proj_lib)
+
+puts "Ruby version #{RUBY_VERSION} - #{RbConfig::CONFIG['RUBY_INSTALL_NAME']}"
+puts "ffi-proj4 version #{Proj4::VERSION}"
+puts "PROJ version #{Proj4.version}"
+if defined?(Proj4::FFIProj4)
+  puts "Using #{Array(Proj4::FFIProj4.proj4_library_path).join(', ')}"
+end
+puts "Using PROJ_LIB #{Proj4.proj_lib}"
+
+module TestHelper
+  TOLERANCE = 0.000001
+
+  PROJ_WGS84 = Proj4::Projection.new('init=epsg:4326')
+  PROJ_GK = Proj4::Projection.new('init=epsg:31467')
+  PROJ_CONAKRY = Proj4::Projection.new('init=epsg:31528')
+  PROJ_ORTEL = Proj4::Projection.new([ 'proj=ortel', 'lon_0=90w' ])
+
+  PROJ_WGS84_HASH = {
+    :init => 'epsg:4326',
+    :proj => 'longlat',
+    :datum => 'WGS84',
+    :no_defs => true,
+    :ellps => 'WGS84',
+    :towgs84 => '0,0,0'
+  }
+
+  PROJ_GK_HASH = {
+    :init => 'epsg:31467',
+    :proj => 'tmerc',
+    :lat_0 => '0',
+    :lon_0 => '9',
+    :k => '1',
+    :x_0 => '3500000',
+    :y_0 => '0',
+    :datum => 'potsdam',
+    :units => 'm',
+    :no_defs => true,
+    :ellps => 'bessel',
+    :towgs84 => '598.1,73.7,418.2,0.202,0.045,-2.455,6.7'
+  }
+
+  PROJ_CONAKRY_HASH = {
+    :init => 'epsg:31528',
+    :proj => 'utm',
+    :zone => '28',
+    :a => '6378249.2',
+    :b => '6356515',
+    :towgs84 => '-23,259,-9,0,0,0,0',
+    :units => 'm',
+    :no_defs => true
+  }
+
+  PROJ_ORTEL_HASH = {
+    :proj => 'ortel',
+    :lon_0 => '90w',
+    :ellps => 'WGS84'
+  }
+end
+
+if RUBY_VERSION >= '1.9'
+  MiniTest::Reporters.use!(MiniTest::Reporters::SpecReporter.new)
+end
+