app/models/geo_record.rb

   1 class GeoRecord < ActiveRecord::Base
   2   before_save :update_tile
   3
   4   # Is this node within -90 >= latitude >= 90 and -180 >= longitude >= 180
   5   # * returns true/false
   6   def in_world?
   7     return false if self.lat < -90 or self.lat > 90
   8     return false if self.lon < -180 or self.lon > 180
   9     return true
  10   end
  11
  12   def self.find_by_area(minlat, minlon, maxlat, maxlon, options)
  13     self.with_scope(:find => {:conditions => OSM.sql_for_area(minlat, minlon, maxlat, maxlon)}) do
  14       return self.find(:all, options)
  15     end
  16   end
  17
  18   def update_tile
  19     self.tile = QuadTile.tile_for_point(lat, lon)
  20   end
  21
  22   def lat=(l)
  23     self.latitude = (l * 10000000).round
  24   end
  25
  26   def lon=(l)
  27     self.longitude = (l * 10000000).round
  28   end
  29
  30   # Return WGS84 latitude
  31   def lat
  32     return self.latitude.to_f / 10000000
  33   end
  34
  35   # Return WGS84 longitude
  36   def lon
  37     return self.longitude.to_f / 10000000
  38   end
  39
  40   # Potlatch projections
  41   def lon_potlatch(baselong,masterscale)
  42     (self.lon-baselong)*masterscale
  43   end
  44
  45   def lat_potlatch(basey,masterscale)
  46     -(lat2y(self.lat)-basey)*masterscale
  47   end
  48
  49   private
  50
  51   def lat2y(a)
  52     180/Math::PI * Math.log(Math.tan(Math::PI/4+a*(Math::PI/180)/2))
  53   end
  54
  55 end
  56