Calculate distance from .kmz files (Python recipe) by Michal Niklas

ActiveState Code (http://code.activestate.com/recipes/576782/)

#!/usr/bin/env python
# -*- coding: utf8 -*-

__version__ = '$Id: kmz2dist.py 609 2009-06-01 05:04:52Z mn $'
# author: Michal Niklas

"""
Calculates distance in kilometers from points saved in .kmz file.
Such files can be created by Google Earth
look at:
http://code.google.com/intl/pl-PL/apis/kml/documentation/kmlreference.html

Recipe by Bartek Górny is used to calculate distance between 2 points:
Recipe 576779: Calculating distance between two geographic points
http://code.activestate.com/recipes/576779/
"""

import sys
import zipfile
import glob
from xml.dom import minidom

# save Recipe 576779: Calculating distance between two geographic points
# http://code.activestate.com/recipes/576779/
# as distance
import distance


def get_distance(coordinates_str):
	"""gets distance of one path from coordinates string in form of:
	14.81363432237944,53.57016581501523,0 14.81411766813742,53.56923005549378,0 14.81880340335202,53.56879451890311 ...
	look at:
		http://code.google.com/intl/pl-PL/apis/kml/documentation/kmlreference.html#coordinates
	"""
	sum_distance = 0.0
	arr = []
	coordinates = []
	if ' ' in coordinates_str:
		arr = coordinates_str.split(' ')
	if len(arr) > 1:
		for s in arr:
			if ',' in s:
				pt = s.split(',')
				pos_latt = (float(pt[0].strip()), 0, 0)
				pos_long = (float(pt[1].strip()), 0, 0)
				position = (pos_latt, pos_long)
				coordinates.append(position)
	if coordinates:
		for i in range(len(coordinates) - 1):
			start = coordinates[i]
			stop = coordinates[i + 1]
			sum_distance += distance.points2distance(start, stop)
	return sum_distance


def show_distance(fname):
	"""calculates distance from points saved in doc.kml which is part of .kmz
	   zip archive file"""
	path_cnt = 0
	zf = zipfile.ZipFile(fname, 'r')
	for fn in zf.namelist():
		if fn.endswith('.kml'):
			content = zf.read(fn)
			xmldoc = minidom.parseString(content)
			placemarks = xmldoc.getElementsByTagName('Placemark')
			for placemark in placemarks:
				name = placemark.getElementsByTagName('name')
				if name:
					name = name[0].firstChild.data.strip()
					coordinates = placemark.getElementsByTagName('coordinates')
					if coordinates:
						coordinates = coordinates[0].firstChild.data.strip()
						if coordinates:
							distance_km = get_distance(coordinates)
							if distance_km > 0.0:
								if path_cnt == 0:
									print '\n%s:' % (fname)
								path_cnt += 1
								print('\t%s\t%5.2f' % (name, distance_km))
	return path_cnt


def main():
	"""show .kmz file name and distance in kilometers"""
	fnames = glob.glob('*.kmz')
	if not fnames:
		print('No *.kmz files found')
	else:
		for fname in fnames:
			show_distance(fname)


if __name__ == '__main__':
	if '--version' in sys.argv:
		print(__version__)
	else:
		main()