geo_swap 0.0.1 → 0.0.2

data/geo_swap.gemspec

@@ -18,4 +18,5 @@ Gem::Specification.new do |gem|
   gem.require_paths = ["lib"]
 
   gem.add_development_dependency 'rake', '10.0.3'
+  gem.add_development_dependency 'rspec', '2.12.0'
 end

data/lib/geo_swap/utilities.rb

@@ -0,0 +1,10 @@
+
+module GeoSwap::Utilities
+
+  def degrees_to_radians(degrees)
+    degrees * (Math::PI / 180)
+  end
+
+  module_function :degrees_to_radians
+
+end

data/lib/geo_swap/version.rb

@@ -1,3 +1,3 @@
 module GeoSwap
-  VERSION = "0.0.1"
+  VERSION = "0.0.2"
 end

data/lib/geo_swap/zone.rb

@@ -0,0 +1,19 @@
+module GeoSwap
+  class Zone
+    def initialize(longitude)
+      @longitude = longitude.to_f
+    end
+
+    def number
+      @number ||= begin
+        number = ((@longitude + 180.0) / 6.0).floor + 1
+        #special case for longitude 180.0
+        (number > 60) ? 60 : number
+      end
+    end
+
+    def origin
+      @origin ||= ((number - 1) * 6) - 177
+    end
+  end
+end

data/lib/geo_swap.rb

@@ -1,5 +1,74 @@
 require "geo_swap/version"
+require 'geo_swap/zone'
+require 'geo_swap/utilities'
 
 module GeoSwap
-  # Your code goes here...
+  extend Utilities
+
+  def lat_long_to_utm(lat, long)
+    validate_range(lat, long)
+    lat_radians = degrees_to_radians(lat)
+    long_radians = degrees_to_radians(long)
+
+    zone = Zone.new(long)
+    origin_radians = degrees_to_radians(zone.origin)
+
+    equator_factor = (EQUATORIAL_RADIUS / Math.sqrt(1 - (ECC_SQUARED * (Math.sin(lat_radians) ** 2)))
+    squared_lat_tangent = Math.tan(lat_radians) ** 2
+    ecc_prime_factor = ECC_PRIME_SQUARED * (Math.cos(lat_radians) ** 2)
+    origin_factor = Math.cos(lat_radians) * (long_radians - origin_radians)
+
+
+    ecc_1 = (1 - ecc(1, 1, 4) - ecc(2, 3, 64) - ecc(3, 5, 256))
+    ecc_2 = ecc(1, 3, 8) + ecc(2, 3, 32) + ecc(3, 45, 1024)
+    ecc_3 = ecc(2, 15, 256) + ecc(3, 45, 1024)
+    ecc_4 = ecc(3, 35, 3072)
+
+    latRad_1 = lat_radians
+    latRad_2 = Math.sin(2 * lat_radians)
+    latRad_3 = Math.sin(4 * lat_radians)
+    latRad_4 = Math.sin(6 * lat_radians)
+
+    northing_factor = EQUATORIAL_RADIUS * (
+      ecc_1 * latRad_1 -
+      ecc_2 * latRad_2 +
+      ecc_3 * latRad_3 -
+      ecc_4 * latRad_4
+    )
+
+  end
+
+  module_function :lat_long_to_utm
+
+
+  private
+
+  MAX_LATITUDE = 90.0
+  MIN_LATITUDE = -90.0
+  MAX_LONGITUDE = 180
+  MIN_LONGITUDE = -180
+  MAX_CONVERSION_LATITUDE = 84.0
+  MIN_CONVERSION_LATITUDE = -80.0
+
+  EQUATORIAL_RADIUS = 6378137.0
+  ECC_SQUARED = 0.006694380023
+  ECC_PRIME_SQUARED = ECC_SQUARED / (1 - ECC_SQUARED)
+
+  def self.validate_range(lat, long)
+    unless lat.between?(MIN_LATITUDE, MAX_LATITUDE) && long.between?(MIN_LONGITUDE, MAX_LONGITUDE)
+      raise InputError, 'Input coordinates are invalid'
+    end
+
+    unless lat.between?(MIN_CONVERSION_LATITUDE, MAX_CONVERSION_LATITUDE)
+      raise InputError, 'Conversions are unreliable close to the polar regions'
+    end
+  end
+
+  def self.ecc(power, numerator, denominator)
+    ecc = ECC_SQUARED ** power
+    (numerator * ecc) / denominator
+  end
+
 end
+
+class GeoSwap::InputError < StandardError; end

data/spec/geo_swap/utilities_spec.rb

@@ -0,0 +1,15 @@
+require_relative '../../lib/geo_swap/utilities'
+
+module GeoSwap
+  describe Utilities do
+    describe 'degrees_to_radians' do
+      DEG_RAD_DATA = [[180, Math::PI], [0, 0], [360, (2 * Math::PI)]]
+
+      it 'converts numbers accurately' do
+        DEG_RAD_DATA.each do |(deg, rad)|
+          Utilities.degrees_to_radians(deg).should == rad
+        end
+      end
+    end
+  end
+end

data/spec/geo_swap/zone_spec.rb

@@ -0,0 +1,45 @@
+require_relative '../../lib/geo_swap/zone'
+
+module GeoSwap
+  describe Zone do
+
+
+    describe 'finding number from lat long' do
+
+
+      ZONE_DATA = [
+        { long: -180.0, zone_number: 1 },
+        { long: -147.0, zone_number: 6 },
+        { long: -139.0, zone_number: 7 },
+        { long: 0.0, zone_number: 31 },
+        { long: 62.0, zone_number: 41 },
+        { long: 71.0, zone_number: 42 },
+        { long: 110.0, zone_number: 49 },
+        { long: 141.0, zone_number: 54 },
+        { long: 180.0, zone_number: 60 },
+      ]
+
+      it 'applys the zone formula correctly' do
+        ZONE_DATA.each do |data|
+          Zone.new(data[:long]).number.should == data[:zone_number]
+        end
+      end
+
+    end
+
+    describe 'determining the zone origin' do
+      it 'can migrate from too far left' do
+        Zone.new(-180).origin.should == -177
+      end
+
+      it 'can migrate from too far right' do
+        Zone.new(180).origin.should == 177
+      end
+
+      it 'doesnt change the value when the long is already the origin' do
+        Zone.new(3).origin.should == 3
+      end
+    end
+
+  end
+end

data/spec/geo_swap_spec.rb

@@ -0,0 +1,49 @@
+require_relative '../lib/geo_swap'
+
+module GeoSwap
+
+  describe GeoSwap do
+    DATA = [
+      {
+        lat: 29.979175,
+        long: 31.1343583,
+        utm: "36R 320010mE 3317942mN",
+      },
+      {
+        lat: 41.8901694,
+        long: 12.4922694,
+        utm: "33T 291952mE 4640623mN",
+      }
+    ]
+
+    describe 'lat_long_to_utm' do
+      def conversion(lat, long)
+        GeoSwap.lat_long_to_utm(lat, long)
+      end
+
+      def check_error(lat, long, message)
+        expect { conversion(lat, long) }.
+          to raise_error(InputError, message)
+      end
+
+      it 'rejects data too close to the north pole' do
+        check_error(84.1, 0.0, 'Conversions are unreliable close to the polar regions')
+      end
+
+      it 'rejects data too close to the south pole' do
+        check_error(-80.1, 0.0, 'Conversions are unreliable close to the polar regions')
+      end
+
+      it 'rejects data outside the possible lat/long range' do
+        check_error(0.0,200.0, 'Input coordinates are invalid')
+      end
+
+      it 'correctly converts several reference points' do
+        DATA.each do |data|
+          GeoSwap.lat_long_to_utm(data[:lat], data[:long]).should == data[:utm]
+        end
+      end
+    end
+  end
+
+end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: geo_swap
 version: !ruby/object:Gem::Version
-  version: 0.0.1
+  version: 0.0.2
   prerelease: 
 platform: ruby
 authors:
@@ -27,6 +27,22 @@ dependencies:
     - - '='
       - !ruby/object:Gem::Version
         version: 10.0.3
+- !ruby/object:Gem::Dependency
+  name: rspec
+  requirement: !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - '='
+      - !ruby/object:Gem::Version
+        version: 2.12.0
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - '='
+      - !ruby/object:Gem::Version
+        version: 2.12.0
 description: Simple utility functions for converting between coordinate systems (Lat/Long,
   UTM, USNG)
 email:
@@ -43,7 +59,12 @@ files:
 - Rakefile
 - geo_swap.gemspec
 - lib/geo_swap.rb
+- lib/geo_swap/utilities.rb
 - lib/geo_swap/version.rb
+- lib/geo_swap/zone.rb
+- spec/geo_swap/utilities_spec.rb
+- spec/geo_swap/zone_spec.rb
+- spec/geo_swap_spec.rb
 homepage: ''
 licenses: []
 post_install_message: 
@@ -58,7 +79,7 @@ required_ruby_version: !ruby/object:Gem::Requirement
       version: '0'
       segments:
       - 0
-      hash: 3296913527661097903
+      hash: -10926259193031771
 required_rubygems_version: !ruby/object:Gem::Requirement
   none: false
   requirements:
@@ -67,7 +88,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
       version: '0'
       segments:
       - 0
-      hash: 3296913527661097903
+      hash: -10926259193031771
 requirements: []
 rubyforge_project: 
 rubygems_version: 1.8.23
@@ -75,4 +96,7 @@ signing_key:
 specification_version: 3
 summary: Simple utility functions for converting between coordinate systems (Lat/Long,
   UTM, USNG)
-test_files: []
+test_files:
+- spec/geo_swap/utilities_spec.rb
+- spec/geo_swap/zone_spec.rb
+- spec/geo_swap_spec.rb