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