cb-ffi-proj4 0.0.3

data/lib/ffi-proj4.rb

@@ -0,0 +1,152 @@
+
+require 'rubygems'
+require 'ffi'
+require 'rbconfig'
+
+module Proj4
+  PROJ4_BASE = File.join(File.dirname(__FILE__), 'ffi-proj4')
+
+  autoload :Projection,
+    File.join(PROJ4_BASE, 'projection')
+  autoload :ProjXY,
+    File.join(PROJ4_BASE, 'projxy')
+  autoload :Tools,
+    File.join(PROJ4_BASE, 'tools')
+  autoload :Point,
+    File.join(PROJ4_BASE, 'point')
+
+  module FFIProj4
+    def self.proj4_library_path
+      return @proj4_library_paths if @proj4_library_paths
+
+      paths = if ENV['PROJ4_LIBRARY_PATH']
+        [ ENV['PROJ4_LIBRARY_PATH'] ]
+      else
+        [ '/usr/local/{lib64,lib}', '/opt/local/{lib64,lib}', '/usr/{lib64,lib}' ]
+      end
+
+      lib = if [
+        Config::CONFIG['arch'],
+        Config::CONFIG['host_os']
+      ].detect { |c| c =~ /darwin/ }
+        'libproj.dylib'
+      else
+        'libproj.so'
+      end
+
+      @proj4_library_path = Dir.glob(paths.collect { |path|
+          "#{path}/#{lib}"
+      }).first
+      
+      if Config::CONFIG['host_os'] =~ /win32/
+        @proj4_library_path = "#{ENV['PROJ4_LIBRARY_PATH'] + '\\' if ENV['PROJ4_LIBRARY_PATH']}proj.dll"
+      end
+      @proj4_library_path
+    end
+
+    extend ::FFI::Library
+
+    ffi_lib(*proj4_library_path)
+
+    FFI_LAYOUT = {
+      :pj_get_release => [
+        :string
+      ],
+
+      :pj_transform => [
+        :int, :pointer, :pointer, :long, :int, :pointer, :pointer, :pointer
+      ],
+
+      :pj_init_plus => [
+        :pointer, :string
+      ],
+
+      :pj_free => [
+        :void, :pointer
+      ],
+
+      :pj_is_latlong => [
+        :int, :pointer
+      ],
+
+      :pj_is_geocent => [
+        :int, :pointer
+      ],
+
+      :pj_get_def => [
+        :string, :pointer, :int
+      ],
+
+      :pj_latlong_from_proj => [
+        :pointer, :pointer
+      ],
+
+      :pj_set_finder => [
+        :void, callback([ :string ], :string)
+      ],
+
+      :pj_set_searchpath => [
+        :void, :int, :pointer
+      ],
+
+      :pj_deallocate_grids => [
+        :void
+      ],
+
+      :pj_strerrno => [
+        :string, :int
+      ],
+
+      :pj_get_errno_ref => [
+        :pointer
+      ],
+
+      :pj_fwd => [
+        Proj4::ProjXY.by_value, Proj4::ProjXY.by_value, :pointer
+      ],
+
+      :pj_inv => [
+        Proj4::ProjXY.by_value, Proj4::ProjXY.by_value, :pointer
+      ],
+
+      :pj_transform => [
+        :int, :pointer, :pointer, :long, :int, :pointer, :pointer, :pointer
+      ],
+
+      :pj_datum_transform => [
+        :int, :pointer, :pointer, :long, :int, :pointer, :pointer, :pointer
+      ]
+    }
+
+    FFI_LAYOUT.each do |fun, ary|
+      ret = ary.shift
+      begin
+        self.class_eval do
+          attach_function(fun, ary, ret)
+        end
+      rescue FFI::NotFoundError
+        # that's okay
+      end
+    end
+  end
+
+  class << self
+    def version
+      FFIProj4.pj_get_release
+    end
+  end
+
+  module Constants
+    VERSION = File.read(File.join(PROJ4_BASE, %w{ .. .. VERSION })).strip
+    PROJ4_VERSION = if Proj4.version =~ /Rel\. (\d+)\.(\d+)\.(\d+)/
+      "#{$1}#{$2}#{$3}".to_f
+    end
+
+    RAD_TO_DEG = 57.29577951308232
+    DEG_TO_RAD = 0.0174532925199432958
+  end
+
+  ENV['PROJ_LIB'] = File.join(File.dirname(PROJ4_BASE), %w{ data }) unless ENV['PROJ_LIB']
+
+  include Constants
+end

data/test/projection_tests.rb

@@ -0,0 +1,119 @@
+
+$: << File.dirname(__FILE__)
+require 'test_helper'
+
+class ProjectionTests < Test::Unit::TestCase
+  include TestHelper
+
+  def definition_sorter(definition)
+    definition.split(/\s+/).sort.join(' ')
+  end
+
+  def test_read_strings
+    tester = lambda { |expected, proj|
+      assert_equal(
+        definition_sorter(expected),
+        definition_sorter(Proj4::Projection.new(proj).definition)
+      )
+    }
+
+    tester[
+      '+init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0',
+      'init=epsg:4326'
+    ]
+
+    tester[
+      '+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',
+      'init=epsg:31467'
+    ]
+
+    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'
+    ]
+
+    tester[
+      '+proj=ortel +lon_0=90w +ellps=WGS84',
+      'proj=ortel lon_0=90w'
+    ]
+
+    assert_raise(Proj4::ProjectionParseError) do
+      tester[
+        '',
+        'gibberish'
+      ]
+    end
+  end
+
+  def test_read_arrays
+    tester = lambda { |expected, proj|
+      assert_equal(
+        definition_sorter(expected),
+        definition_sorter(Proj4::Projection.new(proj).definition)
+      )
+    }
+
+    tester[
+      '+init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0',
+      [ 'init=epsg:4326' ]
+    ]
+
+    tester[
+      '+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',
+      [ 'init=epsg:31467' ]
+    ]
+
+    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' ]
+    ]
+
+    tester[
+      '+proj=ortel +lon_0=90w +ellps=WGS84',
+      [ 'proj=ortel', 'lon_0=90w' ]
+    ]
+
+    assert_raise(Proj4::ProjectionParseError) do
+      tester[
+        '',
+        [ 'gibberish' ]
+      ]
+    end
+  end
+
+  def test_read_hashes
+    tester = lambda { |expected, proj|
+      assert_equal(
+        definition_sorter(expected),
+        definition_sorter(Proj4::Projection.new(proj).definition)
+      )
+    }
+
+    tester[
+      '+init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0',
+      { :init => 'epsg:4326' }
+    ]
+
+    tester[
+      '+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',
+      { :init => 'epsg:31467' }
+    ]
+
+    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' }
+    ]
+
+    tester[
+      '+lon_0=90w +proj=ortel +ellps=WGS84',
+      { :proj => 'ortel', :lon_0 => '90w' }
+    ]
+
+    assert_raise(Proj4::ProjectionParseError) do
+      tester[
+        '',
+        { :proj => 'gibberish' }
+      ]
+    end
+  end
+end

data/test/test_helper.rb

@@ -0,0 +1,29 @@
+
+require 'rubygems'
+require 'test/unit'
+
+if ENV['USE_BINARY_PROJ4']
+  require 'proj4_ruby'
+else
+  require File.join(File.dirname(__FILE__), %w{ .. lib ffi-proj4 })
+end
+
+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 #{ENV['PROJ_LIB']}"
+
+module TestHelper
+  TOLERANCE = 0.00000001
+
+  def self.included(base)
+    base.class_eval do
+    end
+  end
+
+  def setup
+  end
+end

data/test/transformation_tests.rb

@@ -0,0 +1,90 @@
+
+$: << File.dirname(__FILE__)
+require 'test_helper'
+
+class TransformationTests < Test::Unit::TestCase
+  include TestHelper
+
+  def setup
+    @proj_wgs84 = Proj4::Projection.new(["init=epsg:4326"])
+    @proj_gk    = Proj4::Projection.new(["init=epsg:31467"])
+    @proj_merc  = Proj4::Projection.new(["proj=merc"])
+    @lon =  8.4293092923
+    @lat = 48.9896114523
+    @rw = 3458305
+    @hw = 5428192
+    @zw = -5.1790915237
+  end
+
+  def test_gk_to_wgs84
+    to = @proj_gk.transform(@proj_wgs84, @rw, @hw, @zw)
+
+    assert_in_delta @lon, to.x * Proj4::RAD_TO_DEG, TOLERANCE
+    assert_in_delta @lat, to.y * Proj4::RAD_TO_DEG, TOLERANCE
+    assert_in_delta 0, to.z, TOLERANCE
+  end
+
+  def test_wgs84_to_gk
+    point = @proj_wgs84.transform(
+      @proj_gk,
+      @lon * Proj4::DEG_TO_RAD,
+      @lat * Proj4::DEG_TO_RAD,
+      0
+    )
+    assert_equal @rw, point.x.round
+    assert_equal @hw, point.y.round
+    assert_in_delta @zw, point.z, TOLERANCE
+  end
+
+  def test_no_dst_proj
+    assert_raise ArgumentError do
+      @proj_wgs84.transform(
+        nil,
+        @lon * Proj4::DEG_TO_RAD,
+        @lat * Proj4::DEG_TO_RAD,
+        0
+      )
+    end
+  end
+
+  def test_mercator_at_pole_raise
+    assert_raise Proj4::TransformError do
+      @proj_wgs84.transform(
+        @proj_merc,
+        0,
+        90 * Proj4::DEG_TO_RAD,
+        0
+      )
+    end
+  end
+end
+
+class SimpleTransformationTests < Test::Unit::TestCase
+  include TestHelper
+
+  def setup
+    @proj_gk = Proj4::Projection.new(["init=epsg:31467"])
+    @lon =  8.4302123334
+    @lat = 48.9906726079
+    @rw = 3458305
+    @hw = 5428192
+  end
+
+  def test_forward_gk
+    result = @proj_gk.forward_rad(@lon, @lat)
+    assert_in_delta @rw, result.x, 0.1
+    assert_in_delta @hw, result.y, 0.1
+  end
+
+  def test_inverse_gk
+    result = @proj_gk.inverse(@rw, @hw)
+    assert_in_delta @lon, Proj4::Tools.rad_to_deg(result.x), TOLERANCE
+    assert_in_delta @lat, Proj4::Tools.rad_to_deg(result.y), TOLERANCE
+  end
+
+  def test_out_of_bounds
+    assert_raise Proj4::TransformError do
+      @proj_gk.forward_deg(190, 92)
+    end
+  end
+end

metadata

@@ -0,0 +1,100 @@
+--- !ruby/object:Gem::Specification 
+name: cb-ffi-proj4
+version: !ruby/object:Gem::Version 
+  hash: 25
+  prerelease: 
+  segments: 
+  - 0
+  - 0
+  - 3
+  version: 0.0.3
+platform: ruby
+authors: 
+- J Smith
+autorequire: 
+bindir: bin
+cert_chain: []
+
+date: 2011-06-27 00:00:00 -04:00
+default_executable: 
+dependencies: []
+
+description: An ffi wrapper for the PROJ.4 Cartographic Projections library.
+email: dark.panda@gmail.com
+executables: []
+
+extensions: []
+
+extra_rdoc_files: []
+
+files: 
+- MIT-LICENSE
+- Rakefile
+- VERSION
+- 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
+- lib/ffi-proj4.rb
+- lib/ffi-proj4/point.rb
+- lib/ffi-proj4/projection.rb
+- lib/ffi-proj4/projxy.rb
+- lib/ffi-proj4/tools.rb
+- test/projection_tests.rb
+- test/test_helper.rb
+- test/transformation_tests.rb
+has_rdoc: true
+homepage: http://github.com/dark-panda/ffi-proj4
+licenses: []
+
+post_install_message: 
+rdoc_options: []
+
+require_paths: 
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement 
+  none: false
+  requirements: 
+  - - ">="
+    - !ruby/object:Gem::Version 
+      hash: 3
+      segments: 
+      - 0
+      version: "0"
+required_rubygems_version: !ruby/object:Gem::Requirement 
+  none: false
+  requirements: 
+  - - ">="
+    - !ruby/object:Gem::Version 
+      hash: 3
+      segments: 
+      - 0
+      version: "0"
+requirements: []
+
+rubyforge_project: 
+rubygems_version: 1.6.2
+signing_key: 
+specification_version: 3
+summary: An ffi wrapper for the PROJ.4 Cartographic Projections library.
+test_files: []
+