geoprojection 0.1.0 → 0.1.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: b21a5ce188fe4283f22a8ffeae14392942981c101d010f5ae9e9d0f0a67d410e
-  data.tar.gz: 76ac0d07a358785a9505cbcfe4830181b6fe4d7a3b3d923ccf512b13ea906dae
+  metadata.gz: 8560b474f5efac958d7e5296dae0fc6c122d9fc6ea605a8ae406101717a46d06
+  data.tar.gz: 66c955de7126b71538d6d019ad43d767cbed22719929b631bbd0987420070d96
 SHA512:
-  metadata.gz: 9979de234b5ac3661a96561a5e9b836b11351e82d60dd1aece15726a06a1235537324f9baaca026412d7c03712bae235b0ff50951bb0565897e45b717d55f6e1
-  data.tar.gz: 3f572c2025e219ceddb69f6b41e5a87ea84b71ed1ea8a23b3f7c4ab5848a812338af1ee2a34bf3e4b89402bbc3541a6db4461da28002d570bb1b6a44a4e3aea4
+  metadata.gz: efee6e93050df75436a04773199c6721fcfd8b39f979ecbc8092df3717ba5b88086e0351ebf5a8def6b21e0adbaf991002a33748f473134cb7e95361ffea4809
+  data.tar.gz: 720fa6125c7ac803473e0c604c37ca6829ba68ca89bcb539821d5b660e7c81828c72f01726e1c21307b63552bb36e05477c4fc973999583e2e19629fd086631e

data/Gemfile.lock

@@ -1,7 +1,7 @@
 PATH
   remote: .
   specs:
-    geoprojection (0.1.0)
+    geoprojection (0.1.1)
 
 GEM
   remote: https://rubygems.org/

data/doc/example.rb

@@ -0,0 +1,18 @@
+# frozen_string_literal: true
+
+require 'gnuplot'
+require 'geoprojection'
+
+ellipse = Geoprojection::Ellipse.new(latitude: 50.8169, longitude: 0.1367, distance: 10)
+
+Gnuplot.open do |gp|
+  Gnuplot::Plot.new(gp) do |plot|
+    plot.title  'Ellipsis'
+    plot.xlabel 'x'
+    plot.ylabel 'y'
+    plot.data << Gnuplot::DataSet.new(ellipse.points.transpose) do |ds|
+      ds.with = 'points'
+      ds.linewidth = 4
+    end
+  end
+end

data/lib/geoprojection/ellipse.rb

@@ -3,9 +3,9 @@
 module Geoprojection
   # Elliptical projection of a center point (lat, long) with a provided distance from the center.
   class Ellipse
-    EQUATOR_KM_PER_LATITUDE_DEGREE = 110.567
-    EQUATOR_KM_PER_LONGITUDE_DEGREE = 111.321
     DEGREES_PER_RADIAN = Math::PI / 180
+    EQUATOR_KM_PER_LATITUDE_RADIAN = 110.567 * DEGREES_PER_RADIAN
+    EQUATOR_KM_PER_LONGITUDE_RADIAN = 111.321 * DEGREES_PER_RADIAN
 
     def initialize(latitude:, longitude:, distance:, points: 36)
       @latitude = latitude
@@ -33,7 +33,7 @@ module Geoprojection
     end
 
     def derived_y(theta)
-      @latitude + (DEGREES_PER_RADIAN * @distance * latitude_distortion * Math.sin(DEGREES_PER_RADIAN * theta))
+      @latitude + (DEGREES_PER_RADIAN * @distance * Math.sin(DEGREES_PER_RADIAN * theta))
     end
 
     def step
@@ -41,11 +41,7 @@ module Geoprojection
     end
 
     def longitude_distortion
-      @longitude_distortion ||= 1 / (EQUATOR_KM_PER_LONGITUDE_DEGREE * Math.cos(DEGREES_PER_RADIAN * @latitude))
-    end
-
-    def latitude_distortion
-      @latitude_distortion ||= 1 / (EQUATOR_KM_PER_LATITUDE_DEGREE / @latitude)
+      @longitude_distortion ||= 1 / (EQUATOR_KM_PER_LONGITUDE_RADIAN * Math.cos(DEGREES_PER_RADIAN * @latitude))
     end
   end
 end

data/lib/geoprojection/version.rb

@@ -1,5 +1,5 @@
 # frozen_string_literal: true
 
 module Geoprojection
-  VERSION = '0.1.0'
+  VERSION = '0.1.1'
 end

data/test.rb

@@ -0,0 +1,18 @@
+# frozen_string_literal: true
+
+require 'gnuplot'
+require 'geoprojection'
+
+ellipse = Geoprojection::Ellipse.new(latitude: 50.8169, longitude: 0.1367, distance: 10)
+
+Gnuplot.open do |gp|
+  Gnuplot::Plot.new(gp) do |plot|
+    plot.title  'Ellipse'
+    plot.xlabel 'x'
+    plot.ylabel 'y'
+    plot.data << Gnuplot::DataSet.new(ellipse.points.transpose) do |ds|
+      ds.with = 'points'
+      ds.linewidth = 4
+    end
+  end
+end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: geoprojection
 version: !ruby/object:Gem::Version
-  version: 0.1.0
+  version: 0.1.1
 platform: ruby
 authors:
 - Bob Farrell
@@ -25,10 +25,12 @@ files:
 - Makefile
 - README.md
 - Rakefile
+- doc/example.rb
 - lib/geoprojection.rb
 - lib/geoprojection/ellipse.rb
 - lib/geoprojection/version.rb
 - sig/geoprojection.rbs
+- test.rb
 homepage: https://github.com/bobf/geoprojection
 licenses: []
 metadata: