cb-ffi-proj4 0.0.3

data/data/other.extra

@@ -0,0 +1,49 @@
+## 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 <>

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=tmerx 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/lib/ffi-proj4/point.rb

@@ -0,0 +1,5 @@
+
+module Proj4
+  class Point < Struct.new(:x, :y, :z)
+  end
+end

data/lib/ffi-proj4/projection.rb

@@ -0,0 +1,189 @@
+
+module Proj4
+  class ProjectionParseError < RuntimeError; end
+  class TransformError < RuntimeError; end
+
+  class Projection
+    include Tools
+
+    attr_reader :ptr
+
+    def initialize(arg, auto_free = true)
+      args = case arg
+        when Array
+          arg.collect { |a| a.sub(/^\+/, '') }
+        when String
+          arg.strip.split(/ /).collect { |a| a.sub(/^\+/, '') }
+        when Hash
+          arg.collect { |k, v|
+            if v.nil?
+              k.to_s
+            else
+              "#{k.to_s.strip}=#{v.to_s.strip}"
+            end
+          }
+        when Proj4::Projection
+          arg.definition.strip.split(/ /).collect { |a| a.sub(/^\+/, '') }
+        else
+          raise ArgumentError.new("Unknown type #{arg.class} for projection definition")
+      end
+
+      params = args.collect(&:strip).collect { |a|
+        if !(a =~ /^\+/)
+          "+#{a}"
+        else
+          a
+        end
+      }.join(' ')
+
+      ptr = FFIProj4.pj_init_plus(params)
+
+      if ptr.null?
+        errno = FFIProj4.pj_get_errno_ref.read_int
+        raise ProjectionParseError.new(FFIProj4.pj_strerrno(errno))
+      else
+        @ptr = FFI::AutoPointer.new(
+          ptr,
+          auto_free ? self.class.method(:release) : self.class.method(:no_release)
+        )
+      end
+    end
+
+    def self.no_release(ptr) #:nodoc:
+    end
+
+    def self.release(ptr) #:nodoc:
+      FFIProj4.pj_free(ptr)
+    end
+
+    def lat_long?
+      bool_result(FFIProj4.pj_is_latlong(self.ptr))
+    end
+
+    def geocentric?
+      bool_result(FFIProj4.pj_is_geocent(self.ptr))
+    end
+
+    def definition
+      @definition ||= FFIProj4.pj_get_def(self.ptr, 0).strip
+    end
+
+    def definition_as_hash
+      self.definition.split(/ /).inject({}) { |memo, opt|
+        memo.tap {
+          k, v = opt.split(/=/)
+          k.sub!(/^\+/, '')
+          v = true if v.nil?
+          memo[k.to_sym] = v
+        }
+      }
+    end
+
+    #def to_s
+      #"#<Proj4::Projection #{definition}>"
+    #end
+    #alias :inspect :to_s
+
+    def forward(x, y)
+      xy = ProjXY.new(x, y)
+      ret = FFIProj4.pj_fwd(xy, self.ptr)
+      errno = FFIProj4.pj_get_errno_ref.read_int
+      if errno == 0
+        Point.new(ret.x, ret.y)
+      else
+        raise TransformError.new(FFIProj4.pj_strerrno(errno))
+      end
+    end
+    alias :forward_deg :forward
+
+    def forward_rad(x, y)
+      self.forward(deg_to_rad(x), deg_to_rad(y))
+    end
+
+    def inverse(x, y)
+      xy = ProjXY.new(x, y)
+      ret = FFIProj4.pj_inv(xy, self.ptr)
+      errno = FFIProj4.pj_get_errno_ref.read_int
+      if errno == 0
+        Point.new(ret.x, ret.y)
+      else
+        raise TransformError.new(FFIProj4.pj_strerrno(errno))
+      end
+    end
+    alias :inverse_deg :inverse
+
+    def inverse_rad(x, y)
+      self.inverse(deg_to_rad(x), deg_to_rad(y))
+    end
+
+    def transform(proj, x, y, z = nil)
+      if !proj.is_a?(Proj4::Projection)
+        raise ArgumentError.new("Expected a Proj4::Projection")
+      end
+
+      x_ptr = FFI::MemoryPointer.new(:double)
+      y_ptr = FFI::MemoryPointer.new(:double)
+      z_ptr = FFI::MemoryPointer.new(:double) unless z.nil?
+
+      x_ptr.write_double(x)
+      y_ptr.write_double(y)
+      z_ptr.write_double(z) unless z.nil?
+
+      result = FFIProj4.pj_transform(self.ptr, proj.ptr, 1, 1, x_ptr, y_ptr, z_ptr)
+
+      if result >= 0 && !bool_result(result)
+        Point.new(
+          x_ptr.read_double,
+          y_ptr.read_double,
+          z_ptr.read_double
+        )
+      else
+        raise TransformError.new(FFIProj4.pj_strerrno(result))
+      end
+    end
+    alias :transform_rad :transform
+
+    def transform_deg(proj, x, y, z = nil)
+      self.transform(proj, x, y, z).tap { |ret|
+        ret.x = rad_to_deg(ret.x)
+        ret.y = rad_to_deg(ret.y)
+        ret.z = rad_to_deg(ret.z)
+      }
+    end
+
+    def datum_transform(proj, x, y, z = nil)
+      if !proj.is_a?(Proj4::Projection)
+        raise ArgumentError.new("Expected a Proj4::Projection")
+      end
+
+      x_ptr = FFI::MemoryPointer.new(:double)
+      y_ptr = FFI::MemoryPointer.new(:double)
+      z_ptr = FFI::MemoryPointer.new(:double) unless z.nil?
+
+      x_ptr.write_double(x)
+      y_ptr.write_double(y)
+      z_ptr.write_double(z) unless z.nil?
+
+      result = FFIProj4.pj_transform(self.ptr, proj.ptr, 1, 1, x_ptr, y_ptr, z_ptr)
+
+      if result >= 0 && !bool_result(result)
+        Point.new(
+          x_ptr.read_double,
+          y_ptr.read_double,
+          z_ptr.read_double
+        )
+      else
+        raise TransformError.new(FFIProj4.pj_strerrno(result))
+      end
+    end
+    alias :datum_transform_rad :datum_transform
+
+    def datum_transform_deg(proj, x, y, z = nil)
+      self.datum_transform(proj, x, y, z).tap { |ret|
+        ret.x = rad_to_deg(ret.x)
+        ret.y = rad_to_deg(ret.y)
+        ret.z = rad_to_deg(ret.z)
+      }
+    end
+  end
+end

data/lib/ffi-proj4/projxy.rb

@@ -0,0 +1,43 @@
+
+module Proj4
+  class ProjXY < FFI::Struct
+    layout(
+      :x, :double,
+      :y, :double
+    )
+
+    def initialize(*args)
+      if args.first.is_a?(FFI::Pointer)
+        super(*args)
+      else
+        self[:x] = if args.first.respond_to?(:read_double)
+          args.first.read_double
+        else
+          args.first.to_f
+        end
+
+        self[:y] = if args[1].respond_to?(:read_double)
+          args.first.read_double
+        else
+          args[1].to_f
+        end
+      end
+    end
+
+    def x=(v)
+      self[:x] = v
+    end
+
+    def y=(v)
+      self[:y] = v
+    end
+
+    def x
+      self[:x]
+    end
+
+    def y
+      self[:y]
+    end
+  end
+end

data/lib/ffi-proj4/tools.rb

@@ -0,0 +1,27 @@
+
+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)
+      rad * Proj4::RAD_TO_DEG
+    end
+
+    def deg_to_rad(deg)
+      deg * Proj4::DEG_TO_RAD
+    end
+
+    class << self
+      include Proj4::Tools
+    end
+  end
+end