ffi-proj4 0.1.0

data/data/other.extra

@@ -0,0 +1,53 @@
+## NAD83 / BC Albers (this has been superceeded but is kept for compatibility)
+<42102> +proj=aea +ellps=GRS80 +lat_0=45 +lon_0=-126.0 +lat_1=50.0 +lat_2=58.5 +x_0=1000000.0 +y_0=0 +datum=NAD83 +units=m no_defs <>
+
+
+#
+#  OGC-defined extended codes  (41000--41999)
+#  see http://www.digitalearth.gov/wmt/auto.html
+#
+# WGS84 / Simple Mercator
+<41001> +proj=merc +lat_ts=0 +lon_0=0 +k=1.000000 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs  no_defs <>
+#
+#  CubeWerx-defined extended codes  (42100--42199)
+#
+# WGS 84 / LCC Canada
+<42101> +proj=lcc +lat_1=49 +lat_2=77 +lat_0=0 +lon_0=-95 +x_0=0 +y_0=-8000000 +ellps=WGS84 +datum=WGS84 +units=m +no_defs  no_defs <>
+#EPSG:42102,"PROJCS[\"NAD83 / BC Albers\",GEOGCS[\"NAD83\",DATUM[\"North_American_Datum_1983\",SPHEROID[\"GRS_1980\",6378137,298.257222101]],PRIMEM[\"Greenwich\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Albers_conic_equal_area\"],PARAMETER[\"central_meridian\",-126.0],PARAMETER[\"latitude_of_origin\",45],PARAMETER[\"standard_parallel_1\",50.0],PARAMETER[\"standard_parallel_2\",58.5],PARAMETER[\"false_easting\",1000000.0],PARAMETER[\"false_northing\",0],UNIT[\"Meter\",1]]"
+# WGS 84 / LCC USA
+<42103> +proj=lcc +lat_1=33 +lat_2=45 +lat_0=0 +lon_0=-100 +x_0=0 +y_0=0 +ellps=WGS72 +datum=WGS84 +units=m +no_defs  no_defs <>
+# NAD83 / MTM zone 8 Québec
+<42104> +proj=tmerc +lat_0=0 +lon_0=-73.5 +k=0.999900 +x_0=304800 +y_0=0 +ellps=GRS80 +units=m +no_defs  no_defs <>
+# WGS84 / Merc NorthAm
+<42105> +proj=merc +lat_ts=0 +lon_0=-96 +k=1.000000 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs  no_defs <>
+# WGS84 / Lambert Azim Mozambique
+<42106> +proj=laea +lat_0=5 +lon_0=20 +x_0=0 +y_0=0 +a=6370997 +b=6370997 +datum=WGS84 +units=m +no_defs  no_defs <>
+#
+#  CubeWerx-customer definitions  (42300--42399)
+#
+# NAD27 / Polar Stereographic / CM=-98
+<42301> +proj=stere +lat_0=90 +lon_0=-98 +x_0=0 +y_0=0 +ellps=clrk66 +datum=NAD27 +units=m +no_defs  no_defs <>
+# JapanOrtho.09 09
+<42302> +proj=tmerc +lat_0=36 +lon_0=139.833333333333 +k=0.999900 +x_0=0 +y_0=0 +ellps=bessel +units=m +no_defs  no_defs <>
+# NAD83 / Albers NorthAm
+<42303> +proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs  no_defs <>
+# NAD83 / NRCan LCC Canada
+<42304> +proj=lcc +lat_1=49 +lat_2=77 +lat_0=49 +lon_0=-95 +x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs  no_defs <>
+# France_II
+<42305> +proj=lcc +lat_1=45.898918964419 +lat_2=47.696014502038 +lat_0=46.8 +lon_0=2.337229166666667 +x_0=600000 +y_0=2200000 +a=6378249.2 +b=6356514.999904194 +pm=2.337229166666667 +units=m +no_defs  no_defs <>
+# NAD83/QC_LCC
+<42306> +proj=lcc +lat_1=46 +lat_2=60 +lat_0=44 +lon_0=-68.5 +x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs  no_defs <>
+# NAD83 / Texas Central - feet
+<42307> +proj=lcc +lat_1=31.8833333333333 +lat_2=30.1166666666667 +lat_0=29.6666666666667 +lon_0=-100.333333333333 +x_0=700000.0000000001 +y_0=3000000 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192 +no_defs  no_defs <>
+# NAD27 / California Albers
+<42308> +proj=aea +lat_1=34 +lat_2=40.5 +lat_0=0 +lon_0=-120 +x_0=0 +y_0=-4000000 +ellps=clrk66 +datum=NAD27 +units=m +no_defs  no_defs <>
+# NAD 83 / LCC Canada AVHRR-2
+<42309> +proj=lcc +lat_1=49 +lat_2=77 +lat_0=0 +lon_0=-95 +x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs  no_defs <>
+# WGS84+GRS80 / Mercator
+<42310> +proj=merc +lat_ts=0 +lon_0=0 +k=1.000000 +x_0=0 +y_0=0 +ellps=GRS80 +datum=WGS84 +units=m +no_defs  no_defs <>
+# NAD83 / LCC Statcan
+<42311> +proj=lcc +lat_1=49 +lat_2=77 +lat_0=63.390675 +lon_0=-91.86666700000001 +x_0=6200000 +y_0=3000000 +ellps=GRS80 +datum=NAD83 +units=m +no_defs  no_defs <>
+#
+#  Funny epsgish code for google mercator - you should really use EPSG:3857
+#
+<900913> +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs <>

data/data/proj_def.dat

@@ -0,0 +1,17 @@
+# Projection library defaults file
+# SCCSID--- @(#)proj_def.dat	4.3	94/02/23	GIE	REL"
+# very preliminary version
+<general>
+ellps=WGS84 
+<>
+<aea> # Conterminous U.S. map
+lat_1=29.5
+lat_2=45.5
+<>
+<lcc> # Conterminous U.S. map
+lat_1=33
+lat_2=45
+<>
+<lagrng>
+W=2
+<>

data/data/world

@@ -0,0 +1,212 @@
+# SCCSID @(#)world	1.2 95/08/05 GIE REL
+# proj +init files for various non-U.S. coordinate systems.
+#
+<CH1903> # Swiss Coordinate System
+	+proj=somerc +lat_0=46d57'8.660"N +lon_0=7d26'22.500"E
+	+ellps=bessel +x_0=600000 +y_0=200000
+	+k_0=1.  no_defs <>
+<madagascar> # Laborde grid for Madagascar
+	proj=labrd ellps=intl lon_0=46d26'13.95E lat_0=18d54S
+	azi=18d54 k_0=.9995 x_0=400000 y_0=800000
+	no_defs <>
+<new_zealand> # New Zealand Map Grid (NZMG)
+	proj=nzmg  # Projection unique to N.Z. so all factors fixed
+	no_defs <>
+# Secondary grids DMA TM8358.1, p. 4.3
+<bwi> # British West Indies
+	proj=tmerc ellps=clrk80 lon_0=62W
+	x_0=400000 k_0=0.9995
+	no_defs <>
+<costa-n> # Costa Rica Norte
+	proj=lcc ellps=clrk66 lat_1=10d28N lon_0=84d20W
+	x_0=500000 y_0=217820.522 k_0=0.99995696
+	no_defs <>
+<costa-s> # Costa Rica Sud
+	proj=lcc ellps=clrk66 lat_1=9dN lon_0=83d40W
+	x_0=500000 y_0=327987.436 k_0=0.99995696
+	no_defs <>
+<cuba-n> # Cuba Norte
+	proj=lcc ellps=clrk66 lat_1=22d21N lon_0=81dW
+	x_0=500000 y_0=280296.016 k_0=0.99993602
+	no_defs <>
+<cuba-s> # Cuba Sud
+	proj=lcc ellps=clrk66 lat_1=20d43'N lon_0=76d50'W
+	x_0=500000 y_0=229126.939 k_0=0.99994848
+	no_defs <>
+<domin_rep> # Dominican Republic
+	proj=lcc ellps=clrk66 lat_1=18d49'N lon_0=71d30'W
+	x_0=500000 y_0=277063.657 k_0=0.99991102
+	no_defs <>
+<egypt-1> # Egypt
+	proj=tmerc ellps=intl lon_0=25d30'E x_0=300000 k_0=0.99985
+	no_defs <>
+<egypt-2> # Egypt
+	proj=tmerc ellps=intl lon_0=28d30'E x_0=300000 k_0=0.99985
+	no_defs <>
+<egypt-3> # Egypt
+	proj=tmerc ellps=intl lon_0=31d30'E x_0=300000 k_0=0.99985
+	no_defs <>
+<egypt-4> # Egypt
+	proj=tmerc ellps=intl lon_0=34d30'E x_0=300000 k_0=0.99985
+	no_defs <>
+<egypt-5> # Egypt
+	proj=tmerc ellps=intl lon_0=37d30'E x_0=300000 k_0=0.99985
+	no_defs <>
+<el_sal> # El Salvador
+	proj=lcc ellps=clrk66 lat_1=13d47'N lon_0=89dW
+	x_0=500000 y_0=295809.184 k_0=0.99996704
+	no_defs <>
+<guat-n> # Guatemala Norte
+	proj=lcc ellps=clrk66 lat_1=16d49'N lon_0=90d20'W
+	x_0=500000 y_0=292209.579 k_0=0.99992226
+	no_defs <>
+<guat-s> # Guatemala Sud
+	proj=lcc ellps=clrk66 lat_1=14d54'N lon_0=90d20'W
+	x_0=500000 y_0=325992.681 k_0=0.99989906
+	no_defs <>
+<haiti> # Haiti
+	proj=lcc ellps=clrk66 lat_1=18d49'N lon_0=71d30'W
+	x_0=500000 y_0=277063.657 k_0=0.99991102
+	no_defs <>
+<hond-n> # Honduras Norte
+	proj=lcc ellps=clrk66 lat_1=15d30'N lon_0=86d10'W
+	x_0=500000 y_0=296917.439 k_0=0.99993273
+	no_defs <>
+<hond-s> # Honduras Sud
+	proj=lcc ellps=clrk66 lat_1=13d47'N lon_0=87d10'W
+	x_0=500000 y_0=296215.903 k_0=0.99995140
+	no_defs <>
+<levant> # Levant
+	proj=lcc ellps=clrk66 lat_1=34d39'N lon_0=37d21'E
+	x_0=500000 y_0=300000 k_0=0.9996256
+	no_defs <>
+<nica-n> # Nicaragua Norte
+	proj=lcc ellps=clrk66 lat_1=13d52'N lon_0=85d30'W
+	x_0=500000 y_0=359891.816 k_0=0.99990314
+	no_defs <>
+<nica-s> # Nicaragua Sud
+	proj=lcc ellps=clrk66 lat_1=11d40'N lon_0=85d30'W
+	x_0=500000 y_0=288876.327 k_0=0.99992228
+	no_defs <>
+<nw-africa> # Northwest Africa
+	proj=lcc ellps=clrk80 lat_1=34dN lon_0=0dE
+	x_0=1000000 y_0=500000 k_0=0.99908
+	no_defs <>
+<palestine> # Palestine
+	proj=tmerc a=6378300.79 rf=293.488307656
+	lat_0=31d44'2.749"N lon_0=35d12'43.490"E
+	x_0=170251.555 y_0=126867.909 k_0=1
+	no_defs <>
+<panama> # Panama
+	proj=lcc ellps=clrk66 lat_1=8d25'N lon_0=80dW
+	x_0=500000 y_0=294865.303 k_0=0.99989909
+	no_defs <>
+# other grids in DMA TM8358.1
+<bng> # British National Grid
+	proj=tmerc ellps=airy lat_0=49dN lon_0=2dW
+	k_0=0.9996012717 x_0=400000 y_0=-100000
+	no_defs <>
+<malay> # West Malaysian RSO Grid
+	proj=omerc a=6377295.66402 rf=300.8017 alpha=323d01'32.846"
+	no_uoff rot_conv lonc=102d15E lat_0=4dN k_0=0.99984 x_0=804670.240 y_0=0
+	no_defs <>
+<india-I> # India Zone I
+	proj=lcc ellps=everest lon_0=68E lat_1=32d30'N
+	x_0=2743185.69 y_0=914395.23 k_0=.998786408
+	no_defs <>
+<india-IIA> # India Zone IIA
+	proj=lcc ellps=everest lon_0=74E lat_1=26N
+	x_0=2743185.69 y_0=914395.23 k_0=.998786408
+	no_defs <>
+<india-IIB> # India Zone IIB
+	proj=lcc ellps=everest lon_0=90E lat_1=26N
+	x_0=2743185.69 y_0=914395.23 k_0=.998786408
+	no_defs <>
+<india-IIIA> # India Zone IIIA
+	proj=lcc ellps=everest lon_0=80E lat_1=19N
+	x_0=2743185.69 y_0=914395.23 k_0=.998786408
+	no_defs <>
+<india-IIIB> # India Zone IIIB
+	proj=lcc ellps=everest lon_0=100E lat_1=19N
+	x_0=2743185.69 y_0=914395.23 k_0=.998786408
+	no_defs <>
+<india-IVA> # India Zone IVA
+	proj=lcc ellps=everest lon_0=80E lat_1=12N
+	x_0=2743185.69 y_0=914395.23 k_0=.998786408
+	no_defs <>
+<india-IVB> # India Zone IVB
+	proj=lcc ellps=everest lon_0=104E lat_1=12N
+	x_0=2743185.69 y_0=914395.23 k_0=.998786408
+	no_defs <>
+<ceylon> # Ceylon Belt
+	proj=tmerc ellps=everest lon_0=80d46'18.160"E lat_0=7d0'1.729"N
+	x_0=160933.56048 y_0=160933.56048 k_0=1.
+	no_defs <>
+<irish> # Irish Transverse Mercator Grid
+	proj=tmerc ellps=mod_airy lat_0=53d30'N lon_0=8W
+	x_0=200000 y_0=250000 k_0=1.000035
+	no_defs <>
+<neiez> # Netherlands East Indies Equitorial Zone
+	proj=merc ellps=bessel lon_0=110E
+	x_0=3900000 y_0=900000 k_0=0.997
+	no_defs <>
+<n-alger> # Nord Algerie Grid
+	proj=lcc ellps=clrk80 lon_0=2d42E lat_0=36N
+	x_0=500000 y_0=300000 k_0=0.999625544
+	no_defs <>
+<n-maroc> # Nord Maroc Grid
+	proj=lcc ellps=clrk80 lon_0=5d24'W lat_0=33d18'N
+	x_0=500000 y_0=300000 k_0=0.999625769
+	no_defs <>
+<n-tunis> # Nord Tunisie Grid
+	proj=lcc ellps=clrk80 lon_0=9d54E lat_0=36N
+	x_0=500000 y_0=300000 k_0=0.999625544
+	no_defs <>
+<s-alger> # Sud Algerie Grid
+	proj=lcc ellps=clrk80 lon_0=2d42E lat_0=33d18'N
+	x_0=500000 y_0=300000 k_0=0.999625769
+	no_defs <>
+<s-maroc> # Sud Maroc Grid
+	proj=lcc ellps=clrk80 lon_0=5d24W lat_0=29d42'N
+	x_0=500000 y_0=300000 k_0=0.999615596
+	no_defs <>
+<s-tunis> # Sud Tunisie Grid
+	proj=lcc ellps=clrk80 lon_0=9d54'E lat_0=33d18'N
+	x_0=500000 y_0=300000 k_0=0.999625769
+	no_defs <>
+# Gauss Krueger Grid for Germany
+# 
+# The first figure of the easting is lon_0 devided by 3
+# ( 2 for 6d0E, 3 for 9d0E, 4 for 12d0E)
+# For translations you have to remove this first figure
+# and convert northings and eastings from km to meter .
+# The other way round, devide by 1000 and add the figure.
+# I made 3 entrys for the officially used grids in Germany
+# 
+#
+# Und nochmal in deutsch :
+# Die erste Ziffer des Rechtswerts beschreibt den Hauptmeridian
+# und ist dessen Gradzahl geteilt durch 3.
+# Zum Umrechnen in Grad muss daher die erste Ziffer des Rechtswertes
+# entfernt werden und evt. von km auf Metern umgerechnet werden.
+# Zur Umrechnung in Gauss Krueger Koordinaten entsprechend die
+# Ziffer fuer den Hauptmeridian vor dem Rechtswert ergaenzen.
+# Ich hab fuer alle drei in Deutschland ueblichen Hauptmeridiane
+# jeweils einen Eintrag ergaenzt.   
+#
+#
+# added by Michael Goepel <goepel@ffm.plusline.de>
+#
+<gk2-d> # Gauss Krueger Grid for Germany
+        proj=tmerc ellps=bessel lon_0=6d0E lat_0=0
+        x_0=500000
+        no_defs<>
+<gk3-d> # Gauss Krueger Grid for Germany
+        proj=tmerc ellps=bessel lon_0=9d0E lat_0=0
+        x_0=500000
+        no_defs<>
+<gk4-d> # Gauss Krueger Grid for Germany
+        proj=tmerc ellps=bessel lon_0=12d0E lat_0=0
+        x_0=500000
+        no_defs<>
+

data/ffi-proj4.gemspec

@@ -0,0 +1,23 @@
+# -*- encoding: utf-8 -*-
+
+require File.expand_path('../lib/ffi-proj4/version', __FILE__)
+
+Gem::Specification.new do |s|
+  s.name = "ffi-proj4"
+  s.version = Proj4::VERSION
+
+  s.required_rubygems_version = Gem::Requirement.new("> 1.3.1") if s.respond_to? :required_rubygems_version=
+  s.authors = ["J Smith"]
+  s.description = "An ffi wrapper for the PROJ.4 Cartographic Projections library."
+  s.summary = s.description
+  s.email = "dark.panda@gmail.com"
+  s.license = "MIT"
+  s.files = `git ls-files`.split($\)
+  s.executables = s.files.grep(%r{^bin/}).map { |f| File.basename(f) }
+  s.test_files = s.files.grep(%r{^(test|spec|features)/})
+  s.homepage = "http://github.com/dark-panda/ffi-proj4"
+  s.require_paths = ["lib"]
+
+  s.add_dependency("ffi", [">= 1.0.0"])
+end
+

data/lib/ffi-proj4.rb

@@ -0,0 +1,164 @@
+
+require 'rubygems'
+require 'ffi'
+require 'rbconfig'
+require 'ffi-proj4/version'
+require 'ffi-proj4/error'
+
+#ENV['PROJ_LIB'] = File.join(File.dirname(__FILE__), %w{ .. data }) unless ENV['PROJ_LIB']
+#p ENV['PROJ_LIB']
+
+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_path if defined?(@proj4_library_path)
+
+      lib = if FFI::Platform::IS_WINDOWS
+        # For MinGW and the official binaries
+        '{libproj-?,proj}.dll'
+      else
+        "libproj.#{FFI::Platform::LIBSUFFIX}"
+      end
+
+      paths = if ENV['PROJ4_LIBRARY_PATH']
+        [ ENV['PROJ4_LIBRARY_PATH'] ]
+      elsif FFI::Platform::IS_WINDOWS
+        ENV['PATH'].split(File::PATH_SEPARATOR)
+      else
+        [ '/usr/local/{lib64,lib}', '/opt/local/{lib64,lib}', '/usr/{lib64,lib}', '/usr/lib/{x86_64,i386}-linux-gnu' ]
+      end
+
+      @proj4_library_path = Dir.glob(paths.collect { |path|
+        File.expand_path(File.join(path, lib))
+      }).first
+    end
+
+    extend ::FFI::Library
+
+    ffi_lib(*proj4_library_path)
+
+    FFI_LAYOUT = {
+      :pj_get_release => [
+        :string
+      ],
+
+      :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
+      ],
+
+      :setenv => [
+        :int, :string, :string, :int
+      ]
+    }
+
+    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
+    attr_reader :proj_lib
+
+    def version
+      FFIProj4.pj_get_release
+    end
+
+    def proj_lib=(lib)
+      @proj_lib = lib
+      if RUBY_PLATFORM == 'java' && FFIProj4.respond_to?(:setenv)
+        FFIProj4.setenv('PROJ_LIB', lib, 1)
+      else
+        ENV['PROJ_LIB'] = lib
+      end
+    end
+  end
+
+  module Constants
+    PROJ4_VERSION = if Proj4.version =~ /Rel\. (\d+)\.(\d+)\.(\d+)/
+      "#{$1}#{$2}#{$3}".to_f
+    end
+
+    LIBVERSION = PROJ4_VERSION
+
+    RAD_TO_DEG = 57.29577951308232
+    DEG_TO_RAD = 0.0174532925199432958
+  end
+
+  include Constants
+
+  self.proj_lib = ENV['PROJ_LIB'] || File.join(File.dirname(__FILE__), %w{ .. data })
+end