Wrap longitude values explicitly
[rails.git] / lib / osm.rb
1 # The OSM module provides support functions for OSM.
2 module OSM
3
4   require 'time'
5   require 'rexml/parsers/sax2parser'
6   require 'rexml/text'
7   require 'xml/libxml'
8   require 'digest/md5'
9
10   if defined?(SystemTimer)
11     Timer = SystemTimer
12   else
13     require 'timeout'
14     Timer = Timeout
15   end
16
17   # The base class for API Errors.
18   class APIError < RuntimeError
19     def status
20       :internal_server_error
21     end
22
23     def to_s
24       "Generic API Error"
25     end
26   end
27
28   # Raised when an API object is not found.
29   class APINotFoundError < APIError
30     def status
31       :not_found
32     end
33
34     def to_s
35       "Object not found"
36     end
37   end
38
39   # Raised when a precondition to an API action fails sanity check.
40   class APIPreconditionFailedError < APIError
41     def initialize(message = "")
42       @message = message
43     end
44
45     def status
46       :precondition_failed
47     end
48
49     def to_s
50       "Precondition failed: #{@message}"
51     end
52   end
53
54   # Raised when to delete an already-deleted object.
55   class APIAlreadyDeletedError < APIError
56     def initialize(object = "object", object_id = "")
57       @object, @object_id = object, object_id
58     end
59
60     attr_reader :object, :object_id
61
62     def status
63       :gone
64     end
65
66     def to_s
67       "The #{object} with the id #{object_id} has already been deleted"
68     end
69   end
70
71   # Raised when the user logged in isn't the same as the changeset
72   class APIUserChangesetMismatchError < APIError
73     def status
74       :conflict
75     end
76
77     def to_s
78       "The user doesn't own that changeset"
79     end
80   end
81
82   # Raised when the changeset provided is already closed
83   class APIChangesetAlreadyClosedError < APIError
84     def initialize(changeset)
85       @changeset = changeset
86     end
87
88     attr_reader :changeset
89
90     def status
91       :conflict
92     end
93
94     def to_s
95       "The changeset #{@changeset.id} was closed at #{@changeset.closed_at}"
96     end
97   end
98
99   # Raised when a change is expecting a changeset, but the changeset doesn't exist
100   class APIChangesetMissingError < APIError
101     def status
102       :conflict
103     end
104
105     def to_s
106       "You need to supply a changeset to be able to make a change"
107     end
108   end
109
110   # Raised when a diff is uploaded containing many changeset IDs which don't match
111   # the changeset ID that the diff was uploaded to.
112   class APIChangesetMismatchError < APIError
113     def initialize(provided, allowed)
114       @provided, @allowed = provided, allowed
115     end
116
117     def status
118       :conflict
119     end
120
121     def to_s
122       "Changeset mismatch: Provided #{@provided} but only #{@allowed} is allowed"
123     end
124   end
125
126   # Raised when a diff upload has an unknown action. You can only have create,
127   # modify, or delete
128   class APIChangesetActionInvalid < APIError
129     def initialize(provided)
130       @provided = provided
131     end
132
133     def status
134       :bad_request
135     end
136
137     def to_s
138       "Unknown action #{@provided}, choices are create, modify, delete"
139     end
140   end
141
142   # Raised when bad XML is encountered which stops things parsing as
143   # they should.
144   class APIBadXMLError < APIError
145     def initialize(model, xml, message="")
146       @model, @xml, @message = model, xml, message
147     end
148
149     def status
150       :bad_request
151     end
152
153     def to_s
154       "Cannot parse valid #{@model} from xml string #{@xml}. #{@message}"
155     end
156   end
157
158   # Raised when the provided version is not equal to the latest in the db.
159   class APIVersionMismatchError < APIError
160     def initialize(id, type, provided, latest)
161       @id, @type, @provided, @latest = id, type, provided, latest
162     end
163
164     attr_reader :provided, :latest, :id, :type
165
166     def status
167       :conflict
168     end
169
170     def to_s
171       "Version mismatch: Provided #{provided}, server had: #{latest} of #{type} #{id}"
172     end
173   end
174
175   # raised when a two tags have a duplicate key string in an element.
176   # this is now forbidden by the API.
177   class APIDuplicateTagsError < APIError
178     def initialize(type, id, tag_key)
179       @type, @id, @tag_key = type, id, tag_key
180     end
181
182     attr_reader :type, :id, :tag_key
183
184     def status
185       :bad_request
186     end
187
188     def to_s
189       "Element #{@type}/#{@id} has duplicate tags with key #{@tag_key}"
190     end
191   end
192
193   # Raised when a way has more than the configured number of way nodes.
194   # This prevents ways from being to long and difficult to work with
195   class APITooManyWayNodesError < APIError
196     def initialize(id, provided, max)
197       @id, @provided, @max = id, provided, max
198     end
199
200     attr_reader :id, :provided, :max
201
202     def status
203       :bad_request
204     end
205
206     def to_s
207       "You tried to add #{provided} nodes to way #{id}, however only #{max} are allowed"
208     end
209   end
210
211   ##
212   # raised when user input couldn't be parsed
213   class APIBadUserInput < APIError
214     def initialize(message)
215       @message = message
216     end
217
218     def status
219       :bad_request
220     end
221
222     def to_s
223       @message
224     end
225   end
226
227   ##
228   # raised when bounding box is invalid
229   class APIBadBoundingBox < APIError
230     def initialize(message)
231       @message = message
232     end
233
234     def status
235       :bad_request
236     end
237
238     def to_s
239       @message
240     end
241   end
242
243   ##
244   # raised when an API call is made using a method not supported on that URI
245   class APIBadMethodError < APIError
246     def initialize(supported_method)
247       @supported_method = supported_method
248     end
249
250     def status
251       :method_not_allowed
252     end
253
254     def to_s
255       "Only method #{@supported_method} is supported on this URI"
256     end
257   end
258
259   ##
260   # raised when an API call takes too long
261   class APITimeoutError < APIError
262     def status
263       :request_timeout
264     end
265
266     def to_s
267       "Request timed out"
268     end
269   end
270
271   ##
272   # raised when someone tries to redact a current version of
273   # an element - only historical versions can be redacted.
274   class APICannotRedactError < APIError
275     def status
276       :bad_request
277     end
278
279     def to_s
280       "Cannot redact current version of element, only historical versions may be redacted."
281     end
282   end
283
284   # Helper methods for going to/from mercator and lat/lng.
285   class Mercator
286     include Math
287
288     #init me with your bounding box and the size of your image
289     def initialize(min_lat, min_lon, max_lat, max_lon, width, height)
290       xsize = xsheet(max_lon) - xsheet(min_lon)
291       ysize = ysheet(max_lat) - ysheet(min_lat)
292       xscale = xsize / width
293       yscale = ysize / height
294       scale = [xscale, yscale].max
295
296       xpad = width * scale - xsize
297       ypad = height * scale - ysize
298
299       @width = width
300       @height = height
301
302       @tx = xsheet(min_lon) - xpad / 2
303       @ty = ysheet(min_lat) - ypad / 2
304
305       @bx = xsheet(max_lon) + xpad / 2
306       @by = ysheet(max_lat) + ypad / 2
307     end
308
309     #the following two functions will give you the x/y on the entire sheet
310
311     def ysheet(lat)
312       log(tan(PI / 4 + (lat * PI / 180 / 2))) / (PI / 180)
313     end
314
315     def xsheet(lon)
316       lon
317     end
318
319     #and these two will give you the right points on your image. all the constants can be reduced to speed things up. FIXME
320
321     def y(lat)
322       return @height - ((ysheet(lat) - @ty) / (@by - @ty) * @height)
323     end
324
325     def x(lon)
326       return  ((xsheet(lon) - @tx) / (@bx - @tx) * @width)
327     end
328   end
329
330   class GreatCircle
331     include Math
332
333     # initialise with a base position
334     def initialize(lat, lon)
335       @lat = lat * PI / 180
336       @lon = lon * PI / 180
337     end
338
339     # get the distance from the base position to a given position
340     def distance(lat, lon)
341       lat = lat * PI / 180
342       lon = lon * PI / 180
343       return 6372.795 * 2 * asin(sqrt(sin((lat - @lat) / 2) ** 2 + cos(@lat) * cos(lat) * sin((lon - @lon)/2) ** 2))
344     end
345
346     # get the worst case bounds for a given radius from the base position
347     def bounds(radius)
348       latradius = 2 * asin(sqrt(sin(radius / 6372.795 / 2) ** 2))
349
350       begin
351         lonradius = 2 * asin(sqrt(sin(radius / 6372.795 / 2) ** 2 / cos(@lat) ** 2))
352       rescue Errno::EDOM
353         lonradius = PI
354       end
355
356       minlat = (@lat - latradius) * 180 / PI
357       maxlat = (@lat + latradius) * 180 / PI
358       minlon = (@lon - lonradius) * 180 / PI
359       maxlon = (@lon + lonradius) * 180 / PI
360
361       return { :minlat => minlat, :maxlat => maxlat, :minlon => minlon, :maxlon => maxlon }
362     end
363
364     # get the SQL to use to calculate distance
365     def sql_for_distance(lat_field, lon_field)
366       "6372.795 * 2 * asin(sqrt(power(sin((radians(#{lat_field}) - #{@lat}) / 2), 2) + cos(#{@lat}) * cos(radians(#{lat_field})) * power(sin((radians(#{lon_field}) - #{@lon})/2), 2)))"
367     end
368   end
369
370   class GeoRSS
371     def initialize(feed_title='OpenStreetMap GPS Traces', feed_description='OpenStreetMap GPS Traces', feed_url='http://www.openstreetmap.org/traces/')
372       @doc = XML::Document.new
373       @doc.encoding = XML::Encoding::UTF_8
374
375       rss = XML::Node.new 'rss'
376       @doc.root = rss
377       rss['version'] = "2.0"
378       rss['xmlns:geo'] = "http://www.w3.org/2003/01/geo/wgs84_pos#"
379       @channel = XML::Node.new 'channel'
380       rss << @channel
381       title = XML::Node.new 'title'
382       title <<  feed_title
383       @channel << title
384       description_el = XML::Node.new 'description'
385       @channel << description_el
386
387       description_el << feed_description
388       link = XML::Node.new 'link'
389       link << feed_url
390       @channel << link
391       image = XML::Node.new 'image'
392       @channel << image
393       url = XML::Node.new 'url'
394       url << 'http://www.openstreetmap.org/images/mag_map-rss2.0.png'
395       image << url
396       title = XML::Node.new 'title'
397       title << "OpenStreetMap"
398       image << title
399       width = XML::Node.new 'width'
400       width << '100'
401       image << width
402       height = XML::Node.new 'height'
403       height << '100'
404       image << height
405       link = XML::Node.new 'link'
406       link << feed_url
407       image << link
408     end
409
410     def add(latitude=0, longitude=0, title_text='dummy title', author_text='anonymous', url='http://www.example.com/', description_text='dummy description', timestamp=DateTime.now)
411       item = XML::Node.new 'item'
412
413       title = XML::Node.new 'title'
414       item << title
415       title << title_text
416       link = XML::Node.new 'link'
417       link << url
418       item << link
419
420       guid = XML::Node.new 'guid'
421       guid << url
422       item << guid
423
424       description = XML::Node.new 'description'
425       description << description_text
426       item << description
427
428       author = XML::Node.new 'author'
429       author << author_text
430       item << author
431
432       pubDate = XML::Node.new 'pubDate'
433       pubDate << timestamp.to_s(:rfc822)
434       item << pubDate
435
436       if latitude
437         lat_el = XML::Node.new 'geo:lat'
438         lat_el << latitude.to_s
439         item << lat_el
440       end
441
442       if longitude
443         lon_el = XML::Node.new 'geo:long'
444         lon_el << longitude.to_s
445         item << lon_el
446       end
447
448       @channel << item
449     end
450
451     def to_s
452       return @doc.to_s
453     end
454   end
455
456   class API
457     def get_xml_doc
458       doc = XML::Document.new
459       doc.encoding = XML::Encoding::UTF_8
460       root = XML::Node.new 'osm'
461       root['version'] = API_VERSION.to_s
462       root['generator'] = GENERATOR
463       root['copyright'] = COPYRIGHT_OWNER
464       root['attribution'] = ATTRIBUTION_URL
465       root['license'] =  LICENSE_URL
466       doc.root = root
467       return doc
468     end
469   end
470
471   def self.IPToCountry(ip_address)
472     Timer.timeout(4) do
473       ipinfo = Quova::IpInfo.new(ip_address)
474
475       if ipinfo.status == Quova::Success then
476         country = ipinfo.country_code
477       else
478         Net::HTTP.start('api.hostip.info') do |http|
479           country = http.get("/country.php?ip=#{ip_address}").body
480           country = "GB" if country == "UK"
481         end
482       end
483       
484       return country.upcase
485     end
486
487     return nil
488   rescue Exception
489     return nil
490   end
491
492   def self.IPLocation(ip_address)
493     code = OSM.IPToCountry(ip_address)
494
495     if code and country = Country.find_by_code(code)
496       return { :minlon => country.min_lon, :minlat => country.min_lat, :maxlon => country.max_lon, :maxlat => country.max_lat }
497     end
498
499     return nil
500   end
501
502   # Construct a random token of a given length
503   def self.make_token(length = 30)
504     chars = 'abcdefghijklmnopqrtuvwxyzABCDEFGHIJKLMNOPQRSTUVWXYZ0123456789'
505     token = ''
506
507     length.times do
508       token += chars[(rand * chars.length).to_i].chr
509     end
510
511     return token
512   end
513
514   # Return an encrypted version of a password
515   def self.encrypt_password(password, salt)
516     return Digest::MD5.hexdigest(password) if salt.nil?
517     return Digest::MD5.hexdigest(salt + password)
518   end
519
520   # Return an SQL fragment to select a given area of the globe
521   def self.sql_for_area(bbox, prefix = nil)
522     tilesql = QuadTile.sql_for_area(bbox, prefix)
523     bbox = bbox.to_scaled
524
525     return "#{tilesql} AND #{prefix}latitude BETWEEN #{bbox.min_lat} AND #{bbox.max_lat} " +
526                       "AND #{prefix}longitude BETWEEN #{bbox.min_lon} AND #{bbox.max_lon}"
527   end
528
529   def self.legal_text_for_country(country_code)
530     file_name = File.join(Rails.root, "config", "legales", country_code.to_s + ".yml")
531     file_name = File.join(Rails.root, "config", "legales", DEFAULT_LEGALE + ".yml") unless File.exist? file_name
532     YAML::load_file(file_name)
533   end
534 end