tests/steps/terrain.py

   1 from lettuce import *
   2 from nose.tools import *
   3 import logging
   4 import os
   5 import subprocess
   6 import psycopg2
   7 import re
   8 from haversine import haversine
   9 from shapely.wkt import loads as wkt_load
  10 from shapely.ops import linemerge
  11
  12 logger = logging.getLogger(__name__)
  13
  14 class NominatimConfig:
  15
  16     def __init__(self):
  17         # logging setup
  18         loglevel = getattr(logging, os.environ.get('LOGLEVEL','info').upper())
  19         if 'LOGFILE' in os.environ:
  20             logging.basicConfig(filename=os.environ.get('LOGFILE','run.log'),
  21                                 level=loglevel)
  22         else:
  23             logging.basicConfig(level=loglevel)
  24         # Nominatim test setup
  25         self.base_url = os.environ.get('NOMINATIM_SERVER', 'http://localhost/nominatim')
  26         self.source_dir = os.path.abspath(os.environ.get('NOMINATIM_DIR', '../Nominatim'))
  27         self.template_db = os.environ.get('TEMPLATE_DB', 'test_template_nominatim')
  28         self.test_db = os.environ.get('TEST_DB', 'test_nominatim')
  29         self.local_settings_file = os.environ.get('NOMINATIM_SETTINGS', '/tmp/nominatim_settings.php')
  30         self.reuse_template = 'NOMINATIM_REUSE_TEMPLATE' in os.environ
  31         self.keep_scenario_db = 'NOMINATIM_KEEP_SCENARIO_DB' in os.environ
  32         os.environ['NOMINATIM_SETTINGS'] = '/tmp/nominatim_settings.php'
  33
  34         scriptpath = os.path.dirname(os.path.abspath(__file__))
  35         self.scene_path = os.environ.get('SCENE_PATH',
  36                 os.path.join(scriptpath, '..', 'scenes', 'data'))
  37
  38
  39     def __str__(self):
  40         return 'Server URL: %s\nSource dir: %s\n' % (self.base_url, self.source_dir)
  41
  42 world.config = NominatimConfig()
  43
  44 @world.absorb
  45 def write_nominatim_config(dbname):
  46     f = open(world.config.local_settings_file, 'w')
  47     f.write("<?php\n  @define('CONST_Database_DSN', 'pgsql://@/%s');\n" % dbname)
  48     f.close()
  49
  50
  51 @world.absorb
  52 def run_nominatim_script(script, *args):
  53     cmd = [os.path.join(world.config.source_dir, 'utils', '%s.php' % script)]
  54     cmd.extend(['--%s' % x for x in args])
  55     proc = subprocess.Popen(cmd, stdout=subprocess.PIPE, stderr=subprocess.PIPE)
  56     (outp, outerr) = proc.communicate()
  57     assert (proc.returncode == 0), "Script '%s' failed:\n%s\n%s\n" % (script, outp, outerr)
  58
  59 @world.absorb
  60 def make_hash(inp):
  61     return eval('{' + inp + '}')
  62
  63 @world.absorb
  64 def split_id(oid):
  65     """ Splits a unique identifier for places into its components.
  66         As place_ids cannot be used for testing, we use a unique
  67         identifier instead that is of the form <osmtype><osmid>[:class].
  68     """
  69     oid = oid.strip()
  70     if oid == 'None':
  71         return None, None, None
  72     osmtype = oid[0]
  73     assert_in(osmtype, ('R','N','W'))
  74     if ':' in oid:
  75         osmid, cls = oid[1:].split(':')
  76         return (osmtype, int(osmid), cls)
  77     else:
  78         return (osmtype, int(oid[1:]), None)
  79
  80 @world.absorb
  81 def get_placeid(oid):
  82     """ Tries to retrive the place_id for a unique identifier. """
  83     if oid[0].isdigit():
  84         return int(oid)
  85
  86     osmtype, osmid, cls = world.split_id(oid)
  87     if osmtype is None:
  88         return None
  89     cur = world.conn.cursor()
  90     if cls is None:
  91         q = 'SELECT place_id FROM placex where osm_type = %s and osm_id = %s'
  92         params = (osmtype, osmid)
  93     else:
  94         q = 'SELECT place_id FROM placex where osm_type = %s and osm_id = %s and class = %s'
  95         params = (osmtype, osmid, cls)
  96     cur.execute(q, params)
  97     assert_equals (cur.rowcount, 1)
  98     return cur.fetchone()[0]
  99
 100
 101 @world.absorb
 102 def match_geometry(coord, matchstring):
 103     m = re.match(r'([-0-9.]+),\s*([-0-9.]+)\s*(?:\+-([0-9.]+)([a-z]+)?)?', matchstring)
 104     assert_is_not_none(m, "Invalid match string")
 105
 106     logger.debug("Distmatch: %s/%s %s %s" % (m.group(1), m.group(2), m.group(3), m.group(4) ))
 107     dist = haversine(coord, (float(m.group(1)), float(m.group(2))))
 108
 109     if m.group(3) is not None:
 110         expdist = float(m.group(3))
 111         if m.group(4) is not None:
 112             if m.group(4) == 'm':
 113                 expdist = expdist/1000
 114             elif m.group(4) == 'km':
 115                 pass
 116             else:
 117                 raise Exception("Unknown unit '%s' in geometry match" % (m.group(4), ))
 118     else:
 119         expdist = 0
 120
 121     logger.debug("Distances expected: %f, got: %f" % (expdist, dist))
 122     assert dist <= expdist, "Geometry too far away, expected: %f, got: %f" % (expdist, dist)
 123
 124
 125
 126 @world.absorb
 127 def db_dump_table(table):
 128     cur = world.conn.cursor()
 129     cur.execute('SELECT * FROM %s' % table)
 130     print '<<<<<<< BEGIN OF TABLE DUMP %s' % table
 131     for res in cur:
 132             print res
 133     print '<<<<<<< END OF TABLE DUMP %s' % table
 134
 135 @world.absorb
 136 def db_drop_database(name):
 137     conn = psycopg2.connect(database='postgres')
 138     conn.set_isolation_level(0)
 139     cur = conn.cursor()
 140     cur.execute('DROP DATABASE IF EXISTS %s' % (name, ))
 141     conn.close()
 142
 143
 144 world.is_template_set_up = False
 145
 146 @world.absorb
 147 def db_template_setup():
 148     """ Set up a template database, containing all tables
 149         but not yet any functions.
 150     """
 151     if world.is_template_set_up:
 152         return
 153
 154     world.is_template_set_up = True
 155     world.write_nominatim_config(world.config.template_db)
 156     if world.config.reuse_template:
 157         # check that the template is there
 158         conn = psycopg2.connect(database='postgres')
 159         cur = conn.cursor()
 160         cur.execute('select count(*) from pg_database where datname = %s',
 161                      (world.config.template_db,))
 162         if cur.fetchone()[0] == 1:
 163             return
 164     else:
 165         # just in case... make sure a previous table has been dropped
 166         world.db_drop_database(world.config.template_db)
 167     # call the first part of database setup
 168     world.run_nominatim_script('setup', 'create-db', 'setup-db')
 169     # remove external data to speed up indexing for tests
 170     conn = psycopg2.connect(database=world.config.template_db)
 171     psycopg2.extras.register_hstore(conn, globally=False, unicode=True)
 172     cur = conn.cursor()
 173     for table in ('gb_postcode', 'us_postcode', 'us_state', 'us_statecounty'):
 174         cur.execute('TRUNCATE TABLE %s' % (table,))
 175     conn.commit()
 176     conn.close()
 177     # execute osm2pgsql on an empty file to get the right tables
 178     osm2pgsql = os.path.join(world.config.source_dir, 'osm2pgsql', 'osm2pgsql')
 179     proc = subprocess.Popen([osm2pgsql, '-lsc', '-O', 'gazetteer', '-d', world.config.template_db, '-'],
 180     stdin=subprocess.PIPE, stdout=subprocess.PIPE, stderr=subprocess.PIPE)
 181     [outstr, errstr] = proc.communicate(input='<osm version="0.6"></osm>')
 182     world.run_nominatim_script('setup', 'create-functions', 'create-tables', 'create-partition-tables', 'create-partition-functions', 'load-data', 'create-search-indices')
 183
 184
 185 # Leave the table around so it can be reused again after a non-reuse test round.
 186 #@after.all
 187 def db_template_teardown(total):
 188     """ Set up a template database, containing all tables
 189         but not yet any functions.
 190     """
 191     if world.is_template_set_up:
 192         # remove template DB
 193         if not world.config.reuse_template:
 194             world.db_drop_database(world.config.template_db)
 195         try:
 196             os.remove(world.config.local_settings_file)
 197         except OSError:
 198             pass # ignore missing file
 199
 200
 201 ##########################################################################
 202 #
 203 # Data scene handling
 204 #
 205
 206 world.scenes = {}
 207 world.current_scene = None
 208
 209 @world.absorb
 210 def load_scene(name):
 211     if name in world.scenes:
 212         world.current_scene = world.scenes[name]
 213     else:
 214         with open(os.path.join(world.config.scene_path, "%s.wkt" % name), 'r') as fd:
 215             scene = {}
 216             for line in fd:
 217                 if line.strip():
 218                     obj, wkt = line.split('|', 2)
 219                     wkt = wkt.strip()
 220                     scene[obj.strip()] = wkt_load(wkt)
 221             world.scenes[name] = scene
 222             world.current_scene = scene
 223
 224 @world.absorb
 225 def get_scene_geometry(name):
 226     if not ':' in name:
 227         # Not a scene description
 228         return None
 229
 230     geoms = []
 231     for obj in name.split('+'):
 232         oname = obj.strip()
 233         if oname.startswith(':'):
 234             geoms.append(world.current_scene[oname[1:]])
 235         else:
 236             scene, obj = oname.split(':', 2)
 237             oldscene = world.current_scene
 238             world.load_scene(scene)
 239             wkt = world.current_scene[obj]
 240             world.current_scene = oldscene
 241             geoms.append(wkt)
 242
 243     if len(geoms) == 1:
 244         return geoms[0]
 245     else:
 246         return linemerge(geoms)