geodesics 1.0.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: 9fb097bb88b35fc4b47610fbc4781aadebadbab28f6b4c207b144645515d9a19
+  data.tar.gz: 1f1c7ccb9105b0c3dcd8be45fb5bd512028721affd3cf65d59509f358cee92c5
+SHA512:
+  metadata.gz: e102a7f5f1158ae93228e37c1e667ac518f23b47683cf8626ca833f3482c94feb6002bc514c58ac91ebcc2c4e52083f32256c91e1fb4a6b7b615f539c657206b
+  data.tar.gz: 39a076d912a139be79544998462b247e6a425fd4ac7667e2a2dee6bafab64855fbbad0c7d18d21a2d2008b92da8627c6a13b9efc726175cc70b40a15664caf66

data/lib/geodesics.rb

@@ -0,0 +1,34 @@
+# frozen_string_literal: true
+
+require 'geodesics/strategies/lambert'
+require 'geodesics/point'
+require 'forwardable'
+
+class Geodesics
+  class << self
+    extend Forwardable
+    def_delegators :new, :distance, :distance_degree, :distance_radian
+  end
+
+  STRATEGY = Strategies::Lambert.new
+
+  def initialize(strategy: STRATEGY)
+    @strategy = strategy
+  end
+
+  def distance_radian(latitude1, longitude1, latitude2, longitude2)
+    @strategy.distance(
+      Radian.new(latitude1, longitude1),
+      Radian.new(latitude2, longitude2)
+    )
+  end
+
+  def distance_degree(latitude1, longitude1, latitude2, longitude2)
+    @strategy.distance(
+      Degree.new(latitude1, longitude1).to_radian,
+      Degree.new(latitude2, longitude2).to_radian
+    )
+  end
+
+  alias distance distance_degree
+end

data/lib/geodesics/central_angles/haversine.rb

@@ -0,0 +1,18 @@
+# frozen_string_literal: true
+
+class Geodesics
+  module CentralAngles
+    # https://en.wikipedia.org/wiki/Versine#Haversine
+    class Haversine
+      def call(point1, point2)
+        2 * Math.asin(
+          Math.sqrt(
+            Math.sin((point1.latitude - point2.latitude).abs / 2)**2 +
+            Math.cos(point1.latitude) * Math.cos(point2.latitude) *
+            Math.sin((point1.longitude - point2.longitude).abs / 2)**2
+          )
+        )
+      end
+    end
+  end
+end

data/lib/geodesics/central_angles/spherical.rb

@@ -0,0 +1,16 @@
+# frozen_string_literal: true
+
+class Geodesics
+  module CentralAngles
+    # https://en.wikipedia.org/wiki/Spherical_law_of_cosines
+    class Spherical
+      def call(point1, point2)
+        Math.acos(
+          Math.sin(point1.latitude) * Math.sin(point2.latitude) +
+          Math.cos(point1.latitude) * Math.cos(point2.latitude) *
+          Math.cos((point1.longitude - point2.longitude).abs)
+        )
+      end
+    end
+  end
+end

data/lib/geodesics/central_angles/vincenty.rb

@@ -0,0 +1,20 @@
+# frozen_string_literal: true
+
+class Geodesics
+  module CentralAngles
+    # https://en.wikipedia.org/wiki/Vincenty%27s_formulae
+    class Vincenty
+      def call(point1, point2)
+        delta_longitude = (point1.longitude - point2.longitude).abs
+        Math.atan(
+          Math.sqrt(
+            (Math.cos(point2.latitude) * Math.sin(delta_longitude))**2 +
+            (Math.cos(point1.latitude) * Math.sin(point2.latitude) - Math.sin(point1.latitude) * Math.cos(point2.latitude) * Math.cos(delta_longitude))**2
+          ) / (
+            Math.sin(point1.latitude) * Math.sin(point2.latitude) +  Math.cos(point1.latitude) * Math.cos(point2.latitude) * Math.cos(delta_longitude)
+          )
+        )
+      end
+    end
+  end
+end

data/lib/geodesics/point.rb

@@ -0,0 +1,47 @@
+# frozen_string_literal: true
+
+class Geodesics
+  class Point
+    RADIAN = Math::PI / 180
+    DEGREE = 180 / Math::PI
+
+    attr_reader :latitude, :longitude
+
+    def initialize(latitude, longitude)
+      @latitude = latitude
+      @longitude = longitude
+    end
+
+    def ==(other)
+      latitude == other.latitude && longitude == other.longitude
+    end
+
+    def to_radian
+      self
+    end
+
+    def to_degree
+      self
+    end
+  end
+
+  class Radian < Point
+    def to_degree
+      Degree.new(@latitude * DEGREE, @longitude * DEGREE)
+    end
+
+    def ==(other)
+      other.is_a?(Degree) ? super(other.to_radian) : super
+    end
+  end
+
+  class Degree < Point
+    def to_radian
+      Radian.new(@latitude * RADIAN, @longitude * RADIAN)
+    end
+
+    def ==(other)
+      other.is_a?(Radian) ? super(other.to_degree) : super
+    end
+  end
+end

data/lib/geodesics/strategies/lambert.rb

@@ -0,0 +1,51 @@
+# frozen_string_literal: true
+
+require 'geodesics/central_angles/haversine'
+
+class Geodesics
+  module Strategies
+    class Lambert
+      INVERSE_FLATTENING = 298.257223563 # WGS 84
+      FLATTENING = 1.0 / INVERSE_FLATTENING
+      EQUATORIAL_RADIUS = 6_378_137
+      CENTRAL_ANGLE = CentralAngles::Haversine.new
+
+      def initialize(
+        equatorial_radius: EQUATORIAL_RADIUS,
+        flattening: FLATTENING,
+        central_angle: CENTRAL_ANGLE
+      )
+        @equatorial_radius = equatorial_radius
+        @flattening = flattening
+        @central_angle = central_angle
+      end
+
+      def distance(point1, point2)
+        sigma = @central_angle.call(point1, point2)
+        beta1 = parametric_latitude(point1.latitude)
+        beta2 = parametric_latitude(point2.latitude)
+        p = (beta1 + beta2) / 2
+        q = (beta2 - beta1) / 2
+        x = (sigma - sin(sigma)) * sin(p)**2 * cos(q)**2 / cos(sigma / 2)**2
+        y = (sigma + sin(sigma)) * cos(p)**2 * sin(q)**2 / sin(sigma / 2)**2
+        @equatorial_radius * (sigma - @flattening * (x + y) / 2)
+      end
+
+      private
+
+      def parametric_latitude(latitude)
+        Math.atan((1 - @flattening) * Math.tan(latitude))
+      end
+
+      def sin(radian)
+        v = Math.sin(radian)
+        v < 1e-10 ? 1e-10 : v
+      end
+
+      def cos(radian)
+        v = Math.cos(radian)
+        v < 1e-10 ? 1e-10 : v
+      end
+    end
+  end
+end

metadata

@@ -0,0 +1,92 @@
+--- !ruby/object:Gem::Specification
+name: geodesics
+version: !ruby/object:Gem::Version
+  version: 1.0.0
+platform: ruby
+authors:
+- Jian Weihang
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2019-05-11 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: minitest
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 5.11.3
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 5.11.3
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 12.3.2
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 12.3.2
+- !ruby/object:Gem::Dependency
+  name: rubocop
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 0.68.1
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 0.68.1
+description: geodesics calculates the geodesic distance between 2 points with latitude
+  and longitude on ellipsoid Earth using Lambert's formula.
+email: tonytonyjan@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- lib/geodesics.rb
+- lib/geodesics/central_angles/haversine.rb
+- lib/geodesics/central_angles/spherical.rb
+- lib/geodesics/central_angles/vincenty.rb
+- lib/geodesics/point.rb
+- lib/geodesics/strategies/lambert.rb
+homepage: https://github.com/tonytonyjan/geodesics
+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: []
+rubygems_version: 3.0.2
+signing_key: 
+specification_version: 4
+summary: geodesics calculates the geodesic distance between 2 points with latitude
+  and longitude on ellipsoid Earth using Lambert's formula.
+test_files: []