wradmin/controllers/coordtool.py

   1 import re
   2 import xml.dom.minidom as minidom
   3 from xml.parsers.expat import ExpatError
   4 from enum import IntEnum
   5 from flask import request, redirect, url_for, flash, render_template
   6
   7
   8 class MultiGeo(IntEnum):
   9     """Formats multiple coordinates, even in multiple lines to [(latitude, longitude, elevation), ...]
  10     or [(latitude, longitude, None), ...] tuples."""
  11
  12     # Valid for input_format
  13     FORMAT_GUESS = 0  # guesses the input format; default for input_format
  14     FORMAT_NONE = -1  # indicates missing formats
  15
  16     # Valid for input_format and output_format
  17     FORMAT_GEOCACHING = 1  # e.g. "N 47° 13.692 E 011° 25.535"
  18     FORMAT_WINTERRODELN = 2  # e.g. "47.222134 N 11.467211 E"
  19     FORMAT_GMAPPLUGIN = 3  # e.g. "47.232922, 11.452239"
  20     FORMAT_GPX = 4  # e.g. "<trkpt lat="47.181289" lon="11.408827"><ele>1090.57</ele></trkpt>"
  21
  22
  23 def multigeo_from_string(value, input_format):
  24     assert input_format in [MultiGeo.FORMAT_GUESS, MultiGeo.FORMAT_GEOCACHING, MultiGeo.FORMAT_WINTERRODELN,
  25                             MultiGeo.FORMAT_GMAPPLUGIN, MultiGeo.FORMAT_GPX], "input_format %d is not recognized"
  26     lines = [line.strip() for line in value.split("\n") if len(line.strip()) > 0]
  27
  28     result = []
  29     for line in lines:
  30         if input_format == MultiGeo.FORMAT_GUESS or input_format == MultiGeo.FORMAT_GEOCACHING:
  31             r = re.match(r'N ?(\d+)° ?(\d+\.\d+) +E ?(\d+)° ?(\d+\.\d+)', line)
  32             if r is not None:
  33                 g = r.groups()
  34                 result.append((float(g[0]) + float(g[1]) / 60, float(g[2]) + float(g[3]) / 60, None))
  35                 continue
  36
  37         if input_format == MultiGeo.FORMAT_GUESS or input_format == MultiGeo.FORMAT_WINTERRODELN:
  38             r = re.match(r'(\d+\.\d+) N (\d+\.\d+) E', line)
  39             if not r is None:
  40                 result.append((float(r.groups()[0]), float(r.groups()[1]), None))
  41                 continue
  42
  43         if input_format == MultiGeo.FORMAT_GUESS or input_format == MultiGeo.FORMAT_GMAPPLUGIN:
  44             r = re.match(r'(\d+\.\d+), ?(\d+\.\d+)', line)
  45             if r is not None:
  46                 result.append((float(r.groups()[0]), float(r.groups()[1]), None))
  47                 continue
  48
  49         if input_format == MultiGeo.FORMAT_GUESS or input_format == MultiGeo.FORMAT_GPX:
  50             try:
  51                 xml = minidom.parseString(line)
  52                 coord = xml.documentElement
  53                 lat = float(coord.getAttribute('lat'))
  54                 lon = float(coord.getAttribute('lon'))
  55                 try:
  56                     ele = float(coord.childNodes[0].childNodes[0].nodeValue)
  57                 except (IndexError, ValueError):
  58                     ele = None
  59                 result.append((lat, lon, ele))
  60                 continue
  61             except (ExpatError, IndexError, ValueError):
  62                 pass
  63
  64         raise ValueError("Coordinates '%s' have no known format" % line, value)
  65     return result
  66
  67
  68 def multigeo_to_string(value, output_format):
  69     result = []
  70     for latitude, longitude, height in value:
  71         if output_format == MultiGeo.FORMAT_GEOCACHING:
  72             result.append(
  73                 'N %02d° %02.3f E %03d° %02.3f' % (latitude, latitude % 1 * 60, longitude, longitude % 1 * 60))
  74
  75         elif output_format == MultiGeo.FORMAT_WINTERRODELN:
  76             result.append('%.6f N %.6f E' % (latitude, longitude))
  77
  78         elif output_format == MultiGeo.FORMAT_GMAPPLUGIN:
  79             result.append('%.6f, %.6f' % (latitude, longitude))
  80
  81         elif output_format == MultiGeo.FORMAT_GPX:
  82             if height is not None:
  83                 result.append(
  84                     '<trkpt lat="%.6f" lon="%.6f"><ele>%.2f</ele></trkpt>' % (latitude, longitude, height))
  85             else:
  86                 result.append('<trkpt lat="%.6f" lon="%.6f"/>' % (latitude, longitude))
  87
  88         else:
  89             raise ValueError("output_format %d is not recognized" % output_format, value)
  90
  91     return "\n".join(result)
  92
  93
  94 class CoordtoolController:
  95
  96     def index(self):
  97         return render_template('coordtool.html', result=None)
  98
  99     def convert(self):
 100         assert request.method == 'POST'
 101         input = request.form['input']
 102         no_elevation = 'no_elevation' in request.form
 103         simplify = 'simplify' in request.form
 104         swap_latlon = 'swap_latlon' in request.form
 105         no_geoformat = 'no_geoformat' in request.form
 106         no_gpxformat = 'no_gpxformat' in request.form
 107         no_gmapsformat = 'no_gmapsformat' in request.form
 108         no_geocachingformat = 'no_geocachingformat' in request.form
 109
 110         if input is None or len(input.strip()) == 0:
 111             return redirect(url_for('coordtool_index'))
 112
 113         try:
 114             result = multigeo_from_string(input, MultiGeo.FORMAT_GUESS)
 115         except ValueError as e:
 116             flash(str(e))
 117             return redirect(url_for('coordtool_index'))
 118
 119         if swap_latlon:
 120             result = [(latitude, longitude, elevation) for (longitude, latitude, elevation) in result]
 121
 122         if no_elevation:
 123             result = [(longitude, latitude, None) for (longitude, latitude, elevation) in result]
 124
 125         geo_winterrodeln = lambda value: multigeo_to_string(value, MultiGeo.FORMAT_WINTERRODELN)
 126         geo_gmapplugin = lambda value: multigeo_to_string(value, MultiGeo.FORMAT_GMAPPLUGIN)
 127         geo_gpx = lambda value: multigeo_to_string(value, MultiGeo.FORMAT_GPX)
 128         geo_geocaching = lambda value: multigeo_to_string(value, MultiGeo.FORMAT_GEOCACHING)
 129
 130         return render_template('coordtool.html', result=result,
 131                                no_geoformat=no_geoformat, geo_winterrodeln=geo_winterrodeln,
 132                                no_gmapsformat=no_gmapsformat, geo_gmapplugin=geo_gmapplugin,
 133                                no_gpxformat=no_gpxformat, geo_gpx=geo_gpx,
 134                                no_geocachingformat=no_geocachingformat, geo_geocaching=geo_geocaching)