mercator 1.0.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 6cca620d720fe01a33df60ad7be336c101d6bf20
+  data.tar.gz: 58d92e6603f0924c73e8021da333cfde51999bae
+SHA512:
+  metadata.gz: 6e2f8c94e81762e2955b6eb5ac96b30db451b8e4c3062bb8e1c0c806cfa049d4e4852a435196b0d421e831ab5ba66159fc015988a2c71d5a54068872c8152ab0
+  data.tar.gz: e83e145f2b7eae6b125c28afe37e0dd5d2eeb99712074e014fdc1098f2e01be37d0a5b86f9c5be66a2da480ebd31a7daa2ceef809e3fde3d8e375f339bd57633

data/.gitignore

@@ -0,0 +1,15 @@
+/.bundle/
+/.yardoc
+/Gemfile.lock
+/_yardoc/
+/coverage/
+/doc/
+/pkg/
+/spec/reports/
+/tmp/
+*.bundle
+*.so
+*.o
+*.a
+mkmf.log
+*.gem

data/Gemfile

@@ -0,0 +1,4 @@
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in mercator.gemspec
+gemspec

data/LICENSE.txt

@@ -0,0 +1,22 @@
+Copyright (c) 2015 Martin Bergek
+
+MIT License
+
+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/README.md

@@ -0,0 +1,74 @@
+# Mercator
+
+This gem include methods to convert between spherical and grid coordinates. It was designed to handle conversions for the Swedish SWEREF99 projection but can be used for any projection based on Gauss conformal projection (transverse Mercator).
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'mercator'
+```
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install mercator
+
+## Usage
+
+Consider the following code:
+
+```require 'mercator'
+datum = Mercator::WGS84.new
+projection = Mercator::Sweref99TM.new
+Mercator.grid_to_spherical(6097106.672, 356083.438, datum, projection)
+```
+
+Available datums:
+
+* GRS80
+* WGS84
+
+Available projections:
+
+* SWEREF 99 TM
+* SWEREF 99 12 00
+* SWEREF 99 13 30
+* SWEREF 99 15 00
+* SWEREF 99 16 30
+* SWEREF 99 18 00
+* SWEREF 99 14 15
+* SWEREF 99 15 45
+* SWEREF 99 17 15
+* SWEREF 99 18 45
+* SWEREF 99 20 15
+* SWEREF 99 21 45
+* SWEREF 99 23 15
+
+To use other datums and projections, instantiate Mercator::Datum and Mercator::Projection objects and assign suitable parameters, e.g. (for SWEREF 99 TM above):
+
+```datum = Mercator::Datum.new
+datum.semi_major_axis = 6378137.0
+datum.inverse_flattening = 298.257223563
+
+projection = Mercator::Projection.new
+projection.central_meridian = 15.0
+projection.scale = 0.9996
+projection.false_northing = 0.0
+projection.false_easting = 500000.0
+```
+
+For information about the calculations used:  http://www.lantmateriet.se/sv/Kartor-och-geografisk-information/GPS-och-geodetisk-matning/Om-geodesi/Kartprojektioner/
+
+
+## Contributing
+
+1. Fork it ( https://github.com/mbergek/mercator/fork )
+2. Create your feature branch (`git checkout -b my-new-feature`)
+3. Commit your changes (`git commit -am 'Add some feature'`)
+4. Push to the branch (`git push origin my-new-feature`)
+5. Create a new Pull Request

data/Rakefile

@@ -0,0 +1,15 @@
+require "bundler/gem_tasks"
+
+# Test
+require 'rake/testtask'
+Rake::TestTask.new do |t|
+  t.libs << 'test'
+  t.pattern = "test/*_test.rb"
+  t.verbose = true
+end
+
+# Console
+desc "Open an irb session preloaded with this library"
+task :console do
+  sh "irb -rubygems -I lib -r mercator.rb"
+end

data/bin/sweref99tm

@@ -0,0 +1,26 @@
+#!/usr/bin/env ruby
+
+require 'mercator'
+
+if ARGV.length < 2
+  puts <<-END.gsub(/^ {4}/, '')
+    Usage:
+  
+      sweref99tm x y
+  
+    Where x and y are the grid coordinates using SWEREF 99 TM. This will
+    output the latitude and longitude of the coordinates supplied. No
+    checks will be made as to whether the supplied coordinates are valid
+    or not.
+    
+    Example:
+    
+      sweref99tm 6097106.672 356083.438
+
+    END
+    exit
+end
+
+datum = Mercator::WGS84.new
+projection = Mercator::Sweref99TM.new
+p Mercator.grid_to_spherical(ARGV[0].to_f, ARGV[1].to_f, datum, projection)

data/lib/mercator/datum.rb

@@ -0,0 +1,37 @@
+module Mercator
+
+  # Base datum class. Do not use this class in calculations, instead
+  # use one of the subclasses defined below.
+
+  class Datum
+    attr_accessor :semi_major_axis
+    attr_accessor :semi_minor_axis
+    attr_accessor :inverse_flattening
+    
+    def flattening
+      1 / self.inverse_flattening
+    end
+  end
+
+  # The GRS80 datum, or the Geodetic Reference System 1980.
+
+  class GRS80 < Datum
+    def initialize
+      self.semi_major_axis = 6378137.0
+      self.semi_minor_axis = 6356752.314140
+      self.inverse_flattening = 298.257222100882
+    end
+  end
+
+  # The WGS84 datum, World Geodetic System 1984
+
+  class WGS84 < Datum
+    def initialize
+      self.semi_major_axis = 6378137.0
+      self.semi_minor_axis = 6356752.314245
+      self.inverse_flattening = 298.257223563
+    end
+  end
+
+
+end

data/lib/mercator/projection.rb

@@ -0,0 +1,111 @@
+module Mercator
+
+  # Base projection class. Normally you would not use this class in calculations,
+  # instead use one of the subclasses defined below.
+
+  class Projection
+    attr_accessor :central_meridian
+    attr_accessor :scale
+    attr_accessor :false_easting
+    attr_accessor :false_northing
+  end
+
+  # Base SWEREF 99 class that sets some common SWEREF 99 properties.
+  
+  class Sweref99 < Projection
+    def initialize
+      super
+      self.scale = 1.0
+      self.false_northing = 0.0
+      self.false_easting = 150000.0
+    end
+  end
+
+  # SWEREF 99 TM class, one of the most common projection used by
+  # Swedish applications.
+
+  class Sweref99TM < Sweref99
+    def initialize
+      super
+      self.central_meridian = 15.0
+      self.scale = 0.9996
+      self.false_easting = 500000.0
+    end
+  end
+
+  # Various regional SWEREF classes for Sweden
+
+  class Sweref991200 < Sweref99
+    def initialize
+      super
+      self.central_meridian = 12.00
+    end
+  end
+  class Sweref991330 < Sweref99
+    def initialize
+      super
+      self.central_meridian = 13.50
+    end
+  end
+  class Sweref991500 < Sweref99
+    def initialize
+      super
+      self.central_meridian = 15.00
+    end
+  end
+  class Sweref991630 < Sweref99
+    def initialize
+      super
+      self.central_meridian = 16.50
+    end
+  end
+  class Sweref991800 < Sweref99
+    def initialize
+      super
+      self.central_meridian = 18.00
+    end
+  end
+  class Sweref991415 < Sweref99
+    def initialize
+      super
+      self.central_meridian = 14.25
+    end
+  end
+  class Sweref991545 < Sweref99
+    def initialize
+      super
+      self.central_meridian = 15.75
+    end
+  end
+  class Sweref991715 < Sweref99
+    def initialize
+      super
+      self.central_meridian = 17.25
+    end
+  end
+  class Sweref991845 < Sweref99
+    def initialize
+      super
+      self.central_meridian = 18.75
+    end
+  end
+  class Sweref992015 < Sweref99
+    def initialize
+      super
+      self.central_meridian = 20.25
+    end
+  end
+  class Sweref992145 < Sweref99
+    def initialize
+      super
+      self.central_meridian = 21.75
+    end
+  end
+  class Sweref992315 < Sweref99
+    def initialize
+      super
+      self.central_meridian = 23.25
+    end
+  end
+
+end

data/lib/mercator/version.rb

@@ -0,0 +1,3 @@
+module Mercator
+  VERSION = "1.0.0"
+end

data/lib/mercator.rb

@@ -0,0 +1,91 @@
+require "mercator/version"
+require "mercator/datum"
+require "mercator/projection"
+
+include Math
+
+
+# The module contains methods to convert between grid
+# and spherical coordinates. For information about the calculations
+# used: http://bit.ly/1HGZduL
+#
+# Author::      Martin Bergek (mailto:contact@spotwise.com)
+# Copyright::   Copyright (c) 2015 Martin Bergek
+# License::     MIT License
+
+
+module Mercator
+
+  # Convert from spherical coordinates to grid coordinates using the
+  # specified datum and projection.
+
+  def self.spherical_to_grid(latitude, longitude, datum, projection)
+    φ   = latitude * PI / 180.0
+    λ   = longitude * PI / 180.0
+    λ₀  = projection.central_meridian.to_rad
+
+    e²  = datum.flattening * (2 - datum.flattening)
+    n   = datum.flattening / (2 - datum.flattening)
+    â   = datum.semi_major_axis/(1+n) * (1 + n**2/4 + n**4/64)
+
+    a   = e²
+    b   = (5*e²**2 - e²**3)/6
+    c   = (104*e²**3 - 45*e²**4)/120
+    d   = 1237*e²**4/1260
+
+    φ′  = φ - sin(φ)*cos(φ)*(a + b*sin(φ)**2 + c*sin(φ)**4 + d*sin(φ)**6)
+    δ   = λ - λ₀
+
+    ξ   = atan(tan(φ′) / cos(δ))
+    η   = atanh(cos(φ′) * sin(δ))
+
+    β₁  = n/2 - 2*n**2/3 + 5*n**3/16 + 41*n**4/180
+    β₂  = 13*n**2/48 - 3*n**3/5 + 557*n**4/1440
+    β₃  = 61*n**3/240 - 103*n**4/180
+    β₄  = 49561*n**4/161280
+
+    x   = projection.scale * â * (ξ + β₁*sin(2*ξ)*cosh(2*η) + β₂*sin(4*ξ)*cosh(4*η) +
+          β₃*sin(6*ξ)*cosh(6*η) + β₄*sin(8*ξ)*cosh(8*η) ) + projection.false_northing
+    y   = projection.scale * â * (η + β₁*cos(2*ξ)*sinh(2*η) + β₂*cos(4*ξ)*sinh(4*η) +
+          β₃*cos(6*ξ)*sinh(6*η) + β₄*cos(8*ξ)*sinh(8*η) ) + projection.false_easting
+
+    return [x, y]
+  end
+
+  # Convert from grid coordinates to spherical coordinates using the
+  # specified datum and projection.
+
+  def self.grid_to_spherical(north, east, datum, projection)
+    e²  = datum.flattening * (2 - datum.flattening)
+    n   = datum.flattening / (2 - datum.flattening)
+    â   = datum.semi_major_axis/(1+n) * (1 + n**2/4 + n**4/64)
+    λ₀  = projection.central_meridian * PI / 180.0
+
+    ξ   = (north - projection.false_northing) / projection.scale / â
+    η   = (east - projection.false_easting) / projection.scale / â
+
+    δ₁  = n/2 - 2*n**2/3 + 37*n**3/96 - n**4/360
+    δ₂  = n**2/48 + n**3/15 - 437*n**4/1440
+    δ₃  = 17*n**3/480 - 37*n**4/840
+    δ₄  = 4397*n**4/161280
+
+    ξ′  = ξ - δ₁*sin(2*ξ)*cosh(2*η) - δ₂*sin(4*ξ)*cosh(4*η) -
+          δ₃*sin(6*ξ)*cosh(6*η) - δ₄*sin(8*ξ)*cosh(8*η)
+    η′  = η - δ₁*cos(2*ξ)*sinh(2*η) - δ₂*cos(4*ξ)*sinh(4*η) -
+          δ₃*cos(6*ξ)*sinh(6*η) - δ₄*cos(8*ξ)*sinh(8*η)
+
+    φ⋆  = asin(sin(ξ′) / cosh(η′))
+    δ   = atan(sinh(η′) / cos(ξ′))
+
+    a⋆  = e² + e²**2 + e²**3 + e²**4
+    b⋆  = -1.0/6*(7*e²**2 + 17*e²**3 + 30*e²**4)
+    c⋆  = (224*e²**3 + 889*e²**4)/120
+    d⋆  = -1.0/1260*4279*e²**4
+
+    λ   = λ₀ + δ
+    φ   = φ⋆ + sin(φ⋆)*cos(φ⋆) * (a⋆ + b⋆*sin(φ⋆)**2 + c⋆*sin(φ⋆)**4 + d⋆*sin(φ⋆)**6)
+
+    return [λ * 180.0 / PI, φ * 180.0 / PI]
+  end
+
+end

data/mercator.gemspec

@@ -0,0 +1,26 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'mercator/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = "mercator"
+  spec.version       = Mercator::VERSION
+  spec.authors       = ["Martin Bergek"]
+  spec.email         = ["contact@spotwise.com"]
+  spec.summary       = %q{Tools to manage spherical to grid conversions}
+  spec.description   = %q{This library provides classes and methods to perform conversions between spherical and grid coordinates.}
+  spec.homepage      = "https://github.com/mbergek/mercator"
+  spec.license       = "MIT"
+
+  spec.files         = `git ls-files -z`.split("\x0")
+  spec.executables   = spec.files.grep(%r{^bin/}) { |f| File.basename(f) }
+  spec.test_files    = spec.files.grep(%r{^(test|spec|features)/})
+  spec.require_paths = ["lib"]
+
+  spec.add_development_dependency "bundler", "~> 1.7"
+  spec.add_development_dependency "rake", "~> 10.0"
+  spec.add_development_dependency 'minitest', '~> 4.7', '>= 4.7.3'
+  spec.add_runtime_dependency 'tilt', '~> 1.4', '>= 1.4.1'
+  spec.add_runtime_dependency 'redcarpet', '~> 3.3', '>= 3.3.2'
+end

data/test/mercator_test.rb

@@ -0,0 +1,71 @@
+require 'test_helper'
+
+# http://bit.ly/1HGXANA
+
+module MiniTest::Assertions
+  def assert_almost_equal(exp, act)
+    if exp.kind_of?(Array) && act.kind_of?(Array)
+      exp.zip(act).each do |obj|
+        assert_in_delta(obj[0], obj[1])
+      end
+    else
+      assert_in_delta(exp, act)
+    end
+  end
+end
+
+class MercatorTest < Minitest::Unit::TestCase
+
+  def setup
+    @datum = Mercator::WGS84.new
+    @projection = Mercator::Sweref99TM.new    
+  end
+
+  # These are the control points for SWEREF 99 TM as defined on the web page http://bit.ly/1HGXANA
+  def test_grid_to_spherical_55_00_12_45
+    assert_almost_equal [12.75, 55.00], Mercator.grid_to_spherical(6097106.672, 356083.438, @datum, @projection)
+  end
+  def test_grid_to_spherical_55_00_14_15
+    assert_almost_equal [14.25, 55.00], Mercator.grid_to_spherical(6095048.642, 452024.069, @datum, @projection)
+  end
+  def test_grid_to_spherical_57_00_12_45
+    assert_almost_equal [12.75, 57.00], Mercator.grid_to_spherical(6319636.937, 363331.554, @datum, @projection)
+  end
+  def test_grid_to_spherical_57_00_19_30
+    assert_almost_equal [19.50, 57.00], Mercator.grid_to_spherical(6326392.707, 773251.054, @datum, @projection)
+  end
+  def test_grid_to_spherical_59_00_11_15
+    assert_almost_equal [11.25, 59.00], Mercator.grid_to_spherical(6546096.724, 284626.066, @datum, @projection)
+  end
+  def test_grid_to_spherical_59_00_19_30
+    assert_almost_equal [19.50, 59.00], Mercator.grid_to_spherical(6548757.206, 758410.519, @datum, @projection)
+  end
+  def test_grid_to_spherical_61_00_12_45
+    assert_almost_equal [12.75, 61.00], Mercator.grid_to_spherical(6764877.311, 378323.440, @datum, @projection)
+  end
+  def test_grid_to_spherical_61_00_18_45
+    assert_almost_equal [18.75, 61.00], Mercator.grid_to_spherical(6768593.345, 702745.127, @datum, @projection)
+  end
+  def test_grid_to_spherical_63_00_12_00
+    assert_almost_equal [12.00, 63.00], Mercator.grid_to_spherical(6989134.048, 348083.148, @datum, @projection)
+  end
+  def test_grid_to_spherical_63_00_19_30
+    assert_almost_equal [19.50, 63.00], Mercator.grid_to_spherical(6993565.630, 727798.671, @datum, @projection)
+  end
+  def test_grid_to_spherical_65_00_13_30
+    assert_almost_equal [13.50, 65.00], Mercator.grid_to_spherical(7209293.753, 429270.201, @datum, @projection)
+  end
+  def test_grid_to_spherical_65_00_21_45
+    assert_almost_equal [21.75, 65.00], Mercator.grid_to_spherical(7225449.115, 817833.405, @datum, @projection)
+  end
+  def test_grid_to_spherical_67_00_16_30
+    assert_almost_equal [16.50, 67.00], Mercator.grid_to_spherical(7432168.174, 565398.458, @datum, @projection)
+  end
+  def test_grid_to_spherical_67_00_24_00
+    assert_almost_equal [24.00, 67.00], Mercator.grid_to_spherical(7459745.672, 891298.142, @datum, @projection)
+  end
+  def test_grid_to_spherical_69_00_21_00
+    assert_almost_equal [21.00, 69.00], Mercator.grid_to_spherical(7666089.698, 739639.195, @datum, @projection)
+  end
+
+end

data/test/test_helper.rb

@@ -0,0 +1,2 @@
+require 'mercator'
+require 'minitest/autorun'

metadata

@@ -0,0 +1,149 @@
+--- !ruby/object:Gem::Specification
+name: mercator
+version: !ruby/object:Gem::Version
+  version: 1.0.0
+platform: ruby
+authors:
+- Martin Bergek
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2015-07-22 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.7'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.7'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10.0'
+- !ruby/object:Gem::Dependency
+  name: minitest
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '4.7'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 4.7.3
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '4.7'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 4.7.3
+- !ruby/object:Gem::Dependency
+  name: tilt
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.4'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 1.4.1
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '1.4'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 1.4.1
+- !ruby/object:Gem::Dependency
+  name: redcarpet
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3.3'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 3.3.2
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3.3'
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 3.3.2
+description: This library provides classes and methods to perform conversions between
+  spherical and grid coordinates.
+email:
+- contact@spotwise.com
+executables:
+- sweref99tm
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- Gemfile
+- LICENSE.txt
+- README.md
+- Rakefile
+- bin/sweref99tm
+- lib/mercator.rb
+- lib/mercator/datum.rb
+- lib/mercator/projection.rb
+- lib/mercator/version.rb
+- mercator.gemspec
+- test/mercator_test.rb
+- test/test_helper.rb
+homepage: https://github.com/mbergek/mercator
+licenses:
+- MIT
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.4.5
+signing_key: 
+specification_version: 4
+summary: Tools to manage spherical to grid conversions
+test_files:
+- test/mercator_test.rb
+- test/test_helper.rb