proj4rb 0.3.0

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/example/basic.rb

@@ -0,0 +1,18 @@
+#!/usr/bin/ruby -w
+
+$: << File.dirname(__FILE__) + '/../lib/' << File.dirname(__FILE__) + '/../ext/'
+
+require 'proj4'
+include Proj4
+
+# create projection object
+proj = Projection.new(["proj=utm", "zone=21", "units=m", "no.defs"])
+
+# create the point you want to project
+point = Point.new(45000,2700000)
+
+# do the projection
+res = proj.inverse(point)
+
+puts "x=#{point.x}, y=#{point.y} ==> lon=#{res.lon}, lat=#{res.lat}"
+

data/example/list-datums.rb

@@ -0,0 +1,17 @@
+#!/usr/bin/ruby -w
+
+$: << File.dirname(__FILE__) + '/../lib/' << File.dirname(__FILE__) + '/../ext/'
+
+require 'proj4'
+include Proj4
+
+if ARGV[0] == '-v'
+    puts "%-15s %-10s %-50s" % [ 'ID', 'Ellipse ID', 'Comment' ]
+    puts '-----------------------------------------------------------------------------'
+    Datum.list.sort.each do |datum|
+        puts "%-15s %-10s %-50s\n  %s\n " % [ datum.id, datum.ellipse_id, datum.comments, datum.defn ]
+    end
+else
+    puts Datum.list.sort
+end
+

data/example/list-ellipsoids.rb

@@ -0,0 +1,17 @@
+#!/usr/bin/ruby -w
+
+$: << File.dirname(__FILE__) + '/../lib/' << File.dirname(__FILE__) + '/../ext/'
+
+require 'proj4'
+include Proj4
+
+if ARGV[0] == '-v'
+    puts "%-9s %-14s %-18s %-36s" % [ 'ID', 'a', 'b/rf', 'Name' ]
+    puts '-----------------------------------------------------------------------------'
+    Ellipsoid.list.sort.each do |ell|
+        puts "%-9s %-14s %-18s %-36s" % [ ell.id, ell.major, ell.ell, ell.name ]
+    end
+else
+    puts Ellipsoid.list.sort
+end
+

data/example/list-errors.rb

@@ -0,0 +1,11 @@
+#!/usr/bin/ruby -w
+
+$: << File.dirname(__FILE__) + '/../lib/' << File.dirname(__FILE__) + '/../ext/'
+
+require 'proj4'
+include Proj4
+
+Error.list.collect{ |e| e + 'Error' }.each_with_index do |error, index|
+    puts "%2d %s" % [index, error]
+end
+

data/example/list-prime-meridians.rb

@@ -0,0 +1,17 @@
+#!/usr/bin/ruby -w
+
+$: << File.dirname(__FILE__) + '/../lib/' << File.dirname(__FILE__) + '/../ext/'
+
+require 'proj4'
+include Proj4
+
+if ARGV[0] == '-v'
+    puts "%-16s %-50s" % [ 'ID', 'Definition' ]
+    puts '-----------------------------------------------------------------------------'
+    PrimeMeridian.list.sort.each do |pm|
+        puts "%-16s %-50s" % [ pm.id, pm.defn ]
+    end
+else
+    puts PrimeMeridian.list.sort
+end
+

data/example/list-projection-types.rb

@@ -0,0 +1,17 @@
+#!/usr/bin/ruby -w
+
+$: << File.dirname(__FILE__) + '/../lib/' << File.dirname(__FILE__) + '/../ext/'
+
+require 'proj4'
+include Proj4
+
+if ARGV[0] == '-v'
+    puts "%-7s %-70s" % [ 'ID', 'Name/Description' ]
+    puts '-----------------------------------------------------------------------------'
+    ProjectionType.list.sort.each do |pt|
+        puts "%-7s %-70s" % [ pt.id, pt.descr ]
+    end
+else
+    puts ProjectionType.list
+end
+

data/example/list-units.rb

@@ -0,0 +1,17 @@
+#!/usr/bin/ruby -w
+
+$: << File.dirname(__FILE__) + '/../lib/' << File.dirname(__FILE__) + '/../ext/'
+
+require 'proj4'
+include Proj4
+
+if ARGV[0] == '-v'
+    puts "%-8s %-20s %-50s" % [ 'ID', 'To meter', 'Name' ]
+    puts "------------------------------------------------------------------------------"
+    Unit.list.sort.each do |unit|
+        puts "%-8s %-20s %-50s" % [ unit.id, unit.to_meter, unit.name ]
+    end
+else
+    puts Unit.list.sort
+end
+

data/example/version.rb

@@ -0,0 +1,8 @@
+#!/usr/bin/ruby -w
+
+$: << File.dirname(__FILE__) + '/../lib/' << File.dirname(__FILE__) + '/../ext/'
+
+require 'proj4'
+
+puts Proj4::LIBVERSION
+

data/ext/extconf.rb

@@ -0,0 +1,8 @@
+require 'mkmf'
+
+# help with compiling on Mac OS X
+$CFLAGS += " -I/sw/include"
+$LDFLAGS += " -L/sw/lib"
+
+have_library 'proj', 'pj_init'
+create_makefile 'proj4_ruby'

data/ext/mingw/rakefile.rb

@@ -0,0 +1,38 @@
+# We can't use Ruby's standard build procedures
+# on Windows because the Ruby executable is
+# built with VC++ while here we want to build
+# with MingW.  So just roll our own...
+
+require 'rake/clean'
+require 'rbconfig'
+
+RUBY_INCLUDE_DIR = Config::CONFIG["archdir"]
+RUBY_BIN_DIR = Config::CONFIG["bindir"]
+RUBY_LIB_DIR = Config::CONFIG["libdir"]
+RUBY_SHARED_LIB = Config::CONFIG["LIBRUBY"]
+RUBY_SHARED_DLL = RUBY_SHARED_LIB.gsub(/lib$/, 'dll')
+
+EXTENSION_NAME = "proj4_ruby.#{Config::CONFIG["DLEXT"]}"
+    
+CLEAN.include('*.o')
+CLOBBER.include(EXTENSION_NAME)
+
+task :default => "projrb"
+
+SRC = FileList['../*.c']
+OBJ = SRC.collect do |file_name|
+  File.basename(file_name).ext('o')
+end
+
+SRC.each do |srcfile|
+  objfile = File.basename(srcfile).ext('o')
+  file objfile => srcfile do
+    command = "gcc -c -fPIC -O2 -Wall -mms-bitfields -o #{objfile} #{srcfile} -I#{RUBY_INCLUDE_DIR} -I/local/include"
+    sh "sh -c '#{command}'" 
+  end
+end
+
+file "projrb" => OBJ do
+  command = "gcc -shared -o #{EXTENSION_NAME} #{OBJ} -lproj #{RUBY_BIN_DIR}/#{RUBY_SHARED_DLL} -L/local/lib" 
+  sh "sh -c '#{command}'" 
+end