coordinate-converter 0.1.0

data/README.md

@@ -0,0 +1,42 @@
+Converting UTM to Latitude and Longitude
+========================================
+
+Parameters are specified in the following order:
+
+* reference_ellipsoid
+* northing
+* easting
+* zone
+
+Example
+-------
+
+    Coordinates.utm_to_lat_long("WGS-84", 6688940, 219165, "33N")
+    # => {:lat => 60.2399303597183, :long => 9.92496237547609}
+
+Supported ellipsoids
+--------------------
+
+* Airy
+* Australian National
+* Bessel 1841
+* Bessel 1841 (Nambia)
+* Clarke 1866
+* Clarke 1880
+* Everest
+* Fischer 1960 (Mercury)
+* Fischer 1968
+* GRS 1967
+* GRS 1980
+* Helmert 1906
+* Hough
+* International
+* Krassovsky
+* Modified Airy
+* Modified Everest
+* Modified Fischer 1960
+* South American 1969
+* WGS 60
+* WGS 66
+* WGS-72
+* WGS-84

data/Rakefile

@@ -0,0 +1,14 @@
+begin
+  require 'jeweler'
+  Jeweler::Tasks.new do |gemspec|
+    gemspec.name        = "coordinate-converter"
+    gemspec.summary     = "Converting UTM to Latitude and Longitude"
+    gemspec.description = "Converting UTM to Latitude and Longitude"
+    gemspec.email       = "mail@kimjoar.net"
+    gemspec.homepage    = "http://github.com/kjbekkelund/coordinate-converter"
+    gemspec.authors     = "Kim Joar Bekkelund"
+  end
+  Jeweler::GemcutterTasks.new
+rescue LoadError
+  puts "Jeweler not available. Install it with: gem install jeweler"
+end

data/VERSION

@@ -0,0 +1 @@
+0.1.0

data/lib/coordinate-converter.rb

@@ -0,0 +1,88 @@
+module Coordinates
+  def self.utm_to_lat_long(reference_ellipsoid, northing, easting, zone)    
+    a  = equitorial_radius(reference_ellipsoid)
+    ecc_squared = eccentricity_squared(reference_ellipsoid)
+    
+    e1 = (1 - Math.sqrt(1 - ecc_squared)) / (1 + Math.sqrt(1 - ecc_squared))
+    
+    x = easting - 500000.0 
+    y = northing
+    
+    zone_letter = zone[-1].chr
+    zone_number = zone[0..-2].to_i
+    if zone_letter >= 'S'
+      y -= 10000000.0   
+    end
+    
+    long_origin = (zone_number - 1) * 6 - 180 + 3
+    ecc_prime_squared = (ecc_squared) / (1 - ecc_squared)
+    
+    k0 = 0.9996
+    m  = y / k0
+    mu = m / (a * (1 - ecc_squared / 4 - 3 * ecc_squared * ecc_squared / 64 - 5 \
+       * ecc_squared * ecc_squared * ecc_squared / 256))
+    
+    phi1_rad = (mu + (3 * e1 / 2 - 27 * e1 * e1 * e1 / 32) * Math.sin(2 * mu) \
+                   + (21 * e1 * e1 / 16 - 55 * e1 * e1 * e1 * e1 / 32) \
+                   * Math.sin(4 * mu) + (151 * e1 * e1 * e1 / 96) \
+                   * Math.sin(6 * mu))
+    
+    n1 = a / Math.sqrt(1 - ecc_squared * Math.sin(phi1_rad) * Math.sin(phi1_rad))
+    t1 = Math.tan(phi1_rad) * Math.tan(phi1_rad)
+    c1 = ecc_prime_squared * Math.cos(phi1_rad) * Math.cos(phi1_rad)
+    r1 = a * (1 - ecc_squared) / (1 - ecc_squared * Math.sin(phi1_rad) * Math.sin(phi1_rad))**1.5
+    d  = x / (n1 * k0)
+
+    lat = phi1_rad - (n1 * Math.tan(phi1_rad) / r1) * (d * d / 2 - (5 + 3 * t1 \
+        + 10 * c1 - 4 * c1**2 - 9 * ecc_prime_squared) * d**4 / 24 + (61 + 90 \
+        * t1 + 298 * c1 + 45 * t1**2 - 252 * ecc_prime_squared - 3 * c1**2) \
+        * d**6 / 720)
+    lat = lat * RAD_TO_DEG
+
+    long = (d - (1 + 2 * t1 + c1) * d * d * d / 6 + (5 - 2 * c1 + 28 \
+         * t1 - 3 * c1 * c1 + 8 * ecc_prime_squared + 24 * t1 * t1) \
+         * d * d * d * d * d / 120) / Math.cos(phi1_rad)
+
+    long = long_origin + long * RAD_TO_DEG
+
+    { :lat => lat, :long => long }
+  end
+  
+  private 
+  
+  def self.equitorial_radius(ellipsoid)
+    ELLIPSOID[ellipsoid].first
+  end
+  
+  def self.eccentricity_squared(ellipsoid)
+    ELLIPSOID[ellipsoid].last
+  end
+  
+  ELLIPSOID = {
+    "Airy"                   => [6377563, 0.00667054],
+    "Australian National"    => [6378160, 0.006694542],
+    "Bessel 1841"            => [6377397, 0.006674372],
+    "Bessel 1841 (Nambia)"   => [6377484, 0.006674372],
+    "Clarke 1866"            => [6378206, 0.006768658],
+    "Clarke 1880"            => [6378249, 0.006803511],
+    "Everest"                => [6377276, 0.006637847],
+    "Fischer 1960 (Mercury)" => [6378166, 0.006693422],
+    "Fischer 1968"           => [6378150, 0.006693422],
+    "GRS 1967"               => [6378160, 0.006694605],
+    "GRS 1980"               => [6378137, 0.00669438],
+    "Helmert 1906"           => [6378200, 0.006693422],
+    "Hough"                  => [6378270, 0.00672267],
+    "International"          => [6378388, 0.00672267],
+    "Krassovsky"             => [6378245, 0.006693422],
+    "Modified Airy"          => [6377340, 0.00667054],
+    "Modified Everest"       => [6377304, 0.006637847],
+    "Modified Fischer 1960"  => [6378155, 0.006693422],
+    "South American 1969"    => [6378160, 0.006694542],
+    "WGS 60"                 => [6378165, 0.006693422],
+    "WGS 66"                 => [6378145, 0.006694542],
+    "WGS-72"                 => [6378135, 0.006694318],
+    "WGS-84"                 => [6378137, 0.00669438]
+  }
+  
+  RAD_TO_DEG = 180.0 / Math::PI
+end

data/spec/conversion_spec.rb

@@ -0,0 +1,19 @@
+require File.join(File.dirname(__FILE__), '..', 'lib', 'coordinates')
+
+describe Coordinates do
+  describe "#utm_to_lat_long" do
+    it "should return a hash with keys lat and long" do
+      ll = Coordinates.utm_to_lat_long("WGS-84", 6688940, 219165, "33N")
+      ll.should be_an Hash
+      ll.has_key?(:lat).should be_true
+      ll.has_key?(:long).should be_true
+    end
+
+    it "should return a correct lat long for the reference ellipsoid WGS-84 in the zone 33N" do
+      coords = Coordinates.utm_to_lat_long("WGS-84", 6688940, 219165, "33N")
+      
+      coords[:lat].should be_close(60.2399303597183, 0.00000005)
+      coords[:long].should be_close(9.92496237547609, 0.00000005)
+    end
+  end
+end

metadata

@@ -0,0 +1,71 @@
+--- !ruby/object:Gem::Specification 
+name: coordinate-converter
+version: !ruby/object:Gem::Version 
+  hash: 27
+  prerelease: false
+  segments: 
+  - 0
+  - 1
+  - 0
+  version: 0.1.0
+platform: ruby
+authors: 
+- Kim Joar Bekkelund
+autorequire: 
+bindir: bin
+cert_chain: []
+
+date: 2010-08-10 00:00:00 +02:00
+default_executable: 
+dependencies: []
+
+description: Converting UTM to Latitude and Longitude
+email: mail@kimjoar.net
+executables: []
+
+extensions: []
+
+extra_rdoc_files: 
+- README.md
+files: 
+- README.md
+- Rakefile
+- VERSION
+- lib/coordinate-converter.rb
+- spec/conversion_spec.rb
+has_rdoc: true
+homepage: http://github.com/kjbekkelund/coordinate-converter
+licenses: []
+
+post_install_message: 
+rdoc_options: 
+- --charset=UTF-8
+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.3.7
+signing_key: 
+specification_version: 3
+summary: Converting UTM to Latitude and Longitude
+test_files: 
+- spec/conversion_spec.rb