lambert93_ruby 0.0.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 5411da2feee96931bfd61c2067a1162dec075bca
+  data.tar.gz: 72447d8f4a917f95eef569b65c5e4d7e7aff10cd
+SHA512:
+  metadata.gz: 1f3528beaa81795dbe82637f5c3cc04ea7ba309462f70127354e503dd4e97d29b17ad323b8786b588618f167a3a73ca55e9afb79520acec95ac81d3453832e6f
+  data.tar.gz: f5e9a31a6096a130fa1914841e1cea2a3540f4f741fc6d5965c14291a6f9dfd847557df5dafbe5fb696d0eb18669b44ec494b93ee0bffe517fd7ed22532c151c

data/README.md

@@ -0,0 +1,24 @@
+# lambert93_ruby
+Converts GPS coordinates into Lambert 93 coordinates
+
+# Install
+```
+# In  your gemfile
+gem "lambert93_ruby"
+```
+And run
+```
+bundle install
+```
+
+# Usage
+```
+require "lambert93_ruby"
+
+lat  = 44.884566             # gps_latitude
+long = -0.2846110000000408   # gps_longitude
+
+Lambert93Ruby::lambert_proj(lat, long)
+```
+
+It will give you an array with lambert93 coordinates [x, y].

data/lambert93_ruby.gemspec

@@ -0,0 +1,11 @@
+Gem::Specification.new do |s|
+  s.name     = %q{lambert93_ruby}
+  s.version  = "0.0.0"
+  s.date     = %q{2016-03-25}
+  s.platform = Gem::Platform::RUBY
+  s.summary  = %q{Lambert93_ruby calculates lambert 93 coordinates from GPS coordinates.}
+  s.authors       = ["Kulgar, Arpsara, CT2C"]
+  s.email         = 'contact@ct2c.fr'
+  s.files         = `git ls-files`.split("\n")
+  s.require_paths = ["lib"]
+end

data/lib/lambert93_ruby.rb

@@ -0,0 +1,78 @@
+# -*- encoding : utf-8 -*-
+module Lambert93Ruby
+
+  def self.lambert_proj(lat, long)
+    if lat.nil? || long.nil?
+      return ['', '']
+    else
+      lat = deg_to_rad(lat.to_f)
+      long = deg_to_rad(long.to_f)
+      teta = n_lambert * (long - lambert_parameters[:lng_0])
+      rho_lat = rho(lat)
+      rho_lat_0 = rho(lambert_parameters[:lat_0])
+      x = lambert_parameters[:x_0] + rho_lat * Math.sin(teta)
+      y = lambert_parameters[:y_0] + rho_lat_0 - rho_lat * Math.cos(teta)
+      return [x, y]
+    end
+  end
+
+  private
+
+    def self.lambert_parameters
+      return @lambert_parameters if @lambert_parameters
+      a = 6378137
+      f = 1.0 / 298.257222101
+      b = a * (1.0 - f)
+      e = Math.sqrt(a**2 - b**2) / a
+      @lambert_parameters = {
+        # Lambert 93 consMath.tants:
+        a: a,
+        f: f,
+        b: b,
+        e: e,
+        e_2: (e / 2),
+        lng_0: deg_to_rad(3),
+        lat_0: deg_to_rad(46.5), #deg_to_rad(52)
+        x_0: 700000,
+        y_0: 6600000,
+
+        phi_1: deg_to_rad(44), #deg_to_rad(45.89888888888889)
+        phi_2: deg_to_rad(49) #deg_to_rad(47.69583333333333)
+      }
+    end
+
+    def self.deg_to_rad(deg)
+      deg * Math::PI / 180.0
+    end
+
+    def self.rho(phi)
+      first_part = 1.0 / Math.tan( phi / 2.0 + Math::PI / 4.0 )
+      second_part = ( ( 1 + lambert_parameters[:e] * Math.sin(phi) ) / ( 1 - lambert_parameters[:e] * Math.sin(phi) ) )**lambert_parameters[:e_2]
+      rho_phi = rho_0 * ( first_part * second_part )**n_lambert
+      return rho_phi
+    end
+
+    def self.n_lambert
+      return @n_lambert if @n_lambert
+      n_numerator = Math.log(Math.cos(lambert_parameters[:phi_2]) / Math.cos(lambert_parameters[:phi_1]))  +  0.5 * Math.log( (1.0 - (lambert_parameters[:e] * Math.sin(lambert_parameters[:phi_1]))**2.0 ) / (1.0 - (lambert_parameters[:e] * Math.sin(lambert_parameters[:phi_2]))**2.0  )  )
+
+      n_denominator_numerator = Math.tan(lambert_parameters[:phi_1] / 2.0 + Math::PI / 4.0) * (1.0 - lambert_parameters[:e] * Math.sin(lambert_parameters[:phi_1]))**lambert_parameters[:e_2] * (1.0 + lambert_parameters[:e] * Math.sin(lambert_parameters[:phi_2]))**lambert_parameters[:e_2]
+      n_denominator_denominator = Math.tan(lambert_parameters[:phi_2] / 2.0 + Math::PI / 4.0) * (1.0 + lambert_parameters[:e] * Math.sin(lambert_parameters[:phi_1]))**lambert_parameters[:e_2] * (1.0 - lambert_parameters[:e] * Math.sin(lambert_parameters[:phi_2]))**lambert_parameters[:e_2]
+
+      n_denominator = Math.log(n_denominator_numerator / n_denominator_denominator)
+
+      @n_lambert = n_numerator / n_denominator
+
+      return @n_lambert
+    end
+
+    def self.rho_0
+      return @rho_0 if @rho_0
+      rho_0_part_1 = (lambert_parameters[:a] * Math.cos(lambert_parameters[:phi_1])) / (n_lambert * Math.sqrt(1.0 - (lambert_parameters[:e] * Math.sin(lambert_parameters[:phi_1]))**2.0 ))
+      rho_0_part_2 = (Math.tan(lambert_parameters[:phi_1] / 2.0 + Math::PI / 4.0) * ((1.0 - lambert_parameters[:e] * Math.sin(lambert_parameters[:phi_1])) / (1.0 + lambert_parameters[:e] * Math.sin(lambert_parameters[:phi_1])) )**lambert_parameters[:e_2] )**n_lambert
+
+      @rho_0 = rho_0_part_1 * rho_0_part_2
+
+      return @rho_0
+    end
+end

data/test/lib/lambert93_ruby_test.rb

@@ -0,0 +1,33 @@
+class Lambert93RubyTest < Test::Unit::TestCase
+  require  "lambert93_ruby"
+
+  def test_lambert_proj
+    arveyres = Lambert93Ruby::lambert_proj(45.7156058, 4.85958679999999)
+    assert_equal [844649.3015913883, 6514595.073166063], arveyres
+
+    pusignan = Lambert93Ruby::lambert_proj(45.757652, 5.0701733)
+    assert_equal [860903.96, 6519670.35], pusignan.map{|i| i.round(2)}
+
+    ## Ile de France
+    beaumont_sur_oise = Lambert93Ruby::lambert_proj(49.1415510, 2.2846010 )
+    assert_equal [647796.31706025,  6893766.105254374], beaumont_sur_oise
+
+    chateau_landon = Lambert93Ruby::lambert_proj(48.150232, 2.7044260)
+    assert_equal [678018.7566706568, 6783359.741741321], chateau_landon
+
+    provins = Lambert93Ruby::lambert_proj( 48.5601490, 3.2992030)
+    assert_equal [722078.4810329215, 6828922.492311118], provins
+
+    rambouillet = Lambert93Ruby::lambert_proj(48.6438680, 1.8290789)
+    assert_equal [613737.5587523837, 6838827.35755566], rambouillet
+
+    calais = Lambert93Ruby::lambert_proj(50.95129000000001, 1.8586860000000343)
+    assert_equal [619630.8725958228, 7095614.832194619], calais
+
+    nice = Lambert93Ruby::lambert_proj(43.7101728, 7.261953199999994)
+    assert_equal [1043406.5555545087,  6299396.821049856], nice
+
+    andorre = Lambert93Ruby::lambert_proj( 42.506285, 1.5218009999999822)
+    assert_equal [578337.7238319971, 6157413.046731045], andorre
+  end
+end

data/test/run_test.rb

@@ -0,0 +1,9 @@
+base_dir = File.expand_path(File.join(File.dirname(__FILE__), ".."))
+lib_dir  = File.join(base_dir, "lib")
+test_dir = File.join(base_dir, "test")
+
+$LOAD_PATH.unshift(lib_dir)
+
+require 'test/unit'
+
+exit Test::Unit::AutoRunner.run(true, test_dir)

metadata

@@ -0,0 +1,47 @@
+--- !ruby/object:Gem::Specification
+name: lambert93_ruby
+version: !ruby/object:Gem::Version
+  version: 0.0.0
+platform: ruby
+authors:
+- Kulgar, Arpsara, CT2C
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2016-03-25 00:00:00.000000000 Z
+dependencies: []
+description: 
+email: contact@ct2c.fr
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- README.md
+- lambert93_ruby.gemspec
+- lib/lambert93_ruby.rb
+- test/lib/lambert93_ruby_test.rb
+- test/run_test.rb
+homepage: 
+licenses: []
+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.8
+signing_key: 
+specification_version: 4
+summary: Lambert93_ruby calculates lambert 93 coordinates from GPS coordinates.
+test_files: []