# cast42

Cast for two

## Thursday, July 19, 2007

### Calculating the length of Google Maps line segments

The output of the above script is:

distance of fietsweg is 36.35 kilometer

Of course, this calculation is only an approximation since the trajectory is only an approximation (albeit a pretty good one) of the real route I followed and no heigth information is taken into account. Still I guess that its accurate upto a few kilometers. If somebody has a spare Bluetooth GPS receiver, you can always send me one ;-).

The calculation of the distance in kilometer between two points is based on the code in Ruby and the examples from the zips project a sourceforge. The coordinates of the linesegments are obtained from the KML file that Google Maps generates when appending &output=kml to the URL from 'link to this page'. The KML file is parsed using the Python minidom module. I hope that in the future Google will provide a "&output=jsonp&callback=calcfunc" so that a simple webpage with some Javascript can calculate the length. But that's for another programming project.

Dikkie said...

How do I run such a script myself?

cast42 said...
This comment has been removed by the author.
cast42 said...

@dikkie download to a file with extention .py and feed the file to the python interpreter. On osx Python is installed by default, on windoze it's a bit more difficult. But as I explained, I'm working on a more user friendly version that boils down to pasting an URL and pushing a button. In the mean while, you can calculate distances here So stay tuned.

Dikkie said...

Thanks a lot. Btw, I also found this one.

cast42 said...

That's also a nice one, but I would like to create the routes with Google Maps and then be able to calculate the distance. A such I can keep the map stored. The drawback of this mapplet is that you can not store the map after chalking the route on the map.

Glenn said...

cast42,
This tool to save routes on a map seems to do what you are looking for.

cast42 said...

There's now a recipe for this:
Recipe 576782: Calculate distance from .kmz files Anonymous said...

Hi!

I found your app is exactly what I need, however, kml format has probably changed since 2007 or at least I couldn't use it.

I changed your code a little to make it work with the new format

here it is:

#!/usr/bin/python
from xml.dom import minidom
import urllib,sys
import math

# This is the link to the google map
filename=sys.argv;
# Calculate distance based on http://highearthorbit.com/projects/geolocation/gpx.rb
def distance_in_km(lat1, lon1, lat2, lon2):
rtheta = math.radians(lon1 - lon2)
dist = math.sin(rlat1) * math.sin(rlat2) + math.cos(rlat1) * math.cos(rlat2) * math.cos(rtheta)
dist = math.degrees(math.acos(dist))
return (dist * 60.0 * 1.1515 * 1.609344) # distance in km

def getText(nodelist):
rc = ""
for node in nodelist:
if node.nodeType == node.TEXT_NODE:
rc = rc + node.data
return rc

# Retrieve the KML file
#sock =urllib.urlopen(url+"&output=kml")
#sock.close()

f=open(filename);

xmldoc = minidom.parseString(strkml);
placemarks = xmldoc.getElementsByTagName('kml').getElementsByTagName('Document').getElementsByTagName('Placemark')
first=True;
distance=0;
for placemark in placemarks:
name = getText(placemark.getElementsByTagName('name').childNodes)
linestrings = placemark.getElementsByTagName('Point')
for linestring in linestrings:
coordinates = linestring.getElementsByTagName('coordinates').childNodes
listofpoints = getText(coordinates)
points = listofpoints.split(',')
#distance = 0.0
point=points;
values = point;#.split(',')
if len(values) == 3:
(x,y,z) = float(values), float(values), float(values)
if first:
first = False
(x0,y0,z0) = (x,y,z)
else:
dist = distance_in_km(x0,y0,x,y)
# print 'dist='+str(dist);
distance += dist
x0 = x
y0 = y
else:
print 'len points='+str(len(points));
print "Total distance is %.2f kilometer" % distance;

hope the indentation works fine

mehageg at nana10 dot co dot il

cast42 said...

I put the updated code into an pastebin:
http://pastebin.ca/1886961