swedishgrid 0.1.3

data/README.rdoc

@@ -0,0 +1,46 @@
+= swedishgrid
+
+== DESCRIPTION:
+
+Convert coordinates between geodetic WGS84 and Swedish grid RT90 and SWEREF99 systems.
+Read more about Swedish grid systems at http://www.lantmateriet.se/templates/LMV_Page.aspx?id=5197
+
+== SYNOPSIS:
+
+  grid = SwedishGrid.new(:rt90)
+  grid.grid_to_geodetic(7453389.762, 1727060.905)
+  # Outputs => [67.0906813246069, 21.034750437141]
+
+== DEPENDENCIES:
+
+No dependencies.
+
+== INSTALL:
+
+  sudo gem install icehouse-swedishgrid
+
+== LICENSE:
+
+SwedishGrid is released under MIT license with permission from Arnold Andreasson.
+Original Javascript-version written by Arnold Andreasson can be
+found at: http://mellifica.se/geodesi/gausskruger.js
+
+Copyright (c) 2009 ICE House
+
+Permission is hereby granted, free of charge, to any person obtaining a
+copy of this software and associated documentation files (the "Software"),
+to deal in the Software without restriction, including without limitation
+the rights to use, copy, modify, merge, publish, distribute, sublicense,
+and/or sell copies of the Software, and to permit persons to whom the
+Software is furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included
+in all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
+OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
+THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
+

data/Rakefile

@@ -0,0 +1,55 @@
+require 'rubygems'
+require 'rake'
+
+begin
+  require 'jeweler'
+  Jeweler::Tasks.new do |gem|
+    gem.name = "swedishgrid"
+    gem.summary = %Q{TODO}
+    gem.email = "info@entanke.se"
+    gem.homepage = "http://github.com/icehouse/swedishgrid"
+    gem.authors = ["Magnus Enarsson"]
+    # gem is a Gem::Specification... see http://www.rubygems.org/read/chapter/20 for additional settings
+  end
+rescue LoadError
+  puts "Jeweler (or a dependency) not available. Install it with: sudo gem install jeweler"
+end
+
+require 'rake/testtask'
+Rake::TestTask.new(:test) do |test|
+  test.libs << 'lib' << 'test'
+  test.pattern = 'test/**/*_test.rb'
+  test.verbose = true
+end
+
+begin
+  require 'rcov/rcovtask'
+  Rcov::RcovTask.new do |test|
+    test.libs << 'test'
+    test.pattern = 'test/**/*_test.rb'
+    test.verbose = true
+  end
+rescue LoadError
+  task :rcov do
+    abort "RCov is not available. In order to run rcov, you must: sudo gem install spicycode-rcov"
+  end
+end
+
+
+task :default => :test
+
+require 'rake/rdoctask'
+Rake::RDocTask.new do |rdoc|
+  if File.exist?('VERSION.yml')
+    config = YAML.load(File.read('VERSION.yml'))
+    version = "#{config[:major]}.#{config[:minor]}.#{config[:patch]}"
+  else
+    version = ""
+  end
+
+  rdoc.rdoc_dir = 'rdoc'
+  rdoc.title = "swedishgrid #{version}"
+  rdoc.rdoc_files.include('README*')
+  rdoc.rdoc_files.include('lib/**/*.rb')
+end
+

data/VERSION.yml

@@ -0,0 +1,4 @@
+--- 
+:minor: 1
+:patch: 3
+:major: 0

data/lib/swedishgrid.rb

@@ -0,0 +1,112 @@
+class SwedishGrid
+  
+  def initialize(grid)
+    case grid
+    when :rt90
+      # rt90_2.5_gon_v
+      @axis = 6378137.0
+      @flattening = 1.0 / 298.257222101
+      @central_meridian = 15.806284529
+      @lat_of_origin = 0.0
+      @scale = 1.00000561024
+      @false_northing = -667.711
+      @false_easting = 1500064.274 
+    when :sweref99tm
+      @axis = 6378137.0
+      @flattening = 1.0 / 298.257222101
+      @central_meridian = 15.00 
+      @lat_of_origin = 0.0
+      @scale = 0.9996
+      @false_northing = 0.0
+      @false_easting = 500000.0
+    else
+      raise "Unknown grid"
+    end
+  end
+
+  # Conversion from grid coordinates to geodetic coordinates.
+  def grid_to_geodetic(x, y)
+    # Prepare ellipsoid-based stuff.
+    e2 = @flattening * (2.0 - @flattening)
+    n = @flattening / (2.0 - @flattening)
+    a_roof = @axis / (1.0 + n) * (1.0 + n*n/4.0 + n*n*n*n/64.0)
+    delta1 = n/2.0 - 2.0*n*n/3.0 + 37.0*n*n*n/96.0 - n*n*n*n/360.0
+    delta2 = n*n/48.0 + n*n*n/15.0 - 437.0*n*n*n*n/1440.0
+    delta3 = 17.0*n*n*n/480.0 - 37*n*n*n*n/840.0  
+    delta4 = 4397.0*n*n*n*n/161280.0
+
+    astar = e2 + e2*e2 + e2*e2*e2 + e2*e2*e2*e2
+    bstar = -(7.0*e2*e2 + 17.0*e2*e2*e2 + 30.0*e2*e2*e2*e2) / 6.0
+    cstar = (224.0*e2*e2*e2 + 889.0*e2*e2*e2*e2) / 120.0
+    dstar = -(4279.0*e2*e2*e2*e2) / 1260.0
+
+    # Convert.
+    deg_to_rad = Math::PI / 180  
+    lambda_zero = @central_meridian * deg_to_rad
+    xi = (x - @false_northing) / (@scale * a_roof)
+    eta = (y - @false_easting) / (@scale * a_roof)
+    xi_prim = xi - 
+            delta1*Math.sin(2.0*xi) * Math.cosh(2.0*eta) - 
+            delta2*Math.sin(4.0*xi) * Math.cosh(4.0*eta) - 
+            delta3*Math.sin(6.0*xi) * Math.cosh(6.0*eta) - 
+            delta4*Math.sin(8.0*xi) * Math.cosh(8.0*eta)  
+    eta_prim = eta - 
+            delta1*Math.cos(2.0*xi) * Math.sinh(2.0*eta) - 
+            delta2*Math.cos(4.0*xi) * Math.sinh(4.0*eta) - 
+            delta3*Math.cos(6.0*xi) * Math.sinh(6.0*eta) - 
+            delta4*Math.cos(8.0*xi) * Math.sinh(8.0*eta)  
+    phi_star = Math.asin(Math.sin(xi_prim) / Math.cosh(eta_prim))  
+    delta_lambda = Math.atan(Math.sinh(eta_prim) / Math.cos(xi_prim))
+    lng_radian = lambda_zero + delta_lambda  
+    lat_radian = phi_star + Math.sin(phi_star) * Math.cos(phi_star) *
+            (astar + 
+             bstar*(Math.sin(phi_star) ** 2) + 
+             cstar*(Math.sin(phi_star) ** 4) + 
+             dstar*(Math.sin(phi_star) ** 6))
+    return [lat_radian * 180.0 / Math::PI, lng_radian * 180.0 / Math::PI]
+  end
+
+  # Conversion from geodetic coordinates to grid coordinates.
+  def geodetic_to_grid(latitude, longitude)
+    # Prepare ellipsoid-based stuff.
+    e2 = @flattening * (2.0 - @flattening)
+    n = @flattening / (2.0 - @flattening)
+    a_roof = @axis / (1.0 + n) * (1.0 + n*n/4.0 + n*n*n*n/64.0)
+    a = e2
+    b = (5.0*e2*e2 - e2*e2*e2) / 6.0
+    c = (104.0*e2*e2*e2 - 45.0*e2*e2*e2*e2) / 120.0
+    d = (1237.0*e2*e2*e2*e2) / 1260.0
+    beta1 = n/2.0 - 2.0*n*n/3.0 + 5.0*n*n*n/16.0 + 41.0*n*n*n*n/180.0
+    beta2 = 13.0*n*n/48.0 - 3.0*n*n*n/5.0 + 557.0*n*n*n*n/1440.0
+    beta3 = 61.0*n*n*n/240.0 - 103.0*n*n*n*n/140.0
+    beta4 = 49561.0*n*n*n*n/161280.0
+
+    # Convert.
+    deg_to_rad = Math::PI / 180.0
+    phi = latitude * deg_to_rad
+    lambda = longitude * deg_to_rad
+    lambda_zero = @central_meridian * deg_to_rad
+
+    phi_star = phi - Math.sin(phi) * Math.cos(phi) * (a + 
+            b*(Math.sin(phi) ** 2) +
+            c*(Math.sin(phi) ** 4) +
+            d*(Math.sin(phi) ** 6))
+    delta_lambda = lambda - lambda_zero
+    xi_prim = Math.atan(Math.tan(phi_star) / Math.cos(delta_lambda))
+    eta_prim = Math.atanh(Math.cos(phi_star) * Math.sin(delta_lambda))
+    x = @scale * a_roof * (xi_prim +
+            beta1 * Math.sin(2.0*xi_prim) * Math.cosh(2.0*eta_prim) +
+            beta2 * Math.sin(4.0*xi_prim) * Math.cosh(4.0*eta_prim) +
+            beta3 * Math.sin(6.0*xi_prim) * Math.cosh(6.0*eta_prim) +
+            beta4 * Math.sin(8.0*xi_prim) * Math.cosh(8.0*eta_prim)) +
+            @false_northing
+    y = @scale * a_roof * (eta_prim +
+            beta1 * Math.cos(2.0*xi_prim) * Math.sinh(2.0*eta_prim) +
+            beta2 * Math.cos(4.0*xi_prim) * Math.sinh(4.0*eta_prim) +
+            beta3 * Math.cos(6.0*xi_prim) * Math.sinh(6.0*eta_prim) +
+            beta4 * Math.cos(8.0*xi_prim) * Math.sinh(8.0*eta_prim)) +
+            @false_easting
+    return [(x * 1000.0).round / 1000.0, (y * 1000.0).round / 1000.0]
+  end
+  
+end

data/test/swedishgrid_test.rb

@@ -0,0 +1,67 @@
+require "test_helper"
+
+class TestSwedishGrid < Test::Unit::TestCase
+
+  LAT_LNG_TOLERANCE = 1.0e-7
+  GRID_TOLERANCE = 0.001
+
+  def setup
+    # Points from dokument "Kontrollpunkter" at lantmateritet.se
+    # http://www.lantmateriet.se/upload/filer/kartor/geodesi_gps_och_detaljmatning/Transformationer/SWEREF99_RT90_Samband/kontrollpunkter.pdf
+    # http://latlong.mellifica.se/
+    # WGS 84 and SWEREF 99 are, in principle, interchangeable
+    
+    @rt90 =       [ [7453389.762, 1727060.905],  
+                    [7047738.415, 1522128.637],  
+                    [6671665.273, 1441843.186],  
+                    [6249111.351, 1380573.079] ]
+
+    @sweref99tm = [ [7454204.638, 761811.242], 
+                    [7046077.605, 562140.337],
+                    [6669189.376, 486557.055],
+                    [6246136.458, 430374.835] ]
+
+    @lat_lng =    [ [67 +  5/60.0 + 26.452769/3600, 21 +  2/60.0 +  5.101575/3600],
+                    [63 + 32/60.0 + 14.761735/3600, 16 + 14/60.0 + 59.594626/3600],
+                    [60 +  9/60.0 + 33.882413/3600, 14 + 45/60.0 + 28.167152/3600], 
+                    [56 + 21/60.0 + 17.199245/3600, 13 + 52/60.0 + 23.754022/3600] ]
+                    
+  end
+
+
+  def test_rt90_grid_to_geodetic
+    grid = SwedishGrid.new(:rt90)
+    (0..3).each do |i|
+      coord = grid.grid_to_geodetic(*@rt90[i])
+      assert (coord[0] - @lat_lng[i][0]).abs < LAT_LNG_TOLERANCE
+      assert (coord[1] - @lat_lng[i][1]).abs < LAT_LNG_TOLERANCE
+    end
+  end
+  
+  def test_rt90_geodetic_to_grid
+    grid = SwedishGrid.new(:rt90)
+    (0..3).each do |i|
+      coord = grid.geodetic_to_grid(*@lat_lng[i])
+      assert (coord[0] - @rt90[i][0]).abs < GRID_TOLERANCE
+      assert (coord[1] - @rt90[i][1]).abs < GRID_TOLERANCE
+    end
+  end
+  
+  def test_sweref99tm_grid_to_geodetic
+    grid = SwedishGrid.new(:sweref99tm)
+    (0..3).each do |i|
+      coord = grid.grid_to_geodetic(*@sweref99tm[i])
+      assert (coord[0] - @lat_lng[i][0]).abs < LAT_LNG_TOLERANCE
+      assert (coord[1] - @lat_lng[i][1]).abs < LAT_LNG_TOLERANCE
+    end
+  end
+
+  def test_sweref99tm_geodetic_to_grid
+    grid = SwedishGrid.new(:sweref99tm)
+    (0..3).each do |i|
+      coord = grid.geodetic_to_grid(*@lat_lng[i])
+      assert (coord[0] - @sweref99tm[i][0]).abs < GRID_TOLERANCE
+      assert (coord[1] - @sweref99tm[i][1]).abs < GRID_TOLERANCE
+    end
+  end
+end

data/test/test_helper.rb

@@ -0,0 +1,9 @@
+require 'rubygems'
+require 'test/unit'
+
+$LOAD_PATH.unshift(File.dirname(__FILE__))
+$LOAD_PATH.unshift(File.join(File.dirname(__FILE__), '..', 'lib'))
+require 'swedishgrid'
+
+class Test::Unit::TestCase
+end

metadata

@@ -0,0 +1,55 @@
+--- !ruby/object:Gem::Specification
+name: swedishgrid
+version: !ruby/object:Gem::Version
+  version: 0.1.3
+  prerelease: 
+platform: ruby
+authors:
+- Magnus Enarsson
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2009-05-05 00:00:00.000000000 Z
+dependencies: []
+description: 
+email: magnus@icehouse.se
+executables: []
+extensions: []
+extra_rdoc_files:
+- README.rdoc
+files:
+- README.rdoc
+- Rakefile
+- VERSION.yml
+- lib/swedishgrid.rb
+- test/swedishgrid_test.rb
+- test/test_helper.rb
+homepage: 
+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
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  none: false
+  requirements:
+  - - ! '>='
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 1.8.24
+signing_key: 
+specification_version: 2
+summary: Convert coordinates between geodetic WGS84 and Swedish grid RT90 and SWEREF99
+  systems.
+test_files:
+- test/swedishgrid_test.rb
+- test/test_helper.rb