app/controllers/geocoder_controller.rb

   1 # coding: utf-8
   2
   3 class GeocoderController < ApplicationController
   4   require 'uri'
   5   require 'net/http'
   6   require 'rexml/document'
   7
   8   before_filter :authorize_web
   9   before_filter :set_locale
  10   before_filter :require_oauth, :only => [:search]
  11
  12   def search
  13     normalize_params
  14
  15     @sources = []
  16     if params[:lat] && params[:lon]
  17       @sources.push "latlon"
  18       @sources.push "osm_nominatim_reverse"
  19       @sources.push "geonames_reverse" if defined?(GEONAMES_USERNAME)
  20     elsif params[:query].match(/^\d{5}(-\d{4})?$/)
  21       @sources.push "us_postcode"
  22       @sources.push "osm_nominatim"
  23     elsif params[:query].match(/^(GIR 0AA|[A-PR-UWYZ]([0-9]{1,2}|([A-HK-Y][0-9]|[A-HK-Y][0-9]([0-9]|[ABEHMNPRV-Y]))|[0-9][A-HJKS-UW])\s*[0-9][ABD-HJLNP-UW-Z]{2})$/i)
  24       @sources.push "uk_postcode"
  25       @sources.push "osm_nominatim"
  26     elsif params[:query].match(/^[A-Z]\d[A-Z]\s*\d[A-Z]\d$/i)
  27       @sources.push "ca_postcode"
  28       @sources.push "osm_nominatim"
  29     else
  30       @sources.push "osm_nominatim"
  31       @sources.push "geonames" if defined?(GEONAMES_USERNAME)
  32     end
  33
  34     render :layout => map_layout
  35   end
  36
  37   def search_latlon
  38     lat = params[:lat].to_f
  39     lon = params[:lon].to_f
  40     if lat < -90 or lat > 90
  41       @error = "Latitude #{lat} out of range"
  42       render :action => "error"
  43     elsif lon < -180 or lon > 180
  44       @error = "Longitude #{lon} out of range"
  45       render :action => "error"
  46     else
  47       @results = [{:lat => lat, :lon => lon,
  48                    :zoom => params[:zoom],
  49                    :name => "#{lat}, #{lon}"}]
  50
  51       render :action => "results"
  52     end
  53   end
  54
  55   def search_us_postcode
  56     # get query parameters
  57     query = params[:query]
  58
  59     # create result array
  60     @results = Array.new
  61
  62     # ask geocoder.us (they have a non-commercial use api)
  63     response = fetch_text("http://rpc.geocoder.us/service/csv?zip=#{escape_query(query)}")
  64
  65     # parse the response
  66     unless response.match(/couldn't find this zip/)
  67       data = response.split(/\s*,\s+/) # lat,long,town,state,zip
  68       @results.push({:lat => data[0], :lon => data[1],
  69                      :zoom => POSTCODE_ZOOM,
  70                      :prefix => "#{data[2]}, #{data[3]},",
  71                      :name => data[4]})
  72     end
  73
  74     render :action => "results"
  75   rescue Exception => ex
  76     @error = "Error contacting rpc.geocoder.us: #{ex.to_s}"
  77     render :action => "error"
  78   end
  79
  80   def search_uk_postcode
  81     # get query parameters
  82     query = params[:query]
  83
  84     # create result array
  85     @results = Array.new
  86
  87     # ask npemap.org.uk to do a combined npemap + freethepostcode search
  88     response = fetch_text("http://www.npemap.org.uk/cgi/geocoder.fcgi?format=text&postcode=#{escape_query(query)}")
  89
  90     # parse the response
  91     unless response.match(/Error/)
  92       dataline = response.split(/\n/)[1]
  93       data = dataline.split(/,/) # easting,northing,postcode,lat,long
  94       postcode = data[2].gsub(/'/, "")
  95       zoom = POSTCODE_ZOOM - postcode.count("#")
  96       @results.push({:lat => data[3], :lon => data[4], :zoom => zoom,
  97                      :name => postcode})
  98     end
  99
 100     render :action => "results"
 101   rescue Exception => ex
 102     @error = "Error contacting www.npemap.org.uk: #{ex.to_s}"
 103     render :action => "error"
 104   end
 105
 106   def search_ca_postcode
 107     # get query parameters
 108     query = params[:query]
 109     @results = Array.new
 110
 111     # ask geocoder.ca (note - they have a per-day limit)
 112     response = fetch_xml("http://geocoder.ca/?geoit=XML&postal=#{escape_query(query)}")
 113
 114     # parse the response
 115     if response.get_elements("geodata/error").empty?
 116       @results.push({:lat => response.get_text("geodata/latt").to_s,
 117                      :lon => response.get_text("geodata/longt").to_s,
 118                      :zoom => POSTCODE_ZOOM,
 119                      :name => query.upcase})
 120     end
 121
 122     render :action => "results"
 123   rescue Exception => ex
 124     @error = "Error contacting geocoder.ca: #{ex.to_s}"
 125     render :action => "error"
 126   end
 127
 128   def search_osm_nominatim
 129     # get query parameters
 130     query = params[:query]
 131     minlon = params[:minlon]
 132     minlat = params[:minlat]
 133     maxlon = params[:maxlon]
 134     maxlat = params[:maxlat]
 135
 136     # get view box
 137     if minlon && minlat && maxlon && maxlat
 138       viewbox = "&viewbox=#{minlon},#{maxlat},#{maxlon},#{minlat}"
 139     end
 140
 141     # get objects to excude
 142     if params[:exclude]
 143       exclude = "&exclude_place_ids=#{params[:exclude].join(',')}"
 144     end
 145
 146     # ask nominatim
 147     response = fetch_xml("#{NOMINATIM_URL}search?format=xml&q=#{escape_query(query)}#{viewbox}#{exclude}&accept-language=#{http_accept_language.user_preferred_languages.join(',')}")
 148
 149     # create result array
 150     @results = Array.new
 151
 152     # create parameter hash for "more results" link
 153     @more_params = params.reverse_merge({ :exclude => [] })
 154
 155     # extract the results from the response
 156     results =  response.elements["searchresults"]
 157
 158     # parse the response
 159     results.elements.each("place") do |place|
 160       lat = place.attributes["lat"].to_s
 161       lon = place.attributes["lon"].to_s
 162       klass = place.attributes["class"].to_s
 163       type = place.attributes["type"].to_s
 164       name = place.attributes["display_name"].to_s
 165       min_lat,max_lat,min_lon,max_lon = place.attributes["boundingbox"].to_s.split(",")
 166       if type.empty?
 167         prefix_name = ""
 168       else
 169         prefix_name = t "geocoder.search_osm_nominatim.prefix.#{klass}.#{type}", :default => type.gsub("_", " ").capitalize
 170       end
 171       if klass == 'boundary' and type == 'administrative'
 172         rank = (place.attributes["place_rank"].to_i + 1) / 2
 173         prefix_name = t "geocoder.search_osm_nominatim.admin_levels.level#{rank}", :default => prefix_name
 174       end
 175       prefix = t "geocoder.search_osm_nominatim.prefix_format", :name => prefix_name
 176       object_type = place.attributes["osm_type"]
 177       object_id = place.attributes["osm_id"]
 178
 179       @results.push({:lat => lat, :lon => lon,
 180                      :min_lat => min_lat, :max_lat => max_lat,
 181                      :min_lon => min_lon, :max_lon => max_lon,
 182                      :prefix => prefix, :name => name,
 183                      :type => object_type, :id => object_id})
 184       @more_params[:exclude].push(place.attributes["place_id"].to_s)
 185     end
 186
 187     render :action => "results"
 188 #  rescue Exception => ex
 189 #    @error = "Error contacting nominatim.openstreetmap.org: #{ex.to_s}"
 190 #    render :action => "error"
 191   end
 192
 193   def search_geonames
 194     # get query parameters
 195     query = params[:query]
 196
 197     # get preferred language
 198     lang = I18n.locale.to_s.split("-").first
 199
 200     # create result array
 201     @results = Array.new
 202
 203     # ask geonames.org
 204     response = fetch_xml("http://api.geonames.org/search?q=#{escape_query(query)}&lang=#{lang}&maxRows=20&username=#{GEONAMES_USERNAME}")
 205
 206     # parse the response
 207     response.elements.each("geonames/geoname") do |geoname|
 208       lat = geoname.get_text("lat").to_s
 209       lon = geoname.get_text("lng").to_s
 210       name = geoname.get_text("name").to_s
 211       country = geoname.get_text("countryName").to_s
 212       @results.push({:lat => lat, :lon => lon,
 213                      :zoom => GEONAMES_ZOOM,
 214                      :name => name,
 215                      :suffix => ", #{country}"})
 216     end
 217
 218     render :action => "results"
 219   rescue Exception => ex
 220     @error = "Error contacting ws.geonames.org: #{ex.to_s}"
 221     render :action => "error"
 222   end
 223
 224   def search_osm_nominatim_reverse
 225     # get query parameters
 226     lat = params[:lat]
 227     lon = params[:lon]
 228     zoom = params[:zoom]
 229
 230     # create result array
 231     @results = Array.new
 232
 233     # ask nominatim
 234     response = fetch_xml("#{NOMINATIM_URL}reverse?lat=#{lat}&lon=#{lon}&zoom=#{zoom}&accept-language=#{http_accept_language.user_preferred_languages.join(',')}")
 235
 236     # parse the response
 237     response.elements.each("reversegeocode/result") do |result|
 238       lat = result.attributes["lat"].to_s
 239       lon = result.attributes["lon"].to_s
 240       object_type = result.attributes["osm_type"]
 241       object_id = result.attributes["osm_id"]
 242       description = result.get_text.to_s
 243
 244       @results.push({:lat => lat, :lon => lon,
 245                      :zoom => zoom,
 246                      :name => description,
 247                      :type => object_type, :id => object_id})
 248     end
 249
 250     render :action => "results"
 251   rescue Exception => ex
 252     @error = "Error contacting nominatim.openstreetmap.org: #{ex.to_s}"
 253     render :action => "error"
 254   end
 255
 256   def search_geonames_reverse
 257     # get query parameters
 258     lat = params[:lat]
 259     lon = params[:lon]
 260
 261     # get preferred language
 262     lang = I18n.locale.to_s.split("-").first
 263
 264     # create result array
 265     @results = Array.new
 266
 267     # ask geonames.org
 268     response = fetch_xml("http://api.geonames.org/countrySubdivision?lat=#{lat}&lng=#{lon}&lang=#{lang}&username=#{GEONAMES_USERNAME}")
 269
 270     # parse the response
 271     response.elements.each("geonames/countrySubdivision") do |geoname|
 272       name = geoname.get_text("adminName1").to_s
 273       country = geoname.get_text("countryName").to_s
 274       @results.push({:lat => lat, :lon => lon,
 275                      :zoom => GEONAMES_ZOOM,
 276                      :name => name,
 277                      :suffix => ", #{country}"})
 278     end
 279
 280     render :action => "results"
 281   rescue Exception => ex
 282     @error = "Error contacting ws.geonames.org: #{ex.to_s}"
 283     render :action => "error"
 284   end
 285
 286 private
 287
 288   def fetch_text(url)
 289     return Net::HTTP.get(URI.parse(url))
 290   end
 291
 292   def fetch_xml(url)
 293     return REXML::Document.new(fetch_text(url))
 294   end
 295
 296   def format_distance(distance)
 297     return t("geocoder.distance", :count => distance)
 298   end
 299
 300   def format_direction(bearing)
 301     return t("geocoder.direction.south_west") if bearing >= 22.5 and bearing < 67.5
 302     return t("geocoder.direction.south") if bearing >= 67.5 and bearing < 112.5
 303     return t("geocoder.direction.south_east") if bearing >= 112.5 and bearing < 157.5
 304     return t("geocoder.direction.east") if bearing >= 157.5 and bearing < 202.5
 305     return t("geocoder.direction.north_east") if bearing >= 202.5 and bearing < 247.5
 306     return t("geocoder.direction.north") if bearing >= 247.5 and bearing < 292.5
 307     return t("geocoder.direction.north_west") if bearing >= 292.5 and bearing < 337.5
 308     return t("geocoder.direction.west")
 309   end
 310
 311   def format_name(name)
 312     return name.gsub(/( *\[[^\]]*\])*$/, "")
 313   end
 314
 315   def count_results(results)
 316     count = 0
 317
 318     results.each do |source|
 319       count += source[:results].length if source[:results]
 320     end
 321
 322     return count
 323   end
 324
 325   def escape_query(query)
 326     return URI.escape(query, Regexp.new("[^#{URI::PATTERN::UNRESERVED}]", false, 'N'))
 327   end
 328
 329   def normalize_params
 330     query = params[:query]
 331     return unless query
 332
 333     query.strip!
 334
 335     if latlon = query.match(/^([NS])\s*(\d{1,3}(\.\d*)?)\W*([EW])\s*(\d{1,3}(\.\d*)?)$/).try(:captures) # [NSEW] decimal degrees
 336       params.merge!(nsew_to_decdeg(latlon)).delete(:query)
 337     elsif latlon = query.match(/^(\d{1,3}(\.\d*)?)\s*([NS])\W*(\d{1,3}(\.\d*)?)\s*([EW])$/).try(:captures) # decimal degrees [NSEW]
 338       params.merge!(nsew_to_decdeg(latlon)).delete(:query)
 339
 340     elsif latlon = query.match(/^([NS])\s*(\d{1,3})°?\s*(\d{1,3}(\.\d*)?)?['′]?\W*([EW])\s*(\d{1,3})°?\s*(\d{1,3}(\.\d*)?)?['′]?$/).try(:captures) # [NSEW] degrees, decimal minutes
 341       params.merge!(ddm_to_decdeg(latlon)).delete(:query)
 342     elsif latlon = query.match(/^(\d{1,3})°?\s*(\d{1,3}(\.\d*)?)?['′]?\s*([NS])\W*(\d{1,3})°?\s*(\d{1,3}(\.\d*)?)?['′]?\s*([EW])$/).try(:captures) # degrees, decimal minutes [NSEW]
 343       params.merge!(ddm_to_decdeg(latlon)).delete(:query)
 344
 345     elsif latlon = query.match(/^([NS])\s*(\d{1,3})°?\s*(\d{1,2})['′]?\s*(\d{1,3}(\.\d*)?)?["″]?\W*([EW])\s*(\d{1,3})°?\s*(\d{1,2})['′]?\s*(\d{1,3}(\.\d*)?)?["″]?$/).try(:captures) # [NSEW] degrees, minutes, decimal seconds
 346       params.merge!(dms_to_decdeg(latlon)).delete(:query)
 347     elsif latlon = query.match(/^(\d{1,3})°?\s*(\d{1,2})['′]?\s*(\d{1,3}(\.\d*)?)?["″]\s*([NS])\W*(\d{1,3})°?\s*(\d{1,2})['′]?\s*(\d{1,3}(\.\d*)?)?["″]?\s*([EW])$/).try(:captures) # degrees, minutes, decimal seconds [NSEW]
 348       params.merge!(dms_to_decdeg(latlon)).delete(:query)
 349
 350     elsif latlon = query.match(/^\s*([+-]?\d+(\.\d*)?)\s*[\s,]\s*([+-]?\d+(\.\d*)?)\s*$/)
 351       params.merge!({:lat => latlon[1].to_f, :lon => latlon[3].to_f}).delete(:query)
 352     end
 353   end
 354
 355   def nsew_to_decdeg(captures)
 356     begin
 357       Float(captures[0])
 358       captures[2].downcase != 's' ? lat = captures[0].to_f : lat = -(captures[0].to_f)
 359       captures[5].downcase != 'w' ? lon = captures[3].to_f : lon = -(captures[3].to_f)
 360     rescue
 361       captures[0].downcase != 's' ? lat = captures[1].to_f : lat = -(captures[1].to_f)
 362       captures[3].downcase != 'w' ? lon = captures[4].to_f : lon = -(captures[4].to_f)
 363     end
 364     {:lat => lat, :lon => lon}
 365   end
 366
 367   def ddm_to_decdeg(captures)
 368     begin
 369       Float(captures[0])
 370       captures[3].downcase != 's' ? lat = captures[0].to_f + captures[1].to_f/60 : lat = -(captures[0].to_f + captures[1].to_f/60)
 371       captures[7].downcase != 'w' ? lon = captures[4].to_f + captures[5].to_f/60 : lon = -(captures[4].to_f + captures[5].to_f/60)
 372     rescue
 373       captures[0].downcase != 's' ? lat = captures[1].to_f + captures[2].to_f/60 : lat = -(captures[1].to_f + captures[2].to_f/60)
 374       captures[4].downcase != 'w' ? lon = captures[5].to_f + captures[6].to_f/60 : lon = -(captures[5].to_f + captures[6].to_f/60)
 375     end
 376     {:lat => lat, :lon => lon}
 377   end
 378
 379   def dms_to_decdeg(captures)
 380     begin
 381       Float(captures[0])
 382       captures[4].downcase != 's' ? lat = captures[0].to_f + (captures[1].to_f + captures[2].to_f/60)/60 : lat = -(captures[0].to_f + (captures[1].to_f + captures[2].to_f/60)/60)
 383       captures[9].downcase != 'w' ? lon = captures[5].to_f + (captures[6].to_f + captures[7].to_f/60)/60 : lon = -(captures[5].to_f + (captures[6].to_f + captures[7].to_f/60)/60)
 384     rescue
 385       captures[0].downcase != 's' ? lat = captures[1].to_f + (captures[2].to_f + captures[3].to_f/60)/60 : lat = -(captures[1].to_f + (captures[2].to_f + captures[3].to_f/60)/60)
 386       captures[5].downcase != 'w' ? lon = captures[6].to_f + (captures[7].to_f + captures[8].to_f/60)/60 : lon = -(captures[6].to_f + (captures[7].to_f + captures[8].to_f/60)/60)
 387     end
 388     {:lat => lat, :lon => lon}
 389   end
 390
 391 end