#!/usr/bin/python

import string
from LatLongUTMconversion import LLtoUTM, UTMtoLL

#TRK 44.738160 -74.787062
track = map(string.split, open("pots08ne.track").readlines())
track = map(lambda a:(float(a[1]),float(a[2])), track)
track = map(lambda a:LLtoUTM(23, a[0], a[1])[1:3], track)
for x,y in track:
    print x,y

svg = open("pots08ne.svg")
state = 0
path = []
for line in svg:
    fields = line.rstrip().split()
    if not fields: continue   # ignore blank lines
    if state == 0:
        if fields[0] == "C":
            state = 1
    elif state == 1:
        if fields[0] == "C":
            path.append(fields[1])
            state = 2
    elif state == 2:
        if len(fields) == 4:
            path.append(fields[0])
            path.append(fields[1])
            state = 3
        else:
            path.append(fields[0])

path = map(lambda a:a.split(","), path)
path = map(lambda a:(float(a[0]),float(a[1])), path)
for x,y in path:
    print x,y

if len(path) != len(track): raise "lengths different"
for i in xrange(len(track)-1):
    trackslope = (track[i][0] - track[i+1][0]) / (-track[i][1] + track[i+1][1])
    pathslope = (path[i][0] - path[i+1][0]) / (path[i][1] - path[i+1][1])

    print "%12.5f %12.5f %12.5f" % (trackslope, pathslope, trackslope/pathslope)