swissgrid 0.0.1

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 9e5587b567763c2ca15b1d8b2b9b483f98d3d470
+  data.tar.gz: 9cc965ba0e9bf87cde2ff83f86366350c923f613
+SHA512:
+  metadata.gz: d13624a5963cdb45ce43c8dff246732b58c34bacdba02b6f3704504268872d39bef22d197cf2be415d3e674ac8c6b50a0c43dd046cb864ebb776ad2835515d57
+  data.tar.gz: b72d2f7e60cc1dd09b97aa6b33eb1dcb5cd2f0440fb48150dd8b2c5238158cafa56516df6917f0373b547abd137b215585f7edea52bbac1f691ac0f225511ea2

data/.gitignore

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

data/Gemfile

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

data/LICENSE.txt

@@ -0,0 +1,22 @@
+Copyright (c) 2015 Alexander Rueedlinger
+
+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,58 @@
+# Swissgrid
+
+The swissgrid rubygem is conversion library to convert WGS84 (GPS) points into the Swiss coordinate system (CH1903)
+and vice versa.
+
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'swissgrid'
+```
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install swissgrid
+
+##Usage
+Require the swissgrid library as follows:
+```ruby
+require 'swissgrid'
+```
+###From WGS84(GPS) to CH1903
+
+```ruby
+ gps_point = [46.881908, 7.471829] # [lat, lon] Bern
+ swiss_coord = Swissgrid::CH1903(gps_point)
+```
+###From CH1903 to WGS84(GPS)
+```ruby
+ ch1903_point = [600000, 200000] # [y, x] Bern
+ wgs84_coord =  Swissgrid::WGS84(ch1903_point)
+```
+
+###Convert a list of points
+```ruby
+ gps_points = [ [46.881908, 7.471829], [47.220833, 7.028056] ]
+
+# convert to ch1903 points
+ ch1903_points = gps_points.map { |p| Swissgrid::CH1903(p) }
+
+# convert back to gps points
+ new_gps_points = ch1903_points.map { |p| Swissgrid::WGS84(p) }
+```
+
+
+## Contributing
+
+1. Fork it ( https://github.com/[my-github-username]/swissgrid/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,2 @@
+require "bundler/gem_tasks"
+

data/lib/swissgrid/ch1903.rb

@@ -0,0 +1,66 @@
+module Swissgrid
+
+  module Ch1903
+
+    PHI_BERN = 169028.66
+    LAMBDA_BERN = 26782.5
+    SCALE = 10000
+
+    def self.from_wgs84(a_point)
+      lat, lon, z = a_point
+      phi_sec = to_sexagesimal(lat)
+      lambda_sec = to_sexagesimal(lon)
+      lambda_prime = (lambda_sec - LAMBDA_BERN) / SCALE.to_f
+      phi_prime = (phi_sec - PHI_BERN) / SCALE.to_f
+      convert(lambda_prime, phi_prime, z)
+    end
+
+    private
+
+    def self.convert(lambda_prime, phi_prime, h)
+      y = [
+          600072.37,
+          211455.93 * lambda_prime,
+          -10938.51 * lambda_prime * phi_prime,
+          -0.36 * lambda_prime * (phi_prime ** 2),
+          -44.54 * (lambda_prime ** 3)
+      ].reduce(:+)
+
+      x = [
+          200147.07,
+          308807.95 * phi_prime,
+          3745.25 * (lambda_prime ** 2),
+          76.63 * (phi_prime ** 2),
+          -194.56 * (lambda_prime ** 2) * phi_prime,
+          119.79 * (phi_prime ** 3)
+      ].reduce(:+)
+
+      z = if h.nil?
+            h
+          else
+            [
+                h - 49.55,
+                2.73 * lambda_prime,
+                6.94 * phi_prime
+            ].reduce(:+)
+          end
+
+      z.nil? ? [y, x] : [y, x, z]
+    end
+
+
+    def self.to_gps_triple(dec_degress)
+      degrees = dec_degress.to_i
+      minutes = (60*(dec_degress - degrees)).to_i
+      seconds = (60*(dec_degress - degrees) - minutes) * 60
+      [degrees, minutes, seconds]
+    end
+
+    def self.to_sexagesimal(dec_degress)
+      degree, minutes, seconds = to_gps_triple(dec_degress)
+      degree * 3600 + minutes * 60 + seconds
+    end
+
+  end
+
+end

data/lib/swissgrid/version.rb

@@ -0,0 +1,3 @@
+module Swissgrid
+  VERSION = "0.0.1"
+end

data/lib/swissgrid/wgs84.rb

@@ -0,0 +1,52 @@
+module Swissgrid
+
+  module WGS84
+
+    Y_BERN = 600000
+    X_BERN = 200000
+    SCALE = 1000
+
+    def self.from_ch1903(a_point)
+      y, x , z = a_point
+      y_prime = (y - Y_BERN) / SCALE.to_f
+      x_prime = (x - X_BERN) / SCALE.to_f
+      convert(y_prime, x_prime, z)
+    end
+
+    private
+
+    def self.convert(y_prime, x_prime, h_prime)
+      lambda_prime = [
+          2.6779094,
+          4.728982 * y_prime,
+          0.791484 * y_prime * x_prime,
+          0.1306 * y_prime * (x_prime ** 2),
+          -0.0436 * (y_prime ** 3)
+      ].reduce(:+)
+
+      phi_prime = [
+          16.9023892 ,
+          3.238272 * x_prime,
+          -0.270978 * (y_prime ** 2),
+          -0.002528 * (x_prime ** 2),
+          -0.0447 * (y_prime ** 2) * x_prime,
+          -0.0140  * (x_prime ** 3)
+      ].reduce(:+)
+
+      z = if h_prime.nil?
+            h_prime
+          else
+            [
+                h_prime + 49.55,
+                -12.60 * y_prime,
+                -22.64 * x_prime
+            ].reduce(:+)
+          end
+
+      lat, lon = [phi_prime, lambda_prime].map { |v| v * 100/36.to_f }
+      z.nil? ? [lat, lon] : [lat, lon, z]
+    end
+
+  end
+
+end

data/lib/swissgrid.rb

@@ -0,0 +1,15 @@
+require "swissgrid/version"
+require 'swissgrid/ch1903'
+require 'swissgrid/wgs84'
+
+module Swissgrid
+
+  def self.CH1903(a_point)
+    Ch1903.from_wgs84(a_point)
+  end
+
+  def self.WGS84(a_point)
+    WGS84.from_ch1903(a_point)
+  end
+
+end

data/spec/spec_helper.rb

@@ -0,0 +1,2 @@
+require 'swissgrid'
+require 'pp'

data/spec/swissgrid_spec.rb

@@ -0,0 +1,26 @@
+require 'spec_helper'
+
+module Swissgrid
+
+  describe Swissgrid do
+
+    it "test conversion from wgs84 to ch1903" do
+      gps_point = [46.881908, 7.471829]
+      swiss_coord = Swissgrid::CH1903(gps_point)
+      expect(swiss_coord[0]).to be_within(0.5).of(602530.221)
+      expect(swiss_coord[1]).to be_within(0.5).of(192310.331)
+      pp swiss_coord
+    end
+
+
+    it "test conversion from ch1903 to wgs84" do
+      ch1903_point = [600000, 200000]
+      wgs84_coord =  Swissgrid::WGS84(ch1903_point)
+      expect(wgs84_coord[0]).to be_within(0.1/3600).of(46.951082877)
+      expect(wgs84_coord[1]).to be_within(0.1/3600).of(7.438632495)
+      pp wgs84_coord
+    end
+
+  end
+
+end

data/swissgrid.gemspec

@@ -0,0 +1,23 @@
+# coding: utf-8
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+require 'swissgrid/version'
+
+Gem::Specification.new do |spec|
+  spec.name          = "swissgrid"
+  spec.version       = Swissgrid::VERSION
+  spec.authors       = ["Alexander Rueedlinger"]
+  spec.email         = ["a.rueedlinger@gmail.com"]
+  spec.summary       = %q{A library to convert gps points into the Swiss coordinate system (CH1903).}
+  spec.description   = %q{Swissgrid is a library to convert gps points into the Swiss coordinate system (CH1903).}
+  spec.homepage      = ""
+  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"
+end

metadata

@@ -0,0 +1,88 @@
+--- !ruby/object:Gem::Specification
+name: swissgrid
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: ruby
+authors:
+- Alexander Rueedlinger
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2015-01-02 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'
+description: Swissgrid is a library to convert gps points into the Swiss coordinate
+  system (CH1903).
+email:
+- a.rueedlinger@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- ".gitignore"
+- Gemfile
+- LICENSE.txt
+- README.md
+- Rakefile
+- lib/swissgrid.rb
+- lib/swissgrid/ch1903.rb
+- lib/swissgrid/version.rb
+- lib/swissgrid/wgs84.rb
+- spec/spec_helper.rb
+- spec/swissgrid_spec.rb
+- swissgrid.gemspec
+homepage: ''
+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.2.2
+signing_key: 
+specification_version: 4
+summary: A library to convert gps points into the Swiss coordinate system (CH1903).
+test_files:
+- spec/spec_helper.rb
+- spec/swissgrid_spec.rb
+has_rdoc: