blue_geo 0.0.1

data/.gitignore

@@ -0,0 +1,4 @@
+*.gem
+.bundle
+Gemfile.lock
+pkg/*

data/Gemfile

@@ -0,0 +1,4 @@
+source "http://rubygems.org"
+
+# Specify your gem's dependencies in blue_geo.gemspec
+gemspec

data/README.md

@@ -0,0 +1,5 @@
+blue_geo
+
+======================================
+
+Converts easting and northing into latitude and longitude

data/Rakefile

@@ -0,0 +1,8 @@
+require "bundler/gem_tasks"
+
+# Rake::TestTask.new do |t|
+# 	t.libs << 'test'
+# end
+
+# desc "Run tests"
+# task :default => :test

data/blue_geo.gemspec

@@ -0,0 +1,24 @@
+# -*- encoding: utf-8 -*-
+$:.push File.expand_path("../lib", __FILE__)
+require "blue_geo/version"
+
+Gem::Specification.new do |s|
+  s.name        = "blue_geo"
+  s.version     = BlueGeo::VERSION
+  s.authors     = ["Arti Sinani"]
+  s.email       = ["artisinani@gmail.com"]
+  s.homepage    = ""
+  s.summary     = %q{Converts easting and northing into latitude and longitude}
+  s.description = %q{Converts easting and northing into latitude and longitude}
+
+  s.rubyforge_project = "blue_geo"
+
+  s.files         = `git ls-files`.split("\n")
+  s.test_files    = `git ls-files -- {test,spec,features}/*`.split("\n")
+  s.executables   = `git ls-files -- bin/*`.split("\n").map{ |f| File.basename(f) }
+  s.require_paths = ["lib"]
+
+  # specify any dependencies here; for example:
+  # s.add_development_dependency "rspec"
+  # s.add_runtime_dependency "rest-client"
+end

data/lib/blue_geo.rb

@@ -0,0 +1,8 @@
+require_relative "blue_geo/version"
+require_relative "blue_geo/easting_northing_to_latitude_longitude_converter"
+
+module BlueGeo
+  def self.easting_northing_to_lat_lon(easting, northing)
+  	EastingNorthingToLatitudeLongitudeConverter.new.convert(easting, northing)
+  end
+end

data/lib/blue_geo/easting_northing_to_latitude_longitude_converter.rb

@@ -0,0 +1,70 @@
+module BlueGeo
+  class EastingNorthingToLatitudeLongitudeConverter
+  	DEG2RAD  = Math::PI / 180.0
+  	RAD2DEG  = 180.0 / Math::PI
+  	WGS_AXIS = 6378137
+  	LAT_ORIG = 49 * DEG2RAD
+  	LON_ORIG = -2 * DEG2RAD
+  	SCALE    = 0.9996012717
+    FALSE_EASTING 	= 400000
+    FALSE_NORTHING 	= -100000
+    AXIS 				= 6377563.396
+    WGS_ECCENT 	= 0.00669438
+    ECCENT 			= 0.00667054 
+
+
+  	def convert(easting, northing)
+      p       = ECCENT / 8
+      a       = AXIS * (1 - (2 * p) - (3 * p ** 2) - (10 * p ** 3))
+      b       = AXIS * ((6 * p) + (12 * p ** 2) + (45 * p ** 3 )) / 2
+      c       = AXIS * ((15 * p ** 2) + (90 * p ** 3)) / 4
+      mo      = a * LAT_ORIG - b * Math.sin(2 * LAT_ORIG) + c * Math.sin(4 * LAT_ORIG)
+      mp      = mo + ((northing - FALSE_NORTHING) / SCALE)
+      phidash = mp / a
+      phif    = phidash + ((b * Math.sin(2 * phidash)) - (c * Math.sin(4 * phidash))) / (a - (2 * b * Math.cos(2 * phidash)))
+
+      uf = AXIS / Math.sqrt(1 - (ECCENT * (Math.sin(phif) * Math.sin(phif))))
+      h  = (easting - FALSE_EASTING) / (SCALE * uf)
+
+      nsqd = ECCENT * (Math.cos(phif) * Math.cos(phif)) / (1 - ECCENT)
+      tsqd = Math.sin(phif) / Math.cos(phif) ** 2
+
+      lambdap = LON_ORIG + (1 / Math.cos(phif)) * ((h - ((h * h * h) / 6) * (1 + (2 * tsqd) + nsqd)))
+      phip = phif - ((1 + nsqd) * (Math.sin(phif) / Math.cos(phif)) * (((h * h) / 2) - ((h * h * h * h) / 24) * (5 + 3 * tsqd)))
+
+      x, y, z = lat_lon_to_cartesian(phip, lambdap, AXIS, ECCENT, 0)
+
+      lat, lon = cartesian_to_lat_lon(x, y, z, ECCENT, WGS_AXIS)
+
+	    [lat * RAD2DEG, lon * RAD2DEG]
+
+  	end
+
+  	def lat_lon_to_cartesian(lat, lon, a, eccentsq, h)
+  		v = a / (Math.sqrt(1 - eccentsq * (Math.sin(lat) ** 2)))
+      x = (v + h) * Math.cos(lat) * Math.cos(lon)
+      y = (v + h) * Math.cos(lat) * Math.sin(lon)
+      z = ((1 - eccentsq) * v + h) * Math.sin(lat)
+      [x, y, z]
+  	end
+
+  	def cartesian_to_lat_lon(x, y, z, eccentsq, a)
+	    lon      = Math.atan(y / x)
+	    p        = Math.sqrt((x * x) + (y * y))
+	    lat      = Math.atan(z / (p * (1 - eccentsq)))
+	    v        = a / (Math.sqrt(1 - eccentsq * (Math.sin(lat) * Math.sin(lat))))
+	    errvalue = 1.0
+	    lat      = 0
+	    lat0     = 0
+
+	    while errvalue > 0.001
+	    	lat0 = Math.atan((z + eccentsq * v * Math.sin(lat)) / p)
+        errvalue = (lat0 - lat).abs;
+        lat = lat0;
+	    end
+
+	    [lat, lon]
+    end
+
+  end
+end

data/lib/blue_geo/version.rb

@@ -0,0 +1,3 @@
+module BlueGeo
+  VERSION = "0.0.1"
+end

data/test/blue_geo_convertion_test.rb

@@ -0,0 +1,25 @@
+require 'minitest/autorun'
+require File.expand_path("../../lib/blue_geo",  __FILE__)
+
+module BlueGeo
+  class ConvertionTest < MiniTest::Unit::TestCase
+  	# 536861, 176584 -> 51.471150	-0.029269 (easting_north -> lat, lon)
+  	# 	51.47168	-0.03088506	536861	176584
+		# 53.75771	-2.704654	353639	429270
+		# 51.66625	-3.92504	266963	198141
+		# 52.39906	0.2658032	554253	280309
+		# 51.48031	0.1799298	551473	177959
+		# 53.3513	-3.003981	333267	384294
+		# 53.52229	-2.488291	367724	402958
+		# 55.1022	-4.76847	223472	582144
+		# 52.60267	1.733953	652939	307091
+		# 51.88307	-0.5158798	502248	221541
+		# 51.43271	-0.1848878	526273	171974
+  	
+  	def test_converts_correctly
+  		latitude, longitude = BlueGeo.easting_northing_to_lat_lon(536861, 176584)
+  		assert_equal 51.471176444600296, latitude, 'latitude'
+  		assert_equal -0.02940177028571391, longitude, 'longitude'
+  	end
+	end
+end

metadata

@@ -0,0 +1,55 @@
+--- !ruby/object:Gem::Specification
+name: blue_geo
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+  prerelease: 
+platform: ruby
+authors:
+- Arti Sinani
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2013-02-17 00:00:00.000000000 Z
+dependencies: []
+description: Converts easting and northing into latitude and longitude
+email:
+- artisinani@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- .gitignore
+- Gemfile
+- README.md
+- Rakefile
+- blue_geo.gemspec
+- lib/blue_geo.rb
+- lib/blue_geo/easting_northing_to_latitude_longitude_converter.rb
+- lib/blue_geo/version.rb
+- test/blue_geo_convertion_test.rb
+homepage: ''
+licenses: []
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  none: false
+  requirements:
+  - - ! '>='
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  none: false
+  requirements:
+  - - ! '>='
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: blue_geo
+rubygems_version: 1.8.17
+signing_key: 
+specification_version: 3
+summary: Converts easting and northing into latitude and longitude
+test_files:
+- test/blue_geo_convertion_test.rb