simple_mercator_location 1.0.0

data/lib/simple_mercator_location.rb

@@ -0,0 +1,137 @@
+class SimpleMercatorLocation
+
+    attr_accessor :lat_deg, :lon_deg, :zoom
+
+    def initialize(*args)
+        @zoom = 2
+        @lat_deg = 0
+        @lon_deg = 0
+
+        if args.first.is_a?(Hash)
+            @lat_deg = args.first[:lat] if args.first.has_key?(:lat)
+            @lon_deg = args.first[:lon] if args.first.has_key?(:lon)
+            @zoom = args.first[:zoom] if args.first.has_key?(:zoom)
+        end
+    end 
+
+
+    # latitude in radiants
+    def lat_rad
+      shift_to_rad * lat_deg
+    end
+
+
+    # longitude in radiants
+    def lon_rad
+      shift_to_rad * lon_deg
+    end
+
+
+    # factor for scaling to radiants
+    def shift_to_rad
+      Rational(Math::PI, 180)
+    end
+
+
+    # factor for scaling to degrees
+    def shift_to_deg
+      Rational(180.0, Math::PI)
+    end
+
+
+    # size of our tiles
+    def tile_size
+      256
+    end
+
+
+    # earth radius in meters
+    def earth_radius 
+      6378137
+    end
+
+
+    # set the zoom-level and returns self
+    def zoom_at(scale)
+        @zoom = scale
+        self
+    end
+
+
+    # origin of google map lieves on top-left. For getting pixels,
+    # we assume displaying the earh an zoomm-level 0 and a square tile
+    # of tile_size * tile_size pixels. The origin of our wsg84 system
+    # (lat, lon = (0,0) sits and (128,128)px. 
+    def origin_px
+      [tile_size / 2, tile_size / 2]
+    end
+
+
+    # return the pixels per degree
+    def pixel_per_degree
+      Rational(tile_size, 360)
+    end
+
+
+    # return the pixels per radiant
+    def pixel_per_rad
+      tile_size / (2 * Math::PI) # tile_size / (360 * shift_to_rad)
+    end
+
+    
+    # scale the latitude via mercator projections
+    # see http://en.wikipedia.org/wiki/Mercator_projection#Derivation_of_the_Mercator_projection
+    # return the scaled latitude as radiant
+    def lat_scaled_rad
+      Math.log( Math.tan( Rational(Math::PI,4) + Rational(lat_rad,2) )) 
+    end
+
+
+    # returns the scaled latitude as degrees
+    def lat_scaled_deg
+      lat_scaled_rad * shift_to_deg
+    end
+
+
+    # returns the mercotors meters count for the location
+    def to_m
+        mx = earth_radius * lon_rad
+        my = earth_radius * lat_scaled_rad
+
+        my = (my.round(8) == 0)? 0 : my
+
+        [mx.to_f, my.to_f]
+    end
+
+
+    # calculates the google world-coordinates as described here:
+    # https://developers.google.com/maps/documentation/javascript/examples/map-coordinates
+    def to_w
+      px = origin_px.first +  pixel_per_rad * lon_rad
+      py = origin_px.last  -  pixel_per_rad * lat_scaled_rad
+
+      return [px, py]
+    end
+
+
+    # returns the pixel coordinates at the given zoom level
+    def to_px
+      tiles_count = 2**zoom
+      wx,wy = self.to_w
+      px = wx * tiles_count
+      py = wy * tiles_count
+      return [px.to_i, py.to_i]
+    end
+
+
+    # calculates the tile numbers for google maps at the given zoom level
+    def to_tile
+      px,py = self.to_px
+
+      tx = Rational(px, tile_size).to_i
+      ty = Rational(py, tile_size).to_i
+
+      return [tx,ty]
+    end
+
+end

data/spec/lib/simple_mercator_location_spec.rb

@@ -0,0 +1,130 @@
+require "spec_helper"
+
+describe SimpleMercatorLocation do
+
+  describe "accessors" do
+
+    let(:loc) { SimpleMercatorLocation.new }
+
+    attrs = [:lat_deg, :lon_deg, :zoom]
+
+    attrs.each do |attr|
+      it "should have getter and setter for #{attr}" do
+        loc.should respond_to("#{attr}".to_sym)
+        loc.should respond_to("#{attr}=".to_sym)
+        loc.send("#{attr}=".to_sym, "asd")
+        loc.send("#{attr}".to_sym).should eql "asd"
+      end
+    end
+  end
+
+
+  describe "#initialize" do
+    context "without args" do
+
+      let(:loc) { SimpleMercatorLocation.new }
+
+      it "returns a location object" do
+        loc.should be_an_instance_of(SimpleMercatorLocation)
+      end
+    end
+
+
+    context "given a hash with lat, lon and zoom" do
+      let(:loc) { SimpleMercatorLocation.new(lat: 1, lon: 2, zoom: 12) }
+
+      it "returns a location object" do
+        loc.should be_an_instance_of(SimpleMercatorLocation)
+      end
+
+      it "return a location object with lat and lon set" do
+        loc.lat_deg.should eql 1
+        loc.lon_deg.should eql 2
+        loc.zoom.should eql 12
+      end
+    end
+  end
+
+
+  describe "#zoom_at" do
+    let(:loc) { SimpleMercatorLocation.new }
+    it "sets the zoom scale and returns self" do
+      loc.zoom_at(15).should eql loc
+      loc.zoom_at(15).zoom.should eql 15
+    end
+  end
+
+  describe "to_m" do
+    context "given no args (so using lat = 0 and lon = 0)" do
+      let(:loc) { SimpleMercatorLocation.new }
+      it "should return origin" do
+        loc.to_m.should eql  [0.0,0.0]
+      end
+    end
+
+    context "given a location" do
+      places =
+        [
+          { lat: 48.16608541901253, lon: 11.6455078125, mx: 1296371.99971659, my: 6134530.14205511 },
+          { lat: 48.10743118848038, lon: 11.42578125,   mx: 1271912.15066533, my: 6124746.20243460 },
+          { lat: 48.22467264956519, lon: 12.12890625,   mx: 1350183.66762935, my: 6144314.08167561 },
+          { lat: 0,                 lon: 0,             mx: 0.0,              my: 0.0 },
+        ]
+      places.each do |place|
+        it "calculates the mercator projection of (lat: #{place[:lat]}, lon: #{place[:lon]}) to (meters x: #{place[:mx]}, meters y: #{place[:my]})" do
+          meters = SimpleMercatorLocation.new(lon: place[:lon], lat: place[:lat]).to_m
+          meters.map!{|m| m.round(8) }.should eql([place[:mx], place[:my]])
+        end
+      end
+    end
+  end
+
+  describe "to_w" do
+    places =
+        [
+          { lat: 41.850033, lon: -87.65005229999997, wx: 65.67107392000001, wy: 95.1748950436046 },
+          #{ lat: 48.10743118848038, lon: 11.42578125,   wx: 1271912.15066533, wy: 6124746.20243460 },
+          #{ lat: 48.22467264956519, lon: 12.12890625,   wx: 1350183.66762935, wy: 6144314.08167561 },
+          #{ lat: 0,                 lon: 0,             wx: 0.0,              wy: 0.0 },
+        ]
+      places.each do |place|
+        it "calculates the mercator projection of (lat: #{place[:lat]}, lon: #{place[:lon]}) to (world coordinate x: #{place[:wx]}, world coordinate y: #{place[:wy]})" do
+          SimpleMercatorLocation.new(lon: place[:lon], lat: place[:lat]).to_w.should eql([place[:wx], place[:wy]])
+        end
+      end
+
+
+  end
+
+
+  describe "to_px" do
+    places =
+      [
+        { zoom: 11, lat: 49.38237278700955, lon: 8.61328125, px: 274688,  py: 179200  },
+        { zoom: 14, lat: 49.38237278700955, lon: 8.61328125, px: 2197504, py: 1433600 },
+      ]
+
+    places.each do |place|
+      it "calculates the pixels of (lat: #{place[:lat]}, lon: #{place[:lon]}) to (px: #{place[:px]}, py: #{place[:py]})" do
+        SimpleMercatorLocation.new(lon: place[:lon], lat: place[:lat], zoom: place[:zoom]).to_px.should eql([place[:px], place[:py]])
+      end
+    end
+  end
+
+
+  describe "to_tile" do
+    places =
+      [
+        { zoom: 11, lat: 49.412758, lon: 8.671938,   tx: 1073,  ty: 699 },
+        { zoom: 11, lat: 40.689359, lon: -74.045197, tx: 602,   ty: 770  },
+        { zoom: 15, lat: 40.689359, lon: -74.045197, tx: 9644,  ty: 12322 },
+      ]
+
+    places.each do |place|
+      it "calculates the tile of (lat: #{place[:lat]}, lon: #{place[:lon]}) to (tile x: #{place[:tx]}, tile y: #{place[:ty]})" do
+        SimpleMercatorLocation.new(lon: place[:lon], lat: place[:lat], zoom: place[:zoom]).to_tile.should eql([place[:tx], place[:ty]])
+      end
+    end
+  end
+
+end

data/spec/spec_helper.rb

@@ -0,0 +1,19 @@
+require "rspec"
+require "awesome_print"
+
+RSpec.configure do |config|
+  config.color_enabled = true
+  config.filter_run :focus => true
+  config.run_all_when_everything_filtered = true
+  config.treat_symbols_as_metadata_keys_with_true_values = true
+end
+
+def timed(name)
+  start = Time.now
+  puts "\n[STARTED: #{name}]"
+  yield if block_given?
+  finish = Time.now
+  puts "[FINISHED: #{name} in #{(finish - start) * 1000} milliseconds]"
+end
+
+require "simple_mercator_location"

metadata

@@ -0,0 +1,66 @@
+--- !ruby/object:Gem::Specification
+name: simple_mercator_location
+version: !ruby/object:Gem::Version
+  version: 1.0.0
+  prerelease: 
+platform: ruby
+authors:
+- Roman Lehnert
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2013-06-29 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: rspec
+  requirement: !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ~>
+      - !ruby/object:Gem::Version
+        version: '2.5'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    none: false
+    requirements:
+    - - ~>
+      - !ruby/object:Gem::Version
+        version: '2.5'
+description: Converts WSG84 Coordinates via Mercator-projection to meters and tiles
+email: roman.lehnert@googlemail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- lib/simple_mercator_location.rb
+- spec/lib/simple_mercator_location_spec.rb
+- spec/spec_helper.rb
+homepage: https://github.com/romanlehnert/simple_mercator_location
+licenses:
+- MIT
+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: 
+rubygems_version: 1.8.25
+signing_key: 
+specification_version: 3
+summary: A tiny lib for the mercator projecton
+test_files:
+- spec/lib/simple_mercator_location_spec.rb
+- spec/spec_helper.rb