wgs84_distance 0.0.1

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 65ce75ae8c20a7c3728757c84d5fc77ec311e707
+  data.tar.gz: 26bcb444d60e429f5e49d544f4ce538af4a118f3
+SHA512:
+  metadata.gz: fb6af4ace87abe53f2c56277f344d944d6ffac812ddbfcd63f7c571d1322c3a6ad8e7c6f1ea0c60c7a3f3b7c916c787f0ed719c5b5af50035f04d095ee457908
+  data.tar.gz: b04d11e33cfe0da944c8de3736130d63972083373ac6472cda6626009703bc73a7cb825d46166939c7440a6732474284b8fbbb9babf810822e6103dfccc223eb

data/lib/wgs84_distance.rb

@@ -0,0 +1,75 @@
+class WGS84Distance
+  WGS84_MAJOR_AXIS = 6378137.0
+  WGS84_SEMI_MAJOR_AXIS = 6356752.3142
+
+  SQ_CONST1 = (WGS84_MAJOR_AXIS ** 2 - WGS84_SEMI_MAJOR_AXIS ** 2) / WGS84_SEMI_MAJOR_AXIS ** 2
+  SQ_CONST2 = (WGS84_MAJOR_AXIS - WGS84_SEMI_MAJOR_AXIS) / WGS84_MAJOR_AXIS
+
+  def self.between(lat1, lng1, lat2, lng2)
+    lat1, lng1, lat2, lng2 = [lat1, lng1, lat2, lng2].map {|i| i * 3.14159265358979323846 / 180.0}
+
+    iter_a = 0.0
+    u1 = Math.atan((1.0 - SQ_CONST2) * Math.tan(lat1))
+    u2 = Math.atan((1.0 - SQ_CONST2) * Math.tan(lat2))
+
+    cos_u1 = Math.cos(u1)
+    cos_u2 = Math.cos(u2)
+    sin_u1 = Math.sin(u1)
+    sin_u2 = Math.sin(u2)
+    cos_u1_cos_u2 = cos_u1 * cos_u2
+    sin_u1_sin_u2 = sin_u1 * sin_u2
+
+    sigma = 0.0
+    delta_sigma = 0.0
+
+    lamb = lng2 - lng1
+
+    20.times do
+      lambda_orig = lamb
+      coslambda = Math.cos(lamb)
+      sinlambda = Math.sin(lamb)
+
+      sin_sq_sigma = (cos_u2 * sinlambda) ** 2 + (cos_u1 * sin_u2 - sin_u1 * cos_u2 * coslambda) ** 2
+      sinsigma = Math.sqrt(sin_sq_sigma)
+      cossigma = sin_u1_sin_u2 + cos_u1_cos_u2 * coslambda
+      sigma = Math.atan2(sinsigma, cossigma)
+      sinalpha = (sinsigma == 0) ? 0.0 : cos_u1_cos_u2 * sinlambda / sinsigma
+      cos_sq_alpha = 1.0 - sinalpha * sinalpha
+      cos2sm = (cos_sq_alpha == 0) ? 0.0 : cossigma - 2.0 * sin_u1_sin_u2 / cos_sq_alpha
+      u_squared = cos_sq_alpha * SQ_CONST1
+
+      iter_a = 1 + (u_squared / 16384.0) *
+        (4096.0 + u_squared *
+        (-768 + u_squared * (320.0 - 175.0 * u_squared)))
+
+      iter_b = (u_squared / 1024.0) *
+        (256.0 + u_squared *
+        (-128.0 + u_squared * (74.0 - 47.0 * u_squared)))
+
+      iter_c = (SQ_CONST2 / 16.0) *
+        cos_sq_alpha *
+        (4.0 + SQ_CONST2 * (4.0 - 3.0 * cos_sq_alpha))
+
+      cos2smsq = cos2sm * cos2sm
+
+      delta_sigma = iter_b * sinsigma *
+        (cos2sm + (iter_b / 4.0) *
+        (cossigma * (-1.0 + 2.0 * cos2smsq) -
+        (iter_b / 6.0) * cos2sm *
+        (-3.0 + 4.0 * sinsigma * sinsigma) *
+        (-3.0 + 4.0 * cos2smsq)))
+
+      lamb = lng2 - lng1 +
+        (1.0 - iter_c) * SQ_CONST2 * sinalpha *
+        (sigma + iter_c * sinsigma *
+        (cos2sm + iter_c * cossigma *
+        (-1.0 + 2.0 * cos2sm * cos2sm)))
+
+      diff = (lamb - lambda_orig) / lamb
+
+      break if diff.abs < 1.0e-12
+    end
+
+    WGS84_SEMI_MAJOR_AXIS * iter_a * (sigma - delta_sigma)
+  end
+end

metadata

@@ -0,0 +1,44 @@
+--- !ruby/object:Gem::Specification
+name: wgs84_distance
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: ruby
+authors:
+- xyy
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2017-07-25 00:00:00.000000000 Z
+dependencies: []
+description: a wgs84 distance calculator
+email: xyy.zhou@dji.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- lib/wgs84_distance.rb
+homepage: http://rubygems.org/gems/hola
+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.6.10
+signing_key: 
+specification_version: 4
+summary: wgs84 distance
+test_files: []